Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,097,506 G→C 9.6% S137T (AGT→ACT)  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,5060GC 9.6% 158.5 / 7.1 52S137T (AGT→ACT) gpxhydroperoxy fatty acid reductase Gpx1
Reads supporting (aligned to +/- strand):  ref base G (26/21);  new base C (3/2);  total (29/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 = 9.37e-01

GCGTGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGTCG  >  NC_002947/2097359‑2097650
                                                                                                                                                   |                                                                                                                                                
gCGTGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTg                                                                                                                                                <  2:720873/150‑1 (MQ=255)
  gTGAGTTTTCCGCTGTGTCGCAAGGTCGAGCTGAACGGGACGTCAGATGTGTGTCTTTTCGTTGAACTCAAGCATAATTCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGc                                                                                                                                              <  1:1027970/150‑1 (MQ=255)
   tGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGcc                                                                                                                                             <  1:1281000/150‑1 (MQ=255)
       ttttCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGt                                                                                                                                                 <  1:595933/142‑1 (MQ=255)
       ttttCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGt                                                                                                                                                 >  2:595933/1‑142 (MQ=255)
       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                                                                                                  >  1:1408895/1‑133 (MQ=255)
                                                               tcgAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGa                                                                                                  <  2:1408895/133‑1 (MQ=255)
                                                                   aCTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGccc                                                                                                      <  2:117741/125‑1 (MQ=255)
                                                                   aCTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGccc                                                                                                      >  1:117741/1‑125 (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:266159/1‑150 (MQ=255)
                                                                                                      agaagaTCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGc                                          <  2:1369177/150‑1 (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             <  1:256919/150‑1 (MQ=255)
                                                                                                                                   cTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTgg             >  2:901891/1‑150 (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)
                                                                                                                                                   |                                                                                                                                                
GCGTGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGTCG  >  NC_002947/2097359‑2097650

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