Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,523,658	C→T	Q153* (CAA→TAA)	flk →	predicted flagella assembly protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,523,658	0	C	T	100.0%	160.8 / NA	47	Q153* (CAA→TAA)	flk	predicted flagella assembly protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/14);  total (33/14)

TCAGCTCTCACCGGAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGCAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCG  >  minE/1523591‑1523702                                                                    |                                             tCAGCTCTCACCGGAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTaa                                            >  1:603988/1‑70 (MQ=255) tCAGCTCTCACCGGAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAAc                                           >  1:1324525/1‑71 (MQ=255) tCAGCTCTCACCGGAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAAc                                           >  1:858888/1‑71 (MQ=255) tCAGCTCTCACCGGAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAAc                                           >  1:3002425/1‑71 (MQ=255) tCAGCTCTCACCGGAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAAc                                           >  1:136994/1‑71 (MQ=255) tCAGCTCTCACCGGAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAAc                                           >  1:668808/1‑71 (MQ=255) tCAGCTCTCACCGGAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAAc                                           >  1:1511796/1‑71 (MQ=255) tCAGCTCTCACCGGAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAAc                                           >  1:3011083/1‑71 (MQ=255) tCAGCTCTCACCGGAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAAc                                           >  1:3039446/1‑71 (MQ=255) tCAGCTCTCACCGGAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAAc                                           >  1:3509823/1‑71 (MQ=255) tCAGCTCTCACCGGAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAAc                                           >  1:3533599/1‑71 (MQ=255) tCAGCTCTCACCGGAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAAc                                           >  1:2419790/1‑71 (MQ=255)             ggAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCGCAGCAgcgc                                >  1:2834775/1‑70 (MQ=255)             ggAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAgcgc                                >  1:3905262/1‑70 (MQ=255)             ggAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAgcgc                                >  1:4025252/1‑70 (MQ=255)             ggAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAgcgc                                >  1:4180549/1‑70 (MQ=255)             ggAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAgcgc                                >  1:50264/1‑70 (MQ=255)             ggAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAgcgc                                >  1:3245385/1‑70 (MQ=255)             ggAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAgcgc                                >  1:505088/1‑70 (MQ=255)             ggAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAgcgc                                >  1:543675/1‑70 (MQ=255)             ggAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAgcgc                                >  1:831425/1‑70 (MQ=255)             ggAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAgcgc                                >  1:971987/1‑70 (MQ=255)             ggAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAgcgc                                >  1:1354611/1‑70 (MQ=255)             ggAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAgcgc                                >  1:2734988/1‑70 (MQ=255)             ggAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAgcgc                                >  1:1486937/1‑70 (MQ=255)             ggAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAgcgc                                >  1:2211655/1‑70 (MQ=255)             ggAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAgcgc                                >  1:2506037/1‑70 (MQ=255)                     tAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGcgacc                       >  1:1577180/1‑71 (MQ=255)                     tAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGcgacc                       >  1:1582610/1‑71 (MQ=255)                     tAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGcgacc                       >  1:2293793/1‑71 (MQ=255)                     tAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGcgac                        >  1:4354741/1‑70 (MQ=255)                         aaCCATTCTCAGCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCgaccgacc                   <  1:1782009/71‑1 (MQ=255)                         aaCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCgaccgacc                   <  1:1297538/71‑1 (MQ=255)                         aaCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCgaccgacc                   <  1:3313304/71‑1 (MQ=255)                                   aCCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCttt            <  1:3831950/68‑1 (MQ=255)                                   aCCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCttt            <  1:2065395/68‑1 (MQ=255)                                   aCCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCttt            <  1:4194007/68‑1 (MQ=255)                                   aCCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCttt            <  1:2558863/68‑1 (MQ=255)                                   aCCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCttt            <  1:3236230/68‑1 (MQ=255)                                   aCCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCttt            <  1:3135264/68‑1 (MQ=255)                                   aCCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCttt            <  1:2711930/68‑1 (MQ=255)                                   aCCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCttt            <  1:1358063/68‑1 (MQ=255)                                   aCCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCttt            <  1:2940765/68‑1 (MQ=255)                                   aCCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCttt            <  1:2796030/68‑1 (MQ=255)                                   aCCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGAGAGGCGACCGACCGTCCttt            <  1:2007309/68‑1 (MQ=255)                                           gCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcggcg  >  1:2466658/1‑70 (MQ=255)                                           gCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcggcg  >  1:2634766/1‑70 (MQ=255)                                                                    |                                             TCAGCTCTCACCGGAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGCAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCG  >  minE/1523591‑1523702

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