Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 471,714 G→A 12.3% I408I (ATC→ATT rpoD ← RNA polymerase sigma 70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947471,7140GA12.3% 518.7 / 51.7 195I408I (ATC→ATTrpoDRNA polymerase sigma 70 factor
Reads supporting (aligned to +/- strand):  ref base G (91/80);  new base A (11/13);  total (102/93)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 197 total aligned reads displayed.

TGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCG  >  NC_002947/471565‑471861
                                                                                                                                                     |                                                                                                                                                   
tgCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTgg                                                                                                                                                     >  1:389662/1‑150 (MQ=255)
  cACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGAt                                                                                                                                                   >  2:4430935/1‑150 (MQ=255)
    ccGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCa                                                                                                                                                 >  2:2519204/1‑150 (MQ=255)
       gAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGAt                                                                                                                                                   <  2:2635232/145‑1 (MQ=255)
        aaTGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTc                                                                                                                                             <  2:4324989/150‑1 (MQ=255)
         atacGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCg                                                                                                                                            >  1:3240286/4‑150 (MQ=255)
          tacGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGa                                                                                                                                           <  2:2120921/148‑1 (MQ=255)
          tGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGa                                                                                                                                           <  1:5279963/150‑1 (MQ=255)
          tGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGa                                                                                                                                           <  2:1080631/150‑1 (MQ=255)
          tGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGa                                                                                                                                           <  1:4029793/150‑1 (MQ=255)
           gCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAg                                                                                                                                          <  2:3105392/150‑1 (MQ=255)
             ggATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGa                                                                                                                                        >  2:1775890/1‑150 (MQ=255)
              gATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGaa                                                                                                                                       <  2:942900/150‑1 (MQ=255)
               aTGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAAc                                                                                                                                      >  1:4708495/1‑150 (MQ=255)
               aTGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAAc                                                                                                                                      <  2:3507806/150‑1 (MQ=255)
                             gTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGtt                                                                                                                        >  1:1098123/1‑150 (MQ=255)
                                ggCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTTGCCGCGTCGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGAGGTTGCCTTCCTGGATCAGGTCGTGGACCTGCCGGCCTCGGTTGgt                                                                                                                     >  2:1414089/1‑150 (MQ=255)
                                ggCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGgt                                                                                                                     >  1:727540/1‑150 (MQ=255)
                                ggCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGgt                                                                                                                     >  2:5216818/1‑150 (MQ=255)
                                ggCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGgt                                                                                                                     >  1:3769656/1‑150 (MQ=255)
                                   gATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTa                                                                                                                  >  2:5169306/1‑150 (MQ=255)
                                   gATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTa                                                                                                                  <  1:3467875/150‑1 (MQ=255)
                                        aaCGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTActtct                                                                                                             >  2:2439938/1‑150 (MQ=255)
                                         aCGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTActtctt                                                                                                            <  1:4521580/150‑1 (MQ=255)
                                          cGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTg                                                                                                           <  1:2537292/150‑1 (MQ=255)
                                            ggTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGc                                                                                                         >  2:1926148/1‑150 (MQ=255)
                                             gTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGca                                                                                                        >  1:4580067/1‑149 (MQ=255)
                                               gATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGAt                                                                                                      >  1:55200/1‑150 (MQ=255)
                                                aTCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATg                                                                                                     <  2:2335006/150‑1 (MQ=255)
                                                   gCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGag                                                                                                  <  2:920114/150‑2 (MQ=255)
                                                    ccTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGaaa                                                                                                 >  2:4247056/1‑150 (MQ=255)
                                                    ccTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGaga                                                                                                 <  1:2226926/150‑3 (MQ=255)
                                                    ccTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGaga                                                                                                 <  1:1433117/150‑3 (MQ=255)
                                                     cTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAAt                                                                                                <  1:2325116/150‑1 (MQ=255)
                                                       gACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCa                                                                                              >  1:4739836/1‑150 (MQ=255)
                                                       gACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGAtaa                                                                                              >  1:289984/1‑148 (MQ=255)
                                                       gACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGAtaa                                                                                              <  1:2661586/150‑3 (MQ=255)
                                                          ggATCCACCTGGTTGAATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCa                                                                                           <  1:1884146/150‑1 (MQ=255)
                                                           gATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAg                                                                                          >  1:468180/1‑150 (MQ=255)
                                                           gATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAg                                                                                          >  2:529977/1‑150 (MQ=255)
                                                               caccaGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGc                                                                                      <  2:392192/150‑1 (MQ=255)
                                                                accaTGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACgaa                                                                                     >  2:775847/1‑148 (MQ=255)
                                                                accaGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCa                                                                                     >  2:3580962/1‑150 (MQ=255)
                                                                 ccaGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAg                                                                                    <  2:2183214/150‑1 (MQ=255)
                                                                   aGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTATTTCTTGGCGATGGAAATCACCAGACGCAGGt                                                                                  >  2:1285654/1‑150 (MQ=255)
                                                                     gTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTc                                                                                >  2:4080081/1‑150 (MQ=255)
                                                                     gTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTc                                                                                >  2:1328735/1‑150 (MQ=255)
                                                                         cATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCt                                                                            <  1:4825056/150‑1 (MQ=255)
                                                                         cATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCt                                                                            <  2:2850456/150‑1 (MQ=255)
                                                                            aGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTcca                                                                         >  1:1625036/1‑148 (MQ=255)
                                                                             ggTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGAc                                                                        <  1:3598059/150‑1 (MQ=255)
                                                                               tCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGtt                                                                                 <  1:2292477/139‑1 (MQ=255)
                                                                               tCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGtt                                                                                 >  2:2292477/1‑139 (MQ=255)
                                                                               tCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCa                                                                      >  2:298188/1‑150 (MQ=255)
                                                                                 gagaACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAtt                                                                    >  2:3201780/1‑150 (MQ=255)
                                                                                 gagaACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAtt                                                                    <  2:189347/150‑1 (MQ=255)
                                                                                  agaACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAttt                                                                   >  2:3872103/1‑150 (MQ=255)
                                                                                  agaACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAttt                                                                   <  1:4129258/150‑1 (MQ=255)
                                                                                         gTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTtctttctt                                                            <  2:324741/150‑1 (MQ=255)
                                                                                         gTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTg                                                                                                           <  2:4713431/103‑1 (MQ=255)
                                                                                            gCCGCGACGGTATTCGACCGTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGgc                                                         <  2:5049857/150‑1 (MQ=255)
                                                                                            gCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGgc                                                         <  1:2044485/150‑1 (MQ=255)
                                                                                            gCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGgc                                                         <  1:1221674/150‑1 (MQ=255)
                                                                                            gCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGgc                                                         <  1:5200110/150‑1 (MQ=255)
                                                                                                   cGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGcggcgg                                                  <  2:2620143/150‑1 (MQ=255)
                                                                                                    ggTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCggg                                                 >  2:1720034/1‑150 (MQ=255)
                                                                                                       aTTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCt                                              >  2:953899/1‑150 (MQ=255)
                                                                                                       aTTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGACTTCCTGGATCAGGTCGAGGAACTGCAGGCCa                                                                                                                             <  1:1332995/71‑1 (MQ=255)
                                                                                                        ttCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAg                                                                                                                                 >  2:4652760/1‑66 (MQ=255)
                                                                                                        ttCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAg                                                                                                                                 <  1:4652760/66‑1 (MQ=255)
                                                                                                         tCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTg                                            >  1:4404624/1‑150 (MQ=255)
                                                                                                           gAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGc                                          <  1:3393556/150‑1 (MQ=255)
                                                                                                                 gTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCAcc                                    >  2:2662634/1‑150 (MQ=255)
                                                                                                                 gTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCAcc                                    >  1:3846212/1‑150 (MQ=255)
                                                                                                                  tCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCg                                   >  1:757023/1‑150 (MQ=255)
                                                                                                                     aCCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCa                                                                                              <  1:3860942/88‑1 (MQ=255)
                                                                                                                     aCCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCa                                                                                              >  2:3860942/1‑88 (MQ=255)
                                                                                                                     aCCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCt                                              <  2:624953/136‑1 (MQ=255)
                                                                                                                     aCCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATc                                <  2:1474087/150‑1 (MQ=255)
                                                                                                                     aCCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATc                                >  1:4934915/1‑150 (MQ=255)
                                                                                                                     aCCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATc                                >  2:2240723/1‑150 (MQ=255)
                                                                                                                     aCCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATc                                <  2:1546532/150‑1 (MQ=255)
                                                                                                                          tttCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACAt                           <  2:4551397/150‑1 (MQ=255)
                                                                                                                              atcaGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGa                       >  1:1411048/1‑150 (MQ=255)
                                                                                                                              atcaGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGa                       <  1:1962717/150‑1 (MQ=255)
                                                                                                                               tcaGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGAc                      >  2:3068575/1‑150 (MQ=255)
                                                                                                                               tcaGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAGTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGAc                      >  2:3946557/1‑150 (MQ=255)
                                                                                                                                 aGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCgcggc                                                    >  2:1139924/1‑118 (MQ=255)
                                                                                                                                 aGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCgcggc                                                    <  1:1139924/118‑1 (MQ=255)
                                                                                                                                 aGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACgg                    <  2:3468932/150‑1 (MQ=255)
                                                                                                                                 aGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACgg                    >  1:449318/1‑150 (MQ=255)
                                                                                                                                 aGACCGATGTTGCCTTCATGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACgg                    >  2:4154694/1‑150 (MQ=255)
                                                                                                                                   aCCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGtt                  <  1:1742702/150‑1 (MQ=255)
                                                                                                                                     cGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGa                <  2:2687245/150‑1 (MQ=255)
                                                                                                                                           tGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCt          >  2:3452386/1‑150 (MQ=255)
                                                                                                                                               ttCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGAt      <  2:555078/150‑1 (MQ=255)
                                                                                                                                               ttCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGAt      >  1:2437829/1‑150 (MQ=255)
                                                                                                                                                   tGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCg  >  1:955421/1‑150 (MQ=255)
                                                                                                                                                   tGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCg  <  1:5169306/150‑1 (MQ=255)
                                                                                                                                                   tGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCg  >  2:2684002/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCG  >  NC_002947/471565‑471861

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: