Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,927,361	A→T	54.0%	I104I (ATT→ATA)	ygeA ←	predicted racemase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,927,361	0	A	T	54.0%	9.9 / 70.0	63	I104I (ATT→ATA)	ygeA	predicted racemase
Reads supporting (aligned to +/- strand):  ref base A (20/9);  new base T (16/18);  total (36/27)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

TCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGG  >  minE/1927321‑1927425                                         |                                                                 tCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTg                                      >  1:1027327/1‑69 (MQ=255) tCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTg                                      >  1:2133943/1‑69 (MQ=255)   aTTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGa                                     >  1:105073/1‑68 (MQ=255)   aTTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGa                                     >  1:1096904/1‑68 (MQ=255)   aTTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGa                                     >  1:1939033/1‑68 (MQ=255)   aTTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGa                                     >  1:1455019/1‑68 (MQ=255)   aTTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGa                                     >  1:1496640/1‑68 (MQ=255)   aTTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGa                                     >  1:1034574/1‑68 (MQ=255)   aTTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGa                                     >  1:1700243/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:153616/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                          >  1:2306146/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                          >  1:2370343/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                          >  1:2465258/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                          >  1:1712997/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                          >  1:2764270/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                          >  1:2803619/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                          >  1:3014485/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                          >  1:3428885/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                          >  1:571073/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                          >  1:1154449/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                          >  1:1001869/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                          >  1:1126611/1‑69 (MQ=255)             cGGTAATAGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGc                                            >  1:695743/1‑51 (MQ=255)             cGGGAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                          >  1:588080/1‑69 (MQ=255)                       aCGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACtt                 <  1:2758891/68‑1 (MQ=255)                       aCGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACtt                 <  1:1146532/68‑1 (MQ=255)                       aCGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACtt                 <  1:679067/68‑1 (MQ=255)                       aCGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACtt                 <  1:3202289/68‑1 (MQ=255)                       aCGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACtt                 <  1:2989509/68‑1 (MQ=255)                       aCGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACtt                 <  1:2528315/68‑1 (MQ=255)                       aCGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACtt                 <  1:2394380/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGGCACtt                 <  1:2152439/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                 <  1:2693214/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                 <  1:3022010/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                 <  1:1285001/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                 <  1:3108643/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                 <  1:3287947/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                 <  1:1654497/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                 <  1:1683232/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                 <  1:3316804/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                 <  1:2704364/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                 <  1:2301230/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                 <  1:386057/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                 <  1:2420697/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                 <  1:574176/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                 <  1:2370408/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                 <  1:2370169/68‑1 (MQ=255)                        cGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                 <  1:1843503/67‑1 (MQ=255)                        cGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                 <  1:2060992/67‑1 (MQ=255)                                      cGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTgg  >  1:2272222/1‑68 (MQ=255)                                      cGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTgg  >  1:2221063/1‑68 (MQ=255)                                      cGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTgg  >  1:280872/1‑68 (MQ=255)                                      cGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTgg  >  1:340634/1‑68 (MQ=255)                                      cGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTgg  >  1:1168408/1‑68 (MQ=255)                                      cGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTgg  >  1:1819386/1‑68 (MQ=255)                                      cGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTgg  >  1:1167975/1‑68 (MQ=255)                                      cGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTgg  >  1:2399073/1‑68 (MQ=255)                                      cGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTgg  >  1:2185205/1‑68 (MQ=255)                                      cGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTgg  >  1:185625/1‑68 (MQ=255)                                      cGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTg   >  1:2215454/1‑67 (MQ=255)                                         |                                                                 TCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGG  >  minE/1927321‑1927425

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