Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,336,732 C→T 9.1% V35V (GTC→GTT)  gstA → glutathione S‑transferase GstA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,336,7320CT 9.1% 121.3 / 2.7 44V35V (GTC→GTT) gstAglutathione S‑transferase GstA
Reads supporting (aligned to +/- strand):  ref base C (20/20);  new base T (2/2);  total (22/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 = 7.48e-01

ATCAGCCCCATCCCGGCGCTGACCAGGAGCCAAGCCCATGTCGAACCCGATCAAGCTCTACAACTTCCCCAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGCCGACTCCAACGCCATCCTCGTCTACCTGGCCAAAAAGTATGACAACGGTACC  >  NC_002947/1336591‑1336879
                                                                                                                                             |                                                                                                                                                   
aTCAGCCCCATCCCGGCGCTGACCAGGAGCCAAGCCCATGTCGAACCCGATCAAGCTCTACAACTTCCCCAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTgg                                                                                                                                              >  8:275916/1‑149 (MQ=255)
aTCAGCCCCATCCCGGCGCTGACCAGGAGCCAAGCCCATGTCGAACCCGATCAAGCTCTACAACTTCCCCAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTgg                                                                                                                                              <  3:227386/149‑1 (MQ=255)
aTCAGCCCCATCCCGGCGCTGACCAGGAGCCAAGCCCATGTCGAACCCGATCAAGCTCTACAACTTCCCCAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTgg                                                                                                                                              >  2:300087/1‑149 (MQ=255)
     ccccATCCCGGCGCTGACCAGGAGCCAAGCCCATGTCGAACCCGATCAAGCTCTACAACTTCCCCAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTGGTCGATCTGGCCa                                                                                                                                           >  5:113849/1‑147 (MQ=255)
     ccccATCCCGGCGCTGACCAGGAGCCAAGCCCATGTCGAACCCGATCAAGCTCTACAACTTCCCCAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCaaa                                                                                                                                         <  5:148887/149‑1 (MQ=255)
                cgcTGACCAGGAGCCAAGCCCATGTCGAACCCGATCAAGCTCTACAACTTCCCCAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACaa                                                                                                                              <  1:272222/149‑1 (MQ=255)
                      ccAGGAGCCAAGCCCATGTCGAACCCGATCAAGCTCTACAACTTCCCCAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGcc                                                                                                                        <  7:138666/149‑1 (MQ=255)
                       cAGGAGCCAAGCCCATGTCGAACCCGATCAAGCTCTACAACTTCCCCAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCg                                                                                                                       <  8:144886/149‑1 (MQ=255)
                       cAGGAGCCAAGCCCATGTCGAACCCGATCAAGCTCTACAACTTCCCCAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCg                                                                                                                       >  2:172449/1‑149 (MQ=255)
                         ggAGCCAAGCCCATGTCGAACCCGATCAAGCTCTACAACTTCCCCAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGa                                                                                                                     >  5:217443/1‑149 (MQ=255)
                          gAGCCAAGCCCATGTCGAACCCGATCAAGCTCTACAACTTCCCCAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGAc                                                                                                                    <  3:135730/149‑1 (MQ=255)
                             ccAAGCCCATGTCGAACCCGATCAAGCTCTACAACTTCCCCAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTc                                                                                                                 <  6:188056/149‑1 (MQ=255)
                                        tCGAACCCGATCAAGCTCTACAACTTCCCCAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCa                                                                                                       <  1:135295/148‑1 (MQ=255)
                                          gAACCCGATCAAGCTCTACAACTTCCCCTAGGCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAAc                                                                                                     <  5:263357/148‑1 (MQ=255)
                                               cGATCAAGCTCTACAACTTCCCCAAGTCTGGCCACGCCCATCGCAACGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTc                                                                                               >  1:231893/1‑149 (MQ=255)
                                                gATCAAGCTCTACAACTTCCCCAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTTGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACt                                                                                                                   >  6:178497/1‑128 (MQ=255)
                                                gATCAAGCTCTACAACTTCCCCAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTTGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACt                                                                                                                   <  5:178497/128‑1 (MQ=255)
                                                              aCTTCCCCAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGt                                                                                 >  5:78675/1‑148 (MQ=255)
                                                                  ccccAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCAt                                                                              <  1:172449/147‑1 (MQ=255)
                                                                   cccAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGa                                                                           >  4:227134/1‑149 (MQ=255)
                                                                    ccAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTGGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGAt                                                                          <  6:113849/149‑1 (MQ=255)
                                                                    ccAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGAt                                                                          <  2:231893/149‑1 (MQ=255)
                                                                    ccAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGAt                                                                          <  7:275916/149‑1 (MQ=255)
                                                                                  gCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGt                                                            >  7:43988/1‑149 (MQ=255)
                                                                                    ccATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTTGATCTGGCCAAAGGCGCGCACAAGCAGACGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGa                                                          >  5:83825/1‑149 (MQ=255)
                                                                                      atcgcatcgAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGAtt                                                        <  6:78675/149‑1 (MQ=255)
                                                                                       tcgcatcgAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTg                                                       >  7:165689/1‑149 (MQ=255)
                                                                                         gcatcgAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGc                                                      <  1:300087/148‑1 (MQ=255)
                                                                                           atcgAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGCCGa                                                   >  5:1257/1‑149 (MQ=255)
                                                                                                    aTGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGCCGACTCCAACGc                                          >  8:14332/1‑149 (MQ=255)
                                                                                                      gCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGCCGACTCCAACGCCa                                        >  7:294066/1‑149 (MQ=255)
                                                                                                      gCTCTCGCTACTGAACCTCCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAACGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGCCGACTCCAACGcc                                         >  7:193468/1‑148 (MQ=255)
                                                                                                          tcGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGCCGACTCCAACGCCATc                                      <  7:14332/147‑1 (MQ=255)
                                                                                                             cTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGCCGACTCCAACGCCATCCTCGt                                 >  3:83757/1‑149 (MQ=255)
                                                                                                             cTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGCCGACTCCAACGCCATCCTCGt                                 >  2:209698/1‑149 (MQ=255)
                                                                                                                cTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGCCGACTCCAACGCCATCCTCGTCTa                              >  7:278677/1‑149 (MQ=255)
                                                                                                                   aaCCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGCCGACTCCAACGCCATCCTCGTCTAc                             <  6:217443/147‑1 (MQ=255)
                                                                                                                   aaCCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGCCGACTCCAACGCCATCCTCGTCTACCt                           >  8:169334/1‑149 (MQ=255)
                                                                                                                            aCCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGCCGACTCCAACGCCATCCTCGTCTACCTGGCCaaaaa                  >  7:125950/1‑149 (MQ=255)
                                                                                                                                aGCTGGTGTTCGTTGATCTGGCCAAAGGCTCGCACAAGCTGCCGGACTTCCTAGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGCCGACTCCAACGCCATCCTCGTCTACCTGGCCAAAAAGTAt              <  6:83825/149‑1 (MQ=255)
                                                                                                                                aGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGCCGACTCCAACGCCATCCTCGTCTACCTGGCCAAAAAGTAt              <  1:209698/149‑1 (MQ=255)
                                                                                                                                aGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGCCGACTCCAACGCCATCCTCGTCTACCTGGCCAAAAAGTAt              <  7:169334/149‑1 (MQ=255)
                                                                                                                                      tgtTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGCCGACTCCAACGCCATCCTCGTCTACCTGGCCAAAAAGTATGACAAc        <  8:193169/149‑1 (MQ=255)
                                                                                                                                        ttCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGCCGACTCCAACGCCATCCTCGTCTACCTGGCCAAAAAGTATGACAACgg      >  7:80606/1‑149 (MQ=255)
                                                                                                                                           gtcgATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGCCGACTCCAACGCCATCCTCGTCTACCTGGCCAAAAAGTa               <  1:159746/137‑1 (MQ=255)
                                                                                                                                           gtcgATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGCCGACTCCAACGCCATCCTCGTCTACCTGGCCAAAAAGTATGACAACGGTa    <  4:83757/148‑1 (MQ=255)
                                                                                                                                           gtcgATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGACGACTCCAACGCCAACCTCGTCTACCTGCACAAAAAGTa               >  2:159746/1‑137 (MQ=255)
                                                                                                                                            tcgATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGCCGACTCCAACGCCATCCTCGTCTACCTGGCCAAAAAGTATGACAACGGTAcc  >  7:61380/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
ATCAGCCCCATCCCGGCGCTGACCAGGAGCCAAGCCCATGTCGAACCCGATCAAGCTCTACAACTTCCCCAAGTCTGGCCACGCCCATCGCATCGAGTTGATGCTCTCGCTACTGAACCTGCCCACCGAGCTGGTGTTCGTCGATCTGGCCAAAGGCGCGCACAAGCAGCCGGACTTCCTCGCCCTCAACCCATTCGGCCAGGTTCCAGTCATCGATGACAACGGCACCGTGATTGCCGACTCCAACGCCATCCTCGTCTACCTGGCCAAAAAGTATGACAACGGTACC  >  NC_002947/1336591‑1336879

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