Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,976,030 Δ1 bp 5.4% coding (919/1104 nt) pheA → bifunctional chorismate mutase/prephenate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,976,0300C. 5.4% 335.0 / 14.0 112coding (919/1104 nt)pheAbifunctional chorismate mutase/prephenate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (55/48);  new base . (4/2);  total (60/52)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

Only 100 of 112 total aligned reads displayed.

GATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCG  >  NC_002947/1975885‑1976176
                                                                                                                                                  |                                                                                                                                                  
gaTCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGaa                                                                                                                                                >  1:486106/1‑150 (MQ=255)
 aTCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAAc                                                                                                                                               >  1:1638113/1‑150 (MQ=255)
     caGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCa                                                                                                                                           >  2:2971789/2‑150 (MQ=255)
     aaGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGATGGTGCCGTTCTATCAGAACGGCa                                                                                                                                           <  2:89617/150‑1 (MQ=255)
      aGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCAt                                                                                                                                          <  1:2975438/150‑1 (MQ=255)
         ccgccgGGACAAATCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGa                                                                                                                                       <  1:1679733/150‑1 (MQ=255)
         ccGCCCGGACAACTCCACTCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGATGCCGTTCTGATAGAACGGCaccag                                                                                                                                       <  1:615248/150‑5 (MQ=255)
         ccGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGACGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGa                                                                                                                                       >  1:941063/1‑150 (MQ=255)
          cGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGAc                                                                                                                                      <  2:857679/150‑1 (MQ=255)
            cccGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGaa                                                                                                                                                <  2:1169540/138‑1 (MQ=255)
            cccGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGaa                                                                                                                                                >  1:1169540/1‑138 (MQ=255)
            cccGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCt                                                                                                                                    <  2:2740299/150‑1 (MQ=255)
              cGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGa                                                                                                                                  <  1:78395/150‑1 (MQ=255)
                 aCAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGTCGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGAccc                                                                                                                               >  1:1070221/1‑150 (MQ=255)
                     cTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGAc                                                                                                                                      <  2:1494481/139‑1 (MQ=255)
                     cTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTAACGCCAACCGGAGCAGCCAAGACCTCGATTATCGTGTCAATGTGCAACAAGCCATGTGCTTTGCATGTGTTGCTGTTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCAt                                                                                                                           <  2:344563/150‑1 (MQ=255)
                      tCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCTACGTGACCCGCATc                                                                                                                          >  1:632666/1‑150 (MQ=255)
                      tCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGAcc                                                                                                                                >  2:854819/1‑144 (MQ=255)
                      tCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGAcc                                                                                                                                <  1:854819/144‑1 (MQ=255)
                      tCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATc                                                                                                                          >  1:350491/1‑150 (MQ=255)
                      tCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATc                                                                                                                          <  2:941063/150‑1 (MQ=255)
                      tCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATc                                                                                                                          <  1:2614845/150‑1 (MQ=255)
                        cACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTGATAGAACGGc                                                                                                                                            >  1:1681831/1‑130 (MQ=255)
                        cACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTGATAGAACGGc                                                                                                                                            <  2:1681831/130‑1 (MQ=255)
                                tCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGAc                                                                                                                     >  2:973246/1‑145 (MQ=255)
                                tCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGAc                                                                                                                     <  1:973246/145‑1 (MQ=255)
                                   tgatgaTCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGt                                                                                                             >  2:436417/1‑150 (MQ=255)
                                       gatCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTcgcg                                                                                                         >  2:627840/1‑150 (MQ=255)
                                               aCCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCa                                                                                                 >  1:735686/1‑150 (MQ=255)
                                                ccAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCaa                                                                                                <  1:436417/150‑1 (MQ=255)
                                                  aGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGAc                                                                                                                                      <  2:548179/110‑1 (MQ=255)
                                                  aGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGAc                                                                                                                                      >  1:548179/1‑110 (MQ=255)
                                                   ggAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTc                                                                                                            >  2:367767/1‑135 (MQ=255)
                                                   ggAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTc                                                                                                            <  1:367767/135‑1 (MQ=255)
                                                          ccgccgACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGc                                                                                                  <  1:89491/138‑1 (MQ=255)
                                                          ccgccgACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGc                                                                                                  >  2:89491/1‑138 (MQ=255)
                                                          ccgccgACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCACAGCGGCAAGTGGACCTAc                                                                                      <  2:2922255/150‑1 (MQ=255)
                                                           cgccgACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACCAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACg                                                                                     >  1:896037/1‑150 (MQ=255)
                                                           cgccgACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTAc                                                                                      >  2:656317/1‑149 (MQ=255)
                                                           cgccgACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTAc                                                                                      <  1:656317/149‑1 (MQ=255)
                                                             ccgACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACgtg                                                                                   <  1:1017809/150‑1 (MQ=255)
                                                                  cGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGttctt                                                                              >  2:776006/1‑150 (MQ=255)
                                                                      gcctgcAAGACCTTGATCATCGGTTGGATGACCAACAAGCCAGGTGCTTTGCATGTGTTGCTTGTGCCGTTCTATCAGAACGGCATCTACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTCGACCTGCGGGTTCTTCATc                                                                          <  2:212242/145‑1 (MQ=255)
                                                                      gacgacAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATc                                                                          >  2:2564540/1‑150 (MQ=255)
                                                                         gacAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGAc                                                                                                                                      <  2:21478/87‑1 (MQ=255)
                                                                         gacAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGAc                                                                                                                                      >  1:21478/1‑87 (MQ=255)
                                                                         gacAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCaa                                                                                                >  1:2992188/1‑125 (MQ=255)
                                                                         gacAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCaa                                                                                                <  2:2992188/125‑1 (MQ=255)
                                                                          acAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACt                                                                      <  1:627840/150‑1 (MQ=255)
                                                                               aCCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGttc                                                                                <  1:754808/135‑1 (MQ=255)
                                                                               aCCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGttc                                                                                >  2:754808/1‑135 (MQ=255)
                                                                               aCCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTtcgtc                                                                 <  2:1070221/150‑1 (MQ=255)
                                                                                  tCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGc                                                              >  2:2027012/1‑150 (MQ=255)
                                                                                              tCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGc                                                                                                                             <  2:311191/75‑1 (MQ=255)
                                                                                              tCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGc                                                                                                                             >  1:311191/1‑75 (MQ=255)
                                                                                                    aGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCg                                                                                                    <  1:2084538/94‑1 (MQ=255)
                                                                                                    aGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCg                                                                                                    >  2:2084538/1‑94 (MQ=255)
                                                                                                      caacaaGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGAt                                          >  2:512814/1‑150 (MQ=255)
                                                                                                      caacaaGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGAt                                          <  2:632666/150‑1 (MQ=255)
                                                                                                       aacaaGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTc                                                                                                            >  2:1464901/1‑83 (MQ=255)
                                                                                                       aacaaGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTc                                                                                                            <  1:1464901/83‑1 (MQ=255)
                                                                                                       aacaaGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATc                                         <  1:2027012/150‑1 (MQ=255)
                                                                                                       aacaaGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATc                                         >  1:1016381/1‑150 (MQ=255)
                                                                                                        acaaGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCTCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCTCTACCCGCTGATCa                                        >  1:2159551/1‑150 (MQ=255)
                                                                                                         caaGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCaa                                       <  2:1650548/150‑1 (MQ=255)
                                                                                                         caaGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCaa                                       <  2:1638113/150‑1 (MQ=255)
                                                                                                            gCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGc                                    >  2:1805605/1‑150 (MQ=255)
                                                                                                             ccAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTACCAGAACCCCAAAGCCCTCACCCTCACCGAGACCTATCCGTCTCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCg                                   <  2:2647220/150‑1 (MQ=255)
                                                                                                              cAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCAAGAACTGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCgg                                  <  1:2971789/150‑1 (MQ=255)
                                                                                                               aGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGt                                 <  2:1016381/150‑1 (MQ=255)
                                                                                                                  tGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGAc                                                                                                                     <  1:1771453/63‑1 (MQ=255)
                                                                                                                  tGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGAc                                                                                                                     >  2:1771453/1‑63 (MQ=255)
                                                                                                                     gTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGa                           >  2:551804/1‑150 (MQ=255)
                                                                                                                        gCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCCAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGcc                        >  2:867848/1‑149 (MQ=255)
                                                                                                                          aTGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCAt                                                                           <  2:495258/97‑1 (MQ=255)
                                                                                                                           tGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGAc                                                                       <  2:1101277/100‑1 (MQ=255)
                                                                                                                           tGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGAc                                                                       >  1:1101277/1‑100 (MQ=255)
                                                                                                                               ttGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGc                 >  2:220910/1‑150 (MQ=255)
                                                                                                                               ttGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGc                 >  1:2642700/1‑150 (MQ=255)
                                                                                                                                 gtcgatGCCGTTCTGATAGAACGGCACCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCa               >  2:1474534/6‑150 (MQ=255)
                                                                                                                                 gCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCa               <  1:512814/150‑1 (MQ=255)
                                                                                                                                    gatGCCGTTCTGATAGAACGGCACCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCagga            >  1:2577987/3‑150 (MQ=255)
                                                                                                                                     gTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGa                                                                                           <  2:1229424/70‑1 (MQ=255)
                                                                                                                                     gTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGa                                                                                           >  1:1229424/1‑70 (MQ=255)
                                                                                                                                     atGCCGTTCTGATAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAg                       >  2:2713829/2‑138 (MQ=255)
                                                                                                                                     atGCCGTTCTGATAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAg                       <  1:2713829/137‑1 (MQ=255)
                                                                                                                                      tGCCGTTCTGATAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCgg                                  >  2:1420063/1‑126 (MQ=255)
                                                                                                                                      tGCCGTTCTGATAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCgg                                  <  1:1420063/126‑1 (MQ=255)
                                                                                                                                      tGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCaccacc                                                       <  1:2155389/105‑1 (MQ=255)
                                                                                                                                      tGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCaccacc                                                       >  2:2155389/1‑105 (MQ=255)
                                                                                                                                       gCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGccacca                                                         <  2:1917260/102‑1 (MQ=255)
                                                                                                                                       gCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGccacca                                                         >  1:1917260/1‑102 (MQ=255)
                                                                                                                                         cGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCtg       >  2:2592425/1‑150 (MQ=255)
                                                                                                                                          gTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAg                          <  1:484595/130‑1 (MQ=255)
                                                                                                                                          gTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAg                          >  2:484595/1‑130 (MQ=255)
                                                                                                                                            tCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTg                                <  1:1191451/122‑1 (MQ=255)
                                                                                                                                            tCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTg                                >  2:1191451/1‑122 (MQ=255)
                                                                                                                                            tCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTgg    <  1:985025/150‑1 (MQ=255)
                                                                                                                                              tATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGgcg  >  1:1965041/1‑150 (MQ=255)
                                                                                                                                              tATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGgcg  >  1:2748691/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCG  >  NC_002947/1975885‑1976176

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: