Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,927,349	G→A	60.8%	T108T (ACC→ACT)	ygeA ←	predicted racemase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,927,349	0	G	A	60.8%	22.7 / 36.4	56	T108T (ACC→ACT)	ygeA	predicted racemase
Reads supporting (aligned to +/- strand):  ref base G (9/13);  new base A (16/18);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

TGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTT  >  minE/1927288‑1927410                                                              |                                                              tGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCa                                                         <  1:3065140/68‑1 (MQ=255) tGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCa                                                         <  1:3395716/68‑1 (MQ=255) tGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCa                                                         <  1:3338487/68‑1 (MQ=255) tGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCa                                                         <  1:1412306/68‑1 (MQ=255)                                  tCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTg                       >  1:1027327/1‑69 (MQ=255)                                  tCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTg                       >  1:2133943/1‑69 (MQ=255)                                    aTTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGa                      >  1:1496640/1‑68 (MQ=255)                                    aTTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGa                      >  1:1939033/1‑68 (MQ=255)                                    aTTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGa                      >  1:1700243/1‑68 (MQ=255)                                    aTTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGa                      >  1:1455019/1‑68 (MQ=255)                                    aTTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGa                      >  1:1096904/1‑68 (MQ=255)                                    aTTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGa                      >  1:105073/1‑68 (MQ=255)                                    aTTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGa                      >  1:1034574/1‑68 (MQ=255)                                             ccGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCa            <  1:2051182/69‑1 (MQ=255)                                              cGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCa            <  1:1322768/68‑1 (MQ=255)                                              cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGc             >  1:1874026/1‑67 (MQ=255)                                              cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt           >  1:3428885/1‑69 (MQ=255)                                              cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt           >  1:2370343/1‑69 (MQ=255)                                              cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt           >  1:571073/1‑69 (MQ=255)                                              cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt           >  1:3014485/1‑69 (MQ=255)                                              cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt           >  1:2803619/1‑69 (MQ=255)                                              cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt           >  1:2764270/1‑69 (MQ=255)                                              cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt           >  1:2465258/1‑69 (MQ=255)                                              cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt           >  1:1001869/1‑69 (MQ=255)                                              cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt           >  1:2306146/1‑69 (MQ=255)                                              cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt           >  1:1126611/1‑69 (MQ=255)                                              cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt           >  1:1154449/1‑69 (MQ=255)                                              cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt           >  1:153616/1‑69 (MQ=255)                                              cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt           >  1:1712997/1‑69 (MQ=255)                                              cGGTAATAGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGc                             >  1:695743/1‑51 (MQ=255)                                              cGGGAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt           >  1:588080/1‑69 (MQ=255)                                                        aCGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACtt  <  1:679067/68‑1 (MQ=255)                                                        aCGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACtt  <  1:2394380/68‑1 (MQ=255)                                                        aCGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACtt  <  1:1146532/68‑1 (MQ=255)                                                        aCGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACtt  <  1:2528315/68‑1 (MQ=255)                                                        aCGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACtt  <  1:3202289/68‑1 (MQ=255)                                                        aCGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACtt  <  1:2758891/68‑1 (MQ=255)                                                        aCGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACtt  <  1:2989509/68‑1 (MQ=255)                                                        aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGGCACtt  <  1:2152439/68‑1 (MQ=255)                                                        aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt  <  1:1285001/68‑1 (MQ=255)                                                        aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt  <  1:2301230/68‑1 (MQ=255)                                                        aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt  <  1:2370169/68‑1 (MQ=255)                                                        aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt  <  1:574176/68‑1 (MQ=255)                                                        aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt  <  1:2370408/68‑1 (MQ=255)                                                        aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt  <  1:386057/68‑1 (MQ=255)                                                        aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt  <  1:2420697/68‑1 (MQ=255)                                                        aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt  <  1:1683232/68‑1 (MQ=255)                                                        aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt  <  1:3316804/68‑1 (MQ=255)                                                        aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt  <  1:3287947/68‑1 (MQ=255)                                                        aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt  <  1:2693214/68‑1 (MQ=255)                                                        aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt  <  1:3108643/68‑1 (MQ=255)                                                        aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt  <  1:2704364/68‑1 (MQ=255)                                                        aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt  <  1:3022010/68‑1 (MQ=255)                                                        aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt  <  1:1654497/68‑1 (MQ=255)                                                         cGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt  <  1:1843503/67‑1 (MQ=255)                                                         cGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt  <  1:2060992/67‑1 (MQ=255)                                                              |                                                              TGGTGTAACGCGTACCCAGCAGCGCCACACGAGTCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTT  >  minE/1927288‑1927410

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