Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	2,123,260	T→C	R41R (CGT→CGC)	PP_1886 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,123,260	0	T	C	89.7%	82.8 / 0.8	29	R41R (CGT→CGC)	PP_1886	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (14/12);  total (16/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.86e-01

CCTGCGTGTGGCCACATCACCATCCATTGGAGGTTGAGACCATGTCCAATTCACTGGAAGAACTGCTTGAGCTCATGAAAAAGGGCTCAGTAACGCTAGGCTGGGGGGCGGTTGCTGTCTATAGCCGTGACCGGCTCAACCGCTTGCTGGAGCAGCAATACCTGACGCGCTTCAAGGAAAACCGCTACCTCCCGCCTTTCAGTCACACGTTCGCAGGCGTGGGGGTGAAGGAGGCCGTGAGCGTTGAGGCGCTG  >  NC_002947/2123133‑2123386                                                                                                                                |                                                                                                                               ccTGCGTGTGGCCACATCACCATCCATTGGAGGTTGAGACCATGTCCAATTCACTGGAAGAACTGCTTGAGCTCATGAAAAAGGGCTCAGTAACGCTAGGCTGGGGGGCGGTTGCTGTCTATAGCCGCGACCGGCt                                                                                                                        <  2:25654/136‑1 (MQ=255)     cGTGTGGCCACATCACCATCCATTGGAGGTTGAGACCATGTCCAATTCACTGGAAGAACTGCTTGAGCTCATGAAAAAGGGCTCAGTAACGCTAGGCTGGGGGGCGGTTGCTGTCTATAGCCGCGACCGGCTCAAc                                                                                                                    >  2:235940/1‑136 (MQ=255)      gtgtGGCCACATCACCATCCATTGGAGGTTGAGACCATGTCCAATTCACTGGAAGAACTGCTTGAGCTCATGAAAAAGGGCTCAGTAACGCTAGGCGGGGGGGCGGTTGCTGTCTATAGCCGCGACCGGCTCAAcc                                                                                                                   <  2:274890/136‑1 (MQ=255)             cacaTCACCATCCATTGGAGGTTGAGACCATGTCCAATTCACTGGAAGAACTGCTTGAGCTCATGAAAAAGGGCTCAGTAACGCTAGGCTGGGGGGCGGTTGCTGTCTATAGCCGTGACCGGCTCAACCGCTTGCt                                                                                                            >  2:833654/1‑136 (MQ=255)              acaTCACCATCCATTGGAGGTTGAGACCATGTCCAATTCACTGGAAGAACTGCTTGAGCTCATGAAAAAGGGCTCAGTAACGCTAGGCTGGGGGGCGGTTGCTGTCTATAGCCGCGACCGGCTCAACCGCTTGCTg                                                                                                           >  1:26830/1‑136 (MQ=255)              acaTCACCATCCATTGGAGGTTGAGACCATGTCCAATTCACTGGAAGAACTGCTTGAGCTCATGAAAAAGGGCTCAGTAACGCTAGGCTGGGGGGCGGTTGCTGTCTATAGCCGCGACCGGCTCAACCGCTTGCTg                                                                                                           >  1:530554/1‑136 (MQ=255)              acaTCACCATCCATTGGAGGTTGAGACCATGTCCAATTCACTGGAAGAACTGCTTGAGCTCATGAAAAAGGGCTCAGTAACGCTAGGCTGGGGGGCGGTTGCTGTCTATAGCCGCGACCGGCTCAACCGCTTGCTg                                                                                                           <  1:592190/136‑1 (MQ=255)                  cACCATCCATTGGAGGTTGAGACCATGTCCAATTCACTGGAAGAACTGCTTGAGCTCATGAAAAAGGGCTCAGTAACGCTAGGCTGGGGGGCGGTTGCTGTCTATAGCCGCGACCGGCTCAACCGCTTGCTGGagc                                                                                                       >  2:460090/1‑136 (MQ=255)                          atTGGAGGTTGAGACCATGTCCAATTCACTGGAAGAACTGCTTGAGCTCATGAAAAAGGGCTCAGTAACGCTAGGCTGGGGGGCGGTTGCTGTCTATAGCCGCGACCGGCTCAACCGCTTGCTGGAGCAGCAATAc                                                                                               >  2:1290296/1‑136 (MQ=255)                                  ttGAGACCATGTCCAATTCACTGGAAGAACTGCTTGAGCTCATGAAAAAGGGCTCAGTAACGCTAGGCTGGGGGGCGGTTGCTGTCTATAGCCGCGACCGGCTCAACCGCTTGCTGGAGCAGCAATACCTGAcgcg                                                                                       >  1:1109937/1‑136 (MQ=255)                                         cATGTCCAATTCACTGGAAGAACTGCTTGAGCTCATGAAAACGGGCTCAGTAACGCTAGGCTGGGGGGCGGTTGCTGTCTATAGCCGCGACCGGCTCAACCGCTTGCTGGAGCAGCAATACCTGACGCGCTTCAAg                                                                                <  2:626068/136‑1 (MQ=255)                                                          aagaaCTGCTTGAGCTCATGAAAAAGGGCTCAGTAACGCTAGGCTGGGGGGCGGTTGCTGTCTATAGCCGCGACCGGCTCAACCGCTTGCTGGAGCAGCAATACCTGACGCGCTTCAAGGAAAACCGCTACCTccc                                                               >  1:6522/1‑136 (MQ=255)                                                          aagaaCTGCTTGAGCTCATGAAAAAGGGCTCAGTAACGCTAGGCTGGGGGGCGGTTGCTGTCTATAGCCGCGACCGGCTCAACCGCTTGCTGGAGCAGCAATACCTGACGCGCTTCAAGGAAAACCGCTACCTccc                                                               >  1:508961/1‑136 (MQ=255)                                                                 gtttGAGCTCATGAAAAAGGGCTCAGTAACGCTAGGCTGGGGGGCGGTTGCTGTCTATAGCCGCGACCGGCTCAACCGCTTGCTGGAGCAGCAATACCTGACGCGCTTCAAGGAAAACCGCTACCTCCCGCCTTTc                                                        <  2:717949/134‑1 (MQ=255)                                                                            tGAAAATGGGCTCAGTAACGCTAGGGTGGGGGGCGGTTGCTGTCTATAGCCGCGACCGGCTCAACCGCTTGCTGGAGCAGCAATACCTGACGCGCTTCAAGGAAAACCGCTACCTCCCGCCTTTCAGTCACACGtt                                             <  2:1051186/136‑1 (MQ=255)                                                                            tGAAAAAGGGCTCAGTAACGCTAGGCTGGGGGGCGGTTGCTGTCTATAGCCGCGACCGGCTCAACCGCTTGCTGGAGCAGCAATACCTGACGCGCTTCAAGGAAAACCGCTACCTCCCGCCTTTCAGTCACACGtt                                             <  2:1186544/136‑1 (MQ=255)                                                                                  aGGGCTCAGTAACGCTAGGCTGGGGGGCGGTTGCTGTCTATAGCCGCGACCGGCTCAACCGCTTGCTGGAGCAGCAATACCTGACGCGCTTCAAGGAAAACCGCTACCTCCCGCCTTTCAGTCACACGTTCGCAgg                                       <  1:679485/136‑1 (MQ=255)                                                                                         aGTAACGCTAGGCTGGGGGGCGGTTGCTGTCTATAGCCGCGACCGGCTCAACCGCTTGCTGGAGCAGCAATACCTGACGCGCTTCAAGGAAAACCGCTACCTCCCGCCTTTCAGTCACACGTTCGCAGGCGTgggg                                >  2:719587/1‑136 (MQ=255)                                                                                         aGTAACGCTAGGCTGGGGGGCGGTTGCTGTCTATAGCCGCGACCGGCTCAACCGCTTGCTGGAGCAGCAATACCTGACGCGCTCCAAGGAAAACCGCTACCTCCCGCCTTTCAGTCACACGTTCGCAGGCGTgggg                                >  2:764853/1‑136 (MQ=255)                                                                                         aGTAACGCTAGGCTGGGGGGCGGTTGCTGTCTATAGCCGCGACCGGCTCAACCGCTTGCTGGAGCAGCAATACCTGACGCACTTCAAGGAAAACCGCTACCTCCCGCCTTTCAGTCACACGTTCGCAGGCGTgggg                                >  2:802790/1‑136 (MQ=255)                                                                                            aaCGCTAGGCTGGGGGGCGGTTGCTGTCTATAGCCGTGACCGGCTCAACCGCTTGCTGGAGCAGCAATACCTGACGCGCTTCAAGGAAAACCGCTACCTCCCGCCTATCAGTCACACGTTCGCAGGCGTGGGGGTg                             >  1:1339153/1‑136 (MQ=255)                                                                                                    gcgggggggCGGTTGCTGTCTATAGCCGCGACCGGCTCAACCGCTTGCTGGAGCAGCAATACCTGACGCGCTTCAAGGAAAACCGCTACCTCCCGCCTTTCAGTCACACGTTCGCAGGCGTGGGGGTGAAGGAGGc                     <  1:1106895/133‑1 (MQ=255)                                                                                                    gCTGGGGGGCGGTTGCTGTCTATAGCCGCGACCGGCTCAACCGCTTGCTGGAGCAGCAATACCTGACGCGCTTCAAGGAAAACCGCTACCTCCCGCCTTTCAGTCACACGTTCGCAGGCGTGGGGGTGAAGGAGGc                     <  1:1197095/136‑1 (MQ=255)                                                                                                    gCTGGGGGGCGGTTGCTGGCTATAGCCGCGACCGGCTCAACCGCTTGCTGGAGCAGCAATACCTGACGCGCTTCAAGGAAAACCGCTACCTCCCGCCTTTCAGTCACACGTTCGCAGGCGTGGGGGTGAAGGAGGc                     <  1:572167/136‑1 (MQ=255)                                                                                                      tGGGGGGCGGTTGCTGTCTATAGCCGCGACCGGCTCAACCGCTTGCTGGAGCAGCAATACCTGACGCGCTTCAAGGAAAACCGCTACCTCCCGCCTTTCAGTCACACGTTCGCAGGCGTGGGGGTGAAGGAGGCCg                   >  2:1014221/1‑136 (MQ=255)                                                                                                           ggCGGTTGCTGTCTATAGCCGCGACCGGCTCAACCGCTTGCTGGAGCAGCAATACCTGACGCGCTTCAAGGAAAACCGCTACCTCCCGCCTTTCAGTCACACGTTCGCAGGCGTGGGGGTGAAGGAGGCCGTGAgg              >  2:872087/1‑135 (MQ=255)                                                                                                                ttGCTGTCTATAGCCGCGACCGGCTCAACCGCTTGCTGGAGCAGCAATACCTGACGCGCTTCAAGGAAAACCGCTACCTCCCGCCTTTCAGTCACACGTTCGCAGGCGTGGGGGTGAAGGAGGCCGTGAGCGTTGa         >  2:849963/1‑136 (MQ=255)                                                                                                                    tGTCTATCGCCGCGACCGGCTCAACCGCTTGCTGGAGCAGCAATACCTGACGCGCTTCAAGGAAAACCGCTACCTCCCGCCTTTCAGTCACACGTTCGCAGGCGTGGGGGTGAAGGAGGCCGTGAGCGTTGAGgcg     <  1:1087172/136‑1 (MQ=255)                                                                                                                       cTATAGCCGTGACCGGCTCAACCGCTTGCTGGAGCAGCAATACCTGACGCGCTTCAAGGAAAACCGCTACCTCCCGCCTTTCAGTCACACGTTCGCAGGCGTGGGGGTGAAGGAGGCCGTGAGCGTTGAGGCGCTg  <  1:833654/136‑1 (MQ=255)                                                                                                                                |                                                                                                                               CCTGCGTGTGGCCACATCACCATCCATTGGAGGTTGAGACCATGTCCAATTCACTGGAAGAACTGCTTGAGCTCATGAAAAAGGGCTCAGTAACGCTAGGCTGGGGGGCGGTTGCTGTCTATAGCCGTGACCGGCTCAACCGCTTGCTGGAGCAGCAATACCTGACGCGCTTCAAGGAAAACCGCTACCTCCCGCCTTTCAGTCACACGTTCGCAGGCGTGGGGGTGAAGGAGGCCGTGAGCGTTGAGGCGCTG  >  NC_002947/2123133‑2123386

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