Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,346,450 A→G 20.1% intergenic (+55/‑123) pgi‑II → / → acsA‑II glucose 6‑phosphate isomerase/acetyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,346,4500AG20.1% 148.8 / 29.7 73intergenic (+55/‑123)pgi‑II/acsA‑IIglucose 6‑phosphate isomerase/acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base A (34/21);  new base G (7/7);  total (42/31)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.35e-01

TACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATACCCCC  >  NC_002947/5346303‑5346597
                                                                                                                                                      |                                                                                                                                                    
tACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTAcc                                                                                                                                                    <  1:144039/150‑1 (MQ=255)
     gcgcTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTa                                                                                                                                               <  2:5304752/150‑1 (MQ=255)
            gTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCgg                                                                                                                                        <  1:3448399/150‑1 (MQ=255)
              cGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAggcgggt                                                                                                                                        <  1:112124/150‑6 (MQ=25)
              cGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAggcgggt                                                                                                                                        <  1:3492351/150‑6 (MQ=25)
              cGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGt                                                                                                                                      <  2:4575096/150‑1 (MQ=255)
               ggTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCAG‑‑CCTACCAGGCAACGCGTTACCTCGGAGGAGCGGGtt                                                                                                                                     >  2:4183504/1‑150 (MQ=255)
                gTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCgg                                                                                                                                        >  2:2171408/1‑146 (MQ=255)
                gTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCgg                                                                                                                                        <  1:2171408/146‑1 (MQ=255)
                 tAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTa                                                                                                                                   <  2:1979654/150‑1 (MQ=25)
                    ccTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTAccc                                                                                                                                <  2:2217872/150‑1 (MQ=21)
                     cTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTc                                                                                                                                                  >  1:2761463/1‑131 (MQ=255)
                     cTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTc                                                                                                                                                  <  2:2761463/131‑1 (MQ=255)
                      tGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCcgc                                                                                                                              <  1:848547/150‑1 (MQ=255)
                      tGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCcgc                                                                                                                              <  1:4587956/150‑1 (MQ=255)
                         aaGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGaa                                                                                                                           <  1:5144359/150‑1 (MQ=255)
                           gACAGCGCCTAAGACGGTTCCACCCACGGCCTGATCAACTACTTCCGCGGCCGTCAGCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGTTAACGCGTTGCCTGGTAGGAGCGGGTTTACCCGCGaa                                                                                                                           <  2:3375274/150‑1 (MQ=21)
                           gACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACGTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAaga                                                                                                                         <  2:83669/150‑1 (MQ=255)
                            aCAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCACTGGGCGGGGGAAACCCC‑‑CCACCCAGCCACCCCGTTACCACGAGGGGGGGGGTTTGCCCCCGacgg                                                                                                                         >  1:1891544/1‑147 (MQ=19)
                             cAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGTGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAgg                                                                                                                       <  1:1677516/150‑1 (MQ=255)
                                    gAAGCCGGTTCCACCCATGGCCTGATCAACTAGGGCGGGGGCCGTCACCGGGGTTGATCCGTGTTTGTGCCGGCCTCTGCGGGGGTAAACCCG‑‑CCTACCAGGCAAGGGGTTACCTGGTTGGTGGGGGTTTCCCCGCGAAGAGGCCGGccc                                                                                                               <  2:4832859/150‑1 (MQ=255)
                                    gAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCGCGTCGCGGGGGAACCCG‑‑CCGACCCGGCAACGCGTTACCTCGTGGGGGCGGGTTTTCCCGCGGAGGGGGCGGccc                                                                                                               >  1:4823307/1‑150 (MQ=255)
                                     aaGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACCCGTTACCAAGTAGGAGCCGGTTTACCCCCGAAGAGGCCGGCCCa                                                                                                              >  1:3572136/1‑150 (MQ=39)
                                         cGGTTCCACCCATGGCATGATCAACTACTTCCGCGGCCGTCACCGGGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGAAACGCGTTACCTCGTAGGAGCGGTTTTACCCGCGAAGTGGCCGGCCCAGGtt                                                                                                          <  2:4979188/150‑1 (MQ=38)
                                         cGGTTCCACCCACGGCCTTATCAACTACTTCCGCGGCCGTCACGGCGGTTGATCCGAGCTTGTGCCGGCCTCTGCGCGGGTAAACCCG‑‑CCTACCCGGCAACGCGTTCCCTCGTAGGTGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGtt                                                                                                          <  2:5444794/150‑1 (MQ=38)
                                          ggTTCCACCCAAGGCCTTATCAACTACTTCCGCGGCCGTCACCGCGTTTGATCCGAGCTTGTGCCGGCATCTTCGCGGTTAAACCCG‑‑CCTACCAGGAAACGCTTTACCTCGTAGGAGCGGTTTTACCCGCGAAGAGGCCGGCCCAGGTTg                                                                                                         <  2:2244991/150‑1 (MQ=38)
                                                  ccAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGGAAAACCG‑‑CCTACCAGGCCACCCCGTACCACGGAGGAGCGGGGTTTACCGCGAAGAGGCCCGCCCACGTTGAAAAcccc                                                                                                 >  1:886098/1‑148 (MQ=255)
                                                          tGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCgg                                                                                                                                        >  1:4195297/1‑104 (MQ=255)
                                                          tGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCgg                                                                                                                                        <  2:4195297/104‑1 (MQ=255)
                                                          gtaTCAACTCCTCCCGCGGCCGTCACCGGGGTTTATCCGCGCTTGGCCCGGCCACGTCGCGGTAAAACCCG‑‑CCAACCCGGCAACGGGTTACCTCGTAGGCGCGGTTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCCCCCACCt                                                                                         <  2:3572136/148‑1 (MQ=255)
                                                           gATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAggcgggt                                                                                                                                        <  2:351423/105‑6 (MQ=17)
                                                           gATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAggcgggt                                                                                                                                        >  1:351423/1‑100 (MQ=17)
                                                           gATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAAACCC‑‑CCGACCAGGCAACGCGGTTACACGTAGGCGCGGGGTTACCCCCGGAGAGGCCGGCCCAGGGTGAACACCGCGCCCAGcga                                                                                        >  1:2221786/1‑148 (MQ=255)
                                                               aaCTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACt                                                                                    >  2:2233920/1‑150 (MQ=255)
                                                                              cGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGGTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACAGTCTTATCGAATAGc                                                                     >  2:3623106/1‑150 (MQ=255)
                                                                              cGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCGTACCGGGTTAACCCC‑‑CCAACCAGGCCACCCCGTACCCCGTAGGGGCGGGGTTTCCCCCCAAGAGGCCCGGCCCGGGTTAACCCCGCTCCCACCGACACAGTCTTAAAGAATAGc                                                                     >  2:2959779/1‑150 (MQ=255)
                                                                                      cgGTTGATCCGAGCTTGTGCCGGCCTCGTCGCGGGTAAAACCG‑‑CCTACCAGGCAACGCGTTACCTCGGAGGACCGGGGTTACCCCCGAAGAGGCCGGCCCAGGGTGAACACCGCACCCAGCTACACAGTCTTATCGAAGAGCCGGTGCAg                                                             >  2:4475626/1‑150 (MQ=255)
                                                                                             tCCGAGCTTGTGCCGGCCTCGTCGCGGGGAAAACCC‑‑CCGAACACGCCACCCGGTACCACCGAGGGGCGGGTTTAACCCCGCACAGGCCCGCCCCCGGTGTACACCCCACCCCCCAACACAGGCTTATCGAAAAGCCGGTGCAGTACaacc                                                      >  2:696842/1‑147 (MQ=255)
                                                                                                    ttGTGCCGGCCTCTTAGCGGGTAAACCCC‑‑CCAACCAGGCAACGCCTTAACTCGTAGGAGCCGGTTTACCCCCGAAAAGGCCCGCCCAGGTTGAACACCGCTCCCAGCTACACTGTATTAACGAAAAACCGTTGCAGTCCATCGCCCCGTa                                               >  2:1291994/1‑150 (MQ=255)
                                                                                                       ggCCGGCGTCTTCGGGGGTAACCCCG‑‑CCTACCTGGCACGGGGTTCCCTCGTTGGCGGGGGTTTACCGGGGACGGGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACaa                                            <  2:3780820/149‑1 (MQ=255)
                                                                                                            gCCTCTGCGGGGTAAAACCCG‑‑CATACCTGGAAACGCTTTACCTCTTTGTTGGGTTTTTACCCGCTAAGTGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCa                                       <  2:1424298/150‑1 (MQ=255)
                                                                                                             ccTCTTCGGGGGTAAACCCG‑‑CCTACCAGGCAAGGCGTTACGTGGTTGGCGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAg                                      <  2:3255662/150‑1 (MQ=39)
                                                                                                               tctTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAggc                                    >  2:3175721/1‑149 (MQ=255)
                                                                                                                 ttCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGacc                                  >  2:2637770/1‑150 (MQ=255)
                                                                                                                   cgcgGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGaccac                                >  2:478230/1‑150 (MQ=255)
                                                                                                                       ggttAACCCCG‑‑CCTACCCGGCACCGCGTTACGTCGTCGGCGGGGTTTTCCCCGCTATGCGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCg                             <  1:2233920/147‑1 (MQ=255)
                                                                                                                       gggTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCg                             >  1:1368915/1‑150 (MQ=255)
                                                                                                                        ggTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGc                            >  2:3348225/1‑150 (MQ=255)
                                                                                                                         gTAAACCCGCTCCTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGcc                           >  1:4332888/1‑150 (MQ=21)
                                                                                                                         gTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGcc                           >  2:1012599/1‑150 (MQ=255)
                                                                                                                         ctAAACCCG‑‑CCACCCCGGAAACGGTTTACCTCGTTGGAGCGGGTTTACCGGGTATGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGcc                           <  2:4592796/149‑1 (MQ=255)
                                                                                                                             aCCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGt                       >  1:1074599/1‑150 (MQ=255)
                                                                                                                             aCCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGt                       >  1:4016776/1‑150 (MQ=255)
                                                                                                                              cccgctccTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGGGCAGGACCACCCGCCATGtt                      >  1:2888797/7‑150 (MQ=25)
                                                                                                                               ccgctccTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCACAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                     >  2:5310536/6‑150 (MQ=25)
                                                                                                                               ccgctccTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCACAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                     >  2:2437345/6‑150 (MQ=25)
                                                                                                                               ccgctccTACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCACAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                     >  1:5385930/6‑150 (MQ=25)
                                                                                                                               ccG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                     >  1:3693510/1‑150 (MQ=255)
                                                                                                                               ccG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                     >  2:162846/1‑150 (MQ=255)
                                                                                                                               ccG‑‑CCAACCCGGCACCGCGTTCCCGCTTTGGCGCGGGTTCACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTc                     <  2:158979/150‑1 (MQ=255)
                                                                                                                                 g‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGa                   >  2:4991957/1‑150 (MQ=255)
                                                                                                                                      tACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTAc                                              >  1:2460264/1‑118 (MQ=25)
                                                                                                                                      tACGAGGTAACGCGTTGCCTGGTAG‑‑GCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTAc                                              <  2:2460264/118‑1 (MQ=25)
                                                                                                                                       aCCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATc               >  2:411124/1‑150 (MQ=255)
                                                                                                                                         cAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCg             >  1:4805984/1‑150 (MQ=255)
                                                                                                                                           ggCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCa           >  2:768766/1‑150 (MQ=255)
                                                                                                                                           ggCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCa           >  2:768793/1‑150 (MQ=255)
                                                                                                                                           ggCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCa           >  2:768797/1‑150 (MQ=255)
                                                                                                                                               aCGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATa       >  2:4027394/1‑150 (MQ=255)
                                                                                                                                                cgcgTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATAc      >  1:3458223/1‑150 (MQ=255)
                                                                                                                                                 gcgTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATAcc     >  1:3702099/1‑150 (MQ=255)
                                                                                                                                                  cgTTACCTCGTCGGTGGGGGTTTACCCGCGACGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATAccc    <  2:4823307/150‑1 (MQ=255)
                                                                                                                                                    ttACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATAccccc  >  2:2128269/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAACTACTTCCGCGGCCGTCACCGCGGTTGATCCGAGCTTGTGCCGGCCTCTTCGCGGGTAAACCCG‑‑CCTACCAGGCAACGCGTTACCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGTTGAACACCGCTCCCAGCTACACTGTCTTATCGAATAGCCGTTGCAGTCCATCGCCCCTTACAAGAGCAGGACCACCCGCCATGTTCGATATCCGCAAATACCCCC  >  NC_002947/5346303‑5346597

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: