Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,750,092	G→A	*287* (TAG→TAA)  E3K (GAG→AAG)	sgbU → sgbE →	putative L‑xylulose 5‑phosphate 3‑epimerase L‑ribulose‑5‑phosphate 4‑epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,750,092	0	G	A	100.0%	58.0 / NA	21	*287* (TAG→TAA)  E3K (GAG→AAG)	sgbU sgbE	putative L‑xylulose 5‑phosphate 3‑epimerase L‑ribulose‑5‑phosphate 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/10);  total (11/10)

CTATCGCGGTTCTTTCCTGATTGAGATGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGA  >  NC_000913/3749972‑3750225                                                                                                                         |                                                                                                                                      cTATCGCGGTTCTTTCCTGATTGAGATGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACg                                                                                                                     >  1:47351/1‑139 (MQ=255)       cgGTTCTTTCCTGATTGAGATGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTgg                                                                                                               >  1:830896/1‑139 (MQ=255)                              ggACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCt                                                                                        >  2:318160/1‑139 (MQ=255)                                aCCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCa                                                                                      <  2:47351/139‑1 (MQ=255)                                             aaGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACg                                                                         <  2:830896/139‑1 (MQ=255)                                                   cGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCAcc                                                                   <  1:477325/139‑1 (MQ=255)                                                                tattCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGaaa                                                                                                                            >  1:16032/1‑69 (MQ=255)                                                                tattCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGaaa                                                                                                                            <  2:16032/69‑1 (MQ=255)                                                                tattCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCa                                                      <  2:171535/139‑1 (MQ=255)                                                                  ttCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAgc                                                    >  1:564875/1‑139 (MQ=255)                                                                               ttGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGaaa                                       <  2:849308/139‑1 (MQ=255)                                                                                 ggATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAAcg                                     <  1:703192/139‑1 (MQ=255)                                                                                            cgTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAatggat                          >  2:917317/1‑139 (MQ=255)                                                                                                          tGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCtt            >  2:452300/1‑139 (MQ=255)                                                                                                          tGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCtt            >  2:872379/1‑139 (MQ=255)                                                                                                           ggATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTc           >  2:466880/1‑139 (MQ=255)                                                                                                             aTTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCg         <  1:872379/139‑1 (MQ=255)                                                                                                                 atatGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGt     <  1:318160/139‑1 (MQ=255)                                                                                                                    tGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGa  >  1:199648/1‑139 (MQ=255)                                                                                                                      ttAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCt                                                                                               >  2:206815/1‑44 (MQ=255)                                                                                                                      ttAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCt                                                                                               <  1:206815/44‑1 (MQ=255)                                                                                                                         |                                                                                                                                      CTATCGCGGTTCTTTCCTGATTGAGATGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGA  >  NC_000913/3749972‑3750225

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