Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,296,439 A→G 9.8% R16R (CGT→CGC)  PP_2023 ← glutathione S‑transferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,296,4390AG 9.8% 137.4 / 5.5 51R16R (CGT→CGC) PP_2023glutathione S‑transferase family protein
Reads supporting (aligned to +/- strand):  ref base A (22/24);  new base G (2/3);  total (24/27)
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.34e-01

AGGCTCTCTCCCAGTACCTCGGCACCGGTGTCCAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTGAGGCCATTACCCTGGCGCCGCACCTGGATCTGCACCGGGATGCGCTCGTGCATTTCCTGTACGTGGGAGATCACCGCCACTTTGC  >  NC_002947/2296293‑2296586
                                                                                                                                                  |                                                                                                                                                   
aGGCTCTCTCCCAGTACCTCGGCACCGGTGTCCAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGcgcacg                                                                                                                                                   >  1:218826/1‑149 (MQ=255)
        tcCCAGTACCTCGGCACCGGTGTCCAGCACCGGTACCGTGCCCTTGGGCGACGAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGc                                                                                                                                             <  3:223236/147‑1 (MQ=255)
        tcCCAGTACCTCGGCACCGGTGTCCAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGcc                                                                                                                                            >  3:258735/1‑148 (MQ=255)
          ccAGTACCTCGGCACCGGTGTCCAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCAc                                                                                                                                         >  8:93098/1‑149 (MQ=255)
          ccAGTACCTCGGCACCGGTGTCCAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCAc                                                                                                                                         >  7:107349/1‑149 (MQ=255)
           cAGTACCTCGGCACCGGTGTCCAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCATGCGTGCGCGCCTGGccctgc                                                                                                                                        <  7:169505/149‑4 (MQ=255)
                 cTCGGCACCGGTGTCCAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAgc                                                                                                                                  >  8:111086/1‑149 (MQ=255)
                     gCACCGGTGTCCAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAgcggcg                                                                                                                              >  7:238822/1‑149 (MQ=255)
                       aCCGGTGTCCAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGcggcgg                                                                                                                             >  1:222096/1‑148 (MQ=255)
                        ccGGTGTCCAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGaa                                                                                                                           >  6:168367/1‑149 (MQ=255)
                         cGGTGTCCAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGaa                                                                                                                           <  1:150681/148‑1 (MQ=255)
                         cGGTGTCCAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGaa                                                                                                                           >  2:150681/1‑148 (MQ=255)
                         cGGTGTCCAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAAc                                                                                                                          <  7:111086/149‑1 (MQ=255)
                              tCCAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGt                                                                                                                      <  4:4038/148‑1 (MQ=255)
                              tCCAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTa                                                                                                                     <  2:218826/149‑1 (MQ=255)
                                cAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCATGCGTGCGCGCCTGGccctgc                                                                                                                                        <  7:99620/128‑4 (MQ=255)
                                cAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCATGCGTGCGCGCCTGGccctgc                                                                                                                                        >  8:99620/1‑125 (MQ=255)
                                cAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTaga                                                                                                                   <  5:88557/149‑1 (MQ=255)
                                cAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTag                                                                                                                    <  4:258735/148‑1 (MQ=255)
                                cAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTa                                                                                                                     >  4:91620/1‑147 (MQ=255)
                                      cGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGcggcgg                                                                                                                             <  4:127105/133‑1 (MQ=255)
                                      cGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGcggcgg                                                                                                                             >  3:127105/1‑133 (MQ=255)
                                           ccGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGc                                                                                                        <  2:277153/149‑1 (MQ=255)
                                            cGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCATGCGTGCGCGCCTGGccctgc                                                                                                                                        >  5:308684/1‑113 (MQ=255)
                                            cGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCATGCGTGCGCGCCTGGccctgc                                                                                                                                        <  6:308684/116‑4 (MQ=255)
                                             gTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCCGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTATAGGATCACCCGca                                                                                                       >  3:197932/1‑148 (MQ=255)
                                             gTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGca                                                                                                       >  8:1247/1‑148 (MQ=255)
                                                     gggCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTcc                                                                                                <  8:238822/147‑1 (MQ=255)
                                                     gggCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCa                                                                                               <  3:91620/148‑1 (MQ=255)
                                                     gggCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCAc                                                                                              >  6:234827/1‑149 (MQ=255)
                                                         gACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGt                                                                                                                      <  4:236313/121‑1 (MQ=255)
                                                         gACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGt                                                                                                                      >  3:236313/1‑121 (MQ=255)
                                                         gACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCa                                                                                          <  1:126597/149‑1 (MQ=255)
                                                           cAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACt                                                                                        <  2:237093/149‑1 (MQ=255)
                                                              ggCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTGAg                                                                                     <  5:168367/149‑1 (MQ=255)
                                                              ggCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTGAg                                                                                     <  1:162727/149‑1 (MQ=255)
                                                                     agcTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGcac                                                                                                      <  2:111459/125‑1 (MQ=255)
                                                                     agcTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTGAGGCCATTa                                                                              <  2:222096/149‑1 (MQ=255)
                                                                                        tcATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTGAGGCCATTACCCTGGCGCCGCACCTGGa                                                           >  6:50614/1‑149 (MQ=255)
                                                                                          aTCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGTCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTGAGGCCATTACCCTGGCGCCGCACCTGGATc                                                         <  5:52050/149‑1 (MQ=255)
                                                                                          aTCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTGAGGCCATTACCCTGGCGCCGCACCTGGATc                                                         >  6:157568/1‑149 (MQ=255)
                                                                                                    cGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTGAGGCCATTACCCTGGCGCCGCACCTGGATCTGCACCGGGa                                               >  6:331549/1‑149 (MQ=255)
                                                                                                        gATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTGAGGCCATTACCCTGGCGCCGCACCTGGATCTGCACCGGGATgcg                                           <  7:93098/149‑1 (MQ=255)
                                                                                                              cACCTCTCAGCCGGCATAGTGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTGAGGCCATTACCCTGGCGCCGCACCTGGATCTGCACCGGGATGCGCTCGTg                                     >  6:290573/1‑149 (MQ=255)
                                                                                                                ccTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTcc                                                                                                <  7:63415/88‑1 (MQ=255)
                                                                                                                ccTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTcc                                                                                                >  8:63415/1‑88 (MQ=255)
                                                                                                                     cAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTGAGGCCATTACCCTGGCGCCGCACCTGGATCTGCACCGGGATGCGCTCGTGCATTTcc                              <  7:92035/149‑1 (MQ=255)
                                                                                                                                gtgcAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTGAGGCCATTACCCTGGCGCCGCACCTGGATCTGCACCGGGATGCGCTCGTGCATTTCCTGTACGTggg                    <  5:290573/146‑1 (MQ=255)
                                                                                                                                  gcAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTGAGGCCATTACCCTGGCGCCGCACCTGGATCTGCACCGGGATGCGCTCGTGCATTTCCTGTACGTGGgaga                 <  5:331549/149‑1 (MQ=255)
                                                                                                                                        ccAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGCGGATCACCCGCACACCTCCACGTCACTGAGGCCATTACCCTGGCGCCGCACCTGGATCTGCACCGGGATGCGCTCGTGCATTTCCTGTACGTGGGAGATCACCg           >  3:265036/1‑149 (MQ=255)
                                                                                                                                            gcgcgcACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTGAGGCCATTGCCCTGGCGCCGCACCTGGATCTGCACCGGGATGCGCTCGTGCATTTCCTGTACGTGGGAGATCACCGCCAc       >  5:198622/1‑149 (MQ=255)
                                                                                                                                              gcgcACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTGAGGCCATTACCCTGGCGCCGCACCTGGATCTGCACCGGGATGCGCTCGTGCATTTCCTGTACGTGGGAGATCACCGCCACtt     >  7:89544/1‑149 (MQ=255)
                                                                                                                                               cgcacgcaTGGCCCACGGGCAGCGGCCGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTTATTATTATTCATTGGGGGAGCACCTTGATCTTCACCGGGATGCGCTCGTGCATTTCCTGTACGTGGGAGATCACCGCCACttt    <  8:107349/149‑1 (MQ=255)
                                                                                                                                                  acgcaTGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTGAGGCCATTACCCTGGCGCCGCACCTGGATCTGCACCGGGATGCGCTCGTGCATTTCCTGTACGTGGGAGATCACCGCCACTTTGc  >  6:92563/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AGGCTCTCTCCCAGTACCTCGGCACCGGTGTCCAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTGAGGCCATTACCCTGGCGCCGCACCTGGATCTGCACCGGGATGCGCTCGTGCATTTCCTGTACGTGGGAGATCACCGCCACTTTGC  >  NC_002947/2296293‑2296586

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