Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,748,106	G→A	35.3%	D291D (GAC→GAT)	etfA ←	electron transfer flavoprotein subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,748,106	0	G	A	35.3%	13.3 / 11.8	17	D291D (GAC→GAT)	etfA	electron transfer flavoprotein subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (5/6);  new base A (4/2);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

CTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGA  >  NC_002947/4747976‑4748236                                                                                                                                   |                                                                                                                                   cTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGc                                                                                                                               >  1:12233/1‑136 (MQ=255)                                          gggCCGGGTTCTTTATATGTGATGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGTACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCa                                                                                      <  1:409725/136‑1 (MQ=255)                                           ggCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCAc                                                                                     <  2:12233/136‑1 (MQ=255)                                                                  cAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCCCGATCACTTTGGAGTCTTTCATGc                                                              >  2:407941/1‑136 (MQ=255)                                                                  cAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGc                                                              >  2:412680/1‑136 (MQ=255)                                                                    ggCGGATCAGACCAGGTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCTGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTAGATCGCCACGATCACTTTGGAGTCTTTCATGCCg                                                            <  2:390423/136‑1 (MQ=255)                                                                    ggCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCg                                                            <  2:186486/136‑1 (MQ=255)                                                                    ggCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCg                                                            <  2:227756/136‑1 (MQ=255)                                                                                      ttCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATcgcg                                          <  2:291647/136‑1 (MQ=255)                                                                                                ggAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATAc                                >  2:297688/1‑136 (MQ=255)                                                                                                          tgAAGAGATCACCGACGAGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGAt                      >  2:237815/2‑136 (MQ=255)                                                                                                          cGAAGAGATCACCGACGAGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGAt                      <  1:230781/136‑1 (MQ=255)                                                                                                                    accgacgaGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATcg                                            <  2:294148/102‑1 (MQ=255)                                                                                                                    accgacgaGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATcg                                            >  1:294148/3‑104 (MQ=255)                                                                                                                       gACGAGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCtccggt         >  1:417825/1‑131 (MQ=255)                                                                                                                       gACGAGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCtccggt         >  1:414346/1‑131 (MQ=255)                                                                                                                              ccGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGa  >  2:180586/1‑136 (MQ=255)                                                                                                                                   |                                                                                                                                   CTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGA  >  NC_002947/4747976‑4748236

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