Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,027,290	G→A	G57D (GGC→GAC)	ghxQ →	guanine/hypoxanthine transporter GhxQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,027,290	0	G	A	94.7%	40.5 / ‑3.5	19	G57D (GGC→GAC)	ghxQ	guanine/hypoxanthine transporter GhxQ
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (7/11);  total (7/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.68e-01

CAACTAATAATAAAACTTTAACATCCTCGTGAGGACATCATTATGTCTGGAGACATCCTACAAACACCGGACGCACCAAAGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTGCGCGATTTCCTTGACGGCGTTTACCGCATTCAGTCTGGTACTCGGGCAACAAATTAGCGTTCCTGTCGCACTGGGCGCGGTATTTCTGATGGGCGTCATCTTCACCGCCATTTCCGTAACCGGT  >  NC_000913/3027079‑3027522                                                                                                                                                                                                                    |                                                                                                                                                                                                                                         cAACTAATAATAAAACGTTAACATCCTCGTGAGGACATCATTATGTCTGGAGACATCCTACAAACACCGGACGCACCAAAGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTGATTCCGTTATCGTCGTTCCGGTAATGCTGGACAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCAc                                                                                                                                                                                                   <  1:112392/251‑1 (MQ=255)                             gTGAGGACATCATTATGTCTGGAGACATCCTACAAACTCCGGACGCACCAAAGCCGCAGGGCGCGCTGGATAATTATGTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTGCTGGCCATGGTTTATTCCGTTATCGGCGTTCCGGGAAGGCTCGACAAAGCAGGGTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTg                                                                                                                                                                       <  1:245293/251‑1 (MQ=255)                                                       ccttcaatatctgcGGACGCACCAAAGCCACAGTGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGACAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGTGCTTCGGCTCGTTGCTGATGGGATTATGAGCTAATTTGcc                                                                                                                                             <  1:90655/238‑1 (MQ=255)                                                                            ccAAAGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGACAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTGCGCGAt                                                                                                                        >  1:1604/1‑251 (MQ=255)                                                                                             gcTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGACAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTGCGCGATTTCCTTGACGGCGttt                                                                                                        >  1:109700/1‑250 (MQ=255)                                                                                                            tAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGACAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGc                                                                                                                                                                             <  1:7090/166‑1 (MQ=255)                                                                                                                                         gTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGACAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCTGCTCGTTGCTGATGGGATTATgggt                                                                                                                                                      <  1:71844/160‑2 (MQ=255)                                                                                                                                                 gTACTGGCTGGCTTAACGACCTTTCTGGCCCTGGTTTATTACGTTATCGTCGTTCGGGGAATGCTGAACAAAGCAGGTTTTCTTACCGAAGCTGTGTTTGTTGCCACCTGTTTGGTAGAGGGCTTCGCCTCGTTGCTGATGGGATTCTGGGCTAATTTGACAATGGCGATTGGTTGCGCGATTTCCTTGACGGCGTTTACCGCATTCAGTCTGGTACTCGGGCAACAAATTAGCGTTCCTGTCGCACTggg                                                   >  1:177198/1‑251 (MQ=255)                                                                                                                                                  tACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCAGTTATCGTCGTTCCGGGAATGCTGGACAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCAACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTGCGCGATTTCCTTGACGGCGTTTACCGCATTCAGTCTGGTACTCGGGCAACAAATTAGCGTTCCTGTCGc                                                         >  1:247671/1‑244 (MQ=255)                                                                                                                                                   aCTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGCATGCTGGACAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTGCGCGATTTCCTTGACGGCGTTTACCGCATTCAGTCTGGTACTCGGGc                                                                               >  1:141115/1‑221 (MQ=255)                                                                                                                                                   aCTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGACAAAGCAGGTTTTCCGCCCGCAGCTGTGTTTTTTGCCACCTGTCTGGTcg                                                                                                                                                                                        <  1:91453/116‑1 (MQ=255)                                                                                                                                                    ctggctggCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGACAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTGCGCGATTTCCTTGACGGCGTTTACCGCATTCAGTCTGGTACTCGGGCAACAAATTAGCGTTCCTGTCGCACTGGgcgc                                                >  1:18342/1‑251 (MQ=255)                                                                                                                                                          ggctTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGACAAAGCAGGTTTTCCTCCCGCAGCTGTGTttgttg                                                                                                                                                                                                       <  1:200682/94‑1 (MQ=255)                                                                                                                                                               aaCGTCCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGACAAAGCAGGTTTtc                                                                                                                                                                                                                            <  1:95330/68‑1 (MQ=255)                                                                                                                                                                 cGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGtn                                                                                                                                                                                              <  1:76864/96‑2 (MQ=255)                                                                                                                                                                         cTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGACAAAGAAGGTTTTCCTCCCGCAGCTGTCTTTGTTGCCACCTGTCTGGTCCCGGGCTTCGGCTCGTTGATGATGGGATTATGGGCTAATTTGCCAATGGCGGTTGGTTGCGCGATTTCCTTGACGGCGTTTAc                                                                                                      >  1:139815/1‑176 (MQ=255)                                                                                                                                                                                                   cgtTCGGGGAATGCTGGACAAAGCAGGTTTTCCTCCCGCAGCTGGGTTTGTTGCCACCTGTCTGGCCGCGGGTTTCGGCTCGTTGCTGAGGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTGCGCGATTTCCTTGACGGCGTTTACCGCATTCAGTCTGGTACTCGGGCAACAAATGAGCGTTCCTGTCGCACTGGGCGCGGTATTTCTGATGGGCGCCAGCTTCACCGCCATTTCCGTAACCggn  <  1:76449/250‑2 (MQ=255)                                                                                                                                                                                                    gtTCCGGGAATGCTGGACAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTGCGCGATTTCCTTGACGGCGTTTTCCGCag                                                                                                 <  1:73137/154‑2 (MQ=255)                                                                                                                                                                                                              tGCTGGACAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTGCGCGATTTCCTTGACGGCGTTTACCGCATTCAGTCTGGTACTCGGGCAACAAATTAg                                                                     <  1:114381/172‑1 (MQ=255)                                                                                                                                                                                                                    |                                                                                                                                                                                                                                         CAACTAATAATAAAACTTTAACATCCTCGTGAGGACATCATTATGTCTGGAGACATCCTACAAACACCGGACGCACCAAAGCCACAGGGCGCGCTGGATAATTATTTTAAAATTACCGCTCGTGGCAGTACCGTTCGTCAGGAAGTACTGGCTGGCTTAACGACCTTTCTGGCCATGGTTTATTCCGTTATCGTCGTTCCGGGAATGCTGGGCAAAGCAGGTTTTCCTCCCGCAGCTGTGTTTGTTGCCACCTGTCTGGTCGCGGGCTTCGGCTCGTTGCTGATGGGATTATGGGCTAATTTGCCAATGGCGATTGGTTGCGCGATTTCCTTGACGGCGTTTACCGCATTCAGTCTGGTACTCGGGCAACAAATTAGCGTTCCTGTCGCACTGGGCGCGGTATTTCTGATGGGCGTCATCTTCACCGCCATTTCCGTAACCGGT  >  NC_000913/3027079‑3027522

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑