Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,296,432 G→T 9.7% L19M (CTG→ATG)  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,4320GT 9.7% 157.8 / 4.5 51L19M (CTG→ATG) PP_2023glutathione S‑transferase family protein
Reads supporting (aligned to +/- strand):  ref base G (21/25);  new base T (2/3);  total (23/28)
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.96e-01

TGTCCAGGCTCTCTCCCAGTACCTCGGCACCGGTGTCCAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTGAGGCCATTACCCTGGCGCCGCACCTGGATCTGCACCGGGATGCGCTCGTGCATTTCCTGTACGTGGGAGATCACCG  >  NC_002947/2296288‑2296577
                                                                                                                                                |                                                                                                                                                 
tgTCCAGGCTCTCTCCCAGTACCTCGGCACCGGTGTCCAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGgcgc                                                                                                                                               <  2:49324/149‑1 (MQ=255)
     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                                                                                                                               <  6:308684/116‑4 (MQ=255)
                                                 cGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCATGCGTGCGCGCCTGGccctgc                                                                                                                               >  5:308684/1‑113 (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                                                                            <  1:162727/149‑1 (MQ=255)
                                                                   ggCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTGAg                                                                            <  5:168367/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)
                                                                                                                                                |                                                                                                                                                 
TGTCCAGGCTCTCTCCCAGTACCTCGGCACCGGTGTCCAGCACCGGTACCGTGCCCTTGGGCGACAAGGCCAGCAGCTCTGCCGGCTTGTTCTTCATCGCCACCTCGCAGATTTCCACCTCGCAGCCGGCATAGCGCAGGGCCAGGCGCGCACGCATGGCCCACGGGCAGCGGCGGAACGAGTAGAGGATCACCCGCACACCTCCACGTCACTGAGGCCATTACCCTGGCGCCGCACCTGGATCTGCACCGGGATGCGCTCGTGCATTTCCTGTACGTGGGAGATCACCG  >  NC_002947/2296288‑2296577

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