Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,100,317	C→T	49.3%	L39L (CTG→TTG)	ydiN →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,100,317	0	C	T	49.3%	‑1.2 / 68.7	77	L39L (CTG→TTG)	ydiN	predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (30/9);  new base T (6/32);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 6.72e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

TCAGCTACTTCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCGGGCATTGCCTACTTAA  >  minE/1100261‑1100362                                                         |                                              taaGCTATTTCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAACATGACATCTTTGGCAGaaaa                                     >  1:118001/3‑67 (MQ=255) tCAGCTATTTCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAACATGACATCTTTGGCAGaaaa                                     >  1:1848051/1‑67 (MQ=255) tCAGCTATTTCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAACATGACATCTTTGGCAGaaaa                                     >  1:1296371/1‑67 (MQ=255) tCAGCTATTTCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAACATGACATCTTTGGCAGaaaa                                     >  1:3096086/1‑67 (MQ=255) tCAGCTATTTCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAACATGACATCTTTGGCAGaaaa                                     >  1:1684977/1‑67 (MQ=255) tCAGCTATTTCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAACATGACATCTTTGGCAGaaaa                                     >  1:2066949/1‑67 (MQ=255)     cTACTTCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAATATGTCATCTCTGGc                                           >  1:1446425/1‑57 (MQ=255)     cTACTTCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTc                               >  1:1740839/1‑69 (MQ=255)     cTACTTCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTc                               >  1:1173601/1‑69 (MQ=255)     cTACTTCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTc                               >  1:3266185/1‑69 (MQ=255)     cTACTTCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTc                               >  1:3175537/1‑69 (MQ=255)     cTACTTCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTc                               >  1:947411/1‑69 (MQ=255)     cTACTTCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTc                               >  1:2841002/1‑69 (MQ=255)     cTACTTCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTc                               >  1:1644082/1‑69 (MQ=255)          tCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACt                           <  1:669295/68‑1 (MQ=255)          tCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACt                           <  1:1934211/68‑1 (MQ=255)          tCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACt                           <  1:882082/68‑1 (MQ=255)          tCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACt                           <  1:2820970/68‑1 (MQ=255)          tCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACt                           <  1:2130560/68‑1 (MQ=255)          tCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACt                           <  1:293635/68‑1 (MQ=255)          tCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACt                           <  1:612477/68‑1 (MQ=255)          tCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACt                           <  1:497008/68‑1 (MQ=255)                     aTTAGTGTTATTACGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCGGGc               <  1:465426/69‑1 (MQ=255)                     aTTAGTGTTATTACGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCGGGc               <  1:1045506/69‑1 (MQ=255)                      ttAGTGTTATTACGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCGGGCa              <  1:1934381/69‑1 (MQ=255)                        aGTGTTATTACGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCGGGCAtt            >  1:785760/1‑69 (MQ=255)                          ttttattaCGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCGGGCATTg           <  1:1989399/66‑1 (MQ=255)                          tgtTATTACGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCGGGCATTg           <  1:2946709/68‑1 (MQ=255)                          tgtTATTACGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCGGGCATTg           <  1:1527548/68‑1 (MQ=255)                          tgtTATTACGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCGGGCATTg           <  1:2660921/68‑1 (MQ=255)                           gtTATTACGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCGGGCATTg           <  1:802769/67‑1 (MQ=255)                            ttattaCGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCGGGCATTGCCt        >  1:190485/1‑69 (MQ=255)                            ttattaCGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCGGGCATTGCCt        >  1:1295079/1‑69 (MQ=255)                              attaCGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGt                                   >  1:2597411/1‑40 (MQ=255)                              attaCGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTc                               >  1:2465238/1‑44 (MQ=255)                              attaCGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTc                               >  1:1365626/1‑44 (MQ=255)                              attaCGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTc                               >  1:1512318/1‑44 (MQ=255)                              attaCGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTc                               >  1:1534939/1‑44 (MQ=255)                              attaCGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTc                               >  1:3393166/1‑44 (MQ=255)                              attaCGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTc                               >  1:1986176/1‑44 (MQ=255)                              attaCGCTTGCCCAAAATATGTCATCTCTGGCCGAAAAGTTTTCCACTGACAACGCGGGCATTGCCTa       >  1:645117/1‑68 (MQ=255)                              attaCGCTTGCCCAAAATATGTCATATCTGGCGGAAAAGTTTTc                               >  1:3295980/1‑44 (MQ=255)                                taCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:1081464/57‑1 (MQ=255)                                taCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:1100327/57‑1 (MQ=255)                                taCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:3370645/57‑1 (MQ=255)                                taCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:3390352/57‑1 (MQ=255)                                taCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:3399778/57‑1 (MQ=255)                                taCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:1353362/57‑1 (MQ=255)                                taCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:3419994/57‑1 (MQ=255)                                taCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:3122879/57‑1 (MQ=255)                                taCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:3114475/57‑1 (MQ=255)                                taCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:2956316/57‑1 (MQ=255)                                taCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:822916/57‑1 (MQ=255)                                taCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:2201377/57‑1 (MQ=255)                                taCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:2178211/57‑1 (MQ=255)                                taCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:2151983/57‑1 (MQ=255)                                taCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:2119809/57‑1 (MQ=255)                                taCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:1555002/57‑1 (MQ=255)                                taCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:1642379/57‑1 (MQ=255)                                taCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:1953563/57‑1 (MQ=255)                                taCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:1728570/57‑1 (MQ=255)                                taCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:1930863/57‑1 (MQ=255)                                 aCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:974222/56‑1 (MQ=255)                                 aCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:555797/56‑1 (MQ=255)                                 aCGCTTGCCCAAAACATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:468142/56‑1 (MQ=255)                                  cGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCGGGCATTGCCTACTTaa  >  1:388090/1‑69 (MQ=255)                                    cTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCGGGCATTGc          >  1:2919915/1‑59 (MQ=255)                                            aaaCATGACATCTTTGGCAGAAAAGTTTTCCACTGACAACGCggg                <  1:1699533/45‑1 (MQ=255)                                               tatGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCggg                >  1:527931/1‑42 (MQ=255)                                               tatGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCggg                >  1:3266629/1‑42 (MQ=255)                                               tatGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCggg                >  1:772139/1‑42 (MQ=255)                                               tatGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCggg                >  1:2903068/1‑42 (MQ=255)                                               tatGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCggg                >  1:2547176/1‑42 (MQ=255)                                               tatGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCggg                >  1:2533727/1‑42 (MQ=255)                                               tatGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCggg                >  1:1968082/1‑42 (MQ=255)                                               tatGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCggg                >  1:1407903/1‑42 (MQ=255)                                                 tGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCGGGc               <  1:3265569/41‑1 (MQ=255)                                                         |                                              TCAGCTACTTCCTGCACGGCATTAGTGTTATTACGCTTGCCCAAAATATGTCATCTCTGGCGGAAAAGTTTTCCACTGACAACGCGGGCATTGCCTACTTAA  >  minE/1100261‑1100362

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