Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,684,319 G→A 8.5% E286K (GAG→AAG)  PP_1478 → xenobiotic reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,684,3190GA 8.5% 166.9 / 5.8 59E286K (GAG→AAG) PP_1478xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base G (31/23);  new base A (3/2);  total (34/25)
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.18e-01

GCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCAT  >  NC_002947/1684182‑1684468
                                                                                                                                         |                                                                                                                                                      
gCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCa                                                                                                                                            >  1:156092/1‑150 (MQ=255)
 cTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCaa                                                                                                                                           <  2:156092/150‑1 (MQ=255)
   ggAACCCGGATGAAACCGTCGAACTGGCGCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAgg                                                                                                                                         <  1:290917/150‑1 (MQ=255)
    gAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTCTGCAAAGCGCGtc                                                                                                                                                >  2:852795/1‑141 (MQ=255)
    gAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTCTGCAAAGCGCGtc                                                                                                                                                <  1:852795/142‑2 (MQ=255)
      aCCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAg                                                                                                                                      >  1:286712/1‑150 (MQ=255)
        ccGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCa                                                                                                                                    >  2:434049/1‑150 (MQ=255)
                   cGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCAc                                                                                                                        >  1:1588469/1‑150 (MQ=255)
                     tCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCg                                                                                                                      <  1:1216496/150‑1 (MQ=255)
                       gAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGCGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGAAGAGCGCGTGCGCAAGGAAGCaga                                                                                                                                  >  1:1514969/1‑137 (MQ=255)
                       gAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGCGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGAAGAGCGCGTGCGCAAGGAAGCaga                                                                                                                                  <  2:1514969/137‑1 (MQ=255)
                         aCTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCAc                                                                                                                  >  1:979049/1‑150 (MQ=255)
                              cTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTAGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCg                                                                                                             >  2:117305/1‑150 (MQ=255)
                               tcgtcgCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGCGTTTGAGGAAGCAGAGATCG‑CCACCGGCACCGTCgg                                                                                                            <  1:1365475/150‑1 (MQ=255)
                                 gtcgCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCggaa                                                                                                          <  2:1588469/150‑3 (MQ=255)
                                    gCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGa                                                                                                       >  1:567024/1‑150 (MQ=255)
                                     ccTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGAt                                                                                                      >  1:636279/1‑150 (MQ=255)
                                      cTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAg                                                                                                                                          <  2:1350601/114‑1 (MQ=255)
                                      cTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGCCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAg                                                                                                                                          >  1:1350601/1‑114 (MQ=255)
                                           ggTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGa                                                                                                <  2:567024/150‑1 (MQ=255)
                                                          ccTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAAGCCCGGGGGAGCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGa                                                                                 <  2:631835/150‑1 (MQ=255)
                                                             gATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGca                                                                              >  1:363237/1‑150 (MQ=255)
                                                                  aTGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAAc                                                                                              <  1:377084/129‑1 (MQ=255)
                                                                  aTGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAAc                                                                                              >  2:377084/1‑129 (MQ=255)
                                                                       tcctccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGt                                                                    >  1:1462299/1‑150 (MQ=255)
                                                                         ctccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGAACGGAAACAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTAc                                                                  >  1:690816/1‑150 (MQ=255)
                                                                         ctccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACACGCTTTGCAGAACGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATAATCACTGAACCTGCCCAGACCGAGCACATTCTGCGTAc                                                                  >  2:241651/1‑150 (MQ=255)
                                                                           ccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTg                                                                       <  1:952571/143‑1 (MQ=255)
                                                                           ccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTg                                                                       >  2:952571/1‑143 (MQ=255)
                                                                                     cTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCaa                                                                                                                                           >  2:1632932/1‑66 (MQ=255)
                                                                                     cTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCaa                                                                                                                                           <  1:1632932/66‑1 (MQ=255)
                                                                                     cTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCag                                                                                                                                   <  2:792612/74‑1 (MQ=255)
                                                                                     cTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCag                                                                                                                                   >  1:792612/1‑74 (MQ=255)
                                                                                      tCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGAc                                                     <  1:434049/150‑1 (MQ=255)
                                                                                      tCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTACCCAGGCCGa                                                                                 <  2:479532/122‑1 (MQ=255)
                                                                                      tCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTACCCAGGCCGa                                                                                 >  1:479532/1‑122 (MQ=255)
                                                                                      tCGGTCAACGCCGAAATTCCTACAGGCCATGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGGAAAGATCGAACACC‑GCACCGGCGGTATTATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACGGGGCAGGCCGAc                                                     >  1:372207/1‑150 (MQ=255)
                                                                                                ccGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTcc                                           >  2:1576205/1‑150 (MQ=255)
                                                                                                 cGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATc                                                                                                                              >  2:13242/1‑67 (MQ=255)
                                                                                                 cGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATc                                                                                                                              <  1:13242/67‑1 (MQ=255)
                                                                                                  gAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCAc                                                                                                   <  1:1582311/92‑1 (MQ=255)
                                                                                                  gAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCAc                                                                                                   >  2:1582311/1‑92 (MQ=255)
                                                                                                     aTTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCt                                      >  1:1240766/1‑150 (MQ=255)
                                                                                                         cTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGc                                                                      <  1:995414/114‑1 (MQ=255)
                                                                                                         cTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGc                                                                      >  2:995414/1‑114 (MQ=255)
                                                                                                          tACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACg                                                    >  2:300911/1‑131 (MQ=255)
                                                                                                          tACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACg                                                    <  1:300911/131‑1 (MQ=255)
                                                                                                          tACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGa                                 >  2:794607/1‑150 (MQ=255)
                                                                                                                cccTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACCTTCTGCTTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAgctgct                           >  1:351064/1‑150 (MQ=255)
                                                                                                                      cTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTc                                                                         <  2:126694/98‑1 (MQ=255)
                                                                                                                      cTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTc                                                                         >  1:126694/1‑98 (MQ=255)
                                                                                                                         tCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTAc                                                                  >  1:605531/1‑102 (MQ=255)
                                                                                                                         tCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTAc                                                                  <  2:605531/102‑1 (MQ=255)
                                                                                                                            gACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGTAGCTGCTGCGCGATCCTTa               >  2:1076474/1‑150 (MQ=255)
                                                                                                                            caCGCGCTCTGCAAAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGaa                                                                                               <  1:385555/69‑1 (MQ=38)
                                                                                                                            caCGCGCTCTGCAAAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGaa                                                                                               >  2:385555/2‑70 (MQ=38)
                                                                                                                            caCGCGCTCTGCAAAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTa               >  1:481807/2‑150 (MQ=255)
                                                                                                                             aCGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTAc              <  2:979049/150‑1 (MQ=255)
                                                                                                                              cgcgCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTCATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACt             <  1:1329985/150‑1 (MQ=255)
                                                                                                                                         gagCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCGGCGCGATCCTTACTGGCCGCTGCAt  <  1:794607/150‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                      
GCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCG‑CCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCAT  >  NC_002947/1684182‑1684468

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑