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%	281.8 / NA	81	Q153* (CAA→TAA)	flk	predicted flagella assembly protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (57/24);  total (57/24)

ACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGCAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCT  >  minE/1523616‑1523725                                           |                                                                    aCCATTCTCACCCTGTTGCATGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:1004143/71‑1 (MQ=255) aCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:2720954/71‑1 (MQ=255) aCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:2545836/71‑1 (MQ=255) aCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:2252396/71‑1 (MQ=255) aCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:556982/71‑1 (MQ=255) aCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:596444/71‑1 (MQ=255) aCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:220393/71‑1 (MQ=255) aCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:682582/71‑1 (MQ=255) aCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:2019256/71‑1 (MQ=255) aCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:689296/71‑1 (MQ=255) aCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:694198/71‑1 (MQ=255) aCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:702510/71‑1 (MQ=255) aCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:1866715/71‑1 (MQ=255) aCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:1751669/71‑1 (MQ=255) aCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:1480790/71‑1 (MQ=255) aCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:1458672/71‑1 (MQ=255) aCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:813401/71‑1 (MQ=255) aCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:132754/71‑1 (MQ=255) aCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:885624/71‑1 (MQ=255) aCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:940365/71‑1 (MQ=255) aCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:1230233/71‑1 (MQ=255) aCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:1206805/71‑1 (MQ=255)     ttCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:2721284/67‑1 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGcgaccga                                            >  1:3196829/1‑54 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:3085202/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:406460/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:2832293/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:2808048/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:2741296/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:2725524/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:960055/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:314478/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:3176199/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:3285756/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:3299040/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:75911/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:488251/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:95628/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:680451/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:745025/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:1910202/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:2212661/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:2298436/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:2636443/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:2465999/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:1326697/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:2603881/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:2602067/1‑71 (MQ=255)               gTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcgg                           >  1:1323724/1‑71 (MQ=255)                        agggaagATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGcac                   >  1:2880437/1‑70 (MQ=255)                                     aTTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                         <  1:1037335/35‑1 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTgtatta                               >  1:599279/1‑42 (MQ=38)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcggc                          >  1:2540883/1‑47 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:1923873/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:221132/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:2162399/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:2155605/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:1981840/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:1237509/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:695505/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:347860/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:1900791/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:1788372/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:790268/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:1378240/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:543227/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:230167/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:2372020/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:2398434/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:2457358/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:314556/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:253719/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:3067039/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:2963230/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:2926587/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:2834587/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:2562924/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:280914/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:2594822/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:2787177/1‑71 (MQ=255)                                        ccGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCt  >  1:2613102/1‑71 (MQ=255)                                           |                                                                    ACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGCAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCT  >  minE/1523616‑1523725

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: ‑