Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,576,174	G→C	68.6%	A290A (GCC→GCG)	xapB ←	xanthosine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,576,174	0	G	C	68.6%	99.5 / 60.3	86	A290A (GCC→GCG)	xapB	xanthosine transporter
Reads supporting (aligned to +/- strand):  ref base G (18/9);  new base C (37/22);  total (55/31)
Fisher's exact test for biased strand distribution p-value = 8.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCAAATCGCTTTAAAAAGAA  >  minE/1576110‑1576227                                                                 |                                                      gcagcaAAATAAATACGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                                    >  1:936849/1‑68 (MQ=255)  cagcaAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGccac                                                  <  1:2103807/69‑1 (MQ=255)  cagcaAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGccac                                                  <  1:1258314/69‑1 (MQ=255)  cagcaAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGccac                                                  <  1:90322/69‑1 (MQ=255)  cagcaAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGccac                                                  <  1:53191/69‑1 (MQ=255)      aaaaTAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATa                                              <  1:1561827/69‑1 (MQ=255)      aaaaTAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATa                                              <  1:1341646/69‑1 (MQ=255)      aaaaTAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATa                                              <  1:362461/69‑1 (MQ=255)      aaaaTAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATa                                              <  1:1398417/69‑1 (MQ=255)       aaataaatCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                                    >  1:2193407/1‑62 (MQ=255)       aaataaatCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                                    >  1:1572764/1‑62 (MQ=255)       aaataaatCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                                    >  1:3182486/1‑62 (MQ=255)       aaataaatCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                                    >  1:3247490/1‑62 (MQ=255)       aaataaatCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                                    >  1:760793/1‑62 (MQ=255)       aaataaatCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                                    >  1:339684/1‑62 (MQ=255)       aaataaatCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                                    >  1:291942/1‑62 (MQ=255)                        gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:357314/66‑1 (MQ=255)                        gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:1800702/66‑1 (MQ=255)                        gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:1863293/66‑1 (MQ=255)                        gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:186537/66‑1 (MQ=255)                        gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:1596031/66‑1 (MQ=255)                        gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:1515360/66‑1 (MQ=255)                        gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:3171839/66‑1 (MQ=255)                        gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:2580671/66‑1 (MQ=255)                        gACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGccac                                                  <  1:2309410/47‑1 (MQ=255)                                gCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGc                     >  1:1276088/1‑68 (MQ=255)                                gCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGCTTTTAATGc                     >  1:3136816/1‑68 (MQ=255)                                     aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACacatc                                        >  1:996730/1‑44 (MQ=37)                                     aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:1694877/1‑64 (MQ=255)                                     aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:3321130/1‑64 (MQ=255)                                     aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:45313/1‑64 (MQ=255)                                     aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:2092086/1‑64 (MQ=255)                                     aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:460814/1‑64 (MQ=255)                                     aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:471413/1‑64 (MQ=255)                                     aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:559691/1‑64 (MQ=255)                                     aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:766413/1‑64 (MQ=255)                                     aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:1369850/1‑64 (MQ=255)                                     aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:836047/1‑64 (MQ=255)                                     aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:3048902/1‑64 (MQ=255)                                     aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:293713/1‑64 (MQ=255)                                     aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:901355/1‑64 (MQ=255)                                     aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:1123922/1‑64 (MQ=255)                                     aGGCAAAGAAGCCAAAGCGCAGCGTCAACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:2616232/1‑64 (MQ=255)                                      ggCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCAAATCg              >  1:971656/1‑69 (MQ=255)                                      ggCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCAAATCg              >  1:447414/1‑69 (MQ=255)                                              agCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCaaaa      <  1:1292676/69‑1 (MQ=255)                                              agCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCaaaa      <  1:1417020/69‑1 (MQ=255)                                              agCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCaaaa      <  1:1474478/69‑1 (MQ=255)                                              agCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCaaaa      <  1:2457768/69‑1 (MQ=255)                                              agCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCaaaa      <  1:307447/69‑1 (MQ=255)                                              agCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCaaaa      <  1:2436243/69‑1 (MQ=255)                                              agCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCaaaa      <  1:2410078/69‑1 (MQ=255)                                              agCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCaaaa      <  1:3425575/69‑1 (MQ=255)                                                ccAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCaaaa      <  1:1293577/67‑1 (MQ=255)                                                 cAACGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:2079986/44‑1 (MQ=38)                                                 cAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCGAATCGCTTTAAAAAGa   >  1:1796625/1‑69 (MQ=255)                                                 cAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCAAATCGCTTTAAAAAGa   >  1:1130899/1‑69 (MQ=255)                                                 cAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCAAATCGCTTTAAAAAGa   >  1:2576754/1‑69 (MQ=255)                                                 cAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCAAATCGCTTTAAAAAGa   >  1:1642719/1‑69 (MQ=255)                                                 cAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCAAATCGCTTTAAAAAGa   >  1:2745226/1‑69 (MQ=255)                                                 cAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCAAATCGCTTTAAAAAGa   >  1:299675/1‑69 (MQ=255)                                                 cAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:1622974/44‑1 (MQ=255)                                                  aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:585433/1‑69 (MQ=255)                                                  aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:2699936/1‑69 (MQ=255)                                                  aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:2798880/1‑69 (MQ=255)                                                  aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:2530902/1‑69 (MQ=255)                                                  aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:1313507/1‑69 (MQ=255)                                                  aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:83008/1‑69 (MQ=255)                                                  aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:807654/1‑69 (MQ=255)                                                  aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:793250/1‑69 (MQ=255)                                                  aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:709589/1‑69 (MQ=255)                                                  aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:601121/1‑69 (MQ=255)                                                  aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:1979483/1‑69 (MQ=255)                                                  aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:1558652/1‑69 (MQ=255)                                                  aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:533103/1‑69 (MQ=255)                                                  aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:1562002/1‑69 (MQ=255)                                                  aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:483524/1‑69 (MQ=255)                                                  aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:3014563/1‑69 (MQ=255)                                                  aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:3074580/1‑69 (MQ=255)                                                  aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:178617/1‑69 (MQ=255)                                                  aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:3177114/1‑69 (MQ=255)                                                  aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:2254404/1‑69 (MQ=255)                                                  aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:1878450/1‑69 (MQ=255)                                                   aaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:2464888/42‑1 (MQ=255)                                                   aaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:1164526/42‑1 (MQ=255)                                                          gCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCaaaa      <  1:164067/57‑1 (MQ=255)                                                                 |                                                      GCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCAAATCGCTTTAAAAAGAA  >  minE/1576110‑1576227

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