Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,296,443 A→C 9.8% M15R (ATG→AGG)  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,4430AC 9.8% 148.1 / 7.0 51M15R (ATG→AGG) PP_2023glutathione S‑transferase family protein
Reads supporting (aligned to +/- strand):  ref base A (22/24);  new base C (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 = 8.99e-01

TCCCAGTACCTCGGCACCGGTGTCCAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTGAGGCCATTACCCTGGCGCCGCACCTGGATCTGCACCGGGATGCGCTCGTGCATTTCCTGTACGTGGGAGATCACCGCCACTTTGC  >  NC_002947/2296301‑2296586
                                                                                                                                              |                                                                                                                                               
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                                                                                                                                        >  8:99620/1‑125 (MQ=255)
                        cAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCATGCGTGCGCGCCTGGccctgc                                                                                                                                        <  7:99620/128‑4 (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)
                                                                                                                                              |                                                                                                                                               
TCCCAGTACCTCGGCACCGGTGTCCAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTGAGGCCATTACCCTGGCGCCGCACCTGGATCTGCACCGGGATGCGCTCGTGCATTTCCTGTACGTGGGAGATCACCGCCACTTTGC  >  NC_002947/2296301‑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: ‑