Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,576,150	G→A	34.7%	F298F (TTC→TTT)	xapB ←	xanthosine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,576,150	0	G	A	34.7%	50.2 / 37.4	69	F298F (TTC→TTT)	xapB	xanthosine transporter
Reads supporting (aligned to +/- strand):  ref base G (27/18);  new base A (16/8);  total (43/26)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

AAACAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCAAATCG  >  minE/1576092‑1576215                                                           |                                                                  aaaCAATCATCGACAGCAGCAGCAAAATAAATCCTGTTGTTGACGGATCGCCATAGGCGAAGAAGCCaa                                                         >  1:2017163/1‑69 (MQ=255) aaaCAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCaa                                                         >  1:1662496/1‑69 (MQ=255) aaaCAATCATCGACAGCAG‑AGCAAAATAAATCCGGTTGATGACGGATCGCCATAGGCGAAGAAGCCaaa                                                        >  1:1651805/1‑69 (MQ=255)  aaCAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCgaaga                                                               >  1:616987/1‑62 (MQ=255)  aaCAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCaaa                                                        >  1:3083308/1‑69 (MQ=255)  aaCAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCaaa                                                        >  1:2509049/1‑69 (MQ=255)  aaCAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCaaa                                                        >  1:2647190/1‑69 (MQ=255)  aaCAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCaaa                                                        >  1:2418398/1‑69 (MQ=255)  aaCAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCaaa                                                        >  1:1497838/1‑69 (MQ=255)  aaCAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCaaa                                                        >  1:3038282/1‑69 (MQ=255)  aaCAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCaaa                                                        >  1:708529/1‑69 (MQ=255)  aaCAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCaaa                                                        >  1:1258970/1‑69 (MQ=255)  aaCAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCaaa                                                        >  1:2343327/1‑69 (MQ=255)  aaCAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCaa                                                         >  1:2746233/1‑68 (MQ=255)  aaCAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCaa                                                         >  1:926801/1‑68 (MQ=255)    cAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAg                                                       <  1:1789536/68‑1 (MQ=255)    cAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAg                                                       <  1:1075703/68‑1 (MQ=255)    cAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAg                                                       <  1:1857692/68‑1 (MQ=255)    cAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAg                                                       <  1:2315432/68‑1 (MQ=255)    cAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAg                                                       <  1:1632828/68‑1 (MQ=255)    cAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAg                                                       <  1:2996501/68‑1 (MQ=255)    cAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAg                                                       <  1:655019/68‑1 (MQ=255)    cAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAg                                                       <  1:1242963/68‑1 (MQ=255)      atcatcGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAgcgc                                                    <  1:1260915/69‑1 (MQ=255)         atcGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGc                                                 >  1:2592499/1‑69 (MQ=255)                   gcagcaAAATAAATACGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                        >  1:936849/1‑68 (MQ=255)                    cagcaAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGccac                                      <  1:90322/69‑1 (MQ=255)                    cagcaAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGccac                                      <  1:53191/69‑1 (MQ=255)                    cagcaAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGccac                                      <  1:2103807/69‑1 (MQ=255)                    cagcaAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGccac                                      <  1:1258314/69‑1 (MQ=255)                        aaaaTAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATa                                  <  1:1398417/69‑1 (MQ=255)                        aaaaTAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATa                                  <  1:1341646/69‑1 (MQ=255)                        aaaaTAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATa                                  <  1:362461/69‑1 (MQ=255)                        aaaaTAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATa                                  <  1:1561827/69‑1 (MQ=255)                         aaataaatCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                        >  1:2193407/1‑62 (MQ=255)                         aaataaatCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                        >  1:1572764/1‑62 (MQ=255)                         aaataaatCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                        >  1:760793/1‑62 (MQ=255)                         aaataaatCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                        >  1:291942/1‑62 (MQ=255)                         aaataaatCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                        >  1:339684/1‑62 (MQ=255)                         aaataaatCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                        >  1:3247490/1‑62 (MQ=255)                         aaataaatCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                        >  1:3182486/1‑62 (MQ=255)                                          gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                <  1:1596031/66‑1 (MQ=255)                                          gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                <  1:186537/66‑1 (MQ=255)                                          gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                <  1:357314/66‑1 (MQ=255)                                          gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                <  1:3171839/66‑1 (MQ=255)                                          gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                <  1:1863293/66‑1 (MQ=255)                                          gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                <  1:1515360/66‑1 (MQ=255)                                          gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                <  1:2580671/66‑1 (MQ=255)                                          gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                <  1:1800702/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:1123922/1‑64 (MQ=255)                                                       aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga      >  1:766413/1‑64 (MQ=255)                                                       aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga      >  1:836047/1‑64 (MQ=255)                                                       aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga      >  1:901355/1‑64 (MQ=255)                                                       aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga      >  1:559691/1‑64 (MQ=255)                                                       aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga      >  1:471413/1‑64 (MQ=255)                                                       aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga      >  1:460814/1‑64 (MQ=255)                                                       aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga      >  1:45313/1‑64 (MQ=255)                                                       aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga      >  1:3321130/1‑64 (MQ=255)                                                       aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga      >  1:3048902/1‑64 (MQ=255)                                                       aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga      >  1:293713/1‑64 (MQ=255)                                                       aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga      >  1:2092086/1‑64 (MQ=255)                                                       aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga      >  1:1694877/1‑64 (MQ=255)                                                       aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga      >  1:1369850/1‑64 (MQ=255)                                                       aGGCAAAGAAGCCAAAGCGCAGCGTCAACGCCACCATACTCATCAGCATGACGGTTTTAATAccga      >  1:2616232/1‑64 (MQ=255)                                                        ggCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCAAATCg  >  1:447414/1‑69 (MQ=255)                                                        ggCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCAAATCg  >  1:971656/1‑69 (MQ=255)                                                           |                                                                  AAACAATCATCGACAGCAGCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCAAATCG  >  minE/1576092‑1576215

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