Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 93,072 C→G 17.5% L174L (CTG→CTC trpB ← tryptophan synthase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294793,0720CG17.5% 489.1 / 91.2 189L174L (CTG→CTCtrpBtryptophan synthase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (88/68);  new base G (14/19);  total (102/87)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 191 total aligned reads displayed.

AACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGAC  >  NC_002947/92923‑93221
                                                                                                                                                     |                                                                                                                                                     
 aCTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCa                                                                                                                                                      <  1:1473797/150‑1 (MQ=255)
 aCTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCa                                                                                                                                                      <  2:4975677/150‑1 (MQ=255)
   tGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAgg                                                                                                                                                    <  2:2038954/150‑1 (MQ=255)
   tGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAgg                                                                                                                                                    <  1:1268688/150‑1 (MQ=255)
   tGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAgg                                                                                                                                                    >  1:2591777/1‑150 (MQ=255)
              ttcggtGCCGATCACCGACTGGAAATCGCGCACGAGTTCCGGATAAGGAGCAGGACCGGCCGCGGTGCCAATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGcc                                                                                                                                         >  2:4948615/6‑150 (MQ=255)
                 tttGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCcggc                                                                                                                                      >  2:2805018/1‑150 (MQ=255)
                 gttGCCGATGATCGACTGGAAGACGCGAACCATGGCCGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTACAAGGTGTCTACGACGTTGGTAACCCAGTCACGCAGGGCTACGTTCATGGCGTCTTTCAGGGTACAGGTGAcggc                                                                                                                                      >  2:2804953/2‑150 (MQ=255)
                  ttGCCGATGATCGACTGGAAGTCGCGCACCTTGGCGGGGTACGGGTGTGGGCCAGCCTCGGTGCCGATCAGGTAGAAGGTGTCTTCGTCGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCcggcg                                                                                                                                     <  1:2396467/150‑1 (MQ=255)
                  ttGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCcggcg                                                                                                                                     <  1:1601038/150‑1 (MQ=255)
                  ttGCCGATGATCGACTGGAAATCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCcggcg                                                                                                                                     >  1:628558/1‑150 (MQ=255)
                   tGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCcggcgg                                                                                                                                    >  2:4868310/1‑150 (MQ=255)
                       gatgatCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAAc                                                                                                                                >  1:1214260/1‑150 (MQ=255)
                           accGACTGGAAATCGCGCACGAGTTCCGGATAAGGAGCAGGACCGGCCGCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGc                                                                                                                            <  2:4072849/148‑1 (MQ=255)
                             cGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCAc                                                                                                                          <  1:4985217/150‑1 (MQ=255)
                             cGACTGGAAATCGCGCACGAGTTCCGGATAAGGAGCAGGACCGGCCGCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGctt                                                                                                                          >  2:5080934/1‑148 (MQ=255)
                              gACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACg                                                                                                                         >  2:2667725/1‑150 (MQ=255)
                               aCTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGa                                                                                                                        >  1:3774971/1‑150 (MQ=255)
                                 tGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGAtc                                                                                                                      >  1:1398270/1‑150 (MQ=255)
                                 tGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGAtc                                                                                                                      >  1:4222713/1‑150 (MQ=255)
                                    aaGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCg                                                                                                                   >  1:449752/1‑150 (MQ=255)
                                     aGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCgg                                                                                                                  >  2:177969/1‑150 (MQ=255)
                                     aGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCgg                                                                                                                  >  2:1694623/1‑150 (MQ=255)
                                       tCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGgcg                                                                                                                <  1:4045544/150‑1 (MQ=255)
                                         gcgcACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGcc                                                                                                              >  1:4119901/1‑150 (MQ=255)
                                          cgcACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGccc                                                                                                             <  1:606820/150‑1 (MQ=255)
                                             aCCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCcagc                                                                                                          >  2:200142/1‑150 (MQ=255)
                                              ccATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCcagca                                                                                                         <  1:2298603/150‑1 (MQ=255)
                                              ccATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCcagca                                                                                                         <  1:4869592/150‑1 (MQ=255)
                                                     ggggTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCAt                                                                                                  >  2:3632405/1‑150 (MQ=255)
                                                       ggTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGc                                                                                                <  1:3376523/150‑1 (MQ=255)
                                                         tACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCgg                                                                                              >  1:4273282/1‑150 (MQ=255)
                                                               tgtgGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACg                                                                                        >  1:269824/1‑150 (MQ=255)
                                                                gtgGGCCAGCCACGGTGCCGATCAGGTAGATGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATTGCGTCTTTCAGGGTTCCTGTGCCGGCGGTAACTGGCACGATCTCGGCGCCGAGCTGCTTCATGCGGAACACGt                                                                                       <  2:4637928/150‑1 (MQ=255)
                                                                  gggCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTc                                                                                     <  1:3747931/150‑1 (MQ=255)
                                                                  gggCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTc                                                                                     <  1:2480748/150‑1 (MQ=255)
                                                                  ggaccggccgCGGTGCCGATCAGGTAATAGGTGTCGTCACCATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCAGTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGttt                                                                                     <  2:4776785/144‑2 (MQ=255)
                                                                  ggaccggccgCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGttt                                                                                     <  1:4810853/144‑2 (MQ=255)
                                                                  ggaccggccgCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGttt                                                                                     <  2:2423357/144‑2 (MQ=255)
                                                                   ggCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCg                                                                                    >  1:3486436/1‑150 (MQ=255)
                                                                   ggCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCg                                                                                    >  2:611002/1‑150 (MQ=255)
                                                                   ggCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCg                                                                                    >  1:794181/1‑150 (MQ=255)
                                                                      cagccaCGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCct                                                                                 >  1:1960963/1‑150 (MQ=255)
                                                                      cagccaCGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCct                                                                                 >  1:1960886/1‑150 (MQ=255)
                                                                         ccaCGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCctgct                                                                              <  2:1840876/150‑1 (MQ=255)
                                                                            cGGCGCCGATCAGGTAGAAGGTGTATTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGgc                                                                           <  2:3883181/150‑1 (MQ=255)
                                                                                gCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTc                                                                       >  2:2096912/1‑150 (MQ=255)
                                                                                  cGATCAGGTAATAGGTGTCGTCGACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTcaa                                                                     <  2:4854538/150‑3 (MQ=255)
                                                                                  cGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTcaa                                                                     <  1:1669234/150‑3 (MQ=255)
                                                                                     tCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGt                                                                  <  1:202050/150‑1 (MQ=255)
                                                                                       aGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCg                                                                <  1:3056960/150‑1 (MQ=255)
                                                                                         gTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGt                                                              <  2:4519075/150‑1 (MQ=255)
                                                                                         gTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGt                                                              <  1:2343428/150‑1 (MQ=255)
                                                                                          taataGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCCTGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTg                                                             <  2:763162/146‑1 (MQ=255)
                                                                                              aGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCAc                                                         <  1:93912/150‑1 (MQ=255)
                                                                                                   tCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATg                                                    <  2:3785284/150‑1 (MQ=255)
                                                                                                     ttCGACGTTGGTAACCCAGTCTCGCAGGGCTTCGTTCATGGGGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCTCGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGTCACCCTTGTa                                                  <  2:1237864/150‑1 (MQ=255)
                                                                                                     ttCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTa                                                  >  2:2585334/1‑150 (MQ=255)
                                                                                                     ttCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTa                                                  >  1:4028143/1‑150 (MQ=255)
                                                                                                     ttCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTa                                                  >  1:3157957/1‑150 (MQ=255)
                                                                                                     ttCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTa                                                  >  1:3312065/1‑150 (MQ=255)
                                                                                                     ttCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTa                                                  >  2:2488619/1‑150 (MQ=255)
                                                                                                     gtcaacattGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCAtata                                                  >  1:1902883/8‑150 (MQ=255)
                                                                                                      tcaacattGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAg                                                 >  1:436427/7‑150 (MQ=255)
                                                                                                        gACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGAt                                               <  2:2761657/150‑1 (MQ=255)
                                                                                                          cGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCa                                             <  1:3246512/150‑1 (MQ=255)
                                                                                                           gTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCAc                                            <  1:3724062/150‑1 (MQ=255)
                                                                                                            ttGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACg                                           <  1:1004914/150‑1 (MQ=255)
                                                                                                               gtgaCCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCaa                                        <  1:3836165/147‑1 (MQ=255)
                                                                                                               gTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCaa                                        <  2:3822530/150‑1 (MQ=255)
                                                                                                               gTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCaa                                        >  1:4572710/1‑150 (MQ=255)
                                                                                                                tAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAag                                       >  1:2659644/1‑150 (MQ=255)
                                                                                                                 aaCCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAagg                                      <  2:1930346/150‑1 (MQ=255)
                                                                                                                   cccAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAaggca                                    >  2:4341610/1‑150 (MQ=255)
                                                                                                                   cccAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAaggca                                    <  2:1950149/150‑1 (MQ=255)
                                                                                                                    ccAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATg                                                    >  2:2593440/1‑133 (MQ=255)
                                                                                                                    ccAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAaggcag                                   >  2:5118071/1‑150 (MQ=255)
                                                                                                                    ccAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAaggcag                                   <  1:2601599/150‑1 (MQ=255)
                                                                                                                     caatCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAaggcagg                                  >  1:1335683/4‑150 (MQ=255)
                                                                                                                      aatCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAaggcaggc                                 <  2:3335757/148‑1 (MQ=255)
                                                                                                                      aatCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGATTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTAGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAaggcaggc                                 >  2:3074571/3‑150 (MQ=255)
                                                                                                                       gTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGcc                                <  1:3423212/150‑1 (MQ=255)
                                                                                                                       gTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGcc                                >  2:914367/1‑150 (MQ=255)
                                                                                                                          aCGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGaa                             >  2:3653837/1‑150 (MQ=255)
                                                                                                                            gCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAgc                           <  1:5107685/150‑1 (MQ=255)
                                                                                                                            gCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAgc                           >  2:296175/1‑150 (MQ=255)
                                                                                                                            gCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAgc                           >  2:1595189/1‑150 (MQ=255)
                                                                                                                              agcgCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAacgg                         <  2:71402/147‑2 (MQ=255)
                                                                                                                                 gCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGcgcgcg                      >  1:1808243/1‑150 (MQ=255)
                                                                                                                                    tCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCgcggcg                   >  2:2246241/1‑150 (MQ=255)
                                                                                                                                      gTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGAc                 >  1:2699518/1‑150 (MQ=255)
                                                                                                                                      gTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGAc                 >  1:2220529/1‑150 (MQ=255)
                                                                                                                                       gtCATGGCGTCTTTCAGGGGACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCACGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCTTGTACATCACGCAAGGCAGGCCGTAGCGCGCGGCGACa                <  1:1026776/149‑1 (MQ=255)
                                                                                                                                          aTGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtg             >  2:3426551/1‑150 (MQ=255)
                                                                                                                                          aTGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtg             >  2:641747/1‑150 (MQ=255)
                                                                                                                                           tGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtgg            >  1:4865055/1‑150 (MQ=255)
                                                                                                                                             gCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcg          >  1:4646778/1‑150 (MQ=255)
                                                                                                                                                  tttCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcggtggc     >  2:2370387/1‑150 (MQ=255)
                                                                                                                                                  tttCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcggtggc     >  1:4357228/1‑150 (MQ=255)
                                                                                                                                                    tCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGa   <  1:3772555/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGAC  >  NC_002947/92923‑93221

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: