Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,439,105	G→T	intergenic (‑19/+93)	yojL ← / ← ompC	predicted thiamine biosynthesis lipoprotein/outer membrane porin protein C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,439,105	0	G	T	100.0%	295.9 / NA	86	intergenic (‑19/+93)	yojL/ompC	predicted thiamine biosynthesis lipoprotein/outer membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (36/50);  total (36/50)

GCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAA  >  minE/1439049‑1439160                                                         |                                                        gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:2764734/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:1898246/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:1974954/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:2014798/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:2093831/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:2112044/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:2263213/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:2492379/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:2526508/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:2540386/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:2582991/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:262448/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:2703127/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:2756084/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:292865/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:2957703/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:315149/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:3518664/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:3611954/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:648544/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:673519/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:77083/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:860434/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:922247/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:939318/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:971677/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:1192421/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:1354997/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:1167136/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:1158390/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:1137336/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:1663771/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:1664204/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:1711271/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:1319073/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:1809125/62‑1 (MQ=255) gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:1824121/62‑1 (MQ=255)  cagccagcAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCTAACgg                                                    <  1:3434642/61‑1 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTac                                        >  1:404959/1‑41 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTa                    >  1:1105784/1‑61 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:3447756/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:617925/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:874069/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:331267/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:3328871/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:3342161/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:33938/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:3409247/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:558530/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:349478/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:537148/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:304293/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:429935/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:42565/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:3543277/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:1660322/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:1833251/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:2017171/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:2044219/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:1754259/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:2262488/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:175072/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:2322299/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:1860314/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:3150668/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:1631683/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:2770547/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:2809428/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:2817523/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:2918659/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:1486688/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:1479177/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:3052400/1‑62 (MQ=255)                                  tCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAc                   >  1:3140626/1‑62 (MQ=255)                                                     aTGCTAATGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGaaaa  <  1:1201630/60‑1 (MQ=255)                                                     aTGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGaaaa  <  1:714432/60‑1 (MQ=255)                                                     aTGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGaaaa  <  1:1130012/60‑1 (MQ=255)                                                     aTGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGaaaa  <  1:1312952/60‑1 (MQ=255)                                                     aTGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGaaaa  <  1:418428/60‑1 (MQ=255)                                                     aTGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGaaaa  <  1:3526865/60‑1 (MQ=255)                                                     aTGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGaaaa  <  1:1332852/60‑1 (MQ=255)                                                     aTGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGaaaa  <  1:3184440/60‑1 (MQ=255)                                                     aTGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGaaaa  <  1:1560292/60‑1 (MQ=255)                                                     aTGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGaaaa  <  1:1644507/60‑1 (MQ=255)                                                     aTGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGaaaa  <  1:2426204/60‑1 (MQ=255)                                                     aTGCTAACGGTCGCAAGAGTACACCAAAAAACTCAGTTTGTACGCTGAAAACAATGaaaa  <  1:335124/60‑1 (MQ=255)                                                         |                                                        GCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAA  >  minE/1439049‑1439160

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