Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,951,168	T→C	10.8%	S42G (AGC→GGC)	PP_3482 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,951,168	0	T	C	10.8%	92.0 / 4.8	37	S42G (AGC→GGC)	PP_3482	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (14/19);  new base C (2/2);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.30e-01

AGAAGCGGCTGTAACGGGCCGACAGCAGCACCCGCAACTGGGCGCCACGTACCGCATGCTCGGCCAGCAGTGCAGCAAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCAACACGGCGCAGTAACGGTCAGATGAATGTGACACGGTTGATCTCCTCAAGGGTGGTACGGCCATCGCGGAC  >  NC_002947/3951033‑3951298                                                                                                                                        |                                                                                                                                   aGAAGCGGCTGTAACGGGCCGACAGCAGCACCCGCAACTGGGCGCCACGTACCGCATGCTCGGCCAGCAGTGCAGCAAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCt                                                                                                                                    >  1:269491/1‑136 (MQ=255)      cGGCTGTAACGGGCCGACAGCAGCACCCGCAACTGGGCGCCACGTACCGCATGCTCGGCCAGCAGTGCAGCAAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCa                                                                                                                               >  2:314962/1‑136 (MQ=255)                      aCAGCAGCACCCGCAACTGGGCGCCACGTACCGCATGCTCGGCCAGCAGTGCAGCAAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGtt                                                                                                               >  2:231499/1‑136 (MQ=255)                           agcaCCCGCAACTGGGCGCCACGTACCGCATGCTCGGCCAGCAGTGCAGCAAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTtgt                                                                                                          >  2:400838/1‑136 (MQ=255)                             caCCCGCAACTGGGCGCCACGTACCGCATGCTCGGCCAGCAGTGCAGCAAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGc                                                                                                        <  1:126921/136‑1 (MQ=255)                                   cAACTGGGCGCCACGTACCGCATGCTCGGCCAGCAGTGCAGCAAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCgc                                                                                                  >  2:314472/1‑135 (MQ=255)                                             ccACGTACCGCATGCTCGGCCAGCAGTGCAGCAAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCgc                                                                                        >  2:318071/1‑135 (MQ=255)                                             ccACGTACCGCATGCTCGGCCAGCAGTGCAGCAAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGa                                                                                        >  2:309142/1‑136 (MQ=255)                                             ccACGTACCGCATGCTCGGCCAGCAGTGCAGCAAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGa                                                                                        >  1:322787/1‑136 (MQ=255)                                             ccACGTACCGCATGCTCGGCCAGCAGTGCAGCAAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGa                                                                                        >  2:196822/1‑136 (MQ=255)                                                  tACCGCATGCTCGGCCAGCTGTGCAGCAAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGGCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACAccc                                                                                   <  1:50504/136‑1 (MQ=255)                                                  tACCGCATGCTCGGCCAGCAGTGCAGCAAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACAccc                                                                                   <  1:212612/136‑1 (MQ=255)                                                  tACCGCATGCTCGGCCAGCAGTGCAGCAAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACAccc                                                                                   <  1:261074/136‑1 (MQ=255)                                                       cATGCTCGGCCAGCAGTGCAGCAAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGgc                                                                              <  2:68208/136‑1 (MQ=255)                                                            tCGGCCAGCAGTGCAGCAAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCa                                                                         <  2:90003/136‑1 (MQ=255)                                                            tCGGCCAGCAGTGCAGCAAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCa                                                                         <  2:397236/136‑1 (MQ=255)                                                            tCGGCCAGCAGTGCAGCAAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCa                                                                         <  1:314962/136‑1 (MQ=255)                                                              ggCCAGCAGTGCAGCAAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGa                                                                                        >  1:400054/1‑119 (MQ=255)                                                              ggCCAGCAGTGCAGCAAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGa                                                                                        <  2:400054/119‑1 (MQ=255)                                                                 cagcagTGCAGCAAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCAACACg                                                                    <  2:343575/136‑1 (MQ=255)                                                                        gcagcaAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCAACACGGCGCAGt                                                             <  2:269491/136‑1 (MQ=255)                                                                               gTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCAACACGGCGCAGTAACGGTc                                                      <  1:40908/136‑1 (MQ=255)                                                                                    tCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCAACACGGCGCAGTAACGGTCAGATg                                                 >  1:94190/1‑136 (MQ=255)                                                                                       aCGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCAACACGGCGCAGTAACGGTCAGATGAAt                                              >  1:349810/1‑136 (MQ=255)                                                                                        cGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCAACACGGCGCAGTAACGGTCAGATGAAtg                                             <  1:208316/136‑1 (MQ=255)                                                                                           cAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTTGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCAACACGGCGCAGTAACGGTCAGATGAATGTGa                                          <  1:215415/136‑1 (MQ=255)                                                                                              tttCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCCGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCAACACGGCGCAGTAACGGTCAGATGAATGTGacac                                       >  2:103010/1‑136 (MQ=255)                                                                                              tttCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCCGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCAACACGGCGCAGTAACGGTCAGATGAATGTGacac                                       >  2:287915/1‑136 (MQ=255)                                                                                                  ccAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGGCCAACACGGCGCAGTAACGGTCAGATGAATGTGACACGGtt                                   <  1:400838/136‑1 (MQ=255)                                                                                                  ccAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCAACACGGCGCAGTAACGGTCAGATGAATGTGACACGGtt                                   <  1:231499/136‑1 (MQ=255)                                                                                                  ccAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCCGCCCAGCTATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCAACACGGCGCAGTAACGGTCAGATGAATGTGACACGGtt                                   <  1:287915/136‑1 (MQ=255)                                                                                                  ccAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCCGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCAACACGGCGCAGTAACGGTCAGATGAATGTGACACGGtt                                   <  1:103010/136‑1 (MQ=255)                                                                                                                       gCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCAACACGGCGCAGTAACGGTCAGATGAATGTGACACGGTTGATCTCCTCAAGGGTGGTACg              >  2:297527/1‑136 (MQ=255)                                                                                                                       gCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCAACACGGCGCAGTAACGGTCAGATGAATGTGACACGGTTGATCTCCTCAAGGGTGGTACg              >  2:42701/1‑136 (MQ=255)                                                                                                                           ccGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCAACACGGCGCAGTAACGGTCAGATGAATGTGACACGGTTGATCTCCTCAAGGGTGGTACGGCCa          <  2:155767/136‑1 (MQ=255)                                                                                                                                 gcgcgGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCAACACGGCGCAGTAACGGTCAGATGAATGTGACACGGTTGATCTCCTCAAGGGTGGTACGGCCATCGCgg    <  1:314472/136‑1 (MQ=255)                                                                                                                                 gcgcgGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCAACACGGCGCAGTAACGGTCAGATGAATGTGACACGGTTGATCTCCTCAAGGGTGGTACGGCCATCGCgg    <  1:196822/136‑1 (MQ=255)                                                                                                                                   gcgGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCAACACGGCGCAGTAACGGTCAGATGAATGTGACACGGTTGATCTCCTCAAGGGTGGTACGGCCATCGCGGAc  >  2:271346/1‑136 (MQ=255)                                                                                                                                        |                                                                                                                                   AGAAGCGGCTGTAACGGGCCGACAGCAGCACCCGCAACTGGGCGCCACGTACCGCATGCTCGGCCAGCAGTGCAGCAAGTGCATCGACGGCAGTTTCCCAGGCCGGTTGGGGGACGTCGCTGCCGAACGCGCGGCTGCCCAGCCATTGGTGGTGGTTCTTGTGCCAGCGACCCAGGCCGACACCCTCGGCGCCCAACACGGCGCAGTAACGGTCAGATGAATGTGACACGGTTGATCTCCTCAAGGGTGGTACGGCCATCGCGGAC  >  NC_002947/3951033‑3951298

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