Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,097,512 A→G 10.4% Q139R (CAG→CGG)  gpx → hydroperoxy fatty acid reductase Gpx1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,097,5120AG10.4% 130.7 / 5.0 48Q139R (CAG→CGG) gpxhydroperoxy fatty acid reductase Gpx1
Reads supporting (aligned to +/- strand):  ref base A (25/18);  new base G (3/2);  total (28/20)
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.81e-01

TTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGTCGTTCATGCC  >  NC_002947/2097366‑2097658
                                                                                                                                                  |                                                                                                                                                  
ttttCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGt                                                                                                                                                 <  2:1111113/150‑1 (MQ=255)
               cAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGcccc                                                                                                                                  >  2:47474/1‑150 (MQ=255)
               cAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGcccc                                                                                                                                  >  1:571785/1‑150 (MQ=255)
                  ggTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCccac                                                                                                                               <  2:571785/150‑1 (MQ=255)
                     cGAGGTCAACGGCCCGGGCGCTCACCGTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCccaccac                                                                                                                            <  2:108875/150‑1 (MQ=255)
                         gTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAg                                                                                                                        >  1:1257376/1‑150 (MQ=255)
                                                       gtcgAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGt                                                                                          <  2:1257376/150‑1 (MQ=255)
                                                        tcgAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGa                                                                                                          <  2:1408895/133‑1 (MQ=255)
                                                        tcgAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGa                                                                                                          >  1:1408895/1‑133 (MQ=255)
                                                            aCTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGccc                                                                                                              >  1:117741/1‑125 (MQ=255)
                                                            aCTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGccc                                                                                                              <  2:117741/125‑1 (MQ=255)
                                                            aCTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCt                                                                                     >  2:1279781/1‑150 (MQ=255)
                                                             cTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTc                                                                                    >  1:1369177/1‑150 (MQ=255)
                                                              tCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCa                                                                                   >  2:256919/1‑150 (MQ=255)
                                                                aaGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGc                                                                                                       <  1:1168785/128‑1 (MQ=255)
                                                                aaGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGc                                                                                                       >  2:1168785/1‑128 (MQ=255)
                                                                        cgcCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGCTCCAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAg                                                                         <  1:460338/150‑1 (MQ=255)
                                                                              gggTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCAcc                                                                   >  2:906466/1‑150 (MQ=255)
                                                                                 tATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCaa                                                                >  1:1203405/1‑150 (MQ=255)
                                                                                   ttttGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATg                                                              >  2:544327/1‑150 (MQ=255)
                                                                                        ggTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGa                                                         >  2:203186/1‑150 (MQ=255)
                                                                                         gTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGTCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGAc                                                        >  2:603730/1‑150 (MQ=255)
                                                                                         gTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGAc                                                        >  2:1056582/1‑150 (MQ=255)
                                                                                              cAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAggg                                                   <  1:544327/150‑1 (MQ=255)
                                                                                               agaagaTCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGc                                                  <  2:1369177/150‑1 (MQ=255)
                                                                                               agaagaTCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGc                                                  >  2:266159/1‑150 (MQ=255)
                                                                                                        aGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTTGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGctc                                         >  1:706010/1‑150 (MQ=255)
                                                                                                          tGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCtt                                       >  1:179002/1‑150 (MQ=255)
                                                                                                            gAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTcg                                     >  1:1425543/1‑150 (MQ=255)
                                                                                                             aaCTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTcgc                                    >  1:1220749/1‑150 (MQ=255)
                                                                                                                         tttCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCt                        >  1:282501/1‑150 (MQ=255)
                                                                                                                          ttCTGGTCGACCCGGCCAGTGGTCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGTGCACGCAATTGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCtt                       <  1:603730/150‑1 (MQ=255)
                                                                                                                           tCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTg                      <  2:179002/150‑1 (MQ=255)
                                                                                                                            cTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTgg                     >  2:901891/1‑150 (MQ=255)
                                                                                                                            cTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTgg                     <  1:256919/150‑1 (MQ=255)
                                                                                                                                ttgACCTGGCCACTGGCCGGGTCGAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGc                    >  1:380739/3‑147 (MQ=255)
                                                                                                                                ttgACCTGGCCACTGGCCGGGTCGAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGc                    <  2:380739/145‑1 (MQ=255)
                                                                                                                                ttgACCTGGCCACTGGCCGGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAg                 >  1:497321/3‑150 (MQ=255)
                                                                                                                                   acctggCCACTGGCCGGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCa                  <  1:953337/142‑1 (MQ=255)
                                                                                                                                   aCCTGGCCACTGGCCGGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCa                  >  2:953337/1‑146 (MQ=255)
                                                                                                                                   aCCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTa              >  1:215275/1‑150 (MQ=255)
                                                                                                                                     ccGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACAt                                                                                                 >  2:123731/1‑65 (MQ=255)
                                                                                                                                     ccGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACAt                                                                                                 <  1:123731/65‑1 (MQ=255)
                                                                                                                                      cGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGTc           <  1:901891/150‑1 (MQ=255)
                                                                                                                                       ggCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGt                                                            <  2:1252699/100‑1 (MQ=255)
                                                                                                                                       ggCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGt                                                            >  1:1252699/1‑100 (MQ=255)
                                                                                                                                       ggCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGTCg          <  1:1056582/150‑1 (MQ=255)
                                                                                                                                               gCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGTCGTTCATGcc  <  1:47474/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGTCGTTCATGCC  >  NC_002947/2097366‑2097658

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