Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,048,730:1 +C 8.7% coding (572/603 nt) PP_1821 → glutathione S‑transferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,048,7301.C 8.7% 175.1 / 11.6 46coding (572/603 nt)PP_1821glutathione S‑transferase family protein
Reads supporting (aligned to +/- strand):  ref base . (22/20);  new base C (2/2);  total (24/22)
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.51e-01

CAGATGACGCCTTACCTGGTGGGTGAGCAGTACTCGATTGCCGATGTGGCGCTGTATGCCTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCC  >  NC_002947/2048591‑2048877
                                                                                                                                            |                                                                                                                                                   
cAGATGACGCCTTACCTGGTGGGTGAGCAGTACTCGATTGCCGATGTGGCGCTGTATGCCTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCcgc                                                                                                                                             <  6:92467/148‑1 (MQ=255)
cAGATGACGCCTTACCTGGTGGGTGAGCAGTACTCGATTGCCGATGTGGCGCTGTATGCCTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCcgc                                                                                                                                             <  4:285566/148‑1 (MQ=255)
  gATGACGCCTTACCTGGTGGGTGAGCAGTACTCGATTGCCGATGTGGCGCTGTATGCCTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCgg                                                                                                                                           >  4:84141/1‑148 (MQ=255)
        gCCTTACCTGGTGGGTGAGCAGTACTCGATTGCCGATGTGGCGCTGTATGCCTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGc                                                                                                                                    >  7:32131/1‑149 (MQ=255)
             aCCTGGTGGGTGAGCAGTACTCGATTGCCGATGTGGCGCTGTATGCCTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGAtg                                                                                                                               >  1:351503/1‑149 (MQ=255)
               cTGGTGGGTGAGCAGTACTCGATTGCCGATGTGGCGCTGTATGCCTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGAtggt                                                                                                                             >  6:184686/1‑149 (MQ=255)
                    gggtgAGCAGTACTCGATTGCCGATGTGGCGCTGTATGCCTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCAc                                                                                                                                                 >  4:229697/1‑124 (MQ=255)
                    gggtgAGCAGGACTCGATTGCCGATGTGGCGCTGTATGCCTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCAc                                                                                                                                                 <  3:229697/124‑1 (MQ=255)
                      gtgAGCAGTACTCGATTGCCGATGTGGCGCTGTATGCCTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCCACCATCGGCACATGCCGCgggtggctgc                                                                                                                       <  4:292197/149‑5 (MQ=255)
                      gtgAGCAGTACTCGATTGCCGATGTGGCGCTGTATGCCTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCCACCATCGGCACATGCCGCgggtggctgc                                                                                                                       >  3:292197/1‑145 (MQ=255)
                      gtgAGCAGTACTCGATTGCCGATGTGGCGCTGTATGCCTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGa                                                                                                                      <  3:84141/149‑1 (MQ=255)
                         aGCAGTACTCGATTGCCGATGTGGCGCTGTATGCCTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGaaaa                                                                                                                   <  2:17473/149‑1 (MQ=255)
                           cAGTACTCGATTGCCGATGTGGCGCTGTATGCCTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAAcc                                                                                                                 >  1:6405/1‑149 (MQ=255)
                              tACTCGATTGCCGATGTGGCGCTGTATGCCTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGAtg                                                                                                                               >  4:253579/1‑132 (MQ=255)
                               aCTCGATTGCCGATGTGGCGCTGTATGCCTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGAc                                                                                                              <  2:351503/148‑1 (MQ=255)
                                    aTTGCCGATGTGGCGCTGTATGCCTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGttt                                                                                                        >  6:249708/1‑149 (MQ=255)
                                         cGATGTGGCGCTGTATGCCTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAg                                                                                                   >  7:154922/1‑149 (MQ=255)
                                                        tGCCTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGTGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATc                                                                                     <  5:184686/148‑1 (MQ=255)
                                                          ccTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCtt                                                                                  <  5:249708/149‑1 (MQ=255)
                                                           cTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGTCGATGGAGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCtt                                                                                  <  1:89863/148‑1 (MQ=255)
                                                           cTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTtg                                                                                 >  3:109627/1‑149 (MQ=255)
                                                           cTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTtg                                                                                 >  4:149472/1‑149 (MQ=255)
                                                           cTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGAGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCtt                                                                                  >  2:89863/1‑148 (MQ=255)
                                                                      tgGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGc                                                                      <  6:33586/149‑1 (MQ=255)
                                                                        gcgcAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGt                                                                    <  2:147638/149‑1 (MQ=255)
                                                                          gcagcaGGGCGGGTTCGCCCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGc                                                                  >  6:234498/1‑149 (MQ=255)
                                                                                    gggTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAacac                                                        <  3:145606/149‑1 (MQ=255)
                                                                                      gTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCCACCATCGGCACATGCCGCgggtggc                                                                                                                          >  8:16133/1‑81 (MQ=21)
                                                                                      gTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCCACCATCGGCACATGCCGCgggtggc                                                                                                                          <  7:16133/82‑2 (MQ=21)
                                                                                      gTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTgcg                                                                                          <  4:130936/113‑1 (MQ=255)
                                                                                      gTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTgcg                                                                                          >  3:130936/1‑113 (MQ=255)
                                                                                       ttCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACggg                                                     >  5:322019/1‑149 (MQ=255)
                                                                                         cGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGgc                                                    <  8:32131/148‑1 (MQ=255)
                                                                                          gcccTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGAGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCtt                                                                                  >  4:230784/3‑117 (MQ=255)
                                                                                           aCCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGca                                                                     >  5:271097/1‑128 (MQ=255)
                                                                                           aCCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGca                                                                     <  6:271097/129‑2 (MQ=255)
                                                                                             cTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACGTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGc                                               <  6:249469/149‑1 (MQ=255)
                                                                                             cTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGc                                               >  3:323199/1‑149 (MQ=255)
                                                                                                  aGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTg                                          >  2:80189/1‑149 (MQ=255)
                                                                                                         ccGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAacca                                   >  2:2890/1‑149 (MQ=255)
                                                                                                         ccGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAacca                                   >  6:232076/1‑149 (MQ=255)
                                                                                                         ccGGCGGGAAGGGCCGGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAGCACGGGCTGCGCCCGTgttcgcgg                                   <  4:288099/149‑8 (MQ=255)
                                                                                                            gcggTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGAGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACcg                                <  1:80189/149‑1 (MQ=255)
                                                                                                                   gggCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGcc                         <  3:149472/149‑1 (MQ=255)
                                                                                                                       ctggctggAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCt                     >  7:253998/1‑149 (MQ=255)
                                                                                                                           ctggAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTcgc                 <  8:253998/149‑1 (MQ=255)
                                                                                                                                        gcagcCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCt    >  5:311764/1‑149 (MQ=255)
                                                                                                                                         cagcCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTc   <  2:6405/149‑1 (MQ=255)
                                                                                                                                          agcCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTcc  <  8:154922/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
CAGATGACGCCTTACCTGGTGGGTGAGCAGTACTCGATTGCCGATGTGGCGCTGTATGCCTACACCCATGTGGCGCAGCAGGGCGGGTTCGACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCC  >  NC_002947/2048591‑2048877

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: