Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	311,569	G→T	12.5%	P270H (CCC→CAC)	PP_0256 ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	311,569	0	G	T	12.5%	79.1 / 4.2	32	P270H (CCC→CAC)	PP_0256	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (15/13);  new base T (2/2);  total (17/15)
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.39e-01

GGCCTGGGCTTCGCGTTCCCATTGCAGGACGAACTTGGCCATGCCGCGCAGCATCGCCATGTCGCCACCCAGGGCC‑GGGCGGAAGAAGGCGGTGTTGGTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAAT  >  NC_002947/311430‑311705                                                                                                                                             |                                                                                                                                         ggCCTGGGCTTCGCGTTCCCATTGCAGGACGAACTTGGCCATGCCGCGCAGCATCGCCATGTCGCCACCCAGGGCC‑GGGCGGAAGAAGGCGGTGTTGGTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGaa                                                                                                                                <  1:169823/150‑1 (MQ=255)  gCCTGGGCTTCGCGTTCCCATTGCAGGACGAACTTGGCCATGCCGCTCACCATATCCATGTCGCCCCCCAAGTCC‑GGGCGGGAAAAGGGGGGGTTGGGCGGGCGGTCACTGTTGGTGAGCATTTCCCGGTGGTTGTGGGGGCGCTGGaaa                                                                                                                               >  2:3712/1‑149 (MQ=255)       ggCTTCGCGTTCCCATTGCAGGACGAACTTGGCCATGCCGCGCAGCATCGCCATGTCGCCACCCAGGGCC‑GGGCGGAAGAAGGCGGTGTTGGTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTc                                                                                                                          >  2:446374/1‑150 (MQ=255)               gTTCCCATTGCAGGACGAACTTGGCCATGCCGCGCAGCATCGCCATGTCGCCACCCAGGGCC‑GGGCGGAAGAAGGCGGTGTTGGTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCAc                                                                                                                  >  1:1173/1‑150 (MQ=255)                          aGGACGAACTTGGCCATGCCGCGCAGCATCGCCATGTCGCCACCCAGGGCC‑GGGCGGAAGAAGGCGGTGTTGGTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGc                                                                                                       <  1:121686/150‑1 (MQ=255)                                  cTTGGCCATGCCGCGCAGCATCGCCATGTCGCCACCCAGGGCC‑GGGCGGAAGAAGGCGGTGTTGGTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGAt                                                                                               <  1:223226/150‑1 (MQ=255)                                  cTTGGCCATGCCGCGCAGCATCGCCATGTCGCCACCCAGGGCC‑GGGCGGAAGAAGGCGGTGTTGGTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGAt                                                                                               <  1:387090/150‑1 (MQ=255)                                     ggCCATGCCGCGCAGCATCGGCATGTCGCCACCCAGGGCC‑GGGCGGAAGAAGGCGGTGTTGGTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTTGAAGCGCTTCAGACCACGCTCCTTGAGCGGGTTGATGca                                                                                            <  2:557898/150‑1 (MQ=255)                                             cgcgCAGCATCGCCATGTCGCCACCCAGGGCC‑GGGCGGAAGAAGGCGGTGTTGGTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCt                                                                                    >  2:514284/1‑150 (MQ=255)                                                      tCGCCATGTCGCCACCCAGGGCCGGGGCGGAAGAAGGCGGTGTTGGTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCa                                                                            >  2:431170/1‑150 (MQ=255)                                                          cATGTCGCCACCCAGGGCC‑GGGCGGAAGAAGGCGGTGTTGGTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGt                                                                                                   <  1:514461/122‑1 (MQ=255)                                                          cATGTCGCCACCCAGGGCC‑GGGCGGAAGAAGGCGGTGTTGGTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGt                                                                                                   >  2:514461/1‑122 (MQ=255)                                                             gTCGCCACCCAGGGCC‑GGGCGGAAGAAGGCGGTGTTGGTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGTGTTCTTGGGGTGCTGGAAGCGCGCCAGACAAAGCTCCCTGAGCGTGTTGATGCACACCACCTGGGCGCCACGCTTTAc                                                                    <  1:156274/150‑1 (MQ=255)                                                                        gggCC‑GGGCGGAAGAAGGCGGTGTTGGTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCATTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATAACGCa                                                         >  1:356484/1‑150 (MQ=255)                                                                                gCGGAAGAAGGCGGTGTTGGTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTc                                                  <  1:514284/150‑1 (MQ=255)                                                                                  ggAAGAAGGCGGTGTTGGTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGa                                                >  2:356494/1‑150 (MQ=255)                                                                                  ggAAGAAGGCGGTGTTGGTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGa                                                >  1:196729/1‑150 (MQ=255)                                                                                      gaagGCGGTGTTGGTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCAt                                            >  1:192054/1‑150 (MQ=255)                                                                                          gCGGTGTTGGTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCGCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATgcgc                                        >  1:358740/1‑150 (MQ=255)                                                                                            ggTGTTGGTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCgg                                      <  2:557372/150‑1 (MQ=255)                                                                                             gtgtTGGTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCggg                                     <  1:446374/150‑1 (MQ=255)                                                                                                  ggTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGAtt                                <  2:1173/150‑1 (MQ=255)                                                                                                                tGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCt                  >  1:108030/1‑150 (MQ=255)                                                                                                                 gTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCCCACCACCGGGGCGCCAAGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTg                 <  2:250882/150‑1 (MQ=255)                                                                                                                   tGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGc               <  2:196729/150‑1 (MQ=255)                                                                                                                      tCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGcccc            >  2:182284/1‑150 (MQ=255)                                                                                                                         gCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGa         >  2:200051/1‑150 (MQ=255)                                                                                                                                cAGTGGGTTCTGGTGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCa                                                         <  2:11331/95‑1 (MQ=255)                                                                                                                                cAGTGGGTTCTGGTGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCa                                                         >  1:11331/1‑95 (MQ=255)                                                                                                                                cAGTGGGTTCTGGTGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCa                                                         <  2:11305/95‑1 (MQ=255)                                                                                                                                cAGTGGGTTCTGGTGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCa                                                         >  1:11305/1‑95 (MQ=255)                                                                                                                                cAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAAt  >  1:355575/1‑150 (MQ=255)                                                                                                                                             |                                                                                                                                         GGCCTGGGCTTCGCGTTCCCATTGCAGGACGAACTTGGCCATGCCGCGCAGCATCGCCATGTCGCCACCCAGGGCC‑GGGCGGAAGAAGGCGGTGTTGGTCGGGCGGTCACTGTTGGTCAGCATTTCCAGTGGGTTCTGGGGGTGCTGGAAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAAT  >  NC_002947/311430‑311705

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