Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,927,388	T→A	41.9%	S95S (TCA→TCT)	ygeA ←	predicted racemase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,927,388	0	T	A	41.9%	39.2 / 87.6	86	S95S (TCA→TCT)	ygeA	predicted racemase
Reads supporting (aligned to +/- strand):  ref base T (35/15);  new base A (15/21);  total (50/36)
Fisher's exact test for biased strand distribution p-value = 1.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATAGCACAATACCTTCTGCGC  >  minE/1927321‑1927449                                                                    |                                                              tCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTg                                                              >  1:1027327/1‑69 (MQ=255) tCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTg                                                              >  1:2133943/1‑69 (MQ=255)   aTTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGa                                                             >  1:105073/1‑68 (MQ=255)   aTTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGa                                                             >  1:1700243/1‑68 (MQ=255)   aTTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGa                                                             >  1:1496640/1‑68 (MQ=255)   aTTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGa                                                             >  1:1939033/1‑68 (MQ=255)   aTTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGa                                                             >  1:1455019/1‑68 (MQ=255)   aTTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGa                                                             >  1:1096904/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:2465258/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                                                  >  1:2764270/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                                                  >  1:2803619/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                                                  >  1:1001869/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                                                  >  1:2306146/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                                                  >  1:2370343/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                                                  >  1:3014485/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                                                  >  1:1712997/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                                                  >  1:1126611/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                                                  >  1:3428885/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                                                  >  1:153616/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                                                  >  1:571073/1‑69 (MQ=255)             cGGTAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                                                  >  1:1154449/1‑69 (MQ=255)             cGGGAATTGCACGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCAt                                                  >  1:588080/1‑69 (MQ=255)                       aCGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACtt                                         <  1:3202289/68‑1 (MQ=255)                       aCGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACtt                                         <  1:2394380/68‑1 (MQ=255)                       aCGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACtt                                         <  1:1146532/68‑1 (MQ=255)                       aCGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACtt                                         <  1:2989509/68‑1 (MQ=255)                       aCGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACtt                                         <  1:679067/68‑1 (MQ=255)                       aCGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACtt                                         <  1:2758891/68‑1 (MQ=255)                       aCGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACtt                                         <  1:2528315/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGGCACtt                                         <  1:2152439/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                                         <  1:2370408/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                                         <  1:2420697/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                                         <  1:2693214/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                                         <  1:2704364/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                                         <  1:574176/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                                         <  1:386057/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                                         <  1:3022010/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                                         <  1:3108643/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                                         <  1:3316804/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                                         <  1:3287947/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                                         <  1:2370169/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                                         <  1:1654497/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                                         <  1:1683232/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                                         <  1:1285001/68‑1 (MQ=255)                       aCGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                                         <  1:2301230/68‑1 (MQ=255)                        cGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                                         <  1:1843503/67‑1 (MQ=255)                        cGTCCAGTGGCATCCGCTATGTGTAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                                         <  1:2060992/67‑1 (MQ=255)                                      cGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTgg                          >  1:2221063/1‑68 (MQ=255)                                      cGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTgg                          >  1:340634/1‑68 (MQ=255)                                      cGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTgg                          >  1:2399073/1‑68 (MQ=255)                                      cGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTgg                          >  1:1819386/1‑68 (MQ=255)                                      cGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTgg                          >  1:185625/1‑68 (MQ=255)                                      cGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTgg                          >  1:280872/1‑68 (MQ=255)                                      cGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTgg                          >  1:1168408/1‑68 (MQ=255)                                      cGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTgg                          >  1:2185205/1‑68 (MQ=255)                                      cGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTgg                          >  1:1167975/1‑68 (MQ=255)                                      cGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTgg                          >  1:2272222/1‑68 (MQ=255)                                      cGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTg                           >  1:2215454/1‑67 (MQ=255)                                            gtgtAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATa                    >  1:3411245/1‑68 (MQ=255)                                            gtgtAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATAg                   >  1:213746/1‑69 (MQ=255)                                            gtgtAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATAg                   >  1:2974786/1‑69 (MQ=255)                                            gtgtAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATAg                   >  1:433267/1‑69 (MQ=255)                                            gtgtAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATAg                   >  1:267920/1‑69 (MQ=255)                                               tAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATAGcac                <  1:1393746/69‑1 (MQ=255)                                               tAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATAGcac                <  1:1573095/69‑1 (MQ=255)                                               tAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATAGcac                <  1:3092376/69‑1 (MQ=255)                                               tAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATAGcac                <  1:3340978/69‑1 (MQ=255)                                               tAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATAGcac                <  1:3059476/69‑1 (MQ=255)                                               tAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATAGcac                <  1:296663/69‑1 (MQ=255)                                               tAAGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                                         <  1:42484/44‑1 (MQ=38)                                                 aGAAAGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                                         <  1:766524/42‑1 (MQ=38)                                                   aaaGGCAGAGAGCAACGAGACTCAATGGCATCCGCCACtt                                         <  1:215441/40‑1 (MQ=37)                                                            gTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATAGCACAATACCTTCTgcg   >  1:709055/1‑69 (MQ=255)                                                             tGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATAGCACAATACCTTCTgcgc  >  1:1573443/1‑69 (MQ=255)                                                             tGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATAGCACAATACCTTCTgcgc  >  1:321554/1‑69 (MQ=255)                                                             tGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATAGCACAATACCTTCTgcgc  >  1:3406969/1‑69 (MQ=255)                                                             tGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATAGCACAATACCTTCTgcgc  >  1:238400/1‑69 (MQ=255)                                                             tGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATAGCACAATACCTTCTgcgc  >  1:146625/1‑69 (MQ=255)                                                             tGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATAGCACAATACCTTCTgcgc  >  1:2196658/1‑69 (MQ=255)                                                             tGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATAGCACAATACCTTCTgcgc  >  1:95159/1‑69 (MQ=255)                                                             tGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATAGCACAATACCTTCTgcg   >  1:3090650/1‑68 (MQ=255)                                                             tGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATAGCACAATACCTTCTgcg   >  1:3013111/1‑68 (MQ=255)                                                                    |                                                              TCATTCCGGCCCCGGTAATTGCACGTCCGGTGGCATCCGCAATGTGTAAGAAAGGCAGAGTGCAACGTGACTCAATGGCATCCGCCACTTTATGCATCGTATTGGTACATAGCACAATACCTTCTGCGC  >  minE/1927321‑1927449

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