Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,097,463 G→C 8.5% 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,4630GC 8.5% 143.6 / 4.3 47S137T (AGT→ACT) gpxhydroperoxy fatty acid reductase Gpx1
Reads supporting (aligned to +/- strand):  ref base G (22/21);  new base C (2/2);  total (24/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.06e-01

CGTGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGTCGTTC  >  NC_002947/2097317‑2097610
                                                                                                                                                   |                                                                                                                                                    
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                                                                                                                                   <  1:156409/145‑1 (MQ=255)
                    gCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTa                                                                                                                                   >  2:156409/1‑145 (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                                                                                                                                    <  2:16110/96‑1 (MQ=255)
                                                                    tCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCt                                                                                                                                    >  1:16110/1‑96 (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            >  5:90077/1‑148 (MQ=255)
                                                                                                                                        gACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAg            <  6:90077/148‑1 (MQ=255)
                                                                                                                                                  aGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGTCGTTc  <  5:160656/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGTGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACC‑CGGCCAGTGGCCA‑GGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGTCGTTC  >  NC_002947/2097317‑2097610

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