Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,213,715	C→A	14.8%	P791Q (CCA→CAA)	PP_2812 →	transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,213,715	0	C	A	14.8%	80.6 / 3.5	27	P791Q (CCA→CAA)	PP_2812	transporter
Reads supporting (aligned to +/- strand):  ref base C (12/11);  new base A (2/2);  total (14/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 = 8.14e-01

TGGGCGTTCTCGGAGCTCAAGTTCCAGGCCGACATGGGCAAGCTGCTGGCCTTCATGTTCATCGTCAACATGGTCATGGCCATGACCGTGCTGCCGGCCTTCGCGGTCTGGCTGGAGCGGGTATTCCCACGTAAGCGCCCTGTGCGCATGATCGGCGCGCTGGTGCACTGAGGAGAGAAGCATGAAGTTCAGGGTTTCATTCACCGCGGCCTTGCTCGCAGCAGCATTGCCGCTGAT  >  NC_002947/3213588‑3213824                                                                                                                                |                                                                                                              tGGGCGTTCTCGGAGCTCAAGTTCCAGGCCGACATGGGCAAGCTGCTGGCCTTCATGTTCATCGTCAACATGGTCATGGCCATGACCGTGCTGCCGGCCTTCGCGGTCTGGCTGGAGCGGGTATTCCCACGTAAgc                                                                                                       <  2:323693/136‑1 (MQ=255)  gggCGTTCTCGGAGCTCAAGTTCCAGGCCGACATGGGCAAGCTGCTGGCCTTCATGTTCATCGTCAACATGGTCATGGCCATGACCGTGCTGCCGGCCTTCGCGGTCTGGCTGGAGCGGGTATTCCCACGTAAgcg                                                                                                      >  1:664318/1‑136 (MQ=255)  gggCGTTCTCGGAGCTCAAGTTACAGGCCGACATGGGCAAGCAGCTGGCCTTCATGTTCAACGTCAACAAGGTCATGGCCATGACCGTGCAGCCCGCCTAAGCGGTCTGGCTCGAGACGGTATTACCACGAAAgca                                                                                                      >  2:612069/1‑135 (MQ=255)   ggCGTTCTCGGCGCTCAAGTTCCAGGCCGACTTGGGCAAGCTGCTGGCCTTCATGTTCATCGTCAACATGGTCATGGCCATGACCGTGCTGCCGGCCTTCGCTGTCTGGCTGGAGCGGGTATTCCCACGTAAgcgc                                                                                                     <  2:353199/136‑1 (MQ=255)      gTTCTCGGAGCTCAAGTTCCAGGCCGACATGGGCAAGCTGCTGGCCTTCATGTTCATCGTCAACATGGTCATGGCCATGACCGTGCTGCCGGCCTTCGCGGTCTGGCTGGAGCGGGTATTCCCACGTAAGCGCCCt                                                                                                  <  1:573327/136‑1 (MQ=255)        tctcGGAGCTCAAGTTCCAGGCCGACATGGGCAAGCTGCTGGCCTTCATGTTCATCGTCAACATGGTCATGGCCATGACCGTGCTGCCGGCCTTCGCGGTCTGGCTGGAGCGGGTATTCCCACGTAAGCGCCCtgt                                                                                                <  1:662942/136‑1 (MQ=255)         ctcGGAGCTCAAGTTCCAGGCCGACATGGGCAAGCTGCTGGCCTTCATGTTCATCGTCAACATGGTCATGGCCATGACCGTGCTGCCGGCCTTCGCGGTCTGGCTGGAGCGGGTATTCCCACGTAAGCGCCCtgtg                                                                                               >  2:149538/1‑136 (MQ=255)           cGGAGCTCAAGTTCCAGGCCGACATGGGCAAGCTGCTGGCCTTCATGTTCATCGTCAACATGGTCATGGCCATGACCGTGCTGCCGGCCTTCGCGGTCTGGCTGGAGCGGGTATTCCCACGTAAGCGCCCTGTgcg                                                                                             >  2:409550/1‑136 (MQ=255)           cGGAGCTCAAGTTCCAGGCCGACATGGGCAAGCTGCTGGCCTTCATGTTCATCGTCAACATGGTCATGGCCATGACCGTGCTGCCGGCCTTCGCGGTCTGGCTGGAGCGGGTATTCCAACGTAAgc                                                                                                       >  1:323187/1‑126 (MQ=255)           cGGAGCTCAAGTTCCAGGCCGACATGGGCAAGCTGCTGGCCTTCATGTTCATCGTCAACATGGTCATGGCCATGACCGTGCTGCCGGCCTTCGCGGTCTGGCTGGAGCGGGTATTCCAACGTAAgc                                                                                                       <  2:519944/126‑1 (MQ=255)           cGGAGCTCAAGTTCCAGGCCGACATGGGCAAGCTGCTGGCCTTCATGTTCATCGTCAACATGGTCATGGCCATGACCGTGCTGCCGGCCTTCGCGGTCTGGCTGGAGCGGGTATTCCAACGTAAgc                                                                                                       <  2:323187/126‑1 (MQ=255)           cGGAGCTCAAGTTCCAGGCCGACATGGGCAAGCTGCTGGCCTTCATGTTCATCGTCAACATGGTCATGGCCATGACCGTGCTGCCGGCCTTCGCGGTCTGGCTGGAGCGGGTATTCCAACGTAAgc                                                                                                       >  1:519944/1‑126 (MQ=255)                       tCCAGGCCGACATGGGCAAGCTGCTGGCCTTCATGTTCATCGTCAACATGGTCATGGCCATGACCGTGCTGCCGGCCTTCGCGGTCTGGCTGGAGCGGGTATTCCCACGTAAGCGCCCTGTGCGCATGATCGgcgc                                                                                 >  2:528284/1‑136 (MQ=255)                                  aTGGGCAAGCTGCTGGCCTTCATGTTCATCGTCAACATGGTCATGGCCATGACCGTGCTGCCGGCCTTCGCGGTCTGGCTGGAGCGGGTATTCCCACGTAAGCGCCCTGTGCGCATGATCGGCGCGCTGGTGCACt                                                                      <  1:19595/136‑1 (MQ=255)                                    gggCAAGCTGCTGGCCTTCATGTTCATCGTCAACATGGTCATGGCCATGACCGTGCTGCCGGCCTTCGCGGTCTGGCTGGAGCGGGTATTCCCACGTAAGCGCCCTGTGCGCATGATCGGCGCGCTGGTGCACTga                                                                    <  1:435669/136‑1 (MQ=255)                                       cAAGCTGCTGGCCTTCATGTTCATCGTCAACATGGTCATGGCCATGACCGTGCTGCCGGCCTTCGCGGTCTGGCTGGAGCGGGTATTCCCACGTAAGCGCCCTGTGCGCATGATCGGCGCGCTGGTGCACTgagga                                                                 <  2:608205/136‑1 (MQ=255)                                              gtgGCCTTCATGTTCATCGTCAACATGGTCATGGCCATGACCGTGCTGCCGGCCTTCGCGGTCTGGCTGGAGCGGGTATTCCCACGTAAGCGCCCTGTGCGCATGATCGGCGCGCTGGTGCACTGAGGAGAGAAGc                                                          <  1:409550/135‑1 (MQ=255)                                                                              ggCCATGACCGTGCTGCCGGCCTTCGCGGTCTGGCTGGAGCGGTTATTCCCACGTAAGCGCCCTGTGCGCATGATCGGCGCGCTTGTGCACTGAGGAGAGAAGCATGAAGTTCAGGGTTTCATTCACCGCGGCCtt                          <  1:512831/136‑1 (MQ=255)                                                                                  aTGACCGTGCTGCCGGCCTTCGCGGTCTGGCTGGAGCGGGTATTCCCACGTAAGCGCCCTGTGCGCATGATCGGCGCGCTGGTGCACTGAGGAGAGAAGCATGAAGTTCAGGGTTTCATTCACCGCGGCCTTGCTc                      >  1:431573/1‑136 (MQ=255)                                                                                       cGTGCTGCCGGCCTTCGCGGTCTGGCTGGAGCGGGTATTCCCACGTAAGCGCCCTGTGCGCATGATCGGCGCGCTGGTGCACTGAGGAGAGAAGCATGAAGTTCAGGGTTTCATTCACCGCGGCCTTGCTCgcagc                 >  2:41677/1‑136 (MQ=255)                                                                                            tgcCGGCCTTCGCGGTCTGGCTGGAGCGGGTATTCCCACGTAAGCGCCCTGTGCGCATGATCGGCGCGCTGGTGCACTGAGGAGAGAAGCATGAAGTTCAGGGTTTCATTCACCGCGGCCTTGCTCGCAGCAGCAt            >  1:3392/1‑136 (MQ=255)                                                                                            tgcCGGCCTTCGCGGTCTGGCTGGAGCGGGTATTCCCACGTAAGCGCCCTGTGCGCATGATCGGCGCGCTGGTGCACTGAGGAGAGAAGCATGAAGTTCAGGGTTTCATTCACCGCGGCCTTGCTCGCAGCAGCAt            >  2:614171/1‑136 (MQ=255)                                                                                            tgcCGGCCTTCGCGGTCTGGCTGGAGCGGGTATTCCCACGTAAGCGCCCTGTGCGCATGATCGGCGCGCTGGTGCACTGAGGAGAGAAGCATGAAGTTCAGGGTTTCATTCACCGCGGCCTTGCTCGCAGCAGCAt            >  2:655484/1‑136 (MQ=255)                                                                                            tgcCGGCCTTCGCGGTCTGGCTGGAGCGGGTAATCCCACGTAAGCGCCCTGTGCGCATGATCGGCGCGCTGGTGCACAGAGGAGAGAAGCATGAAGTTCAGGCTTTCATTCACCGCGGCCTAGCTCGCAGCAGCAt            >  1:417966/1‑136 (MQ=255)                                                                                                    ttCGCGGTCTGGCTGGAGCGGGTATTCCCACGTAAGCGCCCTGTGCGCATGATCGGCGCGCTGGTGCACTGAGGAGAGAAGCATGAAGTTCAGGGTTTCATTCACCGCGGCCTTGCTCGCAGCAGCATTGCCGCTg    <  2:3392/136‑1 (MQ=255)                                                                                                    ttCGCGGTCTGGCTGGAGCGGGTATTCCCACGTAAGCGCCCTGTGCGCATGATCGGCGCGCTGGTGCACTGAGGAGAGAAGCATGAAGTTCAGGGTTTCATTCACCGCGGCCTTGCTCGCAGCAGCATTGCCGCTg    <  2:519286/136‑1 (MQ=255)                                                                                                      cgcgGTCTGGCTGGAGCGGGTATTCCCACGTAAGCGCCCTGTGCGCATGATCGGCGCGCTGGTGCACTGAGGAGAGAAGCATGAAGTTCAGGGTTTCATTCACCGCGGCCTTGCTCGCAGCAGCATTGCCGCTGAt  >  2:575193/1‑136 (MQ=255)                                                                                                                                |                                                                                                              TGGGCGTTCTCGGAGCTCAAGTTCCAGGCCGACATGGGCAAGCTGCTGGCCTTCATGTTCATCGTCAACATGGTCATGGCCATGACCGTGCTGCCGGCCTTCGCGGTCTGGCTGGAGCGGGTATTCCCACGTAAGCGCCCTGTGCGCATGATCGGCGCGCTGGTGCACTGAGGAGAGAAGCATGAAGTTCAGGGTTTCATTCACCGCGGCCTTGCTCGCAGCAGCATTGCCGCTGAT  >  NC_002947/3213588‑3213824

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