Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	512,180	G→A	51.3%	A741A (GCG→GCA)	ybhJ →	predicted hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	512,180	0	G	A	51.3%	‑4.3 / 69.3	76	A741A (GCG→GCA)	ybhJ	predicted hydratase
Reads supporting (aligned to +/- strand):  ref base G (22/15);  new base A (18/21);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTAACGGAAATTACGCTCTATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCATGTGAATGTAG  >  minE/512117‑512239                                                                |                                                            ggTAACGGAAATTACGCTCTATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAg                                                        <  1:1159316/69‑1 (MQ=255)        gAAATTACGCTCTATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGTTTGAtt                                                 <  1:221896/69‑1 (MQ=255)             tACGCTCTATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGttt                                                     >  1:2510865/1‑60 (MQ=255)             tACGCTCTATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGttt                                                     >  1:660693/1‑60 (MQ=255)             tACGCTCTATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGttt                                                     >  1:2949055/1‑60 (MQ=255)             tACGCTCTATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGttt                                                     >  1:3158493/1‑60 (MQ=255)             tACGCTCTATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGttt                                                     >  1:3143872/1‑60 (MQ=255)             tACGCTCTATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGttt                                                     >  1:1254812/1‑60 (MQ=255)              aCGCTCTATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGTTTGATTAACTTc                                           >  1:2388224/1‑69 (MQ=255)              aCGCTCTATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGTTTGATTAACTTc                                           >  1:2513656/1‑69 (MQ=255)              aCGCTCTATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGTTTGATTAACTTc                                           >  1:3203608/1‑69 (MQ=255)              aCGCTCTATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGTTTGATTAACTTc                                           >  1:1082332/1‑69 (MQ=255)                  tctATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCAGGTAGTTTGATTa                                                <  1:2031068/60‑1 (MQ=255)                  tctATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCAGGTAGTTTGATTa                                                <  1:2098582/60‑1 (MQ=255)                  tctATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCAGGTAGTTTGATTa                                                <  1:1969335/60‑1 (MQ=255)                  tctATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCAGGTAGTTTGATTa                                                <  1:987347/60‑1 (MQ=255)                  tctATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCAGGTAGTTTGATTa                                                <  1:1926402/60‑1 (MQ=255)                  tctATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCAGGTAGTTTGATTa                                                <  1:251232/60‑1 (MQ=255)                  tctATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCAGGTAGTTTGATTa                                                <  1:251670/60‑1 (MQ=255)                  tctATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCAGGTAGTTTGATTa                                                <  1:2536202/60‑1 (MQ=255)                  tctATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCAGGTAGTTTGATTa                                                <  1:3223792/60‑1 (MQ=255)                  tctATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCAGGTAGTTTGATTa                                                <  1:499914/60‑1 (MQ=255)                  tctATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCAGGTAGTTTGATTa                                                <  1:1186611/60‑1 (MQ=255)                  tctATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCAGGTAGTTTGATTa                                                <  1:1709414/60‑1 (MQ=255)                  tctATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCAGGTAGTTTGATTa                                                <  1:1110761/60‑1 (MQ=255)                  tctATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCAGGTAGTTTGATTa                                                <  1:1667611/60‑1 (MQ=255)                  tctATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCAGGTAGTTTGATTa                                                <  1:1251857/60‑1 (MQ=255)                  tctATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCAGGTAGTTTGATTa                                                <  1:1647373/60‑1 (MQ=255)                  tctATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGAAGGTAGTTTGATTa                                                <  1:1379226/60‑1 (MQ=255)                    tataTGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCAGGTAGTTTGATTa                                                <  1:3062677/58‑1 (MQ=255)                       atGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAAc                                  >  1:2960807/1‑69 (MQ=255)                       atGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAAc                                  >  1:1859438/1‑69 (MQ=255)                       atGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAAc                                  >  1:704775/1‑69 (MQ=255)                       atGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAAc                                  >  1:1582536/1‑69 (MQ=255)                            aaGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAACCGTCa                             <  1:1160893/69‑1 (MQ=255)                            aaGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAACCGTCa                             <  1:2780539/69‑1 (MQ=255)                            aaGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAACCGTCa                             <  1:3290977/69‑1 (MQ=255)                            aaGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAACCGTCa                             <  1:3090949/69‑1 (MQ=255)                            aaGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAACCGTCa                             <  1:234689/69‑1 (MQ=255)                            aaGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAACCGTCa                             <  1:2177869/69‑1 (MQ=255)                            aaGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAACCGTCa                             <  1:687816/69‑1 (MQ=255)                            aaGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAACCGTCa                             <  1:1723425/69‑1 (MQ=255)                            aaGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAACCGTCa                             <  1:879798/69‑1 (MQ=255)                            aaGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGAGGGTAGTTTGATTAACTTCAATAAAAACCGTCa                             <  1:2352010/69‑1 (MQ=255)                             aGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCAGGTAGTTTGATTAACTTCAATAAAAACCGTCAg                            <  1:1795296/69‑1 (MQ=255)                             aGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCAGGTAGTTTGATTAACTTCAATAAAAACCGTCAg                            <  1:1931859/69‑1 (MQ=255)                              gTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCAGGTAGTTTGATTAACTTCAATAAAAACCGTCAg                            <  1:1220297/68‑1 (MQ=255)                                            gagCGCGAGATTATCAAGGCATGTAGTTTGATTAACTTCaa                                         >  1:187502/1‑41 (MQ=255)                                            gagCGCGAGATTATCAAGGCATGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGcc              >  1:115312/1‑68 (MQ=255)                                            gagCGCGAGATTATCAAGGCATGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCa             >  1:2885137/1‑69 (MQ=255)                                            gagCGCGAGATTATCAAGGCATGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCa             >  1:2834517/1‑69 (MQ=255)                                            gagCGCGAGATTATCAAGGCATGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCa             >  1:260474/1‑69 (MQ=255)                                            gagCGCGAGATTATCAAGGCATGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCa             >  1:3110980/1‑69 (MQ=255)                                            gagCGCGAGATTATCAAGGCATGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCa             >  1:618352/1‑69 (MQ=255)                                            gagCGCGAGATTATCAAGGCAGGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCa             >  1:3390533/1‑69 (MQ=255)                                            gagCGCGAGATTATCAAGGCAGGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCa             >  1:3218192/1‑69 (MQ=255)                                            gagCGCGAGATTATCAAGGCAGGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCa             >  1:2985425/1‑69 (MQ=255)                                            gagCGCGAGATTATCAAGGCAGGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCa             >  1:2100117/1‑69 (MQ=255)                                            gagCGCGAGATTATCAAGGCAGGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCa             >  1:2217789/1‑69 (MQ=255)                                                gcgAGATTATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCAtgtg         <  1:678796/69‑1 (MQ=255)                                                  gagaTTATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAACCGGCAGATGTAAAAAGCGCCATGTGaa       <  1:1802983/69‑1 (MQ=255)                                                       tATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAAcc                                 >  1:2733862/1‑38 (MQ=255)                                                       tATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCATGTGAATGTa   >  1:1536347/1‑68 (MQ=255)                                                       tATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCATGTGAATGTAg  >  1:1007273/1‑69 (MQ=255)                                                       tATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCATGTGAATGTAg  >  1:3179828/1‑69 (MQ=255)                                                       tATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCATGTGAATGTAg  >  1:692752/1‑69 (MQ=255)                                                       tATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCATGTGAATGTAg  >  1:1381946/1‑69 (MQ=255)                                                       tATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCATGTGAATGTAg  >  1:3278600/1‑69 (MQ=255)                                                       tATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCATGTGAATGTAg  >  1:324086/1‑69 (MQ=255)                                                         tCAAGGCGGGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCAtgtg         <  1:332463/60‑1 (MQ=255)                                                          cAAGGCAGGTAGTTTGATTAACTTCAATAAAAAcc                                 >  1:1633280/1‑35 (MQ=255)                                                          cAAGGCAGGTAGTTTGATTAACTTCAATAAAAAcc                                 >  1:3260450/1‑35 (MQ=255)                                                          cAAGGCAGGTAGTTTGATTAACTTCAATAAAAAcc                                 >  1:1980420/1‑35 (MQ=255)                                                          cAAGGCAGGTAGTTTGATTAACTTCAATAAAAAcc                                 >  1:1281398/1‑35 (MQ=255)                                                          cAAGGCAGGTAGTTTGATTAACTTCAATAAAAAcc                                 >  1:907901/1‑35 (MQ=255)                                                          cAAGGCAGGTAGTTTGATTAACTTCAATAAAAAcc                                 >  1:3147065/1‑35 (MQ=255)                                                                |                                                            GGTAACGGAAATTACGCTCTATATGGAAAGTCTGACTGCTGAAGAGCGCGAGATTATCAAGGCGGGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCATGTGAATGTAG  >  minE/512117‑512239

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