Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,357,719	T→C	12.1%	L256P (CTC→CCC)	PP_2953 →	alcohol dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,357,719	0	T	C	12.1%	81.8 / 3.4	33	L256P (CTC→CCC)	PP_2953	alcohol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (14/15);  new base C (2/2);  total (16/17)
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.67e-01

TCAGCCATGTCGCCAGCCTTACCCATACTGATCAACACTTTGCGCAACTGGTTGAAGTGCTGCGACCACAAGGTCGCCTGGGCGTGATCGACGACCCGCAGACCCTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATGGTGCGCCAGCACGAGCTGCTCGAGCGCGTGGCCGGGCTGATCGATCAAGGCACCCTGCGTACCACCTTGGGCGAG  >  NC_002947/3357599‑3357864                                                                                                                         |                                                                                                                                                  tCAGCCATGTCGCCAGCCTTACCCATACTGATCAACACTTTGCGCAACTGGTTGAAGTGCTGCGACCACAAGGTCGCCTGGGCGTGATCGACGACCCGCAGACCCTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTgg                                                                                                                       <  7:119065/149‑1 (MQ=255) tCAGCCATGTCGCCAGCCTTACCCATACTGATCAACACTTTGCGCAACTGGTTGAAGTGCTGCGACCACAAGGTCGCCTGGGCGTGATCGACGACCCGCAGACCCTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTgg                                                                                                                       <  3:156247/149‑1 (MQ=255)   aGCCATGTCGCCAGCCTTACCCATACTGATCAACACTTTGCGCAACTGGTTGAAGTGCTGCGACCACAAGGTCGCCTGGGCGTGATCGACGACCCGCAGACCCTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGa                                                                                                                     <  8:43426/149‑1 (MQ=255)         gTCGCCAGCCTTACCCATACTGATCAACACTTTGCGCAACTGGTTGAAGTGCTGCGACCACAAGGTCGCCTGGGCGTGATCGACGACCCGCAGACCCTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGAt                                                                                                               >  5:176719/1‑149 (MQ=255)            gCCAGCCTTACCCATACTGATCAACACTTTGCGCAACTGGTTGAAGTGCTGCGACCACAAGGTCGCCTGGGCGTGATCGACGACCCGCAGACCCTGGACGTCATGCCGCTCAAGCGCAAGGCGCTGTCGCTGCACTGGGAGCTGATGtt                                                                                                            >  5:162157/1‑149 (MQ=255)             ccAGCCTTACCCATACTGATCAACACTTTGCGCAACTGGTTGAAGTGCTGCGACCACAAGGTCGCCTGGGCGTGATCGACGACCCGCAGACCCTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGtt                                                                                                            >  6:101019/1‑148 (MQ=255)                            cTGATCAACACTTTGCGCAACTGGTTGAAGTGCTGCGACCACAAGGTCGCCTGGGCGTGATCGACGACCCGCAGACTCTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTAc                                                                                            <  3:245516/149‑1 (MQ=255)                                       tttGCGCAACTGGTTGAAGTGCTGCGACCACAAGGTCGCCTGGGCGTGATCGACGACCCGCAGACCCTGGACGTCATGCCGCCCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACcgacg                                                                                  >  1:173992/1‑148 (MQ=255)                                       tttACGCAACTGGTTGAAGTGCTGCGACCACAAGGTCGCCTGGGCGTGATCGACGACCCGCAGACCCTGGACGTCATGCCGCCCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACcgacg                                                                                  >  1:7270/1‑148 (MQ=255)                                         tacgcAACTGGTTGAAGTGCTGCGACCACAAGGTCGCCTGGGCGTGATCGACGACCCGCAGACCCTGGACGTCATGCCGCCCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACcgacg                                                                                  <  2:7270/144‑1 (MQ=255)                                         tGCGCAACTGGTTGAAGTGCTGCGACCACAAGGTCGCCTGGGCGTGATCGACGACCCGCAGACCCTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACcgacga                                                                                 <  2:80086/147‑1 (MQ=255)                                         tGCGCAACTGGTTGAAGTGCTGCGACCACAAGGTCGCCTGGGCGTGATCGACGACCCGCAGACCCTGGACGTCATGCCGCCCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACcgacg                                                                                  <  2:173992/146‑1 (MQ=255)                                            gcAACTGGTTGAAGTGCTGCGACCACAAGGTCGCCTGGGCGTGATCGACGACCCGCAGACCCTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATgg                                                                            >  4:212334/1‑149 (MQ=255)                                              aaCTGGTTGAAGTGCTGCGACCACAAGGTCGCCTGGGCGTGATCGACGACCCGCAGACCCTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATGGTg                                                                          >  6:93923/1‑149 (MQ=255)                                                cTGGTTGAAGTGCTGCGACCACAAGGTCGCCTGGGCGTGATCGACGACCCGCAGACCCTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATGGTgcg                                                                        >  3:68216/1‑149 (MQ=255)                                                                        ggTCGCCTGGGCGTGATCGACGACCCGCAGACCCTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATGGTGCGCCAGCACGAGCTGCTCGAGCGCGt                                                >  3:138116/1‑149 (MQ=255)                                                                         gTCGCCTGGGCGTGATCGACGACCCGCAGACCCTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATGGTGCGCCAGCACGAGCTGCTCGa                                                      <  3:243492/142‑1 (MQ=255)                                                                         gTCGCCTGGGCGTGATCGACGACCCGCAGACCCTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATGGTGCGCCAGCACGAGCTGCTCGa                                                      >  4:243492/1‑142 (MQ=255)                                                                                   cGTGATCGACGAGCCGCAGACCCTGGACGGCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATGGTGCGCCAGCACGAGCTGCTCGAGCGCGTGGCCGGGCTga                                     <  8:117385/149‑1 (MQ=255)                                                                                    gTGATCGACGACCCGCAGACCCTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATGGTGCGCCAGCACGAGCTGCTCGAGCGCGTGGCCGGGCTgat                                    <  7:139881/149‑1 (MQ=255)                                                                                          gacgacCCGCAGACCCTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATGGTGCGCCAGCACGAGCTGCTCGAGCGCGTGGCCGGGCTgatcgatc                               >  4:8097/1‑148 (MQ=255)                                                                                            cgacCCGCAGACCCTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATGGTGCGCCAGCACGAGCTGCTCGAGCGCGTGGCCGGGCTGATCGATCAAg                            <  3:84928/149‑1 (MQ=255)                                                                                            cgacCCGCAGACCCTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATGGTGCGCCAGCACGAGCTGCTCGAGCGCGTGGCCGGGCTGATCGATCAAg                            <  5:93923/149‑1 (MQ=255)                                                                                             gacCCGCAGACCCTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATGGTGCGCCAGCACGAGCTGCTCGAGCGCGTGGCCGGGCTGATCGATCAAgg                           <  3:8097/149‑1 (MQ=255)                                                                                               cccGCAGACCCTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATGGTGCGCCAGCACGAGCTGCTCGAGCGCGTGGCCGGGCTGATCGATCAAGGCa                         >  7:21782/1‑149 (MQ=255)                                                                                                  gCAGACCCTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATGGTGCGCCAGCACGAGCTGCTCGAGCGCGTGGCCGGGCTGATCGATCAAGGCAccc                      >  4:145266/1‑149 (MQ=255)                                                                                                        ccTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATGGTGCGCCAGCACGAGCTGCTCGAGCGCGTGGCCGGGCTGATCGATCAAGGCACCCTGCGTa                >  7:57445/1‑149 (MQ=255)                                                                                                            gACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATGGTGCGCCAGCACGAGCTGCTCGAGCGCGTGGCCGGGCTGATCGATCAAGGCACCCTGCGTacc              <  3:212334/147‑1 (MQ=255)                                                                                                               gTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATGGTGCGCCAGCACGAGCTGCTCGAGCGCGTGGCCGGGCTGATCGATCAAGGCACCCTGCGTACCACCtt         <  2:206326/149‑1 (MQ=255)                                                                                                                 cATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATGGTGCGCCAGCACGAGCTGCTCGAGCGCGTGGCCGGGCTGATCGATCAAGGCACCCTGCGTACCACCTTgg       >  6:68141/1‑149 (MQ=255)                                                                                                                  aTGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATGGTGCGCCAGCACGAGCTGCTCGAGCGCGTGGCCGGGCTGATCGATCAAGGCACCCTGCGTACCACCTTggg      <  6:162157/149‑1 (MQ=255)                                                                                                                     ccGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATGGTGCGCCAGCACGAGCTGCTCGAGCGCGTGGCCGGGCTGATCGATCAAGGCACCCTGCGTACCACCTTGGGCg    <  3:145266/148‑1 (MQ=255)                                                                                                                       gCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATGGTGCGCCAGCACGAGCTGCTCGAGCGCGTGGCCGGGCTGATCGATCAAGGCACCCTGCGTACCACCTTGGGCGAg  >  7:138527/1‑148 (MQ=255)                                                                                                                         |                                                                                                                                                  TCAGCCATGTCGCCAGCCTTACCCATACTGATCAACACTTTGCGCAACTGGTTGAAGTGCTGCGACCACAAGGTCGCCTGGGCGTGATCGACGACCCGCAGACCCTGGACGTCATGCCGCTCAAGCGCAAGTCGCTGTCGCTGCACTGGGAGCTGATGTTCACCCGCTCGCTTTACGAAACCGACGACATGGTGCGCCAGCACGAGCTGCTCGAGCGCGTGGCCGGGCTGATCGATCAAGGCACCCTGCGTACCACCTTGGGCGAG  >  NC_002947/3357599‑3357864

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