Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 471,678 G→A 20.9% D420D (GAC→GAT rpoD ← RNA polymerase sigma 70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947471,6780GA20.9% 711.1 / 181.5 344D420D (GAC→GATrpoDRNA polymerase sigma 70 factor
Reads supporting (aligned to +/- strand):  ref base G (119/153);  new base A (35/37);  total (154/190)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 344 total aligned reads displayed.

AATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTC  >  NC_002947/471531‑471824
                                                                                                                                                   |                                                                                                                                                  
      gTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGc                                                                                                                                            >  2:6599368/1‑150 (MQ=255)
      aacgatcTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGc                                                                                                                                            >  2:5988848/7‑150 (MQ=255)
       ttGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCt                                                                                                                                           >  1:5704728/1‑150 (MQ=255)
         gAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGATATTCGAACTTGTCCACCGCttt                                                                                                                                         <  1:8360460/150‑1 (MQ=255)
           tcTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTtca                                                                                                                                       <  2:3934412/149‑1 (MQ=255)
           tcTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTtca                                                                                                                                       <  1:2598221/149‑1 (MQ=255)
             ttgttgATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTtcatc                                                                                                                                     <  1:1645725/150‑1 (MQ=255)
             ttgttgATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTtcatc                                                                                                                                     <  1:645605/150‑1 (MQ=255)
              tgttgATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTtcatca                                                                                                                                    <  1:8231085/150‑1 (MQ=255)
               gttgATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAg                                                                                                                                   >  1:3158325/1‑150 (MQ=255)
                ttgATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGa                                                                                                                                  <  1:7361298/150‑1 (MQ=255)
                 tgATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGGACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGAc                                                                                                                                 <  2:2497174/150‑1 (MQ=255)
                  gATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGAcc                                                                                                                                <  1:3475457/150‑1 (MQ=255)
                       tctcGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGt                                                                                                                           >  1:6321182/1‑150 (MQ=255)
                       tctcGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGt                                                                                                                           >  1:6066435/1‑150 (MQ=255)
                        ctcGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGtt                                                                                                                          <  1:2148403/150‑1 (MQ=255)
                        ctcGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGtt                                                                                                                          <  2:5465401/150‑1 (MQ=255)
                           gATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGcc                                                                                                                       >  1:6993948/1‑150 (MQ=255)
                           gATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGcc                                                                                                                       >  2:7897970/1‑150 (MQ=255)
                           gATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGcc                                                                                                                       >  1:5184727/1‑150 (MQ=255)
                           gATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGcc                                                                                                                       >  2:7406015/1‑150 (MQ=255)
                            aTCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCt                                                                                                                      >  1:2226066/1‑150 (MQ=255)
                             tCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCtt                                                                                                                     >  2:8689243/1‑150 (MQ=255)
                               aTGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTcc                                                                                                                   <  1:669493/150‑1 (MQ=255)
                                tgtgCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCt                                                                                                                  <  2:1426131/150‑1 (MQ=255)
                                tgtgCACCGGAATGCGGACGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCt                                                                                                                  <  2:6207121/150‑1 (MQ=255)
                                tgtgCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCt                                                                                                                  <  2:2664455/150‑1 (MQ=255)
                                tgtgCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCt                                                                                                                  <  1:4240555/150‑1 (MQ=255)
                                tgtgCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCt                                                                                                                  <  2:5349137/150‑1 (MQ=255)
                                 gtgCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTg                                                                                                                 <  1:3722815/150‑1 (MQ=255)
                                 gtgCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTg                                                                                                                 >  2:3320847/1‑150 (MQ=255)
                                  tgCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTgg                                                                                                                >  2:1116964/1‑150 (MQ=255)
                                     aCCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATc                                                                                                             >  2:7573972/1‑150 (MQ=255)
                                      ccGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAAtaa                                                                                                            >  2:8380351/1‑148 (MQ=255)
                                       cGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAg                                                                                                           <  2:8112404/150‑1 (MQ=255)
                                            tacGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGa                                                                                                      <  2:1371382/148‑1 (MQ=255)
                                            tGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGa                                                                                                      >  2:6706862/1‑150 (MQ=255)
                                               ggATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGa                                                                                                   >  2:2718739/1‑150 (MQ=255)
                                                gATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGaa                                                                                                  <  1:7590226/150‑1 (MQ=255)
                                                 aTGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAAc                                                                                                 <  2:1569315/150‑1 (MQ=255)
                                                     tgcgtgcCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCa                                                                                             <  1:2462779/150‑1 (MQ=255)
                                                      gcgtgcCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAg                                                                                            <  2:3643321/150‑1 (MQ=255)
                                                       cgtgcCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAgg                                                                                           <  2:5604441/150‑1 (MQ=255)
                                                            cTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACg                                                                                      <  1:5934154/150‑1 (MQ=255)
                                                             tGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACgg                                                                                     <  2:3049106/150‑1 (MQ=255)
                                                                 cGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTgg                                                                                 <  1:1410139/150‑1 (MQ=255)
                                                                  ggCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGgt                                                                                >  1:7982625/1‑150 (MQ=255)
                                                                  ggCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGgt                                                                                <  1:7850945/150‑1 (MQ=255)
                                                                   gCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGgtg                                                                               >  2:203939/1‑150 (MQ=255)
                                                                     gATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTa                                                                             >  1:2161997/1‑150 (MQ=255)
                                                                       tCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTAct                                                                           >  2:2418550/1‑150 (MQ=255)
                                                                          aaCGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTActtct                                                                        <  2:3747449/150‑1 (MQ=255)
                                                                          aaCGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTActtct                                                                        <  2:7708993/150‑1 (MQ=255)
                                                                            cGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTg                                                                      <  1:8223906/150‑1 (MQ=255)
                                                                             gggTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTgg                                                                     >  2:2793728/1‑150 (MQ=255)
                                                                             gggTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTgg                                                                     >  2:6168247/1‑150 (MQ=255)
                                                                             gggTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTgg                                                                     >  1:1045889/1‑150 (MQ=255)
                                                                               gTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCg                                                                   >  2:8567822/1‑150 (MQ=255)
                                                                                  aTCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATg                                                                <  1:950022/150‑1 (MQ=255)
                                                                                   tCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATgg                                                               <  2:1932771/150‑1 (MQ=255)
                                                                                   tCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATgg                                                               <  1:3572355/150‑1 (MQ=255)
                                                                                      ccTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTc                                                                                                                    <  1:8969103/94‑1 (MQ=255)
                                                                                      ccTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTc                                                                                                                    >  2:8969103/1‑94 (MQ=255)
                                                                                        tGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATc                                                          >  2:4070373/1‑150 (MQ=255)
                                                                                         gACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGAtaa                                                         >  2:6827426/1‑148 (MQ=255)
                                                                                          aCGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCAc                                                        <  2:6976816/150‑1 (MQ=255)
                                                                                          aCGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAAc                                                        <  1:4814657/150‑1 (MQ=255)
                                                                                            gaaTCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAAcga                                                      >  1:6575517/3‑148 (MQ=255)
                                                                                            gaaTCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAAcga                                                      <  1:1713724/148‑3 (MQ=255)
                                                                                               tctaccaGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGAc                                                   <  1:4500657/147‑1 (MQ=255)
                                                                                                  accaGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCa                                                >  1:6442923/1‑150 (MQ=255)
                                                                                                  accaGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCa                                                >  2:4675780/1‑150 (MQ=255)
                                                                                                  accaGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCa                                                >  2:1785884/1‑150 (MQ=255)
                                                                                                  accaGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCa                                                >  1:8152833/1‑150 (MQ=255)
                                                                                                   ccaGGTGGCAAAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGAc                                                                                                                                 <  2:1975493/68‑1 (MQ=255)
                                                                                                   ccaGGTGGCAAAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGAc                                                                                                                                 >  1:1975493/1‑68 (MQ=255)
                                                                                                       gTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTc                                           <  1:4136863/150‑1 (MQ=255)
                                                                                                       gTGGCATAGGTAGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTc                                           <  2:5760773/150‑1 (MQ=255)
                                                                                                         ggCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGc                                         <  1:1159698/150‑1 (MQ=255)
                                                                                                            aTAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTc                                      >  1:3251921/1‑150 (MQ=255)
                                                                                                            aTAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTc                                      >  2:7186319/1‑150 (MQ=255)
                                                                                                              aGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGa                                    >  2:7621020/1‑150 (MQ=255)
                                                                                                                gTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGAcc                                  >  2:4848455/1‑150 (MQ=255)
                                                                                                                gTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGAcc                                  >  1:8475640/1‑150 (MQ=255)
                                                                                                                gTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCAcc                                  >  1:1820864/1‑150 (MQ=255)
                                                                                                                  cGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCACCAt                                >  2:4546660/1‑150 (MQ=255)
                                                                                                                   gagaACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAtt                               >  2:2272810/1‑150 (MQ=255)
                                                                                                                   gagaACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAtt                               >  1:5375230/1‑150 (MQ=255)
                                                                                                                        cTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAtttcttt                          >  1:8586511/1‑150 (MQ=255)
                                                                                                                           gTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTtctttctt                       <  2:1244617/150‑1 (MQ=255)
                                                                                                                           gTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTtctttctt                       <  2:511382/150‑1 (MQ=255)
                                                                                                                           gTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTtctttctt                       <  1:2341752/150‑1 (MQ=255)
                                                                                                                           gTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTtctttctt                       <  1:1751585/150‑1 (MQ=255)
                                                                                                                                    aCGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCACCATTTCCTTCTTGgcgatggt              <  1:3880965/150‑6 (MQ=255)
                                                                                                                                     cGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGcggcgg             <  2:2226066/150‑1 (MQ=255)
                                                                                                                                          ttCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCtt        >  1:3403087/1‑150 (MQ=255)
                                                                                                                                           tCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTg       >  2:9061074/1‑150 (MQ=255)
                                                                                                                                           tCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCCCGCCGCGccaag       <  2:6321182/150‑4 (MQ=255)
                                                                                                                                            cGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTgg      <  2:1045889/150‑1 (MQ=255)
                                                                                                                                             gAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGc     <  2:7982625/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTC  >  NC_002947/471531‑471824

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