Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,935,789	C→G	15.4%	A16A (GCC→GCG)	PP_1736 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,935,789	0	C	G	15.4%	62.2 / 4.1	26	A16A (GCC→GCG)	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

TACCTGTGGGCGCACCGCCGGTTCAAGTCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGCGCTGAAAACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCAGTGGTGCGGTGCACTTCAATG  >  NC_002947/1935640‑1935934                                                                                                                                                      |                                                                                                                                                  tACCTGTGGGCGCACCGCCGGTTCAAGTCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGCGCTGAAAACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGcc                                                                                                                                                   >  1:443492/1‑150 (MQ=255)  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)                                                                                                                                                      |                                                                                                                                                  TACCTGTGGGCGCACCGCCGGTTCAAGTCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGCGCTGAAAACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCAGTGGTGCGGTGCACTTCAATG  >  NC_002947/1935640‑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: ‑