Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 256,459 G→A 6.8% intergenic (+24/‑68) gpt → / → frsA xanthine phosphoribosyltransferase; xanthine‑guanine phosphoribosyltransferase/fermentation‑respiration switch protein; PTS Enzyme IIA(Glc)‑binding protein; pNP‑butyrate esterase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913256,4590GA 6.8% 1190.4 / 48.0 369intergenic (+24/‑68)gpt/frsAxanthine phosphoribosyltransferase; xanthine‑guanine phosphoribosyltransferase/fermentation‑respiration switch protein; PTS Enzyme IIA(Glc)‑binding protein; pNP‑butyrate esterase activity
Reads supporting (aligned to +/- strand):  ref base G (150/191);  new base A (15/10);  total (168/201)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 380 total aligned reads displayed.

ATCTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAGTCCGCCGCTTT  >  NC_000913/256310‑256607
                                                                                                                                                      |                                                                                                                                                    
 tCTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCgg                                                                                                                                                     >  2:4431065/1‑150 (MQ=255)
 tCTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCGAATCTTTTCAACGCCTGGCACTGCCgg                                                                                                                                                     >  2:5760774/1‑150 (MQ=255)
  cTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCggg                                                                                                                                                    <  2:1881410/150‑1 (MQ=255)
   ttCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGc                                                                                                                                                   <  1:270911/150‑1 (MQ=255)
   ttCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGc                                                                                                                                                   >  1:2571744/1‑150 (MQ=255)
    tCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCg                                                                                                                                                  >  2:1786374/1‑150 (MQ=255)
    tCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCg                                                                                                                                                  <  2:2142107/150‑1 (MQ=255)
     cGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgt                                                                                                                                                 <  2:4673987/150‑1 (MQ=255)
      gCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgtt                                                                                                                                                <  1:1280424/150‑1 (MQ=255)
      gCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgtt                                                                                                                                                <  1:2380898/150‑1 (MQ=255)
       cAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgttg                                                                                                                                               >  1:1181601/1‑150 (MQ=255)
       cAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgttg                                                                                                                                               <  2:5522234/150‑1 (MQ=255)
         aaaCCGGCTGGTCGTCCGCTGGTTGATGACTAAGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCTTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACTCCAGGCACTACCAGGAgttgtt                                                                                                                                             >  2:2655629/1‑150 (MQ=255)
            ccGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCtt                                                                                                                                          <  1:3121359/150‑1 (MQ=255)
            ccGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCtt                                                                                                                                          <  1:3555031/150‑1 (MQ=255)
             cGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCttt                                                                                                                                         <  1:3775044/150‑1 (MQ=255)
              ggCTGGTCGTCCGCTGGTTGATGACTCTGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCGAATCGTTTCAACGCCGGGCCCTGCCCGGCGGGGGTCtttt                                                                                                                                        >  1:668751/1‑150 (MQ=255)
                 tGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTaa                                                                                                                                     <  1:5232266/150‑1 (MQ=255)
                    tcgtcCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACtt                                                                                                                                  <  2:6509993/150‑1 (MQ=255)
                      gtcCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCa                                                                                                                                <  2:5449925/150‑1 (MQ=255)
                      gtcCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCa                                                                                                                                <  1:471585/150‑1 (MQ=255)
                       tcCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAg                                                                                                                               >  2:1996543/1‑150 (MQ=255)
                          gCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCg                                                                                                                            >  1:4362621/1‑150 (MQ=255)
                          gCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGACCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCg                                                                                                                            >  2:1726070/1‑150 (MQ=255)
                           cTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCgg                                                                                                                           <  2:3212799/150‑1 (MQ=255)
                           cTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCgg                                                                                                                           <  1:6504403/150‑1 (MQ=255)
                           cTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCAGGGCGTTGTTCTTTTTAACTTCAGGCgg                                                                                                                           <  1:4833219/150‑1 (MQ=255)
                             ggTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGt                                                                                                                         <  2:842716/150‑1 (MQ=255)
                                    gaCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATa                                                                                                                  <  1:2803398/150‑1 (MQ=255)
                                           ttgttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCa                                                                                                           <  1:2276828/150‑1 (MQ=255)
                                           ttgttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCa                                                                                                           <  1:5891625/150‑1 (MQ=255)
                                           ttgttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCa                                                                                                           <  1:5778668/150‑1 (MQ=255)
                                           ttggtgATATCCCGCACGATACCTGGATTGAACAGCCGTGGGATAGGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCa                                                                                                           <  1:33596/150‑1 (MQ=255)
                                             gttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagt                                                                                                         >  2:974834/1‑150 (MQ=255)
                                             gttgATATCCCGCAAGATACCTGGAGTGAACTGCCGTGGGATCTGGGCGTCGTA‑TTCGTCCCGCCAGTCTCCGGTCGCTCATCTTTTCAACGCCGGGCACGGCCGGGCGGTGGTGTTGTTAGCTTCAGGCGGGTTACACTCGGTTcccgg                                                                                                         >  2:4092701/1‑147 (MQ=255)
                                                 atatCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaagt                                                                                                     <  1:572709/150‑1 (MQ=255)
                                                  tatCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCaggaagta                                                                                                    >  2:2556893/1‑150 (MQ=255)
                                                      ccGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCt                                                                                                <  1:1021406/150‑1 (MQ=255)
                                                       cGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTATTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTg                                                                                               >  1:1616034/1‑150 (MQ=255)
                                                       cGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCGCATCGTTTCAACGCCGGGCCCAGCCGGGCGGTGTTCTTTTTAACGTCAGGCGGGGTACAATAGTGTCCagtaaggaggcgg                                                                                               >  1:6209812/1‑143 (MQ=255)
                                                         cAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGa                                                                                             <  1:4997757/150‑1 (MQ=255)
                                                          aaGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAg                                                                                            <  1:2659593/150‑1 (MQ=255)
                                                           aGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAgg                                                                                           <  1:290131/150‑1 (MQ=255)
                                                                ccTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCa                                                                                      <  1:4777637/150‑1 (MQ=255)
                                                                 cTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCAt                                                                                     >  2:3374092/1‑150 (MQ=255)
                                                                  tGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGGTGTTCTTTTTAAATTAAGGCGGGGTAAAAAAGTTTCCCGTAAGGATTCTGGGGGCTGCCTc                                                                                    >  2:6234081/1‑150 (MQ=255)
                                                                    gATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCGAATCGTTTCAACGCCTGGCACTGCCGGGCGTGGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGGAGTCTGGAGGCGGCATCCa                                                                                  >  2:5245623/1‑150 (MQ=255)
                                                                      ttGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATATTTTCAAAGCCAGGCACAGCCAGGCGGTGTTATTTTTTAATAAAAGCGGGTTACAAAAAATTACCGTAACTATTAAGTAGGCGGCAtgccat                                                                                >  1:5072626/1‑145 (MQ=255)
                                                                       tGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGa                                                                               <  1:1865760/150‑1 (MQ=255)
                                                                       tGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGa                                                                               <  1:4454586/150‑1 (MQ=255)
                                                                         taCAGCCGTGGGATATGGGAGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCATCGCCTTGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTGTGGAGGCTGCATCCATGaca                                                                             >  2:1002352/2‑150 (MQ=255)
                                                                         aaCAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGaca                                                                             >  2:3465558/1‑150 (MQ=255)
                                                                           cAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGAcaca                                                                           <  2:3976873/150‑1 (MQ=255)
                                                                            aGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAg                                                                          >  2:5019413/1‑150 (MQ=255)
                                                                            aGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCATTTTAACTTCAGGCGGGGTACCATGGGTTCCAGTTAGTAGGGGGGGGGCTGCATCCATGCCCCAg                                                                          >  2:643593/1‑150 (MQ=255)
                                                                             gCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAgg                                                                         <  1:4431065/150‑1 (MQ=255)
                                                                                  gggATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAAc                                                                    <  2:1065290/150‑1 (MQ=255)
                                                                                        tGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGc                                                              <  1:974834/150‑1 (MQ=255)
                                                                                             gtcgtATTTCGTCCCGCCAATCTCCGGTCGCTAATCATTTCAACGCCAGGCACTGCCGGGCGGGGTTAATTTTAAATTCAGGCCGGTTACACAAGTTTACAGTAATTATTCTGGAGGCGGCCAACATGACACAGGCAAACCAGAGCgcca                                                          >  2:3114678/1‑147 (MQ=255)
                                                                                                gtA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCt                                                      <  1:4592942/150‑1 (MQ=255)
                                                                                                gtA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCt                                                      <  2:5709182/150‑1 (MQ=255)
                                                                                                    ttCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGtt                                                   <  1:6391019/150‑1 (MQ=255)
                                                                                                     tCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTATTTTTAAAATAAGGCGGGTTACAATAGTTTCCAGTAAGTATTCGGGGGGCGGCATACAAGACACAGGCAAAACAGATCGAAACCCAGtta                                                  >  2:2849161/1‑149 (MQ=255)
                                                                                                       gTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCaa                                                >  2:269999/1‑150 (MQ=255)
                                                                                                          ccGCCAATCTCCGGTCGCTAATCTTTTCAAAGCCTGGCACTGCCGGGCGTTGTTATTTTTAACTTAATGCGGGTTAAAATATTTTCGAGTTATTATTATTGTGGCTTGATACATTAAACAGGCAAAACTTTGCGAAACCCTGTTaaaaac                                             >  2:2297363/1‑148 (MQ=255)
                                                                                                            gCCAATCTCCGGTCGCTAATCGTTTCAACGCCTGGCACTGCCGGGCGTTGTTATTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACACGCAAACCTGAGCGAAACCCTGTTCAAAcccc                                           >  2:296235/1‑150 (MQ=255)
                                                                                                             ccAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATTACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCg                                          >  2:5873906/1‑150 (MQ=255)
                                                                                                                 tctcCGGTCGCTCAGCTTTTCAACGCCTGGCACTGCCGGGCGGTGTTCTTTTTAACTGCAGGCGGGTTAAAATAGTTTCCCGTAAGTGTTCTGGAGGCTGCATCCATGACACAGGCAAACCTGGGCGAAACCCTGTTCAAACCCCGCttt                                      >  2:147589/1‑150 (MQ=255)
                                                                                                                   tcCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTaa                                    <  2:4366940/150‑1 (MQ=255)
                                                                                                                      ggTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACa                                 >  2:4165871/1‑150 (MQ=255)
                                                                                                                        tctcTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTCTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATc                               <  1:6052866/147‑1 (MQ=255)
                                                                                                                        tCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATc                               <  1:407471/150‑1 (MQ=255)
                                                                                                                        tCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATc                               <  2:681485/150‑1 (MQ=255)
                                                                                                                        tCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATc                               <  2:2430684/150‑1 (MQ=255)
                                                                                                                          gCTAATCTTTTCATCGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCt                             <  2:5349274/150‑1 (MQ=255)
                                                                                                                          gCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATccg                             >  2:6587039/1‑149 (MQ=255)
                                                                                                                          gCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCt                             <  1:6326426/150‑1 (MQ=255)
                                                                                                                            tAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGa                           >  2:5671553/1‑150 (MQ=255)
                                                                                                                            tAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGa                           <  2:369754/150‑1 (MQ=255)
                                                                                                                            tAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGa                           >  2:1238582/1‑150 (MQ=255)
                                                                                                                              aTCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGaaa                         <  2:4463160/150‑1 (MQ=255)
                                                                                                                               tCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAAc                        >  2:1109233/1‑150 (MQ=255)
                                                                                                                                cTTTTCAACGCCTGGCACTGCCGGGCGTTGTTGTTTTTAAATTAAGGCGGGTTACAAAAGTTTGCAGTAAGTATTCTGGAGGCTGCATCCAAGACACAGGCAAAACAGAGCGAAACCCAGTTCAAAACCCCCATTAAAAAACCAGAAAcc                       >  1:5316789/1‑150 (MQ=255)
                                                                                                                                cTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAAcc                       <  2:4922949/150‑1 (MQ=255)
                                                                                                                                 ttttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCt                      <  1:3855251/150‑1 (MQ=255)
                                                                                                                                  tttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTc                     <  1:4211649/150‑1 (MQ=255)
                                                                                                                                  tttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTc                     <  2:6115539/150‑1 (MQ=255)
                                                                                                                                     cAACGCCCGGCAGTGCCAGGCGTTGTTATTTTTAACTTAAGGCGGGTTACAAAAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCGGTTAAAACCCCGCTTTAAACAAACAGAAACCACGAc                  >  1:1800612/1‑150 (MQ=255)
                                                                                                                                      aaCGCCTGGCACTGCCGGGCGTTGTTGTTTTTTAATTCAGGCGGGTTACAATAGGTTACAGTAAGTATTAGGGAGGCTGCATCCAAGAAACACGCAAAACGGAGCGAAACCCAGTTCAAACCCCCCGTTAAAAAGCCGGAAACCacgacg                 >  1:3243142/1‑150 (MQ=255)
                                                                                                                                      aaCGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACTGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACg                 >  2:4749831/1‑150 (MQ=255)
                                                                                                                                      aaCGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAATTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACg                 >  2:4583436/1‑150 (MQ=255)
                                                                                                                                       aCGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGc                <  1:770033/150‑1 (MQ=255)
                                                                                                                                       aCGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGc                <  2:3227200/150‑1 (MQ=255)
                                                                                                                                         gCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTa              >  2:3653477/1‑150 (MQ=255)
                                                                                                                                         gCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTa              <  2:3851407/150‑1 (MQ=255)
                                                                                                                                           ccggCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAGt            >  2:4725471/3‑150 (MQ=255)
                                                                                                                                              gCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAGTccg         >  2:276328/1‑150 (MQ=255)
                                                                                                                                               cACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAGTccgc        <  1:3186722/150‑1 (MQ=255)
                                                                                                                                                aCTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAGTccgcc       <  1:3623637/150‑1 (MQ=255)
                                                                                                                                                     cGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAGTCCGCCGCttt  <  1:141048/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
ATCTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTA‑TTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAGTCCGCCGCTTT  >  NC_000913/256310‑256607

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: ‑