Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,337,036 C→A 7.4% L137I (CTC→ATC)  gstA → glutathione S‑transferase GstA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,337,0360CA 7.4% 170.3 / 3.5 54L137I (CTC→ATC) gstAglutathione S‑transferase GstA
Reads supporting (aligned to +/- strand):  ref base C (24/26);  new base A (2/2);  total (26/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 = 7.87e-01

GAAGAACCTGCCGCCGCCGCCCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGAGCCCTATGCCAATGTGCGCAAC  >  NC_002947/1336889‑1337179
                                                                                                                                                   |                                                                                                                                               
gaagaaCCTGCCGCCGCCGCCCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCt                                                                                                                                                <  1:180405/149‑1 (MQ=255)
gaagaaCCTGCCGCCGCCGCCCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCt                                                                                                                                                >  6:41978/1‑149 (MQ=255)
gaagaaCCTGCCGCCGCCGCCCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCt                                                                                                                                                >  2:180405/1‑149 (MQ=255)
    aaCCTGCCGCCGCCGCCCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAg                                                                                                                                            <  7:68345/149‑1 (MQ=255)
     aCCTGCCGCCGCCGCCCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAgg                                                                                                                                           >  3:130760/1‑149 (MQ=255)
          ccgccgccgccCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGa                                                                                                                                        <  2:321578/147‑1 (MQ=255)
          ccgccgccgccCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGAGCAAGGTgtgg                                                                                                                                      <  5:41292/149‑4 (MQ=255)
            gccgccgccCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGa                                                                                                                                    <  2:292745/149‑1 (MQ=255)
                ccgccCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGcc                                                                                                                                <  3:1036/149‑1 (MQ=255)
                       ggTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAgg                                                                                                                                           >  1:291161/1‑131 (MQ=255)
                       ggTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAgg                                                                                                                                           <  2:291161/131‑1 (MQ=255)
                        gTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGc                                                                                                                        <  8:110452/149‑1 (MQ=255)
                          gCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGcc                                                                                                                       >  2:40987/1‑148 (MQ=255)
                              cGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATAc                                                                                                                  <  5:110889/149‑1 (MQ=255)
                                   gCTATCGGTCGCCGCAGGCCCGGTGGCCTTTGGCCCTGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCa                                                                                                              <  5:281237/148‑1 (MQ=255)
                                     tATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTc                                                                                                           >  3:345267/1‑149 (MQ=255)
                                     tATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCAGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCACACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTc                                                                                                           >  6:188830/1‑149 (MQ=255)
                                          gTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAAGACGCCAATCCTgcc                                                                                                      >  8:333957/1‑147 (MQ=255)
                                             gccaaaGGCCAGCGGGCCTGCGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACggg                                                                                                   >  3:160839/6‑149 (MQ=255)
                                               cgcAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACggg                                                                                                   >  5:228052/1‑147 (MQ=255)
                                                    gCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGcaca                                                                                            <  1:40987/149‑1 (MQ=255)
                                                      ccGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCAcacc                                                                                          >  6:298055/1‑149 (MQ=255)
                                                         cTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGTGGATTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCAcacccac                                                                                       <  4:66224/149‑1 (MQ=255)
                                                          tGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCAATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCAcacccacc                                                                                      <  4:130760/149‑1 (MQ=255)
                                                            gCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACccacca                                                                                     >  5:145562/1‑148 (MQ=255)
                                                                                          aCGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGGACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGc                                                      >  5:327296/1‑149 (MQ=255)
                                                                                          aCGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCg                                                       >  4:73009/1‑148 (MQ=255)
                                                                                                       cTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCaa                                                                     <  3:62938/121‑1 (MQ=255)
                                                                                                       cTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCaa                                                                     >  4:62938/1‑121 (MQ=255)
                                                                                                       cTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCg                                           >  8:151174/1‑147 (MQ=255)
                                                                                                              cAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAAc                                   <  5:188830/148‑1 (MQ=255)
                                                                                                              cAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAAc                                   <  7:333957/148‑1 (MQ=255)
                                                                                                               aaTACAGAAGAAGTGATTTCCCGCGACCACACATTGATCAAGGTGATCGATGCCTAGTTGGCCAATACGCAATTCCTGAAGTGAAGCAGACCCACAAGAGCGGACATTGCCAACTACTCTTACATCGCTCACGCAAAGGAAGGAAACgt                                 >  4:284848/1‑149 (MQ=255)
                                                                                                                aTACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGTCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACgtg                                >  2:64257/1‑149 (MQ=255)
                                                                                                                aTACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCTTAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACgtg                                >  1:310894/1‑149 (MQ=255)
                                                                                                                aTACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACgtg                                <  6:228052/149‑1 (MQ=255)
                                                                                                                 tACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACgtgt                               >  3:236402/1‑149 (MQ=255)
                                                                                                                   cagacgaACTGATTTCCCGAGCCCACACATTGCTCAAGGTTATCCATGCCGAGTTGGCCAATACGCCATTCATCACTGGACTCCCACCCACCAAACCGGACATTCGCAACTAAACATACATCCCACACGCCACATAAGTCAGCGTATCg                             >  8:355334/3‑149 (MQ=255)
                                                                                                                    cgacgaAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACg                                                 >  4:92210/1‑128 (MQ=255)
                                                                                                                    cgacgaAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACg                                                 <  3:92210/128‑1 (MQ=255)
                                                                                                                        gaAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCAtt                                                                                                            >  2:14566/1‑65 (MQ=255)
                                                                                                                        gaAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCAtt                                                                                                            <  1:14566/65‑1 (MQ=255)
                                                                                                                        gaAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGa                        <  2:310894/149‑1 (MQ=255)
                                                                                                                         aaGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGAg                       <  5:41978/149‑1 (MQ=255)
                                                                                                                           gTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCACCCGAAGGCAACGTGTCGCTGGAGcc                     >  3:221616/1‑149 (MQ=255)
                                                                                                                                   ccGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGAGCCCTATGCCa             <  3:73009/149‑1 (MQ=255)
                                                                                                                                    cgcgCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCAAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGAGCCCTATGCCaa            <  3:284848/149‑1 (MQ=255)
                                                                                                                                      atcgatcacCTTGAGCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGAGCCCTATGCCAAtg          >  3:162834/7‑149 (MQ=255)
                                                                                                                                       tcgatcacCTTGAGCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGAGCCCTATGCCAAtgt         <  4:162834/144‑1 (MQ=255)
                                                                                                                                         ccACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAAc                                   <  4:40801/121‑1 (MQ=255)
                                                                                                                                         ccACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGAGCCCTATGCCAATGTgc       <  6:327296/149‑1 (MQ=255)
                                                                                                                                         ccACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGAGCCCTATGCCAATGTgc       >  3:16992/1‑149 (MQ=255)
                                                                                                                                           acacCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGAGCCCTATGCCAATGTgcgc     <  4:236402/149‑1 (MQ=255)
                                                                                                                                              ccTTGCTCAAGGTGATCGATGCGGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGAGCCCTATGCCAATGTGCGCAAc  <  4:16992/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GAAGAACCTGCCGCCGCCGCCCGGGTGCAACGCTGGCTATCGGTCGCCGCAGGCCCGCTGGCCTTTGGCCCGGCAGCCGCTCGCCTGGTCACGGTGTTCGGGGCTTCGTTCAATACCGACGAAGTGATTTCCCGCGCCCACACCTTGCTCAAGGTGATCGATGCCGAGTTGGCCAATACGCCATTCCTGACGGGCAGCACACCCACCATAGCGGACATTGCCAACTACTCGTACATCGCTCACGCCCCCGAAGGCAACGTGTCGCTGGAGCCCTATGCCAATGTGCGCAAC  >  NC_002947/1336889‑1337179

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