Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,638,300	G→A	G220D (GGC→GAC)	pitA →	metal phosphate:H(+) symporter PitA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,638,300	0	G	A	94.7%	45.1 / NA	20	G220D (GGC→GAC)	pitA	metal phosphate:H(+) symporter PitA
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (8/10);  minor base C (1/0);  total (9/11)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

GACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGAAGCTGGCGCAACGC  >  NC_000913/3638055‑3638539                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                gaCCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTAATAAGTATTGTCGGTTCTCTGATCGTTTCCCCTAGTGTCGGCCTGGTGTGTGCTGGCGGTCTGATTTTCTTGCAGCGTCGCTACGGGAGCTGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTCCTGGACGCGTATTGCGCGGATCCTTTCCGCTATCgacgtgg                                                                                                                                                                                                                                            <  1:329566/251‑1 (MQ=255)         cGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCg                                                                                                                                                                                                                                    >  1:10503/1‑251 (MQ=255)                        cACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTCTTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTAATGGAGAGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCCGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGACGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGAAGTGGCGTTTTCGCACGGGGCGAACGAt                                                                                                                                                                                                                     >  1:382711/1‑251 (MQ=255)                                                        ggTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTAGCTACTGGAGCGGCAACCAGACACGCGCCCGTATCCACCCGGCCCCAGCGGAGCGTGAAAAGAAAGACGGCCAGCAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTggg                                                                                                                                                                                     >  1:206523/1‑250 (MQ=255)                                                                          ccTATTGTCGGCCTGGTGTTTGCTGGCGGTATGAATTTCTTGCTGCGTCACTACTGGAGCGGCACCACGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCCAGAAAAAGCCGCCGTTCTGGAAGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCAAGGCGCGGACGATGGTAAGAAAGGAATTGGTCTGGTTATGTTGGTATTGATTAGCGTCGCGcc                                                                                                                                                                   >  1:33894/1‑251 (MQ=255)                                                                                                         tGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCCCCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACTCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTTGCGACGCGCCAGCAGGTTTCGTGGTGAACATGAATg                                                                                                                                         <  1:380286/246‑1 (MQ=255)                                                                                                             tttCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCCCCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGGAAGAAAATGCCGCCGTTCTGGACGCGGAGTGGTCTGATTCGTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAGCGGTGGGTAGTAATGCAGGGGTCTGTTTATTTTTGTATTGATTGGCTTCTCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGc                                                                                                                                 <  1:224966/250‑1 (MQ=255)                                                                                                                          cGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGCCCCCAGCGGAGCGTGAAATGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCaccng                                                                                                                   <  1:348174/251‑1 (MQ=255)                                                                                                                                gggAGCTGCACCAAGAAACGCGCCCGTATCCACCTTACCCCAGCGGAGCGTTAAAAGAAAGACTGCAAGAAAAAGCCGCCATTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACgg                                                                                                                                                                                                                               <  1:409801/136‑1 (MQ=38)                                                                                                                                    gCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGACAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGcc                                                                                                                                                                   >  1:228762/1‑193 (MQ=255)                                                                                                                                      ggCACCAAGAAACGCGCCCGTATCCCCCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAgg                                                                                                                                                                                                          <  1:1866/152‑1 (MQ=255)                                                                                                                                              naaaCGCGCACGTATAAACCTGACCCCCGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGc                                                                                                                                                                           >  1:111568/2‑175 (MQ=39)                                                                                                                                                         gTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCAGCCGTTCTGGACGCGTATTGCGATGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCt                                                                                    >  1:372117/1‑251 (MQ=255)                                                                                                                                                                 gggataccAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTg                                                                             <  1:99372/244‑1 (MQ=255)                                                                                                                                                                          gcggaccnggANAAGANAGNCGGCAAGAAAAAGCCNCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTNNCNNNTNCGCACGgc                                                                                                                                                                                                                              <  2:313221/87‑1 (MQ=37)                                                                                                                                                                          gCGGAGCNTGANAAGANAGNCGGCAAGAAAAAGCCNCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTNNCNNNTNCGCACGgc                                                                                                                                                                                                                              <  2:313231/96‑1 (MQ=38)                                                                                                                                                                          gCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGCCGTGGCGTTTTCGCAAGgc                                                                                                                                                                                                                              >  1:313231/1‑96 (MQ=255)                                                                                                                                                                                                    gAAAAAGCCGCCGGTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTGTTCGCACGGCGCGAACGAGGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTattgatt                                                                                                                                                                              <  1:152730/118‑1 (MQ=255)                                                                                                                                                                                                                  tCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTTGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAATATGAATGCCACTGGCTGCGAGATCGCCCGTACTCGTGATGCCATCGACAACGTCTATGCTGACTTGGAGCAGCATCCTGCGCTGCTCAAACAgn                                                  <  1:49124/228‑2 (MQ=255)                                                                                                                                                                                                                  tCTGGACCNGTATNGNNNTGATCNTTTCCGCTANCGACGTGGCGTTTTCGCACNGCGCGan                                                                                                                                                                                                                         >  2:49124/1‑60 (MQ=17)                                                                                                                                                                                                                                            tCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGGCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGAAGCGGTcgcaacgn  <  1:358867/250‑2 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                GACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGAAGCTGGCGCAACGC  >  NC_000913/3638055‑3638539

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑