Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,469,680	G→A	26.6%	R202K (AGG→AAG)	queF →	NADPH‑dependent 7‑cyano‑7‑deazaguanine reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,469,680	0	G	A	26.6%	22.0 / 4.2	15	R202K (AGG→AAG)	queF	NADPH‑dependent 7‑cyano‑7‑deazaguanine reductase
Reads supporting (aligned to +/- strand):  ref base G (2/9);  new base A (2/2);  total (4/11)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTGCGCTGCAACCCGGAGCGGGTGGTGGAAGAAACCCTGCACAGCCATCTGCTGAAGTCCAACTGCCCGGTCACCGGCCAGCCGGATTGGGGCAGCGTGGTGGTCCAGTACAAGGGCAGGGCGCTGGACCATGCCAGCCTGCTGACCTATCTGATCAGCTTCCGCCAGCATGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCA  >  NC_002947/2469560‑2469800                                                                                                                         |                                                                                                                         tGTTGCGCTGCAACCCGGAGCGGGTGGTGGAAGAAACCCTGCACAGCCATCTGCTGAAGTCCAACTGCCCGGTCACCGGCCAGCCGGATTGGGGCAGCGTGGTGGTCCAGTACAAGGGCAGGGCGCTGGACCATGc                                                                                                           >  2:97198/1‑136 (MQ=255)              cccGGAGCGGGTGGTGGAAGAAACCCTGCACAGCCATCTGCTGAAGTCCAACTGCCCGGTCACCGGCCAGCCGGATTGGGGCAGCGTGGTGGTCCAGTACAAGGGCAGGGCGCTGGACCATGCCAGCCTGCTGAcc                                                                                              <  2:374926/136‑1 (MQ=255)                          ggtggAAGAAACCCTGCACAGCCATCTGCTGAAGTCCAACTGCCCGGTCACCGGCCAGCCGGATTGGGGCAGCGTGGTGGTCCAGTACAAGGGCAAGGCGCTGGACCATGCCAGCCTGCTGACCTATCTGATCAGc                                                                                  >  2:189225/1‑136 (MQ=255)                          ggtggAAGAAACCCTGCACAGCCATCTGCTGAAGTCCAACTGCCCGGTCACCGGCCAGCCGGATTGGGGCAGCGTGGTGGTCCAGTACAAGGGCAAGGCGCTGGACCATCCCAGCCTGCAGACCTATCTCATCAGc                                                                                  >  2:189227/1‑136 (MQ=255)                                      ccTGCACAGCCATCTGCTGAAGTCCTACTGCCCGGTCACCGGCCAGCCGGATTGGGGCAGCGTGGTGGTCCAGTACAAGGGCAAGGCGCTGGACCATGCCAGCCTGCTGACCTATCTGATCAGCTTCCGCCAGCAt                                                                      <  1:189227/136‑1 (MQ=255)                                      ccTGCACAGCCATCTGCTGAAGTCCAACTGCCCGGTCACCGGCCAGCCGGATTGGGGCAGCGTGGTGGTCCAGTACAAGGGCAAGGCGCTGGACCATGCCAGCCTGCTGACCTATCTGATCAGCTTCCGCCAGCAt                                                                      <  1:189225/136‑1 (MQ=255)                                        ggCACAGCCATCTGCTGAAGTCCAACTGCCCGGTCACCGGCCAGCCGGATTGGGGCAGCGTGGTGGTCCAGTACAAGGGCAGGGCGCTGGACCATGCCAGCCTGCTGACCTATCTGATCAGCTTCCGCCAGCATGc                                                                    <  1:354890/135‑1 (MQ=255)                                                 aTCTGCTGAAGTCCAACTGCCCGGTCACCGGCCAGCCGGATTGTGGCAGCGTGGTGGTCCAGTACAAGGGCAGGGCGCTGGACCATGCCAGCCTGCTGACCTATCTGATCAGCTTCCGCCAGCATGCCGACTTCCa                                                           <  1:363349/136‑1 (MQ=255)                                                 aTCTGCTGAAGTCCAACTGCCCGGTCACCGGCCAGCCGGATTGGGGCAGCGTGGTGGTCCAGTACAAGGGCAGGGCGCTGGACCATGCCAGCCTGCTGACCTATCTGATCAGCTTCCGCCAGCATGCCGACTTCCa                                                           <  1:451718/136‑1 (MQ=255)                                                 aTCTGCTGAAGTACAACTGCCCGGTCACCGGCCAGCCGGATTGGGGCAGCGTGGTGGTCCAGTACAAGGGCAGGGCGCTGGACCATGCCAGCCTGCTGACCTATCTGATCAGCTTCCGCCAGCATGCCGACTTCCa                                                           <  1:355321/136‑1 (MQ=255)                                                    tgctgAAGTCCAACTGCGCTGTCACCGTCCAGCCGGATTGGGGCAGCGTGGTGGTCCAGTACAAGGGCAGGGCGCTGGACCATGCCAGCCTGCTGACCTATCTGATCAGCTTCCGACAGCATGCCGACTTCCATGa                                                        <  1:50413/136‑1 (MQ=255)                                                    tgctgAAGTCCAACTGCCCGGTCACCGGCCAGCCGGATTGGGGCAGCGTGGTGGTCCAGTACAAGGGCAGGGCGCTGGACCATGCCAGCCTGCTGACCTATCTGATCAGCTTCCGCCAGCATGCCGACTTCCATGa                                                        <  2:471696/136‑1 (MQ=255)                                                                        gTCACCGGCCAGCCGGATTGGGGCAGCGTGGTGGTCCAGTACAAGGGCAGGGCGCTGGACCATGCCAGCCTGCTGACCTATCTGATCAGCTTCCGCCAGCATGCCGACTTCCATGAGCAGTGTGTGGAGCGGAttt                                    <  1:95273/136‑1 (MQ=255)                                                                                                    tggtggTCCAGTACAAGGGCAGGGCGCTGGACCATGCCAGCCTGCTGACCTATCTGATCAGCTTCCGCCAGCATGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGcc        >  2:455914/1‑136 (MQ=255)                                                                                                          tCCAGTACAAGGGCAGGGCGCTGGACCATGCCAGCCTGCTGACCTATCTGATCAGCTTCCGCCAGCATGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCa  <  1:300326/136‑1 (MQ=255)                                                                                                                         |                                                                                                                         TGTTGCGCTGCAACCCGGAGCGGGTGGTGGAAGAAACCCTGCACAGCCATCTGCTGAAGTCCAACTGCCCGGTCACCGGCCAGCCGGATTGGGGCAGCGTGGTGGTCCAGTACAAGGGCAGGGCGCTGGACCATGCCAGCCTGCTGACCTATCTGATCAGCTTCCGCCAGCATGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCA  >  NC_002947/2469560‑2469800

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