Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,784,296	G→A	19.6%	K456K (AAG→AAA)	yjbB →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,784,296	0	G	A	19.6%	114.7 / 16.7	76	K456K (AAG→AAA)	yjbB	predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (46/15);  new base A (10/5);  total (56/20)
Fisher's exact test for biased strand distribution p-value = 5.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

ATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAAC  >  minE/2784233‑2784360                                                                |                                                                 aTTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGc                                                             <  1:2683517/69‑1 (MQ=255)  ttttCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGc                                                             <  1:214791/68‑1 (MQ=255)             tGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGtt                                                 >  1:2224311/1‑69 (MQ=255)             tGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGtt                                                 >  1:2377066/1‑69 (MQ=255)             tGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGtt                                                 >  1:3323871/1‑69 (MQ=255)               ggTTGAAGGATCTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGttct                                               >  1:1076280/1‑69 (MQ=255)               ggTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGttct                                               >  1:317347/1‑69 (MQ=255)               ggTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGttct                                               >  1:234528/1‑69 (MQ=255)               ggTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGttct                                               >  1:538232/1‑69 (MQ=255)               ggTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGttct                                               >  1:68426/1‑69 (MQ=255)               ggTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGttct                                               >  1:711762/1‑69 (MQ=255)               ggTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGttct                                               >  1:1185251/1‑69 (MQ=255)               ggTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGttct                                               >  1:746232/1‑69 (MQ=255)                       gAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGc                                       <  1:641842/69‑1 (MQ=255)                         aCTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAACTGGCAATGTCGGTTTTCTTCTCTGGCg                                      >  1:374408/1‑68 (MQ=255)                         aCTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAACTGGCAATGTCGGTTTTCTTCTCTGGCGa                                     >  1:2670836/1‑69 (MQ=255)                         aCTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAACTGGCAATGTCGGTTTTCTTCTCTGGCGa                                     >  1:2579267/1‑69 (MQ=255)                         aCTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAACTGGCAATGTCGGTTTTCTTCTCTGGCGa                                     >  1:2717133/1‑69 (MQ=255)                         aCTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAACTGGCAATGTCGGTTTTCTTCTCTGGCGa                                     >  1:2222518/1‑69 (MQ=255)                         aCTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAACTGGCAATGTCGGTTTTCTTCTCTGGCGa                                     >  1:716479/1‑69 (MQ=255)                         aCTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAACTGGCAATGTCGGTTTTCTTCTCTGGCGa                                     >  1:3285642/1‑69 (MQ=255)                         aCTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAACTGGCAATGTCGGTTTTCTTCTCTGGCGa                                     >  1:1411484/1‑69 (MQ=255)                         aCTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAACTGGCAATGTCGGTTTTCTTCTCTGGCGa                                     >  1:822832/1‑69 (MQ=255)                         aCTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAACTGGCAATGTCGGTTTTCTTCTCTGGCGa                                     >  1:2176185/1‑69 (MQ=255)                          cTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGAt                                    <  1:2039834/69‑1 (MQ=255)                          cTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGAt                                    <  1:1757952/69‑1 (MQ=255)                             gATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTc                                 >  1:433400/1‑69 (MQ=255)                               tGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGc                                                             <  1:2178683/39‑1 (MQ=255)                                        tGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCa                                >  1:56505/1‑59 (MQ=255)                                        tGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcg                      >  1:3364262/1‑69 (MQ=255)                                        tGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcg                      >  1:2359785/1‑69 (MQ=255)                                        tGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcg                      >  1:943835/1‑69 (MQ=255)                                        tGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcg                      >  1:2294727/1‑69 (MQ=255)                                        tGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcg                      >  1:1614032/1‑69 (MQ=255)                                        tGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcg                      >  1:703565/1‑69 (MQ=255)                                        tGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcg                      >  1:1725251/1‑69 (MQ=255)                                        tGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcg                      >  1:495164/1‑69 (MQ=255)                                          aGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCAGCAGCGCtcgtc                    >  1:1258488/1‑69 (MQ=255)                                          aGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcgtc                    >  1:1876776/1‑69 (MQ=255)                                          aGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcgtc                    >  1:1754225/1‑69 (MQ=255)                                          aGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcgtc                    >  1:1411459/1‑69 (MQ=255)                                          aGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcgtc                    >  1:3323815/1‑69 (MQ=255)                                              aTTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCAc                               <  1:295525/54‑1 (MQ=255)                                              aTTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCAc                               <  1:2525383/54‑1 (MQ=255)                                              aTTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCAc                               <  1:2226463/54‑1 (MQ=255)                                              aTTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCAc                               <  1:2319599/54‑1 (MQ=255)                                              aTTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCAc                               <  1:1607761/54‑1 (MQ=255)                                              aTTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCAc                               <  1:1334527/54‑1 (MQ=255)                                              aTTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtc                       >  1:2461561/1‑62 (MQ=255)                                              aTTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGttt                >  1:3351931/1‑69 (MQ=255)                                              aTTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGttt                >  1:498922/1‑69 (MQ=255)                                              aTTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGttt                >  1:3232537/1‑69 (MQ=255)                                              aTTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGttt                >  1:2480567/1‑69 (MQ=255)                                              aTTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGttt                >  1:2700889/1‑69 (MQ=255)                                              aTTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGttt                >  1:1383914/1‑69 (MQ=255)                                              aTTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGttt                >  1:1982197/1‑69 (MQ=255)                                              aTTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGttt                >  1:2023649/1‑69 (MQ=255)                                              aTTCCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTgg                                        >  1:2669794/1‑45 (MQ=255)                                               ttGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCAc                               <  1:2958638/53‑1 (MQ=255)                                                         aTTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcgtc                    >  1:2728922/1‑54 (MQ=255)                                                         aTTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCgcagca     >  1:2969832/1‑69 (MQ=255)                                                         aTTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCgcagca     >  1:1304763/1‑69 (MQ=255)                                                         aTTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCgcagca     >  1:2529477/1‑69 (MQ=255)                                                         aTTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCgcagca     <  1:1617347/69‑1 (MQ=255)                                                         aTTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCgcagca     >  1:3035734/1‑69 (MQ=255)                                                         aTTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCgcagca     <  1:207578/69‑1 (MQ=255)                                                         aTTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCgcagc      >  1:2190521/1‑68 (MQ=255)                                                         aTTTAAAACTGGCAATGTCGGTTTTCTTCTCTGGCGATGTTACCAGCGCCCGTCGTTTGCGCCgcagca     <  1:2189716/69‑1 (MQ=255)                                                         aTTTAAAACTGGCAATGTCGGTTTTCTTCTCTGGCGATGTTACCAGCGCCCGTCGTTTGCGCCgcagca     <  1:2393913/69‑1 (MQ=255)                                                         aTTTAAAACTGGCAATGTCGGTTTTCTTCTCTGGCGATGTTACCAGCGCCCGTCGTTTGCGCCgcagca     <  1:2501763/69‑1 (MQ=255)                                                         aTTTAAAACTGGCAATGTCGGTTTTCTTCTCTGGCGATGTTACCAGCGCCCGTCGTTTGCGCCgcagca     <  1:33266/69‑1 (MQ=255)                                                          tttAAAACTGGCAATGTCGGTTTTCTTCTCTGGCGATGTTACCAGCGCCCGTCGTTTGCGCCgaagca     <  1:1142739/68‑1 (MQ=255)                                                            tAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCaaa   >  1:1717158/1‑68 (MQ=255)                                                            tAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAAc  >  1:2622039/1‑69 (MQ=255)                                                            tAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAAc  >  1:1330718/1‑69 (MQ=255)                                                            tAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAAc  >  1:128765/1‑69 (MQ=255)                                                                |                                                                 ATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTGCGTCGCAGCAAAC  >  minE/2784233‑2784360

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