Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,097,457	C→T	8.3%	P135L (CCG→CTG)	gpx →	hydroperoxy fatty acid reductase Gpx1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,097,457	0	C	T	8.3%	131.8 / 2.5	48	P135L (CCG→CTG)	gpx	hydroperoxy fatty acid reductase Gpx1
Reads supporting (aligned to +/- strand):  ref base C (23/21);  new base T (2/2);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.76e-01

ATTTCGGCGTGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAG  >  NC_002947/2097310‑2097600                                                                                                                                                     |                                                                                                                                                 aTTTCGGCGTGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑Cg                                                                                                                                                 >  1:118430/1‑149 (MQ=255) aTTGCGGCGTGAGTTTGCCGCTGTTTCGCAAGGTCGAGGGCAACGGCCCGGGCGCTCACCCTCTGTGCGTCGAACTCAAGCACCGCGCCCCGGGTATTTTGGGTAGCCAGAAGAGCAAGTGGAACTTCACCAAGTTTCTGGTCGAcc‑c                                                                                                                                                  <  5:92686/148‑1 (MQ=255)    tCGGCGTGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGc                                                                                                                                               <  1:103204/148‑1 (MQ=255)        cGTGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCTTCGAACTCAAGCAACGCGGCCGGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTg                                                                                                                                          <  7:18656/149‑1 (MQ=255)              ttttCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑gg                                                                                                                                   >  2:61081/1‑149 (MQ=255)                ttCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGt                                                                                                                                  <  3:111705/148‑1 (MQ=255)                   ctcTGTTTCGCAAGTTCGAGGTCAACGTCCCGGTCGCTCACCCTCTGTTCGTCTAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAATATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAg                                                                                                                              <  8:285568/147‑1 (MQ=255)                   cGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAg                                                                                                                              >  2:138338/1‑149 (MQ=255)                    gCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAgc                                                                                                                             >  2:88957/1‑149 (MQ=255)                     cTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑TGGCCACTGGCCG‑GGTcgacca                                                                                                                            <  2:265362/149‑6 (MQ=255)                     cTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAgcg                                                                                                                            <  1:61081/149‑1 (MQ=255)                            gCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTa                                                                                                                         >  2:156409/1‑145 (MQ=255)                            gCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTa                                                                                                                         <  1:156409/145‑1 (MQ=255)                                       gTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCaa                                                                                                          <  6:38230/149‑1 (MQ=255)                                                     cgcTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTgg                                                                                            <  4:18941/149‑1 (MQ=255)                                                       cTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCACGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGa                                                                                           <  1:88957/149‑1 (MQ=255)                                                             ctctGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGAc                                                                                    <  1:204172/149‑1 (MQ=255)                                                                tGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATc                                                                                 >  7:92807/1‑149 (MQ=255)                                                                     gtcgAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCg                                                                            <  1:184206/149‑1 (MQ=255)                                                                        gAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGttt                                                                         >  8:81886/1‑149 (MQ=255)                                                                          aCTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGc                                                                       >  3:181784/1‑149 (MQ=255)                                                                            tCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCt                                                                                                                          >  1:16110/1‑96 (MQ=255)                                                                            tCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCt                                                                                                                          <  2:16110/96‑1 (MQ=255)                                                                              aaGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAg                                                                   >  3:114775/1‑149 (MQ=255)                                                                                    cgcgCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCAGACATCGAGCGTTTGCTCAGCCGCt                                                              >  8:279830/1‑148 (MQ=255)                                                                                           cGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCa                                                      >  5:139596/1‑149 (MQ=255)                                                                                           cGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCAGACATCGAGCGTTTGCTCAGCCGCTGAACGGCa                                                      <  7:279830/149‑1 (MQ=255)                                                                                            gggTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCAc                                                     >  4:138372/1‑149 (MQ=255)                                                                                                     gggTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGt                                             <  2:118430/148‑1 (MQ=255)                                                                                                            cAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTg                                                                        <  2:251477/114‑1 (MQ=255)                                                                                                            cAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTg                                                                        >  1:251477/1‑114 (MQ=255)                                                                                                                     aaGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCATGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAg                             >  1:66690/1‑149 (MQ=255)                                                                                                                                cACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCgcggcg                 >  6:191805/1‑149 (MQ=255)                                                                                                                                cACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCgcggcg                 >  6:187304/1‑149 (MQ=255)                                                                                                                                   cAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGaa              >  5:72421/1‑149 (MQ=255)                                                                                                                                   cAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGaa              <  1:138338/149‑1 (MQ=255)                                                                                                                                     aGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATg            >  6:160656/1‑149 (MQ=255)                                                                                                                                      gTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATgg           <  7:81886/149‑1 (MQ=255)                                                                                                                                      gTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCACTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATgg           >  2:84451/1‑149 (MQ=255)                                                                                                                                       tttCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGc          <  3:138372/149‑1 (MQ=255)                                                                                                                                             gTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAg                             <  7:298888/124‑1 (MQ=255)                                                                                                                                             gTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTg       <  6:72421/146‑1 (MQ=255)                                                                                                                                             gTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGc     >  1:38647/1‑148 (MQ=255)                                                                                                                                             gTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGACGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAg                             >  8:298888/1‑124 (MQ=255)                                                                                                                                             cttgACC‑TGGCCACTGGCCG‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGt                                             <  7:298849/105‑1 (MQ=255)                                                                                                                                             cttgACC‑TGGCCACTGGCCG‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGCGCACCCAATGGt                                             >  8:298849/4‑108 (MQ=255)                                                                                                                                              ttgACC‑TGGCCACTGGCCG‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGc     >  5:123259/3‑147 (MQ=255)                                                                                                                                                gACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAg  <  6:90077/148‑1 (MQ=255)                                                                                                                                                gACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAg  >  5:90077/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                 ATTTCGGCGTGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAG  >  NC_002947/2097310‑2097600

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