Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,976,271 T→A 13.1% intergenic (+56/‑109) pheA → / → PP_1770 bifunctional chorismate mutase/prephenate dehydratase/3‑phosphoshikimate 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,976,2710TA13.1% 711.7 / 77.4 256intergenic (+56/‑109)pheA/PP_1770bifunctional chorismate mutase/prephenate dehydratase/3‑phosphoshikimate 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base T (127/94);  new base A (20/14);  total (148/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 261 total aligned reads displayed.

CCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGA  >  NC_002947/1976122‑1976420
                                                                                                                                                     |                                                                                                                                                     
ccGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGt                                                                                                                                                       <  2:2423523/150‑1 (MQ=255)
ccGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCACTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGt                                                                                                                                                       >  1:3677007/1‑150 (MQ=255)
  gcgACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGa                                                                                                                                                     <  1:8570543/150‑1 (MQ=255)
      cccGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCAt                                                                                                                                                 <  1:1342096/150‑1 (MQ=255)
      cccGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCAt                                                                                                                                                 <  1:8746195/150‑1 (MQ=255)
       ccGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCAtg                                                                                                                                                >  1:2127926/1‑150 (MQ=255)
       ccGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCAtg                                                                                                                                                <  2:1640620/150‑1 (MQ=255)
        cGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCAtgt                                                                                                                                               <  1:6254535/150‑1 (MQ=255)
          cTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGa                                                                                                                                             <  2:4888280/150‑1 (MQ=255)
           tGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGag                                                                                                                                            >  1:6912572/1‑149 (MQ=255)
             aTCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGc                                                                                                                                          <  1:1594553/150‑1 (MQ=255)
             aTCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGc                                                                                                                                          <  1:316984/150‑1 (MQ=255)
                 aGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCTAGCCCGCTCCCACAGAGAtcggcatcacatgct                                                                                                                                      <  1:6001669‑M1/150‑16 (MQ=255)
                  ggCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGAtc                                                                                                                                     >  2:7643275/1‑150 (MQ=255)
                       tGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCtgt                                                                                                                                <  2:3941993/150‑1 (MQ=255)
                        gCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTAGCCGATGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCtgtg                                                                                                                               <  1:8856349/150‑1 (MQ=255)
                         cTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGAGAtcggcatcacatgctcactgtgg                                                                                                                              <  1:2373324‑M1/150‑24 (MQ=255)
                           cGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGa                                                                                                                            <  2:2956530/150‑1 (MQ=255)
                             aGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGAGAtcggcatcacatgctcactgtgggagc                                                                                                                          <  1:3634486‑M1/150‑28 (MQ=255)
                              gCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGACCGCTCCCACAGTGAGCATGTGATGCCGATCTGTGTTGGAGCg                                                                                                                         >  2:2054400/1‑150 (MQ=255)
                              gCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCAGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCg                                                                                                                         >  2:8624607/1‑150 (MQ=255)
                               cAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCTAGCCCGCTCCCACAGAGAtcggcatcacatgctcactgtgggagcgg                                                                                                                        <  1:7976213‑M1/150‑30 (MQ=255)
                               cAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCTAGCCCGCTCCCACAGAGAtcggcatcacatgctcactgtgggagcgg                                                                                                                        <  2:6387719‑M1/150‑30 (MQ=255)
                               cAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGAGAtcggcatcacatgctcactgtgggagcgg                                                                                                                        <  1:4735437‑M1/150‑30 (MQ=255)
                                  aTCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCt                                                                                                                     >  2:2152053/1‑150 (MQ=255)
                                  aTCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCt                                                                                                                     >  1:1817051/1‑150 (MQ=255)
                                   tCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTa                                                                                                                    <  1:4236117/150‑1 (MQ=255)
                                   tCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTa                                                                                                                    <  1:5502846/150‑1 (MQ=255)
                                     aGCCAGGAGGCTGTGGCGCTGAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACCGGGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGcc                                                                                                                 >  2:1491442/1‑150 (MQ=255)
                                       ccAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGccc                                                                                                                >  2:3273325/1‑150 (MQ=255)
                                         aggaggCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCcgc                                                                                                              >  2:8116998/1‑150 (MQ=255)
                                         aggaggCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCcgc                                                                                                              >  2:3503750/1‑150 (MQ=255)
                                                 gtgGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGgct                                                                                                      <  1:5771580/150‑2 (MQ=255)
                                                   ggCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGgccgc                                                                                                    >  2:2665825/1‑150 (MQ=255)
                                                   ggCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGgccgc                                                                                                    >  2:7284831/1‑150 (MQ=255)
                                                    gcgcTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGgccgcg                                                                                                   >  1:8678085/1‑150 (MQ=255)
                                                     cgcTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTTTTTCAGGGGCCTTGGGGCTGCCTTGCGGCCGTTCGCGGGCTAGCCCGCTCCCACAGAGAtcggcatcacatgctcactgtgggagcgggcatgcccgcgaagggctgcaa                                                                                                  <  1:5622673‑M1/150‑52 (MQ=255)
                                                        tGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGgccgcgccgc                                                                                               <  2:5870387/150‑1 (MQ=255)
                                                           aGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCCGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGc                                                                                            <  2:214672/150‑1 (MQ=255)
                                                                atGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGAGAGCATGAGATGCCGATATCTGTGGGAGCGGGCTACCCAGCGAAGGGCCGCGCCGCCGCCCCaa                                                                                       >  1:1067538/2‑150 (MQ=39)
                                                                 tGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCCAt                                                                                      >  1:52582/1‑150 (MQ=255)
                                                                 tGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAAt                                                                                      >  1:6440440/1‑150 (MQ=255)
                                                                 tGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAAt                                                                                      >  2:7552980/1‑150 (MQ=255)
                                                                  gggCTCTTATCCGACGGCGGTGCTTTGATGCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGAGGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATg                                                                                     <  2:7038374/150‑1 (MQ=255)
                                                                    gCTCTTATCCGACGGCGGTGCTTTGATTCCGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGGGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGc                                                                                   <  1:4815337/150‑1 (MQ=255)
                                                                    gCTCTTATCCGACGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGc                                                                                   <  2:3159105/150‑1 (MQ=255)
                                                                      tctTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGccc                                                                                 >  2:5614172/1‑150 (MQ=255)
                                                                        ttATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTc                                                                               >  2:29012/1‑150 (MQ=255)
                                                                        ttATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGAGAtcggcatcacatgctcactgtgg                                                                                                                              <  1:2110457‑M1/103‑24 (MQ=255)
                                                                          aTCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTTCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTcac                                                                             >  1:6327677/1‑150 (MQ=255)
                                                                          aTCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTcac                                                                             <  2:1199666/150‑1 (MQ=255)
                                                                          aTCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTcac                                                                             >  1:8361078/1‑150 (MQ=255)
                                                                           tCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTcaca                                                                            <  2:8265141/150‑1 (MQ=255)
                                                                            ccGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAg                                                                           <  1:6599069/150‑1 (MQ=255)
                                                                              gAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGc                                                                         >  1:7220727/1‑150 (MQ=255)
                                                                                     gTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGATCAAGTGATGACGATCTCTGTGGGAGCGGGCAAGCCCCCGAAGGGCCCCGCCGCCGCCCCAAAGGCCCACAAAGGCTACCAAg                                                                  >  1:2523633/1‑150 (MQ=37)
                                                                                     gTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAg                                                                  >  2:829192/1‑150 (MQ=255)
                                                                                         tttGATTCCGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGAttt                                                              <  1:6499005/150‑1 (MQ=255)
                                                                                         tttGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGAttt                                                              >  2:5101247/1‑150 (MQ=255)
                                                                                         tttGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCGAGAttt                                                              >  2:62357/1‑150 (MQ=255)
                                                                                         tttGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAAAGGCCCACACACGCAACCGAGAttt                                                              >  1:5010623/1‑150 (MQ=39)
                                                                                            gATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGTGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTg                                                           <  2:7626512/150‑1 (MQ=255)
                                                                                            gATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTg                                                           <  2:1407948/150‑1 (MQ=255)
                                                                                            gATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCCCGCCGCGGCCCCACTGGCCCGCACAGGCTACCGAGATTTCTg                                                           >  1:1173811/1‑150 (MQ=255)
                                                                                              ttCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGaa                                                         >  2:4360742/1‑150 (MQ=255)
                                                                                              ttCAGGGGCCTTGGGGCTGCCTGGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGCGCTTGTGTTGCCGATCTCTGTGGGAGCGGGCTTGCCCGCGATGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGaa                                                         <  2:7551135/150‑1 (MQ=255)
                                                                                                  ggggCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAgg                                                     <  2:6923364/150‑1 (MQ=255)
                                                                                                  ggggCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAgg                                                     <  2:6251271/150‑1 (MQ=255)
                                                                                                   gggCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTATTAACAggg                                                    >  1:8610393/1‑150 (MQ=255)
                                                                                                     gCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGtt                                                  >  2:4066393/1‑150 (MQ=255)
                                                                                                       cTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCg                                                >  2:6766827/1‑150 (MQ=255)
                                                                                                        ttGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGGTCGc                                               >  1:6780212/1‑150 (MQ=255)
                                                                                                         tGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGCTGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCa                                              >  2:413468/1‑150 (MQ=255)
                                                                                                         aagggCTGCAAGGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCa                                              >  2:1852254/3‑150 (MQ=255)
                                                                                                               tGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTg                                        >  1:3765872/1‑150 (MQ=255)
                                                                                                               tGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTg                                        >  1:3516677/1‑150 (MQ=255)
                                                                                                                       aGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATgc                                >  2:6556717/1‑150 (MQ=255)
                                                                                                                               gcgGGCATGCCCGCTCCCACCGTGAGCCTGGGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCCATGGCCCTCACAGGCTACCTCGAGTTCTGAACAGGGTTCGCACCAGTGGTCAATgcagcccaac                        >  2:5444211/1‑145 (MQ=255)
                                                                                                                                       gCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGtt                                                  <  2:5226192/116‑1 (MQ=255)
                                                                                                                                         ccGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCcgc              >  2:8746426/1‑150 (MQ=255)
                                                                                                                                         ccGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCGCACAGGCTACCGAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCcgc              >  1:5043815/1‑150 (MQ=255)
                                                                                                                                         ccGCTCCCACAGAGAtcggcatcacatgctcactgtgggagcgggctagcccgcgaagggccgcgccgcggccccaatggccctcacaggctacctagatttctgaacagggttcgcaccagtggtaaatgcagcaaaaaccaaacccgc              >  1:7129527‑M1/1‑15 (MQ=255)
                                                                                                                                         ccGCTCCCACAGAGAtcggcatcacatgctcactgtgggagcgggctagcccgcgaagggccgcgccgcggccccaatggccctcacaggctacctagatttctgaacagggttcgcaccagtggtaaatgcagcaaaaaccaaacccgc              >  2:6084730‑M1/1‑15 (MQ=255)
                                                                                                                                         ccGCTCACACAGTGAGCATTTGATGACGATATCTGTGGGATCTGGCTAGACCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGCGTTCGCACCAGAGGAAAAAGCAGCAAAAACCAAAACcgc              >  2:7444665/1‑150 (MQ=255)
                                                                                                                                             tCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGccc          >  2:8995892/1‑150 (MQ=255)
                                                                                                                                              cccACAGAGAtcggcatcacatgctcactgtgggagcgggctagcccgcgaagggccgcgccgcggccccaatggccctcacaggctacctagatttctgaacagggttcgcaccagtggtaaatgcagcaaaaaccaaacccgcgccca         >  2:3730120‑M1/1‑10 (MQ=255)
                                                                                                                                               ccACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCcat        >  2:6859729/1‑150 (MQ=255)
                                                                                                                                               ccACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAAACAAACCCGCGCCcat        >  1:8938632/1‑150 (MQ=255)
                                                                                                                                               ccACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGGGGTAAATGCAGCAAAAACCAAACCCGCGCCcat        >  1:1389990/1‑150 (MQ=255)
                                                                                                                                               ccACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGGGGTAAATGCAGCAAAAAACAAACCCGCGCCcat        >  1:1920872/1‑150 (MQ=255)
                                                                                                                                               ccACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTAACTAGATTTGTTTACATGGTTCTCACCAGTGGTTAATGCATCAAAAAAAAAAAACGCGCCcat        >  1:48797/1‑150 (MQ=255)
                                                                                                                                               ccACAGAGAtcggcatcacatgctcactgtgggagcgggcatgcccgcgaagggccgcgccgcggccccaatggccctcacaggctacctagatttctgaacagggttcgcaccagtggtaaatgca                               <  1:3978360‑M1/127‑119 (MQ=255)
                                                                                                                                               ccACAGAGAtcggcatcacatgctcactgtgggagcgggcatgcccgcgaagggccgcgccgcggccccaatggccctcacaggctacctagatttctgaacagggttcgcaccagtggtaaatgca                               >  2:3978360‑M1/1‑9 (MQ=255)
                                                                                                                                                cacaGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCcatc       >  2:8419208/1‑150 (MQ=255)
                                                                                                                                                cacaGAGAtcggcatcacatgctcactgtgggagcgggctagcccgcgaagggccgcgccgcggccccaatggccctcacaggctacctagatttctgaacagggttcgcaccagtggtaaatgcagcaaaaaccaaacccgcgcccatc       >  2:6010711‑M1/1‑8 (MQ=255)
                                                                                                                                                cacaGAGAtcggcatcacatgctcactgtgggagcgggctagcccgcgaagggccgcgccgcggccccaatggccctcacaggctaccgagatttatgaacagggttcgcaccagtggtaaatgcagcaaaaaccaaacccgcgcccatc       >  1:2703022‑M1/1‑8 (MQ=255)
                                                                                                                                                 acaGAGAtcggcatcacatgctcactgtgggagcgggctagcccgcgaagggccgcgccgcggccccaatggccctcacaggctacctagatttctgaacagggttcgcaccagtggtaaatgcagcaaaaaccaaacccgcgcccatca      >  1:9094277‑M1/1‑7 (MQ=255)
                                                                                                                                                  caGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCcatcat     <  1:8926820/150‑1 (MQ=255)
                                                                                                                                                   aGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAAACAAACCCGCGCCCatcatc    >  1:6397790/1‑150 (MQ=255)
                                                                                                                                                     tGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGa  <  1:4047683/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGA  >  NC_002947/1976122‑1976420

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