Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	335,632	G→A	17.4%	T51I (ACC→ATC)	PP_5437 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	335,632	0	G	A	17.4%	50.9 / 4.6	23	T51I (ACC→ATC)	PP_5437	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (14/5);  new base A (2/2);  total (16/7)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

GATGATACACCGTTGCGGGGGTGGAGCATTAGCCTTGTGTTGTCTGGTATGGCCTCTTCGCGGGTGAACCCGCTCCTACACGGGGGCGGGTAAGCCGGGTTATCTGTGTAAGGGAATTCTGAACCATGGGGTCAAAAGCTTCAGGGTTTGGATGCTTCGGAAGTATCCTACGCGCGCCGCCGAAGTAACGCTGTTGTCGGGGAAGTGGCTCACGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAACCC  >  NC_002947/335502‑335758                                                                                                                                   |                                                                                                                               gatgatACACCGTTGCGGGGGTGGAGCATTAGCCTTGTGTTGTCTGGTATGGCCTCTTCGCGGGTGAACCCGCTCCTACACGGGGGCGGGTAAGCCGGGTTATCTGTGTAAGGGAATTCTGAACCATGGGGTCaaa                                                                                                                           >  1:384463/1‑136 (MQ=255)                 gggggTGGAGCATTAGCCTTGTGTTGTCTGGTATGGCCTCTTCGCGGGTGAACCCGCTCCTACACGGGGGCGGGTAAGCCGGGTTATCTGTGTAAGGGAATTCTGAACCATGGGGTCAAAAGCTTCAGGGTTTGGa                                                                                                           >  1:514287/1‑136 (MQ=255)                   gggTGGAGCATTAGCCTTGTGTTGTCTGGTATGGCCTCTTCGCGGGTGAACCCGCTCCTACACGGGGGCGGGTAAGCCGGGTTATCTGTGTAAGGGAATTCTGAACCATGGGGTCAAAAGCTTCAGGGTTTGGAtt                                                                                                         >  1:242371/1‑135 (MQ=255)                   gggTGGAGCATTAGCCTTGTGTTGTCTGGTATGGCCTCTTCGCGGGTGAACCCGCTCCTACACGGGGGCGGGTAAGCCGGGTTATCTGTGTAAGGGAATTCTGAACCATGGGGTCAAAAGCTTCAGGGTTTGGATg                                                                                                         >  1:540016/1‑136 (MQ=255)                   gggTGGAGCATTAGCCTTGTGTTGTCTGGTATGGCCTCTTCGCGGGTGAACCCGCTCCTACACGGGGGCGGGTAAGCCGGGTTATCTGTGTAAGGGAATTCTGAACCATGGGGTCAAAAGCTTCAGGGTTTGGATg                                                                                                         >  2:350519/1‑136 (MQ=255)                                gCCTTGTGTTGTCTGGTATGGCCTCTTCGCGGGTGAACCCGCTCCTACACGGGGGCGGGTAAGCCGGGTTATCTGTGTAAGGGAATTCTGAACCATGGGGTCAAAAGCTTCAGGGTTTGGATGCTTCGGAAGTATc                                                                                            >  1:21240/1‑136 (MQ=255)                                        ttgtCTGGTATGGCCTCTTCGCGGGTGAACCCGCTCCTACACGGGGGCGGGTAAGCCGGGTTATCTGTGTAAGGGAATTCTGAACCATGGGGTCAAAAGCTTCAGGGTTTGGATGCTTCGGAAGTATCCTAcgcgc                                                                                    >  2:487822/1‑136 (MQ=255)                                        ttgtCTGGTATGGCCTCTTCGCGGGTGAACCCGCTCCTACACGGGGGCGGGTAAGCCGGGTTATCTGTGTAAGGGAATTCTGAACCATGGGGTCAAAAGCTTCAGGGTTTGGATGCTTCGGAAGTATCCTAcgcgc                                                                                    >  2:27060/1‑136 (MQ=255)                                        ttgtCTGGTATGGCCTCTTCGCGGGTGAACCCGCTCCTACACGGGGGCGGGTAAGCCGGGTTATCTGTGTAAGGGAATACTGAACCATGGGGTCAAAAGCTTCAGGGTTTGGATGCTTCGGAAGTATCCTAcgcgc                                                                                    >  1:223363/1‑136 (MQ=255)                                        ttgtCTGGTATGGCCTCTTCGCGGGTGAACCCGCTCCTACACGGGGGCGGGTAAGCCGGCTTATCTGTGTAAGGTAATTCTGAACCATGCGGTCAAAAGCATCAGCGTTTGGATGCTTCGGAAGTATCCTActcgc                                                                                    >  1:74548/1‑136 (MQ=39)                                                    gCCTCTTCGCGGGTGAACCCGCTCCTACACGGGGGCGGGTAAGCCGGGTTATCTGTGTAAGGGAATTCTGAACCATGGGGTCAAAAGCTTCAGGGTTTGGATGCTTCGGAAGTATCCTACGCGCGCCGCCGAAGTa                                                                        <  1:10905/136‑1 (MQ=255)                                                        ctTCGCGGGTGAACCCGCTCCTACACGGGGGCGGGTAAGCCGGGTTATCTGTGTAAGGGAATTCTGAACCATGGGATCAAAAGCTTCAGGGTTTGGATGCTTCGGAAGTATCCTACGCGCGCCGCCGAAGTAACGc                                                                    >  1:410418/1‑136 (MQ=255)                                                        ctTCGCGGGTGAACCCGCTCCTACACGGGGGCGGGTAAGCCGGGTTATCTGTGTAAGGGAATTCTGAACCATGGGATCAAAAGCTTCAGGGTTTGGATGCTTCGGAAGTATCCTACGCGCGCCGCCGAAGTAACGc                                                                    >  1:410419/1‑136 (MQ=255)                                                                      ccGCTCCTACACGGGGGCGGGTAAGCCGGGTTATCTGTGTAAGGGAATTCTGAACAATGGGGTCATAAGCTTCAGGGTTTGGATGATTCGGAAGTATCCTACGCGCGCCGCCGAAGTAACGCTGTTGTCGGGGAAg                                                      >  1:117648/1‑136 (MQ=255)                                                                        gCTCCTACACGGCGGCGGGTAAGCCGGGTTATCTGTGTAAGGGAATTCTGAACCATGGGGTCAAAAGCTTCAGGGTTTGGATGCTTCGGAAGTATCCTACGCGCGCCGCCGAAGTAACGCTGTTGTCGGGGAAGTg                                                    <  2:54602/136‑1 (MQ=255)                                                                                       ccggTAAACCGGGTTATATTTGTAATGGAAATAAGAACCATTGGGTAAAAAGCTTCAGGGTTTGGATTCTTAGGAAGTATACTACGCGCGCCGCCGAAGTAAAGCAGTTGTCCGGGAAGAGGCTAACGGATACCGt                                     >  2:307365/3‑136 (MQ=255)                                                                                            aaGCCGGGTTATCTGTGTAAGGGAATTCTGAACCATGGGATCAAAAGCTTCAGGGTTTGGATGCTTCGGAAGTATCCTACGCGCGCCGCCGAAGTAACGCTGTTGTCGGGGAAGTGGCTCACGGATACCGTTTTCg                                <  2:410418/136‑1 (MQ=255)                                                                                            aaGCCGGGTTATCTGTGTAAGGGAATTCTGAACCATGGGATCAAAAGCTTCAGGGTTTGGATGCTTCGGAAGTATCCTACGCGCGCCGCCGAAGTAACGCTGTTGTCGGGGAAGTGGCTCACGGATACCGTTTTCg                                <  2:410419/136‑1 (MQ=255)                                                                                                                        tGAACCATGGGGTCAAAAGCTTCAGGGTTTGGATGCTTCGGAAGTATCCTACGCGCGCGGCCGAAGTAACGCTGTTGTCGGGGAAGTGGCTCACGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAAc    >  1:91132/1‑136 (MQ=255)                                                                                                                        tGAACCATGGGGTCAAAAGCTTCAGGGTTTGGATGCTTCGGAAGTATCCTACGCGCGCCGCCGAAGTAACGCTGTTGTCGGGGAAGTGGCTCACGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAAc    >  1:110682/1‑136 (MQ=255)                                                                                                                          aaCCATGGGGTCAAAAGCTTCAGGGTTTGGATGCTTCGGAAGTATCCTACGCGCGCCGCCGAAGTAACGCTGTTGTCGGGGAAGTGGCTCACGTATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAAccc  <  2:21240/136‑1 (MQ=255)                                                                                                                          aaCCATGGGGTCAAAAGCTTCAGGGTTTGGATGCTTCGGAAGTATCCTACGCGCGCCGCCGAAGTAACGCTGTTGTCGGGGAAGTGGCTCACGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAAccc  <  2:384463/136‑1 (MQ=255)                                                                                                                          aaCCATGGGGTCAAAAGCTTCAGGGTTTGGATGCTTCGGAAGTATCCTACGCGCGCCGCCGAAGTAACGCTGTTGTCGGGGAAGTGGCTCACGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAAccc  <  2:540016/136‑1 (MQ=255)                                                                                                                                   |                                                                                                                               GATGATACACCGTTGCGGGGGTGGAGCATTAGCCTTGTGTTGTCTGGTATGGCCTCTTCGCGGGTGAACCCGCTCCTACACGGGGGCGGGTAAGCCGGGTTATCTGTGTAAGGGAATTCTGAACCATGGGGTCAAAAGCTTCAGGGTTTGGATGCTTCGGAAGTATCCTACGCGCGCCGCCGAAGTAACGCTGTTGTCGGGGAAGTGGCTCACGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAACCC  >  NC_002947/335502‑335758

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