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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947471,7650AG 6.6% 566.5 / 17.3 182I391I (ATT→ATCrpoDRNA polymerase sigma 70 factor
Reads supporting (aligned to +/- strand):  ref base A (89/80);  new base G (7/5);  total (97/85)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 188 total aligned reads displayed.

CCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTGGCAAGC  >  NC_002947/471617‑471914
                                                                                                                                                    |                                                                                                                                                     
 cTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAAt                                                                                                                                                     <  1:2325116/150‑1 (MQ=255)
   gACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCa                                                                                                                                                   >  1:2820951/1‑150 (MQ=255)
   gACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCa                                                                                                                                                   >  1:4739836/1‑150 (MQ=255)
     cGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCAcc                                                                                                                                                 <  2:4297909/150‑1 (MQ=255)
       gATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAg                                                                                                                                               >  1:4513983/1‑150 (MQ=255)
       gATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAg                                                                                                                                               >  1:468180/1‑150 (MQ=255)
                 gTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTc                                                                                                                                     >  2:1328735/1‑150 (MQ=255)
                 gTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTc                                                                                                                                     >  2:4080081/1‑150 (MQ=255)
                    gCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGcc                                                                                                                                  <  2:5189811/150‑1 (MQ=255)
                    gCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGct                                                                                                                                  >  2:951743/1‑149 (MQ=255)
                     cATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCt                                                                                                                                 <  1:4825056/150‑1 (MQ=255)
                     cATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCt                                                                                                                                 <  2:2850456/150‑1 (MQ=255)
                        aGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTcca                                                                                                                              >  1:1625036/1‑148 (MQ=255)
                         ggTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGAc                                                                                                                             >  2:3920083/1‑150 (MQ=255)
                         ggTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGAc                                                                                                                             <  1:3598059/150‑1 (MQ=255)
                           tCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGtt                                                                                                                                      <  1:2292477/139‑1 (MQ=255)
                           tCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGtt                                                                                                                                      >  2:2292477/1‑139 (MQ=255)
                           tCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCa                                                                                                                           <  1:1902352/150‑1 (MQ=255)
                           tCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCa                                                                                                                           >  2:298188/1‑150 (MQ=255)
                            cGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAt                                                                                                                          >  1:5176727/1‑150 (MQ=255)
                            cGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCACCAt                                                                                                                          <  1:45901/150‑1 (MQ=255)
                             gagaACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAtt                                                                                                                         <  2:189347/150‑1 (MQ=255)
                             gagaACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAtt                                                                                                                         >  2:3201780/1‑150 (MQ=255)
                             gagaACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAtt                                                                                                                         >  2:5247540/1‑150 (MQ=255)
                              agaACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAttt                                                                                                                        >  2:3872103/1‑150 (MQ=255)
                              agaACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAttt                                                                                                                        <  1:4129258/150‑1 (MQ=255)
                                    tGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATttctttct                                                                                                                  <  2:2188615/150‑1 (MQ=255)
                                     gTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTtctttctt                                                                                                                 <  2:324741/150‑1 (MQ=255)
                                     gTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTtctttctt                                                                                                                 <  2:382040/150‑1 (MQ=255)
                                      tAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTg                                                                                                                <  2:4361761/150‑1 (MQ=255)
                                       aaccacgacgaTATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCACCATTTCCTTCTTgg                                                                                                               >  2:1089107/6‑150 (MQ=255)
                                        gCCGCGACGGTATTCGACCGTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGgc                                                                                                              <  2:5049857/150‑1 (MQ=255)
                                        gCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGgc                                                                                                              <  1:2472355/150‑1 (MQ=255)
                                               cGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGcggcgg                                                                                                       <  2:2620143/150‑1 (MQ=255)
                                                ggTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCggg                                                                                                      >  2:1720034/1‑150 (MQ=255)
                                                ggTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCggg                                                                                                      >  2:2154950/1‑150 (MQ=255)
                                                    ttCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCtt                                                                                                  <  1:1998446/150‑1 (MQ=255)
                                                    ttCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCtt                                                                                                  >  2:3122064/1‑150 (MQ=255)
                                                       gAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGc                                                                                               <  1:3393556/150‑1 (MQ=255)
                                                           ttGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCa                                                                                           >  2:4163542/1‑150 (MQ=255)
                                                           ttGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCa                                                                                           >  1:936122/1‑150 (MQ=255)
                                                             gTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCAcc                                                                                         >  1:3846212/1‑150 (MQ=255)
                                                             gTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCAcc                                                                                         >  2:2662692/1‑150 (MQ=255)
                                                                 aCCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCa                                                                                                                                                   <  1:3860942/88‑1 (MQ=255)
                                                                 aCCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCa                                                                                                                                                   >  2:3860942/1‑88 (MQ=255)
                                                                 aCCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCt                                                                                                   >  1:624953/1‑136 (MQ=255)
                                                                 aCCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATc                                                                                     >  2:2240723/1‑150 (MQ=255)
                                                                 aCCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATc                                                                                     >  1:4934915/1‑150 (MQ=255)
                                                                 aCCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATc                                                                                     <  2:1546532/150‑1 (MQ=255)
                                                                 aCCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATc                                                                                     <  2:1474087/150‑1 (MQ=255)
                                                                   cGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGa                                                                                   >  2:1437816/1‑150 (MQ=255)
                                                                      tttCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACAt                                                                                >  2:3537845/1‑150 (MQ=255)
                                                                          atcaGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGa                                                                            >  1:1411048/1‑150 (MQ=255)
                                                                           tcaGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGAc                                                                           >  2:3068575/1‑150 (MQ=255)
                                                                             aGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCgcggc                                                                                                         <  1:1139924/118‑1 (MQ=255)
                                                                             aGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCgcggc                                                                                                         >  2:1139924/1‑118 (MQ=255)
                                                                             aGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACgg                                                                         <  2:3468932/150‑1 (MQ=255)
                                                                                 cGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGa                                                                     >  1:5320899/1‑150 (MQ=255)
                                                                                 cGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGa                                                                     <  2:2687245/150‑1 (MQ=255)
                                                                                   aTGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGAtt                                                                   <  2:571170/150‑1 (MQ=255)
                                                                                     gTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTc                                                                 >  1:1206258/1‑150 (MQ=255)
                                                                                       tGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCt                                                               >  2:3452386/1‑150 (MQ=255)
                                                                                            tCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGAtc                                                          >  2:1067088/1‑150 (MQ=255)
                                                                                            tCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACAATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGAta                                                          >  2:4135994/1‑149 (MQ=255)
                                                                                              cTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGAtctc                                                        <  2:4708495/150‑1 (MQ=255)
                                                                                               tGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCg                                                       <  1:1269849/150‑1 (MQ=255)
                                                                                               tGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCAACAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGAGTTCCTTGATCTCg                                                       >  2:4490909/1‑150 (MQ=255)
                                                                                                      ggTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACgg                                                >  2:3980416/1‑150 (MQ=255)
                                                                                                      ggTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACgg                                                >  2:3980350/1‑150 (MQ=255)
                                                                                                        tCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTc                                              <  1:185680/150‑1 (MQ=255)
                                                                                                           aGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAgg                                           <  1:5247540/150‑1 (MQ=255)
                                                                                                             gAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGcc                                         <  1:1775890/150‑1 (MQ=255)
                                                                                                             gAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGcc                                         >  1:2941525/1‑150 (MQ=255)
                                                                                                              aaCTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCg                                        <  2:2820951/150‑1 (MQ=255)
                                                                                                                 tGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGtt                                     >  2:3271648/1‑150 (MQ=255)
                                                                                                                    aGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCg                                  <  1:5398720/150‑1 (MQ=255)
                                                                                                                     ggCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCgg                                 <  2:2378438/150‑1 (MQ=255)
                                                                                                                      gCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCAc                                                                                          >  2:3939546/1‑92 (MQ=255)
                                                                                                                      gCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCAc                                                                                          <  1:3939546/92‑1 (MQ=255)
                                                                                                                      gCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGt                                <  1:2519252/150‑1 (MQ=255)
                                                                                                                        cACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGtct                              >  2:181284/1‑150 (MQ=255)
                                                                                                                         aCGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGtctc                             <  1:1926148/150‑1 (MQ=255)
                                                                                                                           ggTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTACTTGGCGCGGCGGGCCTTGGCATCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCaa                           >  2:1911128/1‑150 (MQ=255)
                                                                                                                            gTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAg                          <  2:727540/150‑1 (MQ=255)
                                                                                                                                 tgtACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGa                     <  1:2406207/150‑1 (MQ=255)
                                                                                                                                 tgtACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGa                     <  2:55200/150‑1 (MQ=255)
                                                                                                                                 tgtACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGa                     >  1:4374629/1‑150 (MQ=255)
                                                                                                                                   tACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATc                   <  1:4430935/150‑1 (MQ=255)
                                                                                                                                     cttcttGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAg                 <  1:3455413/150‑1 (MQ=255)
                                                                                                                                      ttcttGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGc                <  1:2439938/150‑1 (MQ=255)
                                                                                                                                      ttcttGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGc                >  1:2373754/1‑150 (MQ=255)
                                                                                                                                            gCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGt          >  2:1234851/1‑150 (MQ=255)
                                                                                                                                             cGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGtt         <  2:4739836/150‑1 (MQ=255)
                                                                                                                                             cGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGtt         <  1:1720034/150‑1 (MQ=255)
                                                                                                                                             cGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGtt         <  1:2045661/150‑1 (MQ=255)
                                                                                                                                             cGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGtt         >  1:4756278/1‑150 (MQ=255)
                                                                                                                                              gATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTg        >  1:3114128/1‑150 (MQ=255)
                                                                                                                                              gATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTg        >  2:2657248/1‑150 (MQ=255)
                                                                                                                                                   aaaTCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTGGCAAg   >  2:3663392/1‑150 (MQ=255)
                                                                                                                                                    aaTCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTCTCGGTCTCAAGGTCGATCAGCTTCTGTTGGCAAGc  >  1:854316/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTGGCAAGC  >  NC_002947/471617‑471914

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: