Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,346,556 A→C 9.8% V179G (GTC→GGC)  gor ← glutathione reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,346,5560AC 9.8% 133.6 / 7.0 51V179G (GTC→GGC) gorglutathione reductase
Reads supporting (aligned to +/- strand):  ref base A (19/27);  new base C (2/3);  total (21/30)
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.97e-01

GGCGTTGAACTGCAGGTCGAGGCCGCGTTTCTCCAGCTCTTGCTTCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCTTGAGGTAGAAGGCTTCGTTGGAGGTGATGGCCAGTTCCTTGCCGGGAATGTCCGGTACCTGCGGCCAGCCACCGGTAGCGATGAGGA  >  NC_002947/4346408‑4346691
                                                                                                                                                    |                                                                                                                                       
ggCGTTGAACTGCAGGTCGAGGCCGCGTTTCTCCAGCTCTTGCTTCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGa                                                                                                                                         >  6:116898/1‑149 (MQ=255)
ggCGTTGAACTGCAGGTCGAGGCCGCGTTTCTCCAGCTCTTGCTTCAGGTGGGTGCGCACCCAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGa                                                                                                                                         >  2:101533/1‑149 (MQ=255)
   gTTGAACTGCAGGTCGAGGCCGCGTTTCTCCAGCTCTTGCTTCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACgg                                                                                                                                      >  8:106569/1‑149 (MQ=255)
       aaCTGCAGGTCGAGGCCGCGTTTCTCCAGCTCTTGCTTCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGTCGAATTCGCCGGcat                                                                                                                                   <  8:175442/148‑3 (MQ=255)
             aGGTCGAGGCCGCGTTTCTCCAGCTCTTGCTTCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAAcc                                                                                                                            >  5:245096/1‑149 (MQ=255)
              ggTCGAGGCCGCGTTTCTCCAGCTCTTGCTTCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGat                                                                                                                                  <  7:166973/142‑1 (MQ=255)
              ggTCGAGGCCGCGTTTCTCCAGCTCTTGCTTCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGat                                                                                                                                  >  8:166973/1‑142 (MQ=255)
                  gctgCCGCGTTTATCCAGCTCTTGGTTCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTTCAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCcccac                                                                                                                       <  6:44096/146‑1 (MQ=255)
                  gAGGCCGCGTTTCTCCAGCTCTTGCTTCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTTTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCccca                                                                                                                        <  5:179161/148‑1 (MQ=255)
                  gAGGCCGCGTTTCTCCAGCTCTTGCTTCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGAACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCccca                                                                                                                        <  5:174913/148‑1 (MQ=255)
                    ggCCGCGTTTCTCCAGCTCTTGCTTCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCcccaccc                                                                                                                     <  8:212048/149‑1 (MQ=255)
                       cgcgTTTCTCCAGCTCTTGCTTCAGGTGGGTGCGCACCGAGCCGTCTAAGCCGCGCATGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCcacccacc                                                                                                                  <  1:211899/149‑1 (MQ=255)
                       cgcgTTTCTCCAGCTCTTGCTTCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCcacccacc                                                                                                                  >  4:219319/1‑149 (MQ=255)
                          gTTTCTCCAGCTCTTGCTTCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACccaccacc                                                                                                               <  1:75069/149‑1 (MQ=255)
                          gTTTCTCCAGCTCTTGCTTCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACccaccac                                                                                                                >  4:278240/1‑148 (MQ=255)
                           tttCTCCAGCTCTTGCTTCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACccaccacca                                                                                                              <  3:104450/149‑1 (MQ=255)
                                 cAGCTCTTGCTTCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCAcgc                                                                                                        >  4:148984/1‑149 (MQ=255)
                                   gCTCTTGCTTCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCAcgcgg                                                                                                      >  8:30469/1‑149 (MQ=255)
                                     tctTGCTTCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGTCGAATTCGCCGGcat                                                                                                                                   >  6:134548/1‑116 (MQ=255)
                                     tctTGCTTCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGTCGAATTCGCCGGcat                                                                                                                                   <  5:134548/118‑3 (MQ=255)
                                       ttGCTTCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCAcgcggcgcg                                                                                                  >  7:38794/1‑149 (MQ=255)
                                         gCTTCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCAcgcggcgcgg                                                                                                 >  6:154590/1‑148 (MQ=255)
                                           ttCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGc                                                                                                            <  4:114764/135‑1 (MQ=255)
                                           ttCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGc                                                                                                            >  3:114764/1‑135 (MQ=255)
                                            tCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCaaa                                                                                             >  8:129657/1‑149 (MQ=255)
                                               ggtgggtgCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGACGAATTCGa                                                                                                                                         >  2:229340/1‑102 (MQ=255)
                                               ggtgggtgCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGACGAATTCGa                                                                                                                                         <  1:229340/102‑1 (MQ=255)
                                                 tgggtgCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCATCACGCGGCGCGGCAAATCCtt                                                                                        <  6:245096/149‑1 (MQ=255)
                                                 tgggtgCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCtt                                                                                        <  8:63807/149‑1 (MQ=255)
                                                 tgggtgCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCtt                                                                                        <  1:101533/149‑1 (MQ=255)
                                                 tgggtgCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCtt                                                                                        <  3:219319/149‑1 (MQ=255)
                                                 tgggtgCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCtt                                                                                        <  3:278240/149‑1 (MQ=255)
                                                   ggtgCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCTTGa                                                                                      <  7:30469/149‑1 (MQ=255)
                                                   ggtgCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCTTGa                                                                                      <  7:106569/149‑1 (MQ=255)
                                                       cgcACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCTTGAGGt                                                                                   <  1:16591/148‑1 (MQ=255)
                                                               gCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCTTGAGGTAGaa                                                                               <  5:1998/144‑1 (MQ=255)
                                                               gCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCTTGAGGTAGaa                                                                               >  6:1998/1‑144 (MQ=255)
                                                                  gTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGTCGAATTCGCCGGcat                                                                                                                                   <  2:123818/89‑3 (MQ=255)
                                                                       aGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCTTGAGGTAGAAGGCTTCGTTGGAg                                                                  <  3:148984/149‑1 (MQ=255)
                                                                        gCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCTTGAGGTAGAAGGCTTCGTTGGAgg                                                                 <  7:129657/149‑1 (MQ=255)
                                                                            cgcAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCTTGAGGTAGAAGGCTTCGTTGGAGGTGAt                                                             <  5:116898/149‑1 (MQ=255)
                                                                                            cACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCCCCCGCACGCGGCGCGGCAAATCCTTGAGGTAGAAGGGTTCGTGGGAGGTGATGGCCAGTTCCTTGCCg                                             >  7:82186/1‑149 (MQ=255)
                                                                                            cACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCTTGAGGTAGAAGGCTTCGTTGGAGGTGATGGCCAGTTCCTTGCCg                                             >  7:178348/1‑149 (MQ=255)
                                                                                                 tACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCTTGAGGTAGa                                                                                <  7:163313/109‑1 (MQ=255)
                                                                                                 tACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCTTGAGGTAGa                                                                                >  8:163313/1‑109 (MQ=255)
                                                                                                                tCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCTTGAgg                                                                                    <  1:241710/90‑1 (MQ=255)
                                                                                                                tCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCTTGAgg                                                                                    >  2:241710/1‑90 (MQ=255)
                                                                                                                  ggCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCTTGAGGTAGAAGGCTTCGTTGGAGGTGATGGCCAGTTCCTTGCCGGGAATGTCCGGTACCTGCGgcc                       <  8:38794/149‑1 (MQ=255)
                                                                                                                   gCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCTTGAGGTAGAAGGCTTCGTTGGAGGTGCTGGCCAGTTCCTTGCCGGGAATGTCCGGTACCTGCGgcca                      <  8:178348/149‑1 (MQ=255)
                                                                                                                               tGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCTTGAGGTAGAAGGCTTCGTTGGAGGTGATGGCCAGTTCCTTGCCGGGAATGTCCGGTACCTGCGGCCAGCCACCGGTAGc          >  1:96699/1‑149 (MQ=255)
                                                                                                                                 gAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCTTGAGGTAGAAGGCTTCGTTGGAGGTGATGGCCAGTTcc                                                   <  2:263459/106‑1 (MQ=255)
                                                                                                                                 gAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCTTGAGGTAGAAGGCTTCGTTGGAGGTGATGGCCAGTTcc                                                   >  1:263459/1‑106 (MQ=255)
                                                                                                                                     tcgCCGTCGAATTCGCCGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCTTGAGGTAGAAGGCTTCGTTGGAGGTGATGGCCAGTTCCTTGCCGGGAATGTCCGGTACCTGCGGCCAGCCACCGGTAGCGATGAg    >  3:143233/3‑149 (MQ=255)
                                                                                                                                        ccGGCGAATTCGACGGCGATATAACCCCCACCCCCCACCAGCACGCGGCGCGGCAAATCCTTGAGGTAGAAGGCTTCGTTGGAGGTGATGGCCAGTTCCTTGCCGGGAATGTCCGGTACCTGCGGCCAGCCTCCGGTAGCGATGAGGa  >  6:9567/1‑148 (MQ=255)
                                                                                                                                        ccGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCTTGAGGTAGAAGGTTTCGTTGGAGGTGATGGCCAGTTCCTTGCCGGGAATGTCCGGTACCTGCGGCCAGCCACCGGTAGCGATGAGGa  <  5:154590/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                       
GGCGTTGAACTGCAGGTCGAGGCCGCGTTTCTCCAGCTCTTGCTTCAGGTGGGTGCGCACCGAGCCGTCGAAGCCGCGCAGGAACAGGTCGCCACGGTACAGCAGTGTGGTGTCGGCCCCCAGGCCCTGGAAGATGCCGGCGAATTCGACGGCGATATAACCCCCACCCACCACCAGCACGCGGCGCGGCAAATCCTTGAGGTAGAAGGCTTCGTTGGAGGTGATGGCCAGTTCCTTGCCGGGAATGTCCGGTACCTGCGGCCAGCCACCGGTAGCGATGAGGA  >  NC_002947/4346408‑4346691

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