Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 471,687 T→C 22.3% K417K (AAA→AAG rpoD ← RNA polymerase sigma 70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947471,6870TC22.3% 355.3 / 98.0 179K417K (AAA→AAGrpoDRNA polymerase sigma 70 factor
Reads supporting (aligned to +/- strand):  ref base T (80/59);  new base C (23/17);  total (103/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

Only 100 of 182 total aligned reads displayed.

AGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACAT  >  NC_002947/471541‑471836
                                                                                                                                                  |                                                                                                                                                     
aGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTtc                                                                                                                                                    <  1:3733230/150‑1 (MQ=255)
   ttgttgATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTtcatc                                                                                                                                                 >  2:3099487/1‑150 (MQ=255)
     gttgATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAg                                                                                                                                               >  2:2467625/1‑150 (MQ=255)
       tgATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGAc                                                                                                                                             <  1:3671756/150‑1 (MQ=255)
       tgATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGAc                                                                                                                                             <  2:4317414/150‑1 (MQ=255)
        gATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGAcc                                                                                                                                            >  2:4999756/1‑150 (MQ=255)
           cGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGAt                                                                                                                                         <  1:2995531/150‑1 (MQ=255)
           cGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGAt                                                                                                                                         <  1:2500638/150‑1 (MQ=255)
             tctcGATCATGTGCACTGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGt                                                                                                                                       >  1:4458993/1‑150 (MQ=255)
             tctcGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGt                                                                                                                                       >  2:3514725/1‑150 (MQ=255)
               tcGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTg                                                                                                                                     >  1:2332366/1‑150 (MQ=255)
                 gatgatGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGcc                                                                                                                                   >  1:2905263/1‑150 (MQ=255)
                  aTCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCt                                                                                                                                  <  1:1514464/150‑1 (MQ=255)
                   tCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCtt                                                                                                                                 >  1:551081/1‑150 (MQ=255)
                   tCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCtt                                                                                                                                 >  1:2160513/1‑150 (MQ=255)
                   tCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCAACCGCTTTCATCAGACCGATGTTGCCtt                                                                                                                                 >  1:1233412/1‑150 (MQ=255)
                      tgtgCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCt                                                                                                                              >  2:2595459/1‑150 (MQ=255)
                      tgtgCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCt                                                                                                                              <  1:489692/150‑1 (MQ=255)
                      tgtgCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCt                                                                                                                              <  2:3847685/150‑1 (MQ=255)
                        tgCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTgg                                                                                                                            >  1:2150717/1‑150 (MQ=255)
                          cACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAAt                                                                                                                          <  1:4303089/150‑1 (MQ=255)
                           aCCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGACTTCCTGGATc                                                                                                                         >  1:1343175/1‑150 (MQ=255)
                             cGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAAGTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAg                                                                                                                       >  2:37302/1‑150 (MQ=255)
                             cGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAg                                                                                                                       >  2:1647313/1‑150 (MQ=255)
                               gAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAggn                                                                                                                     <  1:2952862/150‑2 (MQ=255)
                               gAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGt                                                                                                                     <  2:1059115/150‑1 (MQ=255)
                                  tGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGa                                                                                                                  <  2:2177654/150‑1 (MQ=255)
                                    cGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAgg                                                                                                                >  1:3349699/1‑150 (MQ=255)
                                     ggATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGTGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGa                                                                                                               >  2:1215201/1‑150 (MQ=255)
                                     ggATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGa                                                                                                               >  2:4129842/1‑150 (MQ=255)
                                     ggATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGa                                                                                                               <  2:2420982/150‑1 (MQ=255)
                                          gtgcgtgcCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGc                                                                                                          <  2:1351746/150‑1 (MQ=255)
                                           tgcgtgcCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCa                                                                                                         <  1:823075/150‑1 (MQ=255)
                                                   tGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACgg                                                                                                 <  2:3714058/150‑1 (MQ=255)
                                                    ggTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGt                                                                                                <  1:4336995/150‑1 (MQ=255)
                                                      tCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTg                                                                                              <  2:4950682/150‑1 (MQ=255)
                                                         gcaaTCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGgtg                                                                                           >  2:1442884/4‑150 (MQ=255)
                                                           gATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTa                                                                                         <  2:4537645/150‑1 (MQ=255)
                                                           gATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTa                                                                                         >  1:2662900/1‑150 (MQ=255)
                                                                aaCGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTActtct                                                                                    >  1:4606971/1‑150 (MQ=255)
                                                                 aCGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGATAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTActtctt                                                                                   <  2:4726154/150‑1 (MQ=255)
                                                                 aCGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTActtctt                                                                                   >  1:1898858/1‑150 (MQ=255)
                                                                 aCGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTActtctt                                                                                   >  1:1899101/1‑150 (MQ=255)
                                                                   gggTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTAcagatcgg                                                                                 >  1:4707737/1‑143 (MQ=255)
                                                                   gggTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTgg                                                                                 <  1:178222/150‑1 (MQ=255)
                                                                      tGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGa                                                                              <  1:4100305/150‑1 (MQ=255)
                                                                       gATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGAt                                                                             <  2:5103271/150‑1 (MQ=255)
                                                                       gATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGAt                                                                             >  1:1481083/1‑150 (MQ=255)
                                                                         tCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATgg                                                                           >  2:605239/1‑150 (MQ=255)
                                                                           gCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGaa                                                                         <  1:2871128/150‑1 (MQ=255)
                                                                             cTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGAt                                                                       >  2:643332/1‑150 (MQ=255)
                                                                               gACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATTGAGAtaa                                                                     >  2:858280/1‑148 (MQ=255)
                                                                                aCGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCAc                                                                    <  2:1693023/150‑1 (MQ=255)
                                                                                  ggATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCa                                                                  >  1:2447889/1‑150 (MQ=255)
                                                                                  ggATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCa                                                                  >  2:1877912/1‑150 (MQ=255)
                                                                                   gATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAg                                                                 >  2:401614/1‑150 (MQ=255)
                                                                                     tCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGAc                                                               <  2:4566644/150‑1 (MQ=255)
                                                                                     tCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGAc                                                               <  2:381974/150‑1 (MQ=255)
                                                                                      ccaccaTGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACg                                                              >  2:4448516/1‑150 (MQ=255)
                                                                                         ccaTGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAg                                                           >  1:1359548/1‑150 (MQ=255)
                                                                                          caGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAgg                                                          <  1:2101830/150‑1 (MQ=255)
                                                                                                  aTAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTc                                                  >  1:1044320/1‑150 (MQ=255)
                                                                                                    aGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGa                                                <  2:787751/150‑1 (MQ=255)
                                                                                                     ggTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGAc                                               >  1:4348688/1‑150 (MQ=255)
                                                                                                      gTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTACATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCAcc                                              <  1:558717/150‑1 (MQ=255)
                                                                                                         gagaACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAtt                                           >  1:4602312/1‑150 (MQ=255)
                                                                                                         gagaACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAtt                                           >  1:3332406/1‑150 (MQ=255)
                                                                                                         gagaACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCACCAtt                                           <  1:36940/150‑1 (MQ=255)
                                                                                                          agaACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAttt                                          <  2:2024891/150‑1 (MQ=255)
                                                                                                          agaACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAttt                                          >  1:1166122/1‑150 (MQ=255)
                                                                                                            aaCTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTAGCCTCGACCAtttct                                        <  1:2457623/150‑1 (MQ=255)
                                                                                                             aCTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAtttctt                                       <  1:3620737/150‑1 (MQ=255)
                                                                                                             aCTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAtttctt                                       >  2:2595812/1‑150 (MQ=255)
                                                                                                              cTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAtttcttt                                      <  1:2595459/150‑1 (MQ=255)
                                                                                                                 gTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTtctttctt                                   <  2:3007675/150‑1 (MQ=255)
                                                                                                                 gTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTtctttctt                                   >  1:385023/1‑150 (MQ=255)
                                                                                                                     ccGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGgcg                               >  1:3228019/1‑150 (MQ=255)
                                                                                                                      cacgacgaTATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCACCATTTCCTTCTTGgcga                              >  1:1739489/3‑149 (MQ=255)
                                                                                                                       acgacgaTATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCACCATTTCCTTCTTGgcgat                             >  1:513995/2‑148 (MQ=255)
                                                                                                                         gACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCgcggc                           <  1:1998942/150‑1 (MQ=255)
                                                                                                                         gACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCACCATTTCCTTCTTGgcgatgg                           >  1:3979959/1‑146 (MQ=255)
                                                                                                                          aCGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCgcggcg                          >  1:3090263/1‑150 (MQ=255)
                                                                                                                              tATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGcc                      <  2:143161/150‑1 (MQ=255)
                                                                                                                                ttCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCtt                    <  1:2880364/150‑1 (MQ=255)
                                                                                                                                ttCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCtt                    >  1:1978682/1‑150 (MQ=255)
                                                                                                                                ttCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCtt                    >  1:3409258/1‑150 (MQ=255)
                                                                                                                                 tCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTg                   >  2:4351170/1‑150 (MQ=255)
                                                                                                                                   gAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGtt                                                                                               >  1:2263055/1‑72 (MQ=255)
                                                                                                                                   gAACTTGTCCACCGCTGTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGc                 >  2:2181473/1‑150 (MQ=255)
                                                                                                                                    aaCTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCt                >  2:1765385/1‑150 (MQ=255)
                                                                                                                                    aaCTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCt                >  2:3825144/1‑150 (MQ=255)
                                                                                                                                       ttGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCa             >  2:513965/1‑150 (MQ=255)
                                                                                                                                         gTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCAcc           >  2:4116452/1‑150 (MQ=255)
                                                                                                                                            cACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGAt        >  1:1014137/1‑150 (MQ=255)
                                                                                                                                             aCCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATc       >  2:5046623/1‑150 (MQ=255)
                                                                                                                                             aCCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATc       >  1:4574120/1‑150 (MQ=255)
                                                                                                                                             aCCGCTTTCATCAGACCGATGGTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATc       <  1:2707349/150‑1 (MQ=255)
                                                                                                                                              ccGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCg      >  1:5028225/1‑150 (MQ=255)
                                                                                                                                              ccGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAAcga                                                                  >  1:128071/1‑88 (MQ=255)
                                                                                                                                                  tttCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACAt  >  1:3210524/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACAT  >  NC_002947/471541‑471836

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: