Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 471,783 C→G 6.1% A385A (GCG→GCC rpoD ← RNA polymerase sigma 70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947471,7830CG 6.1% 819.5 / 31.4 245A385A (GCG→GCCrpoDRNA polymerase sigma 70 factor
Reads supporting (aligned to +/- strand):  ref base C (119/111);  new base G (8/7);  total (127/118)
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.81e-01

Only 100 of 253 total aligned reads displayed.

GTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTGGCAAGCGACGATGGCAGCGTTCTT  >  NC_002947/471634‑471932
                                                                                                                                                     |                                                                                                                                                     
gTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTc                                                                                                                                                       >  2:8165708/1‑150 (MQ=255)
  ggCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGc                                                                                                                                                     <  1:1159698/150‑1 (MQ=255)
  cgCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGc                                                                                                                                                     <  1:6963939/149‑1 (MQ=255)
   gCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGct                                                                                                                                                    >  2:7627271/1‑149 (MQ=255)
     aTAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTc                                                                                                                                                  >  1:2207332/1‑150 (MQ=255)
     aTAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTc                                                                                                                                                  >  2:8047080/1‑150 (MQ=255)
     aTAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTc                                                                                                                                                  >  2:7186219/1‑150 (MQ=255)
       aGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGa                                                                                                                                                <  1:5201875/150‑1 (MQ=255)
       aGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGa                                                                                                                                                <  1:160213/150‑1 (MQ=255)
         gTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGACCTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGAcc                                                                                                                                              >  2:9011421/1‑150 (MQ=255)
         gTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGAcc                                                                                                                                              <  1:2585927/150‑1 (MQ=255)
          tCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTccacca                                                                                                                                             >  1:7409048/1‑150 (MQ=255)
          tCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTccacca                                                                                                                                             >  1:3418519/1‑150 (MQ=255)
            gagaACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAtt                                                                                                                                           >  2:278847/1‑150 (MQ=255)
               aaCTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAtttct                                                                                                                                        <  2:1362237/150‑1 (MQ=255)
                 cTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAtttcttt                                                                                                                                      <  2:8639532/150‑1 (MQ=255)
                  ttGTAGCCGAGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAtttctttc                                                                                                                                     >  2:1118682/1‑150 (MQ=255)
                    gTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTtctttctt                                                                                                                                   <  1:2341752/150‑1 (MQ=255)
                    gTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTtctttctt                                                                                                                                   <  2:1484983/150‑1 (MQ=255)
                    gTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTtctttctt                                                                                                                                   <  1:1834601/150‑1 (MQ=255)
                       gCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGgc                                                                                                                                >  1:995187/1‑150 (MQ=255)
                        ccGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGgcg                                                                                                                               >  1:8179191/1‑150 (MQ=255)
                           cgacgaTATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCACCATTTCCTTCTTGgcgatg                                                                                                                            <  2:5688073/150‑4 (MQ=255)
                            gACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCgcggc                                                                                                                           >  1:8597092/1‑150 (MQ=255)
                             aCGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCACCATTTCCTTCTTGgcgatggt                                                                                                                          <  1:3880965/150‑6 (MQ=255)
                              cGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGcggcgg                                                                                                                         <  2:2226066/150‑1 (MQ=255)
                                gTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGc                                                                                                                       >  2:4868612/1‑150 (MQ=255)
                                gTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGc                                                                                                                       >  1:1735249/1‑150 (MQ=255)
                                   ttCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCtt                                                                                                                    >  2:9127231/1‑150 (MQ=255)
                                   gtCGAACTTGTCCACCGCTTTCATCATACCGATGTTCCCTTCCTGGATCAGGACGAGGAACTGCAGGCCACGGTTGGTCTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCTACCATTACATTCTTGGCGCGGCGGGCCtt                                                                                                                    >  2:4527193/2‑150 (MQ=255)
                                    tCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTg                                                                                                                   <  2:2397588/150‑1 (MQ=255)
                                      gAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGAGGCGGGCCTTGGc                                                                                                                 >  2:8094239/1‑150 (MQ=255)
                                       aaCTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCt                                                                                                                <  2:7681949/150‑1 (MQ=255)
                                       aaCTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCt                                                                                                                <  2:1902412/150‑1 (MQ=255)
                                        aCTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCtt                                                                                                               >  1:8859934/1‑150 (MQ=255)
                                              ccACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCgt                                                                                                         >  2:6029876/1‑149 (MQ=255)
                                              ccACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGa                                                                                                         <  1:6408312/150‑1 (MQ=255)
                                               cACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGAt                                                                                                        >  1:5297700/1‑150 (MQ=255)
                                                aCCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATc                                                                                                       >  1:3248969/1‑150 (MQ=255)
                                                aCCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATc                                                                                                       >  1:92829/1‑150 (MQ=255)
                                                   gCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGAc                                                                                                    <  1:1116964/150‑1 (MQ=255)
                                                    cTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACa                                                                                                   <  2:3622020/150‑1 (MQ=255)
                                                       tcatcaGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGc                                                                                                >  2:6621056/1‑150 (MQ=255)
                                                         atcaGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGa                                                                                              >  1:3449130/1‑150 (MQ=255)
                                                         atcaGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGa                                                                                              <  1:5086466/150‑1 (MQ=255)
                                                         atcaGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGa                                                                                              >  2:6928767/1‑150 (MQ=255)
                                                          tcaGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGAc                                                                                             >  1:79472/1‑150 (MQ=255)
                                                             gACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGt                                                                                          <  2:1704136/150‑1 (MQ=255)
                                                             gACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGt                                                                                          >  1:294222/1‑150 (MQ=255)
                                                              aCCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGtt                                                                                         >  2:845196/1‑150 (MQ=255)
                                                                cGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGa                                                                                       >  1:2212874/1‑150 (MQ=255)
                                                                cGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGa                                                                                       <  2:8170612/150‑1 (MQ=255)
                                                                     ttGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTcc                                                                                  <  1:88024/150‑1 (MQ=255)
                                                                      tGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCt                                                                                 <  2:3068770/150‑1 (MQ=255)
                                                                       gCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCtt                                                                                <  1:2418550/150‑1 (MQ=255)
                                                                       gCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCtt                                                                                <  2:8274256/150‑1 (MQ=255)
                                                                       gCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCtt                                                                                <  2:2539685/150‑1 (MQ=255)
                                                                          ttCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGAt                                                                             <  1:8472431/150‑1 (MQ=255)
                                                                             cTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGAtctc                                                                          <  1:3320847/150‑1 (MQ=255)
                                                                               ggATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCgg                                                                        <  1:8567822/150‑1 (MQ=255)
                                                                                      gTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGt                                                                 >  2:5147546/1‑150 (MQ=255)
                                                                                       tCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTc                                                                >  2:5042274/1‑150 (MQ=255)
                                                                                           ggAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGc                                                            >  1:5512624/1‑150 (MQ=255)
                                                                                           ggAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGc                                                            >  2:7807266/1‑150 (MQ=255)
                                                                                               cTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGt                                                        >  1:806554/1‑150 (MQ=255)
                                                                                                tGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGtt                                                       >  1:8397743/1‑150 (MQ=255)
                                                                                                   aGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCg                                                    <  2:7583336/150‑1 (MQ=255)
                                                                                                   aGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCg                                                    <  1:4070373/150‑1 (MQ=255)
                                                                                                   aGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCg                                                    >  1:1765945/1‑150 (MQ=255)
                                                                                                   aGGCCACGGTTGGTGTACTTCTTGGCAATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCg                                                    <  2:633557/150‑1 (MQ=255)
                                                                                                    ggCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCgg                                                   <  2:2161667/150‑1 (MQ=255)
                                                                                                     gCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGTGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGt                                                  >  2:114489/1‑150 (MQ=255)
                                                                                                     gCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGt                                                  <  2:7868630/150‑1 (MQ=255)
                                                                                                     gCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGt                                                  <  2:3252648/150‑1 (MQ=255)
                                                                                                       cACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGtct                                                >  2:4692256/1‑150 (MQ=255)
                                                                                                        aCGGTTGGTGTACTTCTTGGTGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGtctc                                               <  1:6706862/150‑1 (MQ=255)
                                                                                                                tgtACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGa                                       <  1:785085/150‑1 (MQ=255)
                                                                                                                tgtACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGa                                       >  2:6291732/1‑150 (MQ=255)
                                                                                                                 gtACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGAt                                      <  1:5990788/150‑1 (MQ=255)
                                                                                                                  tACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATc                                     >  2:2819023/1‑150 (MQ=255)
                                                                                                                       cttGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCtt                                <  1:7937076/150‑1 (MQ=255)
                                                                                                                         tGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCt                              <  1:7621020/150‑1 (MQ=255)
                                                                                                                          ggCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTg                             <  1:6103240/150‑1 (MQ=255)
                                                                                                                           gCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGt                            >  1:2029730/1‑150 (MQ=255)
                                                                                                                                ggAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTGGCa                       >  1:193533/1‑150 (MQ=255)
                                                                                                                                 gAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTGTCTGGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTGGCaa                      <  2:5870716/150‑1 (MQ=255)
                                                                                                                                  aaaTCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTGGCAAg                     <  2:6993948/150‑1 (MQ=255)
                                                                                                                                  aaaTCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTGGCAAg                     >  2:2361016/1‑150 (MQ=255)
                                                                                                                                    aTCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTGGCAAGcg                   <  2:1399122/150‑1 (MQ=255)
                                                                                                                                    aTCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTGGCAAGcg                   >  2:2024583/1‑150 (MQ=255)
                                                                                                                                       aCCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTGGCAAGcgacg                >  2:5452417/1‑150 (MQ=255)
                                                                                                                                         cAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTGGCAAGCGACGAt              >  1:6369084/1‑150 (MQ=255)
                                                                                                                                          aGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCTACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTGGCAAGCGACGATg             >  1:2592167/1‑150 (MQ=255)
                                                                                                                                          aGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTGGCAAGCGACGATg             <  1:6585608/150‑1 (MQ=255)
                                                                                                                                          aGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTGGCAAGCGACGATg             <  1:7104384/150‑1 (MQ=255)
                                                                                                                                           gACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTGGCAAGCGACGATgg            <  1:7573972/150‑1 (MQ=255)
                                                                                                                                              gCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTGGCAAGCGACGATGGCAg         >  2:4189469/1‑150 (MQ=255)
                                                                                                                                              gCAGGTTCGCCTCGACCATTTCTTTCGTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTGGCAAGCGACGATGGCAg         <  1:4712473/150‑1 (MQ=255)
                                                                                                                                                  gTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTGGCAAGCGACGATGGCAGCGtt     >  1:5779654/1‑150 (MQ=255)
                                                                                                                                                     cGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTGGCAAGCGACGATGGCAGCGTTCtt  <  1:7186319/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCGATCGACATGCGACGGTTGATTTCCTTGATCTCGGCAACGGTCAGGCCGGTTTCGGTCTCAAGGTCGATCAGCTTCTGTTGGCAAGCGACGATGGCAGCGTTCTT  >  NC_002947/471634‑471932

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: