Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 471,759 G→A 14.7% I393I (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,7590GA14.7% 422.0 / 55.7 170I393I (ATC→ATTrpoDRNA polymerase sigma 70 factor
Reads supporting (aligned to +/- strand):  ref base G (84/61);  new base A (16/9);  total (100/70)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 174 total aligned reads displayed.

GTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTG  >  NC_002947/471610‑471908
                                                                                                                                                     |                                                                                                                                                     
gTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCg                                                                                                                                                       >  1:2914260/1‑150 (MQ=255)
  gATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGAt                                                                                                                                                     <  1:2307150/150‑1 (MQ=255)
      gCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGaa                                                                                                                                                 >  1:809878/1‑150 (MQ=255)
        cTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGAt                                                                                                                                               >  2:643332/1‑150 (MQ=255)
             ggATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCa                                                                                                                                          >  1:2447891/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)
                   accaTGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACgaa                                                                                                                                    >  2:5025047/1‑148 (MQ=255)
                   accaGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCa                                                                                                                                    >  2:1377953/1‑150 (MQ=255)
                    ccaTGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAg                                                                                                                                   >  1:1359548/1‑150 (MQ=255)
                     caGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAgg                                                                                                                                  <  1:2101830/150‑1 (MQ=255)
                       tgtcgCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGtt                                                                                                                                <  1:4936859/146‑1 (MQ=255)
                             aTAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTc                                                                                                                          >  1:1044320/1‑150 (MQ=255)
                               aGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTcca                                                                                                                        >  2:8517/1‑148 (MQ=255)
                                ggTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGAc                                                                                                                       >  1:4348688/1‑150 (MQ=255)
                                 gTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTACATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCAcc                                                                                                                      <  1:558717/150‑1 (MQ=255)
                                  tCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCa                                                                                                                     >  1:3397010/1‑150 (MQ=255)
                                  tCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTccacca                                                                                                                     >  1:1717221/1‑150 (MQ=255)
                                    gagaACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAtt                                                                                                                   >  2:1152978/1‑150 (MQ=255)
                                    gagaACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAtt                                                                                                                   >  2:1741335/1‑150 (MQ=255)
                                    gagaACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAtt                                                                                                                   >  1:3332406/1‑150 (MQ=255)
                                    gagaACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCACCAtt                                                                                                                   <  1:36940/150‑1 (MQ=255)
                                     agaACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAttt                                                                                                                  <  1:3410817/150‑1 (MQ=255)
                                     agaACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAttt                                                                                                                  <  2:2024891/150‑1 (MQ=255)
                                       aaCTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCACCAtttcc                                                                                                                >  2:617735/1‑149 (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)
                                                 cacgacgaTATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCACCATTTCCTTCTTGgcga                                                                                                      <  2:586251/148‑2 (MQ=255)
                                                 cacgacgaTATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCACCATTTCCTTCTTGgcga                                                                                                      >  2:694541/3‑149 (MQ=255)
                                                  gcgACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGcgcg                                                                                                     >  1:268462/1‑150 (MQ=255)
                                                  gcgACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGcgcg                                                                                                     >  2:1498789/1‑150 (MQ=255)
                                                  gcgACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCGTGGcgcg                                                                                                     <  2:4348688/150‑1 (MQ=255)
                                                  acgacgaTATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCACCATTTCCTTCTTGgcgat                                                                                                     >  1:513995/2‑148 (MQ=255)
                                                    gACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCgcggc                                                                                                   <  1:1998942/150‑1 (MQ=255)
                                                    gACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCACCATTTCCTTCTTGgcgatgg                                                                                                   >  1:3979959/1‑146 (MQ=255)
                                                    gACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCACCATTTCCTTCTTGgcgatgg                                                                                                   <  1:4252037/150‑5 (MQ=255)
                                                     aCGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCgcggcg                                                                                                  >  1:3090263/1‑150 (MQ=255)
                                                      cGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGcggcgg                                                                                                 <  2:4628271/150‑1 (MQ=255)
                                                         tATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGcc                                                                                              <  2:143161/150‑1 (MQ=255)
                                                         tATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGcc                                                                                              >  2:1847287/1‑150 (MQ=255)
                                                           ttCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCtt                                                                                            >  1:1978682/1‑150 (MQ=255)
                                                           ttCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCtt                                                                                            <  1:2880364/150‑1 (MQ=255)
                                                              gAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGc                                                                                         >  2:4920223/1‑150 (MQ=255)
                                                               aaCTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCt                                                                                        >  2:3825144/1‑150 (MQ=255)
                                                                  ttGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCa                                                                                     >  2:513965/1‑150 (MQ=255)
                                                                  ttGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCa                                                                                     >  2:3476254/1‑150 (MQ=255)
                                                                   tGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCAc                                                                                    <  2:5079537/150‑1 (MQ=255)
                                                                    gTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCAcc                                                                                   >  2:4116452/1‑150 (MQ=255)
                                                                       cACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGAt                                                                                >  1:1014137/1‑150 (MQ=255)
                                                                        aCCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATc                                                                               >  1:4574120/1‑150 (MQ=255)
                                                                        aCCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATc                                                                               >  1:458001/1‑150 (MQ=255)
                                                                        aCCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATc                                                                               <  2:3833862/150‑1 (MQ=255)
                                                                         ccGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAAcga                                                                                                                                          <  2:128071/90‑3 (MQ=255)
                                                                         ccGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAAcga                                                                                                                                          >  1:128071/1‑88 (MQ=255)
                                                                             tttCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACAt                                                                          >  1:3210524/1‑150 (MQ=255)
                                                                                 atcaGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGa                                                                      >  1:4830320/1‑150 (MQ=255)
                                                                                  tcaGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGAc                                                                     >  2:1465471/1‑150 (MQ=255)
                                                                                  tcaGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGAc                                                                     >  2:4574938/1‑150 (MQ=255)
                                                                                    aGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACgg                                                                   >  1:4645806/1‑150 (MQ=255)
                                                                                      aCCGATGTTGTCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGtt                                                                 <  1:4999756/150‑1 (MQ=255)
                                                                                      aCCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGtt                                                                 <  2:2494303/150‑1 (MQ=255)
                                                                                      aCCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGACTTCACCGATCGACATGCGACGGtt                                                                 >  1:1522678/1‑150 (MQ=255)
                                                                                       ccGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTg                                                                >  2:4744462/1‑150 (MQ=255)
                                                                                          aTGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGAtt                                                             >  2:4797262/1‑150 (MQ=255)
                                                                                             ttGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTcc                                                          <  1:4227975/150‑1 (MQ=255)
                                                                                               gCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCtt                                                        >  2:1729479/1‑150 (MQ=255)
                                                                                               gCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCtt                                                        <  1:1647313/150‑1 (MQ=255)
                                                                                               gCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCtt                                                        <  2:2570065/150‑1 (MQ=255)
                                                                                               gCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCtt                                                        <  1:1074437/150‑1 (MQ=255)
                                                                                                 cTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGa                                                      >  1:4071083/1‑150 (MQ=255)
                                                                                                 cTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGa                                                      >  1:2415476/1‑150 (MQ=255)
                                                                                                  ttCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGAt                                                     >  1:582718/1‑150 (MQ=255)
                                                                                                    ccTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGAtct                                                   >  1:693392/1‑150 (MQ=255)
                                                                                                     cTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGAtctc                                                  >  2:5057992/1‑150 (MQ=255)
                                                                                                     cTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGAtctc                                                  >  2:2100679/1‑150 (MQ=255)
                                                                                                      tGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCg                                                 <  1:991892/150‑1 (MQ=255)
                                                                                                      tGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCg                                                 <  2:2150717/150‑1 (MQ=255)
                                                                                                            aGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACg                                           <  2:3349699/150‑1 (MQ=255)
                                                                                                              gTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGt                                         >  1:2716893/1‑150 (MQ=255)
                                                                                                                    gAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGcc                                   >  2:3700999/1‑150 (MQ=255)
                                                                                                                          cAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTc                             >  2:3008341/1‑150 (MQ=255)
                                                                                                                           aGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCg                            <  1:2467625/150‑1 (MQ=255)
                                                                                                                           aGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCg                            <  2:4606874/150‑1 (MQ=255)
                                                                                                                               cACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGtct                        >  2:1185244/1‑150 (MQ=255)
                                                                                                                               cACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAAAGGTCAGGCCGGTTTCGGtct                        >  2:2855958/1‑150 (MQ=255)
                                                                                                                                         gtACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGAt              >  2:424127/1‑150 (MQ=255)
                                                                                                                                         gtACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGAt              >  2:2819171/1‑150 (MQ=255)
                                                                                                                                         gtACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGAt              <  1:4129842/150‑1 (MQ=255)
                                                                                                                                          tACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATc             <  2:2180232/150‑1 (MQ=255)
                                                                                                                                           aCTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCa            <  2:3170611/150‑1 (MQ=255)
                                                                                                                                             ttcttGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGc          >  1:3610320/1‑150 (MQ=255)
                                                                                                                                             ttcttGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGc          >  1:444557/1‑150 (MQ=255)
                                                                                                                                               cttGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCtt        <  2:4435242/150‑1 (MQ=255)
                                                                                                                                                     gATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTg  >  1:3508799/1‑150 (MQ=255)
                                                                                                                                                     gATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTg  <  2:2349120/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTG  >  NC_002947/471610‑471908

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: