Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,903,222	A→G	23.6%	Q1244R (CAG→CGG)	PP_1703 →	assimilatory nitrate reductase/sulfite reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,903,222	0	A	G	23.6%	19.8 / 5.0	17	Q1244R (CAG→CGG)	PP_1703	assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base A (6/7);  new base G (2/2);  total (8/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.70e-01

TGGGCCCTGGTACTGGCATTGCACCGTTCCGTGCGTTTCTCGAGGAGCGAGAGGCACGGGGGGCAAAAGGCGACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTACCAGGAGCAGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGT  >  NC_002947/1903093‑1903371                                                                                                                                  |                                                                                                                                                      tgagCCCTGGTACAGGCATTGCACCGTTCCGTGCGTTTCTCGAGGAGCCAGAAGCACGGGGGGCAAAAGCCGACAAGTGGCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTATCAGGAGCAGTTGCAGGCCTGGCAAccgg                                                                                                                                   >  2:17504/4‑150 (MQ=255)        tGGTACTGGCATTGCACCGTTCCGTGCGTTTCTCGAGGAGCGAGAGGCACGGGGGGCAAAAGGCGACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTACCAGGAGCAGTTGCAGGCCTGGCAAGCGGCGGGGCa                                                                                                                            >  1:89100/1‑150 (MQ=255)                gCATTGCACCGTTCCGTGCGTTTCTCGAGGAGCGAGAGGCACGGGGGGCAAAAGGCGACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTACCAGGAGCGGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGc                                                                                                                    >  1:127393/1‑150 (MQ=255)                gCATTGCACCGTTCCGTGCGTTTCTCGAGGAGCGAGAGGCACGGGGGGCAAAAGGCGACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTACCAGGAGCGGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGc                                                                                                                    >  1:127364/1‑150 (MQ=255)                       aCCGTTCCGTGCGTTTCTCGAGGAGCGAGAGGCACGGGGGGCAAAAGGCGACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTACCAGGAGCAGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGacac                                                                                                             <  1:360064/150‑1 (MQ=255)                                    tttCTCGAGGAGCGAGAGGCACGGGGGGCAAAAGGCGACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTACCAGGAGCGGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCgg                                                                                                <  2:127393/150‑1 (MQ=255)                                    tttCTCGAGGAGCGAGAGGCACGGGGGGCAAAAGGCGACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTACCAGGAGCGGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCgg                                                                                                <  2:127364/150‑1 (MQ=255)                                        acGAGGAGCGAGAGGCACGGGGGGCAAAGCGGGACAACTGGCTGGTCAACGGCGAGCAGTACGCAGGGACCGATGTCTATTACGAGGAGCAGTTGCAGGCCTGGGAAGCGGCGGGACAACTGAGACTGGAAACGGCGTTCTCGAGGGAca                                                                                            >  2:184409/2‑149 (MQ=255)                                                        aCGGGGGGCAAAAGGCGACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTACCAGGAGCAGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTa                                                                            <  2:89100/150‑1 (MQ=255)                                                        aCGGGGGGCAAAAGGCGACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGAACGATTTCTATTACCAGGAGCAGTTg                                                                                                                                                   <  2:169756/79‑1 (MQ=255)                                                        aCGGGGGGCAAAAGGCGACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGAACGATTTCTATTACCAGGAGCAGTTg                                                                                                                                                   >  1:169756/1‑79 (MQ=255)                                                                        gACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTACCAGGAGCAGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATg                                                            <  2:373233/150‑1 (MQ=255)                                                                                gCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTACCAGGAGCAGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCa                                                    <  2:270507/150‑1 (MQ=255)                                                                                                      gCAGCGACCGATTTCTATTACCAGGAGCAGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGagaaga                                                                                <  1:355971/100‑1 (MQ=255)                                                                                                      gCAGCGACCGATTTCTATTACCAGGAGCAGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGagaaga                                                                                >  2:355971/1‑100 (MQ=255)                                                                                                                          ctaGGAGCAGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTAttt          >  2:368608/3‑150 (MQ=255)                                                                                                                              gAGCAGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTAt      <  1:368608/150‑1 (MQ=255)                                                                                                                                  aGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTAtgtgt  <  1:275338/150‑1 (MQ=255)                                                                                                                                  |                                                                                                                                                      TGGGCCCTGGTACTGGCATTGCACCGTTCCGTGCGTTTCTCGAGGAGCGAGAGGCACGGGGGGCAAAAGGCGACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTACCAGGAGCAGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGT  >  NC_002947/1903093‑1903371

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