Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 980,284 C→A 6.8% R20S (CGC→AGC)  iscA → copper(I) binding iron‑sulfur cluster assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947980,2840CA 6.8% 228.8 / 6.1 73R20S (CGC→AGC) iscAcopper(I) binding iron‑sulfur cluster assembly protein
Reads supporting (aligned to +/- strand):  ref base C (36/31);  new base A (3/2);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.21e-01

CCGAAGATGCCATCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGATCATCGATCCC  >  NC_002947/980140‑980424
                                                                                                                                                  |                                                                                                                                            
ccGAAGATGCCATCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCGAACCACATTCGCCGCTCCCTCGACGGGCGCgg                                                                                                                                          <  1:235769/149‑1 (MQ=255)
ccGAAGATGCCATCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCgg                                                                                                                                          <  7:222943/149‑1 (MQ=255)
     gATGCCATCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCGCCCGTCGAGGGAGCGGCgaat                                                                                                                                     <  2:302017/149‑2 (MQ=255)
       tGCCATCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGt                                                                                                                                   >  8:196673/1‑149 (MQ=255)
          cATCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAg                                                                                                                                <  3:39343/149‑1 (MQ=255)
           aTCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGcgcg                                                                                                                                           >  2:359255/1‑137 (MQ=255)
           aTCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGcgcg                                                                                                                                           <  1:359255/137‑1 (MQ=255)
           aTCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAgg                                                                                                                               >  3:37407/1‑149 (MQ=255)
           aTCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAgg                                                                                                                               <  7:123104/149‑1 (MQ=255)
           aTCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAgg                                                                                                                               <  2:88738/149‑1 (MQ=255)
            tCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAggg                                                                                                                              >  6:73464/1‑149 (MQ=255)
             cAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGcgcg                                                                                                                                           >  4:48727/1‑135 (MQ=255)
             cAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGcgcg                                                                                                                                           <  3:48727/135‑1 (MQ=255)
                 gCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTc                                                                                                                         >  5:76065/1‑149 (MQ=255)
                      cttaggcgAGTACTAGCAGAAGATAGGCTTGATCTAAGTCGCCTGTTGCAAGGTATGGAGTTCTGTTGGCTCTCAGGATGACAGAAGCCGACGCCAACCACATTCGCCTCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCt                                                                                                                     <  3:280101/143‑1 (MQ=255)
                       gTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTgg                                                                                                                   <  7:290984/149‑1 (MQ=255)
                       gTACGCGACTACAAGCAGAAGAAAGCCTTGATCTACGTCGCCTGTTGCAAGGTAAGGAGGTCTGATGGCTACCAGCATGACAGAAGCCGCCCCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTgg                                                                                                                   <  4:107131/149‑1 (MQ=255)
                           gcgACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGt                                                                                                               <  1:56767/149‑1 (MQ=255)
                             gACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTgc                                                                                                             <  7:202360/149‑1 (MQ=255)
                                   aaGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGcacca                                                                                                       >  5:191883/1‑149 (MQ=255)
                                      cAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCg                                                                                                    >  3:288505/1‑149 (MQ=255)
                                        gaagaaAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCGCCCGTCGAGGGAGCGGCgaat                                                                                                                                     >  2:239378/1‑113 (MQ=38)
                                        gaagaaAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCGCCCGTCGAGGGAGCGGCgaat                                                                                                                                     <  1:239378/114‑2 (MQ=38)
                                           gaaAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGc                                                                                               >  6:165552/1‑149 (MQ=255)
                                                 cTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGt                                                                                         >  6:327636/1‑149 (MQ=255)
                                                    gATGTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGGGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCTCCACCGGTTGCTCgg                                                                                           <  3:174779/144‑1 (MQ=255)
                                                    caTCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCCCCGCTCACTCGACGGGCGAGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTAGCTCgggagaa                                                                                      >  4:174779/2‑145 (MQ=255)
                                                     aTCTAAGTCGCCTGTTTCAAGGTAAGGCGTTCTGTTGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTg                                                                                      <  4:372943/148‑1 (MQ=255)
                                                     aTCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGAAAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTgg                                                                                     >  7:125724/1‑149 (MQ=255)
                                                      tCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGc                                                                                    >  5:113219/1‑149 (MQ=255)
                                                          aGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGc                                                                                    <  4:37407/145‑1 (MQ=255)
                                                          aGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTa                                                                                 <  6:76065/148‑1 (MQ=255)
                                                          aGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTAc                                                                                <  2:353731/149‑1 (MQ=255)
                                                            tCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGt                                                                              >  1:46232/1‑149 (MQ=255)
                                                               ccTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCt                                                                           <  5:288273/149‑1 (MQ=255)
                                                                 tGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTg                                                                          >  8:341140/1‑148 (MQ=255)
                                                                 tGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTg                                                                          >  2:327412/1‑148 (MQ=255)
                                                                   ttGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGaa                                                                       >  4:200123/1‑149 (MQ=255)
                                                                     gCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAAtt                                                                     >  7:18821/1‑149 (MQ=255)
                                                                      cAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATtc                                                                    <  7:196673/149‑1 (MQ=255)
                                                                      cAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATtc                                                                    >  3:112204/1‑149 (MQ=255)
                                                                         ggTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATtcgtc                                                                 <  7:341140/149‑1 (MQ=255)
                                                                          gTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATtcgtc                                                                 >  6:90676/1‑148 (MQ=255)
                                                                           tAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTcga                                                               >  5:33100/1‑149 (MQ=255)
                                                                                 gTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACt                                                         >  7:168018/1‑149 (MQ=255)
                                                                                  ttCTGATGGCTATCAGCATGACAGAAGCCGCCGCGCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACt                                                         >  1:247201/1‑149 (MQ=255)
                                                                                     tGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCg                                                     >  4:181407/1‑149 (MQ=255)
                                                                                     tGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCg                                                     >  4:175468/1‑149 (MQ=255)
                                                                                       aTGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGc                                                      <  8:168018/146‑1 (MQ=255)
                                                                                       aTGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGa                                                   >  8:244869/1‑149 (MQ=255)
                                                                                             aTCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCgg                                                    <  8:140114/142‑1 (MQ=255)
                                                                                             aTCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCgg                                                    >  7:140114/1‑142 (MQ=255)
                                                                                              tCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGAc                                            <  3:200123/149‑1 (MQ=255)
                                                                                                 gCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAg                                         <  2:46232/149‑1 (MQ=255)
                                                                                                  cATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAgg                                        >  6:32334/1‑149 (MQ=255)
                                                                                                   aTGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGgt                                       >  5:342015/1‑149 (MQ=255)
                                                                                                    tGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGgtg                                      <  1:327412/149‑1 (MQ=255)
                                                                                                       cAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTc                                   >  4:151504/1‑149 (MQ=255)
                                                                                                                   ccAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGa                                             >  3:230421/1‑127 (MQ=255)
                                                                                                                   ccAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGa                                             <  4:230421/127‑1 (MQ=255)
                                                                                                                      aaCCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGc                                                                                    <  7:49287/86‑1 (MQ=255)
                                                                                                                      aaCCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGc                                                                                    >  8:49287/1‑86 (MQ=255)
                                                                                                                       ctcaccTTTGCCGCGCCCGTCGAGGGAGCGGCGAATGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGt                     >  3:274118/3‑148 (MQ=37)
                                                                                                                        tcaccTTTGCCGCGCCCGTCGAGGGAGCGGCGAATGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTc                    >  7:111572/2‑148 (MQ=38)
                                                                                                                         cacaTTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTc                    >  1:295022/1‑147 (MQ=255)
                                                                                                                           caTTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAgg                <  6:191883/149‑1 (MQ=255)
                                                                                                                           caTTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAgg                <  5:90676/149‑1 (MQ=255)
                                                                                                                               cgccgcTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGa             <  6:113219/148‑1 (MQ=255)
                                                                                                                                gccgcTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGat            <  8:125724/148‑1 (MQ=255)
                                                                                                                                gccgcTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGat            >  6:314165/1‑148 (MQ=255)
                                                                                                                                   gcTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGatcatc        <  5:314165/149‑1 (MQ=255)
                                                                                                                                    cTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCCGGTGTTCGAAAACCATGGCGTCAAGGTGATCATCg       >  8:283206/1‑149 (MQ=255)
                                                                                                                                       ccTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGATCATCGATc    >  4:15932/1‑149 (MQ=255)
                                                                                                                                        cTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGATCATCGATcc   >  4:21867/1‑149 (MQ=255)
                                                                                                                                          tCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGATCATCGATcc   >  5:380766/1‑148 (MQ=255)
                                                                                                                                            gACGGGCGCGTCAAAGGTGAGGGCATTCTCCTGTGCGTGCTCACCACCGGTTGCTCTGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGATCATCGATccc  <  3:281800/147‑1 (MQ=255)
                                                                                                                                            gACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGATCATCGATccc  <  5:165552/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
CCGAAGATGCCATCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGATCATCGATCCC  >  NC_002947/980140‑980424

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: