Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,658,742	G→T	11.1%	P1006P (CCC→CCA)	PP_1449 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,658,742	0	G	T	11.1%	117.5 / 3.0	36	P1006P (CCC→CCA)	PP_1449	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (21/11);  new base T (2/2);  total (23/13)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.10e-01

CCTGGGGTAGGGCGAGCGACCCTGCACTGTGCATGACCACGCTGCCTTCGCCGCCATCGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGACCCCCCTTTGCCGGCGGCGCGTACATTCAGGTCTTCGCCGGTGCTGGTGTCCAGGCGCGCATCGCCACCATAGGCCAGGCGTTGAACG  >  NC_002947/1658612‑1658868                                                                                                                                   |                                                                                                                               ccTGGGGTAGGGCGAGCGACCCTGCACTGTGCATGACCACGCTGCCTTCGCCGCCATCGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCAc                                                                                                                           <  1:621736/136‑1 (MQ=255)              gagcgaCCCTGCACTGTGCATGACCACGCTGCCTTCGCCGCCATCGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCg                                                                                                              >  2:418318/1‑136 (MQ=255)               agcgaCCCTGCACTGTGCATGACCACGCTGCCTTCGCCGCCATCGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCgg                                                                                                             >  1:639004/1‑136 (MQ=255)               agcgaCCCTGCACTGTGCATGACCACGCTGCCTTCGCCGCCATCGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCgg                                                                                                             <  1:333042/136‑1 (MQ=255)                                 aTGACCACGCTGCCTTCGCCGCCATCGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCg                                                                                           >  1:543364/1‑136 (MQ=255)                                 aTGACCACGCTGCCTTCGCCGCCATCGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCg                                                                                           >  1:517216/1‑136 (MQ=255)                                 aTGACCACGCTGCCTTCGCCGCCATCGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCg                                                                                           >  1:516696/1‑136 (MQ=255)                                 aTGACCACGCTGCCTTCGCCGCCATCGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCg                                                                                           >  1:429651/1‑136 (MQ=255)                                 aTGACCACGCTGCCTTCGCCGCCATCGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCg                                                                                           >  1:392340/1‑136 (MQ=255)                                 aTGACCACGCTGCCTTCGCCGCCATCGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCg                                                                                           >  2:475510/1‑136 (MQ=255)                                 aTGACCACGCTGCCTTCGCCGCCATCGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCg                                                                                           >  1:348163/1‑136 (MQ=255)                                 aTGACCACGCTGCCTTCGCCGCCATCGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCg                                                                                           >  2:56430/1‑136 (MQ=255)                                 aTGACCACGCTGCCTTCGCCGCCATCGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCg                                                                                           >  1:181718/1‑136 (MQ=255)                                 aTGACAACGCTGCCTTCGCCGCCATCGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACACTGTTGACCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGcc                                                                                           >  1:545783/1‑135 (MQ=255)                                      cACGCTGCCTTCGCCGCCATCGATGTGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCTGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGAcccc                                                                                      >  1:15963/1‑136 (MQ=255)                                      cACGCTGCCTTCGCCGCCATCGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCTGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGAcccc                                                                                      >  1:523633/1‑136 (MQ=255)                                      cACGCTGCCTTCGCCGCCATCGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGAcccc                                                                                      >  1:666719/1‑136 (MQ=255)                                      cACGCTGCCTTCGCCGCCATCGAAGCGGGCCCCTTCGTAGCGGCCATCGCAGCCCAGCTGTACCTGGATACCCTGTTGCCAGAAAAGGCTGCCGGGCACGCCGCTGCTGGCGGTATCCTTTTTATCCAGCGAcccc                                                                                      >  1:419613/1‑136 (MQ=255)                                        cGCTGCCTTCGCCGCCATCGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGAcccccc                                                                                    <  2:392340/136‑1 (MQ=255)                                           tGCCTTCGCCGCCATCGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGACCCCCCttt                                                                                 >  2:473288/1‑136 (MQ=255)                                           tGCCTTCGCCGCCATCGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGACCCCCCttt                                                                                 >  1:138126/1‑136 (MQ=255)                                                    cgccATCGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGACCCCCCTTTGCCggcggc                                                                        <  1:475510/136‑1 (MQ=255)                                                          cGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGA‑CCCCCTTTGCCGGCGGCGCGTACa                                                                 <  2:314834/136‑1 (MQ=255)                                                           gATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGACCCCCCTTTGCCGGCGGCGCGTACa                                                                 <  1:56430/136‑1 (MQ=255)                                                                                 ccATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGACCCCCCTTTGCCGGCGGCGCGTACATTCAGGTCTTCGCCGGTGCTgg                                           >  2:622532/1‑136 (MQ=255)                                                                                           ccAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGACCCCCCTTTGCCGGCGGCGCGTACATTCAGGTCTTCGCCGGTGCTGGTGTCCAGgcg                                 <  2:138126/136‑1 (MQ=255)                                                                                               cTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGACCCCCCTTTGCCGGCGGCGCGTACATTCAGGTCTTCGCCGGTGCTGGTGTCCAGGCGCGCa                             <  2:516696/136‑1 (MQ=255)                                                                                                tGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCTGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGACCCCCCTTTGCCGGCGGCGCGTACATTCAGGTCTTCGCCGGTGCTGGTGTCCAGGCGCGCAt                            <  2:15963/136‑1 (MQ=255)                                                                                                tGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCTGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGACCCCCCTTTGCCGGCGGCGCGTACATTCAGGTCTTCGCCGGTGCTGGTGTCCAGGCGCGCAt                            <  2:523633/136‑1 (MQ=255)                                                                                                           aCCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGACCCCCCTTTGCCGGCGGCGCGTACATTCAGGTCTTCGCCGGTGCTGGTGTCCAGGCGCGCATCGCCACCATAg                 <  1:418318/136‑1 (MQ=255)                                                                                                               tGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGACCCCCCTTTGCCGGCGGCGCGTACATTCAGGTCTTCGCCGGTGCTGGTGTCCAGGCGCGCATCGCCACCATAGGCCa             >  1:370028/1‑136 (MQ=255)                                                                                                                  tGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGACCCCCCTTTGCCGGCGGCGCGTACATTCAGGTCTTCGCCGGTGCTGGTGTCCAGGCGCGCATCGCCACCATAGGCCAGGc          >  1:133332/1‑136 (MQ=255)                                                                                                                     ccGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGACCCCCCTTTGCcgg                                                                            >  2:232716/1‑67 (MQ=255)                                                                                                                     ccGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGACCCCCCTTTGCcgg                                                                            <  1:232716/67‑1 (MQ=255)                                                                                                                         aaaGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGACCCCCCTTTGCCGGCGGCGCGTACATTCAGGTCTTCGCCGGTGCTGGTGTCCAGGCGCGCATCGCCACCATAGGCCAGGCGTTGAAc   >  2:35589/1‑136 (MQ=255)                                                                                                                          aaGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGACCCCCCTTTGCCGGCGGCGCGTACATTCAGGTCTTCGCCGGTGCTGGTGTCCAGGCGCGCATCGCCACCATAGGCCAGGCGTTGAACg  <  2:348163/136‑1 (MQ=255)                                                                                                                                   |                                                                                                                               CCTGGGGTAGGGCGAGCGACCCTGCACTGTGCATGACCACGCTGCCTTCGCCGCCATCGATGCGGGCCCCTTCGTAGCGGCCATCGCTGCCCAGCTGTACCTGGATACCCTGTTGCCCGTAAAGGCTGCCGGGCACGGCGGTGCTGGCGGTATCCTTTTTATCCAGCGACCCCCCTTTGCCGGCGGCGCGTACATTCAGGTCTTCGCCGGTGCTGGTGTCCAGGCGCGCATCGCCACCATAGGCCAGGCGTTGAACG  >  NC_002947/1658612‑1658868

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: ‑