Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,680,754	T→C	Y1407C (TAC→TGC)	yfhM ←	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,680,754	0	T	C	100.0%	226.0 / NA	61	Y1407C (TAC→TGC)	yfhM	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/25);  total (36/25)

GCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGTACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAAC  >  minE/1680711‑1680804                                            |                                                   gCGTTCACCAAACGCCTGCTGGGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:986917/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGtt                               >  1:1698278/1‑65 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:2072394/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:2469428/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:2139974/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:2301543/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:2352389/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:2431646/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:2509165/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:2610076/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:2722364/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:2762975/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:2800289/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:3042134/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:55409/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:598912/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:644618/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:8982/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:90189/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:1355936/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:1391829/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:1808625/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:1433375/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:146047/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:1376381/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:172864/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:1688152/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                         >  1:1129139/1‑71 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAttt                          >  1:779765/1‑70 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAttt                          >  1:2841587/1‑70 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAttt                          >  1:1675553/1‑70 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAttt                          >  1:2547642/1‑70 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAttt                          >  1:966781/1‑70 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAttt                          >  1:744546/1‑70 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAttt                          >  1:1345594/1‑70 (MQ=255) gCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAttt                          >  1:2247229/1‑70 (MQ=255)                      cGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTa    <  1:3055721/71‑1 (MQ=255)                      cGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTa    <  1:2853973/71‑1 (MQ=255)                      cGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTa    <  1:309326/71‑1 (MQ=255)                      cGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTa    <  1:318386/71‑1 (MQ=255)                      cGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTa    <  1:2647914/71‑1 (MQ=255)                      cGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTa    <  1:982618/71‑1 (MQ=255)                      cGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTa    <  1:2305914/71‑1 (MQ=255)                      cGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTa    <  1:178493/71‑1 (MQ=255)                      cGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTa    <  1:2265777/71‑1 (MQ=255)                      cGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTa    <  1:2003175/71‑1 (MQ=255)                       gAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTa    <  1:1779589/70‑1 (MQ=255)                       gAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTa    <  1:192164/70‑1 (MQ=255)                        aaaGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAAc  <  1:1214475/71‑1 (MQ=255)                        aaaGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAAc  <  1:900772/71‑1 (MQ=255)                        aaaGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAAc  <  1:1105816/71‑1 (MQ=255)                        aaaGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAAc  <  1:336981/71‑1 (MQ=255)                        aaaGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAAc  <  1:2989232/71‑1 (MQ=255)                        aaaGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAAc  <  1:2813229/71‑1 (MQ=255)                        aaaGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAAc  <  1:280696/71‑1 (MQ=255)                        aaaGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAAc  <  1:1068959/71‑1 (MQ=255)                        aaaGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAAc  <  1:2454437/71‑1 (MQ=255)                        aaaGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAAc  <  1:2433827/71‑1 (MQ=255)                        aaaGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAAc  <  1:2113702/71‑1 (MQ=255)                        aaaGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCATCCAGCAAGGAGAGCATTAAc  <  1:2792957/71‑1 (MQ=255)                          aGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGATAGCATTAAc  <  1:1345738/69‑1 (MQ=255)                                         gtgcACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAAc  <  1:310090/54‑1 (MQ=255)                                            |                                                   GCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGTACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAAC  >  minE/1680711‑1680804

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