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%	113.5 / NA	36	*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 (20/16);  total (20/16)

ACTATCGCGGTTCTTTCCTGATTGAGATGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAAC  >  NC_000913/3749971‑3750212                                                                                                                          |                                                                                                                         aCTATCGCGGTTCTTTCCTGATTGAGATGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGAc                                                                                                         >  2:650507/1‑139 (MQ=255)      cgcgGTTCTTTCCTGATTGAGATGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCt                                                                                                    >  2:1062169/1‑139 (MQ=255)                         agaTGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCtt                                                                                 >  2:862349/1‑139 (MQ=255)                         agaTGCGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTGATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCtt                                                                                 >  1:681731/1‑139 (MQ=255)                          gaTGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGc                                                                                                                       <  2:884695/100‑1 (MQ=255)                          gaTGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGc                                                                                                                       >  1:884695/1‑100 (MQ=255)                          gaTGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCt                                                                                  <  1:528798/137‑1 (MQ=255)                          gaTGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCt                                                                                  >  2:528798/1‑137 (MQ=255)                           aTGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTcc                                                                               >  1:389622/1‑139 (MQ=255)                           aTGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTcc                                                                               <  1:650507/139‑1 (MQ=255)                                gACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGCCGGCGAATCTGGCGCTTCCCGCTc                                                                          >  1:600196/1‑139 (MQ=255)                                 aCCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCa                                                                         <  1:862349/139‑1 (MQ=255)                                     aaaaaGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCAt                                                                     >  1:175465/1‑139 (MQ=255)                                         aGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTgg                                                                 <  2:389622/139‑1 (MQ=255)                                           ccAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTg                                                               >  1:344264/1‑139 (MQ=255)                                                                             gcgTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACg                             >  2:1049263/1‑139 (MQ=255)                                                                                ttGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGaaa                          <  1:257986/139‑1 (MQ=255)                                                                                    aTTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAAcgcg                      <  2:175465/139‑1 (MQ=255)                                                                                        aaGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCaa                  >  2:267699/1‑139 (MQ=255)                                                                                        aaGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCa                   >  2:686907/1‑138 (MQ=255)                                                                                        aaGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCa                   <  1:686907/138‑1 (MQ=255)                                                                                         aGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAat                 <  1:1049263/139‑1 (MQ=255)                                                                                         aGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAat                 <  2:344264/139‑1 (MQ=255)                                                                                             cgTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCa                                                                         >  1:786324/1‑79 (MQ=255)                                                                                             cgTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCa                                                                         <  2:786324/79‑1 (MQ=255)                                                                                             cgTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACg                                                            >  2:155157/1‑92 (MQ=255)                                                                                             cgTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACg                                                            <  1:155157/92‑1 (MQ=255)                                                                                              gTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAatggatg            <  1:267699/139‑1 (MQ=255)                                                                                                    aggaggCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATc      <  2:113753/139‑1 (MQ=255)                                                                                                      gaggCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCaa    >  2:486990/1‑139 (MQ=255)                                                                                                        ggCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCa                                                                         <  1:572371/68‑1 (MQ=255)                                                                                                        ggCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCa                                                                         >  2:572371/1‑68 (MQ=255)                                                                                                        ggCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAAc  >  2:208253/1‑139 (MQ=255)                                                                                                                  atatGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGaaa                          >  2:741702/1‑105 (MQ=255)                                                                                                                  atatGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGaaa                          <  1:741702/105‑1 (MQ=255)                                                                                                                       ttAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAatg                >  1:857087/1‑110 (MQ=255)                                                                                                                          |                                                                                                                         ACTATCGCGGTTCTTTCCTGATTGAGATGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAAC  >  NC_000913/3749971‑3750212

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