Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,935,790	C→G	15.4%	R17G (CGC→GGC)	PP_1736 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,935,790	0	C	G	15.4%	62.7 / 5.8	26	R17G (CGC→GGC)	PP_1736	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (11/11);  new base G (2/2);  total (13/13)
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.46e-01

ACCTGTGGGCGCACCGCCGGTTCAAGTCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGCGCTGAAAACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCAGTGGTGCGGTGCACTTCAATG  >  NC_002947/1935641‑1935934                                                                                                                                                      |                                                                                                                                                 aCCTGTGGGCGCACCGCCGGTTCAAGTCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGCGCTGAAAACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGccc                                                                                                                                                  >  2:146729/1‑150 (MQ=255)                                  ccGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGCGCTGAAAACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCa                                                                                                                 >  1:592795/1‑150 (MQ=255)                                  ccGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGCGCTGAAAACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTCGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCa                                                                                                                 >  1:476802/1‑150 (MQ=255)                                   cGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGCGCTGAAAACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCAt                                                                                                                <  1:382352/150‑1 (MQ=255)                                    gAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGCGCTGAAAACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCGTACCAGGTCGGTGTGCTGGCGGGCATc                                                                                                               <  1:579871/150‑1 (MQ=255)                                               gcgGCTGTATGACAAGAAGAAGCGCTGAAAACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGAGGCGGCGCGGGCgccgccgcc                                                                                                                                       <  2:291860/115‑4 (MQ=255)                                               gcgGCTGTATGACAAGAAGAAGCGCTGAAAACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGAGGCGGCGCGGGCgccgccgcc                                                                                                                                       >  1:291860/1‑112 (MQ=255)                                                         gACAAGAAGAAGCGCTGAAAACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCGGGCgccgccgcca                                                                                                                                      >  2:463545/1‑102 (MQ=255)                                                         gACAAGAAGAAGCGCTGAAAACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCGGGCgccgccgcca                                                                                                                                      <  1:463545/106‑5 (MQ=255)                                                                        tGAAAACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCtttt                                                                           >  2:756461/1‑150 (MQ=255)                                                                             acacATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGt                                                                      <  1:146729/150‑1 (MQ=255)                                                                               acaTCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCgg                                                                                      >  1:290209/1‑132 (MQ=255)                                                                               acaTCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCgg                                                                                      <  2:290209/132‑1 (MQ=255)                                                                                    tACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCAtcgtc                                                               <  2:443492/150‑1 (MQ=255)                                                                                        ccAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCa                                                           >  1:842240/1‑150 (MQ=255)                                                                                          aggaaggaCGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCAcc                                                         >  1:616878/1‑150 (MQ=255)                                                                                                gaCGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTccgcc                                                   >  1:212333/1‑150 (MQ=255)                                                                                                    ccATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCgggg                                               <  2:774393/150‑1 (MQ=255)                                                                                                      aTGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGcc                                             <  2:842240/150‑1 (MQ=255)                                                                                                         gCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATc                                          >  1:522090/1‑150 (MQ=255)                                                                                                           accaccCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCa                                            <  2:476922/146‑1 (MQ=255)                                                                                                           accaccCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCa                                            >  1:476922/1‑146 (MQ=255)                                                                                                               cccccGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGcc                                    <  2:592795/150‑1 (MQ=255)                                                                                                                                 aTCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCaa                                        <  2:186647/128‑1 (MQ=255)                                                                                                                                 aTCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCaa                                        >  1:186647/1‑128 (MQ=255)                                                                                                                                       aatgccggcgtctCCCGCGCTGCCTACTAGTTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCAGTggtgcggtg            <  2:476802/145‑1 (MQ=255)                                                                                                                                                 gcgcCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCAGTGGTGCGGTGCACTTCAATg  >  1:593852/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                 ACCTGTGGGCGCACCGCCGGTTCAAGTCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGCGCTGAAAACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCAGTGGTGCGGTGCACTTCAATG  >  NC_002947/1935641‑1935934

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