Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 2,738,435 C→T Y216Y (TAC→TAT)  RPPX_12165 → benzene 1,2‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099742,738,4350CT93.2% 15.5 / 0.3 44Y216Y (TAC→TAT) RPPX_12165benzene 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base T (25/16);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 7.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

ACCACTGGTCGAGCACCTGGGCGAGTCGGCGAAGATCATCGACATGATCGTCGATCAGTCGCCTGAAGGCCTGGAAGTGCTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTCGACCATGGACACCTGCTG  >  CP009974/2738285‑2738582
                                                                                                                                                      |                                                                                                                                                   
gccgcTGGTCGAGGACCTGGGCGAGTCGGCGAAGATCATCGACATGATCGTTGATCAGTCGCCTGCAGGCCTGGAAGTGCTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTAc                                                                                                                                                     <  2:3787806/147‑1 (MQ=0)
gccgcTGGTCGAGCACCTGGGCGAGTCGGCGAAGATCATCGACATGATCGTCGATCAGTCGCCTGAAGGCCTGGAAGTGCTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTAc                                                                                                                                                     >  1:1976348/5‑151 (MQ=1)
 ccgcTGGTCGAGCACCTGGGCGAGTCGGCGAAGATCATCGACATGATCGTCGATCAGTCGCCTGAAGGCCTGGAAGTGCTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTAcc                                                                                                                                                    >  2:1675479/4‑151 (MQ=1)
     tgctCGAGCGCCTGGGGGAGTCGGCTAAGATCATCGACGTGATCGTCGATCAGTCGCCTGAAGGGCTGGAAGTGCTGCGCGGTCCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCGCCGCCGAAAACGGCGCCGACGGCTACCACGt                                                                                                                                                <  2:2200548/148‑1 (MQ=1)
           aGCACCTGGGCGAGTCGGCGAAGATCATCGACATGATCGTCGACCAGTCGCCGGACGGTCTGGAAGTGCTGCGTGGAGCTAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTc                                                                                                                                          <  2:5382427/151‑1 (MQ=1)
                                agctcatcGACATGAGCGTCGATCAGTCGCCTGAAGGCCTGGAAGTGCTGCGCGGTGCCAGCTCGTACATCTACGAAGACAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGc                                                                                                                     <  2:3817212/148‑1 (MQ=0)
                                                    gaccAGTCGCCGGACGGTCTGGAAGTGCTGCGTGGAGCTAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAagc                                                                                                 <  1:950757/148‑1 (MQ=0)
                                                    gaccAGTCGCCGGACGGTCTGGAAGTGCTGCGTGGAGCTAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAagc                                                                                                 >  2:3350667/4‑151 (MQ=1)
                                                                                    tggagctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGACCTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATg                                                                 <  1:1420902/144‑1 (MQ=1)
                                                                                    tggagctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATg                                                                 <  1:635532/144‑1 (MQ=1)
                                                                                    tggagctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATg                                                                 <  1:2678076/144‑1 (MQ=1)
                                                                                      gagctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAg                                                               >  1:423224/6‑151 (MQ=1)
                                                                                       agctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgc                                                              >  2:4290040/5‑150 (MQ=1)
                                                                                        gctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgcg                                                             >  1:3967119/4‑149 (MQ=1)
                                                                                        gctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgcg                                                             <  2:4719174/148‑3 (MQ=1)
                                                                                              tcctACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGgct                                                       >  1:4831024/4‑150 (MQ=0)
                                                                                              tcctACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGgcc                                                       >  2:2158023/4‑150 (MQ=1)
                                                                                                       accggaGGCAACTGGAAGCTCCCCGCCGAAACCGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGAAACTACGCTGCCACCCAGAACCAGCGAAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCa                                              <  2:1457271/146‑1 (MQ=1)
                                                                                                             ggtaaCTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAgg                                        >  2:1779345/4‑151 (MQ=0)
                                                                                                             ggtaaCTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAgg                                        >  2:4119968/4‑151 (MQ=0)
                                                                                                             ggtaaCTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAgg                                        >  2:324255/4‑151 (MQ=0)
                                                                                                                      aaGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGtt                               >  2:2054432/1‑151 (MQ=0)
                                                                                                                         ctgaCCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCt                            <  1:387176/148‑1 (MQ=0)
                                                                                                                         ctgaCCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCt                            <  2:2876788/148‑1 (MQ=0)
                                                                                                                                  gagaaCGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCg                   >  2:2858554/4‑151 (MQ=0)
                                                                                                                                       cGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCAt              >  1:2642600/1‑151 (MQ=2)
                                                                                                                                       cGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCAt              <  1:4204967/151‑1 (MQ=2)
                                                                                                                                        ggCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATg             <  1:3496694/151‑1 (MQ=2)
                                                                                                                                         gcgcCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATgg            >  1:5226025/1‑151 (MQ=2)
                                                                                                                                         gcgcCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATgg            >  2:4214548/1‑151 (MQ=11)
                                                                                                                                         gcgcCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATgg            <  2:3370886/151‑1 (MQ=2)
                                                                                                                                           gcCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGac          <  1:5318288/151‑1 (MQ=2)
                                                                                                                                           gcCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGac          >  2:3806900/1‑151 (MQ=2)
                                                                                                                                           gcCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGac          >  1:5016622/1‑151 (MQ=2)
                                                                                                                                            ccGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGaca         <  1:743505/151‑1 (MQ=2)
                                                                                                                                               aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACct      >  2:2646067/1‑151 (MQ=2)
                                                                                                                                               aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACct      >  2:2644653/1‑151 (MQ=2)
                                                                                                                                               aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACct      <  2:4236934/151‑1 (MQ=11)
                                                                                                                                               aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACct      >  1:2596495/1‑151 (MQ=2)
                                                                                                                                               aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACct      >  2:668703/1‑151 (MQ=2)
                                                                                                                                               aCGGCTATCATGTTAGCGCGGTGCACTGGAACTACGCCGCTAATCGCAACCAGCGCAAGCAGCGCGAGGCAGGCGACTAAGTTAAATCCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTTGTTGTTTCTATTCCTTCGAGCATGGACTCct      <  2:3041600/151‑1 (MQ=2)
                                                                                                                                                cGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACctg     <  2:2365483/151‑1 (MQ=2)
                                                                                                                                                cGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACctg     >  1:2506923/1‑151 (MQ=11)
                                                                                                                                                cGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGGTTCTATTCCTTCCACCATGGACACctg     >  2:4326824/1‑151 (MQ=2)
                                                                                                                                                 ggCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACctgc    >  2:493223/1‑151 (MQ=11)
                                                                                                                                                  gCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACctgct   >  1:4515264/1‑151 (MQ=11)
                                                                                                                                                   cTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTCCTATTCCTTCGACCATGGACACctgctg  >  2:3228592/1‑151 (MQ=11)
                                                                                                                                                      |                                                                                                                                                   
ACCACTGGTCGAGCACCTGGGCGAGTCGGCGAAGATCATCGACATGATCGTCGATCAGTCGCCTGAAGGCCTGGAAGTGCTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTCGACCATGGACACCTGCTG  >  CP009974/2738285‑2738582

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