Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,561,929	T→C	intergenic (‑106/‑84)	glpE ← / → glpD	thiosulfate sulfurtransferase GlpE/aerobic glycerol 3‑phosphate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,561,929	0	T	C	100.0%	124.4 / NA	35	intergenic (‑106/‑84)	glpE/glpD	thiosulfate sulfurtransferase GlpE/aerobic glycerol 3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/13);  total (22/13)

TACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAG  >  NC_000913/3561859‑3562006                                                                       |                                                                              tACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTCCAATGTTAcc                                                                     >  1:677136/1‑81 (MQ=255)  aCCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTCCAATGTTAcc                                                                     >  1:677131/1‑80 (MQ=255)    cGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTCCAATGTTACCTaa                                                                  <  2:677096/81‑1 (MQ=255)     gCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTCCAATGTTACCTa                                                                   >  2:677138/1‑79 (MQ=255)     gCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTCCAATGTTACCTa                                                                   >  1:677137/1‑79 (MQ=255)             tACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGAtt                                                         >  2:677144/1‑81 (MQ=255)              aCGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTc                                                        >  2:677133/1‑81 (MQ=255)               cGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCt                                                       <  2:677097/81‑1 (MQ=255)                gTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCtt                                                      >  2:677139/1‑81 (MQ=255)                 tCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCttt                                                     <  1:677098/81‑1 (MQ=255)                   gAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGc                                                   >  1:677149/1‑81 (MQ=255)                   gAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGc                                                   >  2:677156/1‑81 (MQ=255)                     aGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCt                                                  >  2:677162/1‑80 (MQ=255)                             aTAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTc                                         <  2:677099/81‑1 (MQ=255)                                aaCGCCATAATGTTATACATATCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGAt                                      <  2:677100/81‑1 (MQ=255)                                   gCCATAATGTTATACATATCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAAc                                   >  1:677129/1‑81 (MQ=255)                                     cATAATGTTATACATATCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGa                                 <  2:677101/81‑1 (MQ=255)                                      aTAATGTTATACATATCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCTAGTATGTTCGATAACGaa                                >  2:677152/1‑81 (MQ=255)                                      aTAATGTTATACATATCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGaa                                >  2:677150/1‑81 (MQ=255)                                         aTGTTATACATATCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAt                             <  2:677102/81‑1 (MQ=255)                                             tataCATATCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTAt                         >  1:677160/1‑81 (MQ=255)                                              ataCATATCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTAt                         >  2:677145/1‑80 (MQ=255)                                               taCATATCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGa                       <  2:677103/81‑1 (MQ=255)                                                 cATATCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAg                      <  2:677104/80‑1 (MQ=255)                                                  atatCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCt                    >  2:677146/1‑81 (MQ=255)                                                  atatCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCt                    <  1:677105/81‑1 (MQ=255)                                                     tCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCttt                  >  1:677159/1‑80 (MQ=255)                                                     tCACTCTAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTa                 <  1:677106/81‑1 (MQ=255)                                                           tAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGa             >  1:677171/1‑79 (MQ=255)                                                           tAAAATGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGa             >  2:677153/1‑79 (MQ=255)                                                            aaaaTGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAg          >  2:677158/1‑81 (MQ=255)                                                             aaaTGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGaaa           <  1:677108/79‑1 (MQ=255)                                                             aaaTGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGt         >  1:677172/1‑81 (MQ=255)                                                             aaaTGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGt         <  1:677107/81‑1 (MQ=255)                                                              aaTGTTTTTCCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGt         >  1:677142/1‑80 (MQ=255)                                                                    tttcCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAg  >  2:677181/1‑81 (MQ=255)                                                                       |                                                                              TACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAG  >  NC_000913/3561859‑3562006

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑