Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,796,598	T→G	42.7%	D571A (GAC→GCC)	PP_3354 ←	acyl‑CoA dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,796,598	0	T	G	42.7%	‑0.3 / 4.5	21	D571A (GAC→GCC)	PP_3354	acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (7/5);  new base G (5/4);  total (12/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 = 9.92e-01

TGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCAACGCCTGGATGCAGCGCGCGCTGCG  >  NC_002947/3796471‑3796722                                                                                                                                |                                                                                                                             tgtgTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGGCGATgggc                                                                                                                      >  2:193518/1‑134 (MQ=255)      tCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGATGGCCGCCAc                                                                                                                 >  2:256108/1‑136 (MQ=255)           ggACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGGCGCCGCGCTGGCCCCCCCCCGcc                                                                                                            >  1:174160/1‑136 (MQ=255)           ggACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCCCGTCGCAGGCCGCCCCCAGcc                                                                                                            >  1:276935/1‑136 (MQ=255)                gCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGGCGCGGCGAGGGCCGCCACCAGcccggcg                                                                                                       >  2:7237/1‑131 (MQ=255)                gCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGATGGCCGCCACCTGCCGGtgg                                                                                                       >  1:475916/1‑134 (MQ=255)                gCTGCATGCGCTACCGCCCCGCGGCACGCCAGCCGCAGCCGGTGTTCATGCCGACCTCGCACGGGGGTTTCAGAACGATTCGCCAATGAATGGCAACGGCGTGGCGGCCGCGTCTATGGCCGCCACCTGCCGGTCg                                                                                                       <  1:445544/136‑1 (MQ=255)                       gcgcTACCGCCCCTCGGCACGCCAGCCGCAGCCGGGGTTCATGACGACCGCGCCCGGGGGTTTCAGAACGATGCGGCAATGAATGGCAACGGCGTGGCGGCCGCGGCGAGGGCCGCCACCTGCCGGTCGAGGTCgg                                                                                                <  1:89789/136‑1 (MQ=255)                         gcTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGATGGCCGCCACCAGCCGGTCGAGGTCgggg                                                                                              >  1:309076/1‑135 (MQ=255)                               gCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGGGGCTGGGCGCCCCCGGCCGGTCGCGGGCGGGCCGCCgg                                                                                        >  1:385691/1‑136 (MQ=255)                                       gCAGGCCAGCCGCTGCCGGTGTTCTTGACGACCTCGCAGGTGGGTTTCAGAACGATGCGGCATTGAATGGAAACGGGGTGGCGGCCGGGTCGTTGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAg                                                                                <  2:443751/136‑1 (MQ=255)                                       gCACGCCAGCCGCCGCCGGTGTGCATGCCGACCGCGCAGGTGGGTTTCAGAAGGATTCGGCAATTATTGGCACGGGGGGGGGGGCGGCGTCGTGGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAg                                                                                <  2:342287/136‑1 (MQ=255)                                               gccgccgcggGTGTTAATTCCTCCCTCGCAGGTGGGTTCAAGAACGATCCGGCAATTATTGCCACGGGCGGGGCGGCCGCGGCGAGGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAAcc                                                                        <  2:85455/136‑1 (MQ=255)                                                ccgcagccgGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGATGGCCCCCCCCTGCCGGTCGAGGTCGGGCAGCACGTAGCCGAGGAAGAAccc                                                                       >  2:53992/1‑136 (MQ=255)                                                      ccgGTGTTCATGACGCCCTCGCACGTGGGTTTCAGACCTTTTCGGAAATGAATGGCACCGGCGGGGGGGCGGCGGCTATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCt                                                                 <  1:256108/136‑1 (MQ=255)                                                                         cGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGATTGCCGCCCCCAGCCCGTCGCGGGCGGGCCGCAGGTGGCCGCGGAAGCACCCCACGCGTTCCTGCGTTTACCCgcgc                                              >  1:108381/1‑136 (MQ=255)                                                                         cGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCCCGTCGATGGCCCCCCCCAGCCGGTCGAGGTCGGGCCGCAGGTAGCCCAGGAAGAACCCAACGCATTACGGCTTTTACCCgcgc                                              >  1:132706/1‑136 (MQ=255)                                                                                 ggTTTCAGAACGATTCGGCAATCAATGGCAACGGCGGGGCGGCCGCGTCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGcgccgc                                      <  2:194984/136‑1 (MQ=255)                                                                                        caccGATTCGGCAATTAATGGCAACGGCGGGGCGGCCGCGGCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACt                               <  1:111903/133‑1 (MQ=255)                                                                                          aCGATTCGGCATTCAATGCCAACGGGGGGGCGGCCGCGTCTTTGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTATCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTgg                             <  1:47826/136‑1 (MQ=255)                                                                                                                     tGGCGGCCGCGGCGATGGCCGCCCCCTGCCGGTCGCGGTCGGGCAGCAGGTAGCCGCGGAAGAACCCAACGCGTTGCTGCTTTTCCCGGCGCTGCGCCGCCCCCCCTGGCAACGCCGGGATGCAGCGCGCGCTGCg  >  1:464614/1‑136 (MQ=255)                                                                                                                                |                                                                                                                             TGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCAACGCCTGGATGCAGCGCGCGCTGCG  >  NC_002947/3796471‑3796722

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