Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,839,384	T→A	intergenic (‑120/‑140)	ascG ← / → ascF	asc operon transcriptional repressor, prpBC operon repressor/cellobiose/arbutin/salicin‑specific PTS enzymes, IIB and IC components

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,839,384	0	T	A	100.0%	168.8 / NA	52	intergenic (‑120/‑140)	ascG/ascF	asc operon transcriptional repressor, prpBC operon repressor/cellobiose/arbutin/salicin‑specific PTS enzymes, IIB and IC components
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (25/27);  total (25/27)

TCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGATAACAGGATAA  >  NC_000913/2839255‑2839517                                                                                                                                  |                                                                                                                                      tcgtcgTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCActctct                                                                                                                              <  1:146268/139‑1 (MQ=255) tcgtcgTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCActctct                                                                                                                              <  2:352107/139‑1 (MQ=255)   gtcgtcATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTg                                                                                                                            >  1:794965/1‑139 (MQ=255)   gtcgtcATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTg                                                                                                                            >  1:33998/1‑139 (MQ=255)    tcgtcATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGc                                                                                                                           >  2:544672/1‑139 (MQ=255)     cgtcATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCt                                                                                                                          >  1:760022/1‑139 (MQ=255)         atcatTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTAtt                                                                                                                      >  2:1866/1‑139 (MQ=255)                tCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGAc                                                                                                               <  2:317931/139‑1 (MQ=255)                  aTCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACaa                                                                                                             <  2:532067/139‑1 (MQ=255)                  aTCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACaa                                                                                                             <  2:672733/139‑1 (MQ=255)                        ttACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTc                                                                                                       >  1:638551/1‑139 (MQ=255)                                          gtgtACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCgg                                                                                     <  2:760022/139‑1 (MQ=255)                                          gtgtACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCgg                                                                                     >  1:867598/1‑139 (MQ=255)                                          gtgtACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATCATTTTGTGAAACCgg                                                                                     <  2:705165/139‑1 (MQ=255)                                                     cGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCa                                                                                                                                    <  2:585660/81‑1 (MQ=255)                                                     cGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCa                                                                                                                                    >  1:585660/1‑81 (MQ=255)                                                     cGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTc                                                                          <  1:1866/139‑1 (MQ=255)                                                     cGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTc                                                                          <  1:66330/139‑1 (MQ=255)                                                       tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAGTGAACAAATCACTCTCTTGc                                                                                                                           <  2:621253/88‑1 (MQ=255)                                                       tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATc                                                                                                                                     >  1:1037818/1‑78 (MQ=255)                                                       tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATc                                                                                                                                     <  2:1037818/78‑1 (MQ=255)                                                       tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCAct                                                                                                                                  >  2:179248/1‑81 (MQ=255)                                                       tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCAct                                                                                                                                  <  2:816394/81‑1 (MQ=255)                                                       tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCAct                                                                                                                                  >  1:816394/1‑81 (MQ=255)                                                       tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCAct                                                                                                                                  <  1:179248/81‑1 (MQ=255)                                                       tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGc                                                                                                                           >  1:621253/1‑88 (MQ=255)                                                       tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCg                                                                        >  2:941860/1‑139 (MQ=255)                                                       tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCg                                                                        >  2:243870/1‑139 (MQ=255)                                                       tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCg                                                                        >  1:519415/1‑139 (MQ=255)                                                       tttGCCGATTGTCCTTGCACAATCGGCGGAAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACt                                                                                                           <  1:200073/104‑1 (MQ=255)                                                       tttGCCGATTGTCCTTGCACAATCGGCGGAAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACt                                                                                                           >  2:200073/1‑104 (MQ=255)                                                                             tCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCg                                                  >  1:176732/1‑139 (MQ=255)                                                                             tCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCg                                                  >  2:492788/1‑139 (MQ=255)                                                                              cggcggGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGt                                                 <  2:781975/139‑1 (MQ=255)                                                                                gcggGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCgg                                                                                     <  2:453899/101‑1 (MQ=255)                                                                                gcggGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCgg                                                                                     >  1:453899/1‑101 (MQ=255)                                                                                   ggAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGt                                            <  1:941860/139‑1 (MQ=255)                                                                                      aaaaTATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGa                                         <  1:29599/139‑1 (MQ=255)                                                                                       aaaTATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGAc                                        <  2:867598/139‑1 (MQ=255)                                                                                           atTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGAt                                    <  2:33998/139‑1 (MQ=255)                                                                                           atTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGAt                                    <  2:519415/139‑1 (MQ=255)                                                                                           atTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGAt                                    <  2:619342/139‑1 (MQ=255)                                                                                               aGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATaaca                                <  1:997753/139‑1 (MQ=255)                                                                                                         tttCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGa                                         <  1:264687/120‑1 (MQ=255)                                                                                                         tttCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGa                                         >  2:264687/1‑120 (MQ=255)                                                                                                         tttCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTa                      >  1:313240/1‑139 (MQ=255)                                                                                                         tttCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTa                      <  2:794965/139‑1 (MQ=255)                                                                                                           tCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTAcc                    >  2:498532/1‑139 (MQ=255)                                                                                                                   ataAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTcc                                                   >  1:375701/1‑100 (MQ=255)                                                                                                                   ataAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTcc                                                   <  2:375701/100‑1 (MQ=255)                                                                                                                            gAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGATAACAGGATa   >  1:455406/1‑139 (MQ=255)                                                                                                                             aaCAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGATAACAGGATaa  >  1:969367/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                      TCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGATAACAGGATAA  >  NC_000913/2839255‑2839517

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑