Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 93,087 G→A 16.8% A169A (GCC→GCT trpB ← tryptophan synthase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294793,0870GA16.8% 400.9 / 66.7 172A169A (GCC→GCTtrpBtryptophan synthase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (81/62);  new base A (11/18);  total (92/80)
Fisher's exact test for biased strand distribution p-value = 7.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 185 total aligned reads displayed.

TCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGG  >  NC_002947/92938‑93235
                                                                                                                                                     |                                                                                                                                                    
tctTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCcg                                                                                                                                                      >  1:4908149/1‑150 (MQ=255)
  tttGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCcggc                                                                                                                                                    >  2:2804877/1‑150 (MQ=255)
  tttGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCcggc                                                                                                                                                    <  2:3514490/150‑1 (MQ=255)
  tttGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGCCGTCTTTCAGGGTACCGGTGCcggc                                                                                                                                                    >  2:2804983/1‑150 (MQ=255)
   ttGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCcggcg                                                                                                                                                   <  1:1601038/150‑1 (MQ=255)
   ttGCCGATGATCGACTGGAAATCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCcggcg                                                                                                                                                   >  1:628558/1‑150 (MQ=255)
    tGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCcggcgg                                                                                                                                                  >  2:4868310/1‑150 (MQ=255)
      ccGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTa                                                                                                                                                >  1:1342423/1‑150 (MQ=255)
        gatgatCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAAc                                                                                                                                              >  1:1214260/1‑150 (MQ=255)
        gatgatCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAAc                                                                                                                                              >  1:3802423/1‑150 (MQ=255)
           gatCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTgg                                                                                                                                           <  1:1667596/150‑1 (MQ=255)
            accGACTGGAAATCGCGCACGAGTTCCGGATAAGGAGCAGGACCGGCCGCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGc                                                                                                                                          <  2:4072849/148‑1 (MQ=255)
              cGACTGGAAATCGCGCACGAGTTCCGGATAAGGAGCAGGACCGGCCGCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGctt                                                                                                                                        >  2:5080934/1‑148 (MQ=255)
                  tGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGAtc                                                                                                                                    >  1:4222713/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)
                                          tACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCgg                                                                                                            >  1:4273282/1‑150 (MQ=255)
                                                tgtgGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACg                                                                                                      >  1:269824/1‑150 (MQ=255)
                                                   gggCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTc                                                                                                   <  1:3747931/150‑1 (MQ=255)
                                                   ggaccggccgCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGttt                                                                                                   <  2:2423357/144‑2 (MQ=255)
                                                   ggaccggccgCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGttt                                                                                                   <  2:178660/144‑2 (MQ=255)
                                                   ggaccggccgCGGTGCCGATCAGGTAATAGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGttt                                                                                                   <  1:4810853/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)
                                                     gccagccaCGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGc                                                                                                 >  2:45192/1‑150 (MQ=255)
                                                       cagccaCGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCct                                                                                               >  1:1960886/1‑150 (MQ=255)
                                                             cGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGgc                                                                                         <  1:45192/150‑1 (MQ=255)
                                                             cGGCGCCGATCAGGTAGAAGGTGTATTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGgc                                                                                         <  2:3883181/150‑1 (MQ=255)
                                                               gTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGgcgc                                                                                       <  1:2633743/150‑1 (MQ=255)
                                                                 gCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTc                                                                                     >  2:2096912/1‑150 (MQ=255)
                                                                   cGATCAGGTAATAGGTGTCGTCGACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTcaa                                                                                   <  2:4854538/150‑3 (MQ=255)
                                                                        aGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCg                                                                              <  1:3056960/150‑1 (MQ=255)
                                                                          gTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGt                                                                            <  2:1859716/150‑1 (MQ=255)
                                                                          gTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGt                                                                            <  1:2343428/150‑1 (MQ=255)
                                                                           taataGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCCTGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTg                                                                           <  2:763162/146‑1 (MQ=255)
                                                                           taataGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTg                                                                           >  1:2240965/5‑150 (MQ=255)
                                                                               aGGTGTCGTCAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCAc                                                                       <  1:93912/150‑1 (MQ=255)
                                                                                   gtcgtcAACATTGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCAt                                                                   <  2:1706642/150‑1 (MQ=255)
                                                                                    tCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATg                                                                  <  2:3785284/150‑1 (MQ=255)
                                                                                      ttCGACGTTGGTAACCCAGTCTCGCAGGGCTTCGTTCATGGGGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCTCGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGTCACCCTTGTa                                                                <  2:1237864/150‑1 (MQ=255)
                                                                                      ttCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTa                                                                >  1:4028143/1‑150 (MQ=255)
                                                                                      ttCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTa                                                                >  1:1406406/1‑150 (MQ=255)
                                                                                      ttCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTa                                                                >  1:3157957/1‑150 (MQ=255)
                                                                                      ttCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTa                                                                >  1:3312065/1‑150 (MQ=255)
                                                                                      gtcaacattGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCAtata                                                                >  1:1902883/8‑150 (MQ=255)
                                                                                       tcaacattGGTGACCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAg                                                               >  1:436069/7‑150 (MQ=255)
                                                                                         gACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGAt                                                             <  2:2761657/150‑1 (MQ=255)
                                                                                          aCGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATc                                                            <  2:1610536/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)
                                                                                                gtgaCCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCAcc                                                                      >  1:2463331/4‑134 (MQ=255)
                                                                                                gtgaCCCAATCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCAcc                                                                      <  2:2463331/131‑1 (MQ=255)
                                                                                                gTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCaa                                                      <  2:3822530/150‑1 (MQ=255)
                                                                                                  aaCCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAagg                                                    <  2:1930346/150‑1 (MQ=255)
                                                                                                    cccAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAaggca                                                  >  2:4341610/1‑150 (MQ=255)
                                                                                                    cccAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAaggca                                                  >  2:4341829/1‑150 (MQ=255)
                                                                                                    cccAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAaggca                                                  <  2:1950149/150‑1 (MQ=255)
                                                                                                     ccAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATg                                                                  <  1:2593440/133‑1 (MQ=255)
                                                                                                     ccAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAaggcag                                                 >  2:5118071/1‑150 (MQ=255)
                                                                                                      caatCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAaggcagg                                                >  1:3915333/4‑150 (MQ=255)
                                                                                                      caatCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAaggcagg                                                >  1:1335678/4‑150 (MQ=255)
                                                                                                       aatCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAaggcaggc                                               <  2:3335757/148‑1 (MQ=255)
                                                                                                       aatCACGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGATTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTAGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAaggcaggc                                               >  2:3074571/3‑150 (MQ=255)
                                                                                                       aatCACGAAGCGATTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGTCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGATTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAaggcaggc                                               <  1:1973850/148‑1 (MQ=255)
                                                                                                        gTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGcc                                              <  1:3423204/150‑1 (MQ=255)
                                                                                                        gTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGcc                                              <  1:3423212/150‑1 (MQ=255)
                                                                                                        gTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGcc                                              >  2:914367/1‑150 (MQ=255)
                                                                                                           tcGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGaa                                           <  1:4439776/149‑1 (MQ=255)
                                                                                                           aCGAAGCGCTTCGTTCATGGCGTCTTTGAGCGTGCCATTGCCAGCCGTTACCGGCTTTACTTCGGCACCCAGCAGCTTCATGCGGAAGACGTTTGGCTTCTGACGTTCAACGTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGaa                                           >  2:3653837/1‑150 (MQ=255)
                                                                                                             gCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAgc                                         <  1:5107685/150‑1 (MQ=255)
                                                                                                             gCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAgc                                         >  2:296175/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:2032626/1‑150 (MQ=255)
                                                                                                                       gTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGAc                               <  1:4868310/150‑1 (MQ=255)
                                                                                                                       gTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGAc                               >  1:2699518/1‑150 (MQ=255)
                                                                                                                        gtCATGGCGTCTTTCAGGGGACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCACGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCTTGTACATCACGCAAGGCAGGCCGTAGCGCGCGGCGACa                              <  1:1026776/149‑1 (MQ=255)
                                                                                                                           aTGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtg                           >  2:641747/1‑150 (MQ=255)
                                                                                                                            tGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtgg                          >  1:4865055/1‑150 (MQ=255)
                                                                                                                              gCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcg                        >  1:4646778/1‑150 (MQ=255)
                                                                                                                               cGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcgg                       >  2:4740864/1‑150 (MQ=255)
                                                                                                                                gTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGGTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcggt                      <  2:2591777/150‑1 (MQ=255)
                                                                                                                                 tCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcggtg                     <  1:5085746/150‑1 (MQ=255)
                                                                                                                                  cTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcggtgg                    <  2:2659644/150‑1 (MQ=255)
                                                                                                                                   tttCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAgtggcggtggc                   >  1:4357228/1‑150 (MQ=255)
                                                                                                                                     tCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGa                 <  1:3772555/150‑1 (MQ=255)
                                                                                                                                     tCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGa                 <  2:4966688/150‑1 (MQ=255)
                                                                                                                                     tCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGa                 >  2:3644992/1‑150 (MQ=255)
                                                                                                                                        gggTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGacac              >  1:2380035/1‑150 (MQ=255)
                                                                                                                                           tACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGt           >  1:3841668/1‑150 (MQ=255)
                                                                                                                                           tACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGt           <  2:2525099/150‑1 (MQ=255)
                                                                                                                                           tACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGt           <  2:2581484/150‑1 (MQ=255)
                                                                                                                                           tACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGt           <  2:2667805/150‑1 (MQ=255)
                                                                                                                                                    cggcggTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCgg  <  1:860023/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGGGCTTCGTTCATGGCGTCTTTCAGGGTACCGGTGCCGGCGGTAACTGGCACGATCTCGGCGCCCAGCAGCTTCATGCGGAACACGTTCGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGG  >  NC_002947/92938‑93235

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: