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,4350CT97.7% 17.8 / NA 44Y216Y (TAC→TAT) RPPX_12165benzene 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (26/17);  total (27/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCACTGGTCGAGCACCTGGGCGAGTCGGCGAAGATCATCGACATGATCGTCGATCAGTCGCCTGAAGGCCTGGAAGTGCTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTCGACCATGGACACCTG  >  CP009974/2738285‑2738579
                                                                                                                                                      |                                                                                                                                                
gccgcTGGTCGAGCACCTGGGCGAGTCGGCGAAGATCATCGACATGATCGTCGATCAGTCGCCTGAAGGCCTGGAAGTGCTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTAc                                                                                                                                                  <  1:760952/147‑1 (MQ=1)
 ccgctcgTCGAGCACCTGGGCGAGTCGGCGAAGCTCATCCACATGATCGTCGATCAGTCGCCTGAAGGCCTGGAAGTGCTGCGCGGTGCCAGCGCGTACATCTACGACGGCACCTGGCCGCTCACCGCCCGAAACGGCGCCGACGGCTccc                                                                                                                                                 >  2:4346206/7‑151 (MQ=1)
 ccgcTGGTCGAGCACCTGGGCGAGTCGGCGAAGATCATCGACATGATCGTCGATCAGTCGCCTGAAGGCCTGGAAGTGCTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTAcc                                                                                                                                                 >  2:5350450/4‑151 (MQ=1)
 ccgcTGGTCGAGCACCTGGGCGAGTCGGCGAAGATCATCGACATGATCGTCGATCAGTCGCCTGAAGGCCTGGAAGTGCTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTAcc                                                                                                                                                 <  2:4126059/148‑1 (MQ=1)
                                                    gaccAGTCGCCGGACGGTCTGGAAGTGCTGCGTGGAGCTAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAagc                                                                                              >  1:2415216/4‑151 (MQ=0)
                                                    gaccAGTCGCCGGACGGTCTGGAAGTGCTGCGTGGAGCTAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAagc                                                                                              <  1:5751820/148‑1 (MQ=0)
                                                                      cTGGAAGTGCTGCGTGGAGCTAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACg                                                                            >  1:7084630/1‑151 (MQ=1)
                                                                                    tggagctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATg                                                              >  1:2906196/8‑151 (MQ=1)
                                                                                    tggagctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATg                                                              >  1:359946/8‑151 (MQ=1)
                                                                                    tggagctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATg                                                              <  1:5002263/144‑1 (MQ=1)
                                                                                     ggagctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAGTCAAAACCATGa                                                             >  2:1360147/7‑151 (MQ=1)
                                                                                      gagctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAg                                                            <  1:2083604/146‑1 (MQ=1)
                                                                                      gagctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTACCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGaa                                                            >  2:5660016/6‑150 (MQ=1)
                                                                                      gagctaGCTCCTACATCGACGAAGGTGACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAg                                                            <  2:110436/146‑1 (MQ=1)
                                                                                       agctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgc                                                           <  1:6709857/147‑2 (MQ=1)
                                                                                        gctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgcg                                                          >  1:2869055/4‑149 (MQ=1)
                                                                                        gctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgcg                                                          >  2:3047589/4‑149 (MQ=1)
                                                                                              tcctACATCTACGAAGGTAACTGGACGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCCCTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGgct                                                    <  2:1706468/148‑2 (MQ=0)
                                                                                                                         ctgaCCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCt                         >  1:523479/4‑151 (MQ=0)
                                                                                                                                  gagaaCGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCg                >  1:3622923/4‑151 (MQ=0)
                                                                                                                                  gagaaCGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCg                >  1:3695646/4‑151 (MQ=0)
                                                                                                                                  gagaaCGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCg                >  1:5920433/4‑151 (MQ=0)
                                                                                                                                  gagaaCGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCg                >  1:6030034/4‑151 (MQ=0)
                                                                                                                                  gagaaCGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCg                <  2:1672564/148‑1 (MQ=1)
                                                                                                                                       cGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCAt           >  2:3785492/1‑151 (MQ=2)
                                                                                                                                       cGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCAt           <  2:3744259/151‑1 (MQ=11)
                                                                                                                                        ggCGCTGACGGCTATCTTGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACTAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGTGCTAAGCAGGGTGGTGGTTTCTATTTCTTCTACCAtt          >  2:1718050/1‑150 (MQ=2)
                                                                                                                                        ggCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCAtt          >  1:3317212/1‑150 (MQ=2)
                                                                                                                                        ggCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATg          <  2:2684838/151‑1 (MQ=2)
                                                                                                                                        ggCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATg          <  2:2477312/151‑1 (MQ=2)
                                                                                                                                        ggCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATg          <  2:7021841/151‑1 (MQ=2)
                                                                                                                                        ggCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATg          >  2:4213504/1‑151 (MQ=11)
                                                                                                                                        ggCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATg          >  2:1717807/1‑151 (MQ=2)
                                                                                                                                        ggCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATg          >  2:1419084/1‑151 (MQ=2)
                                                                                                                                        ggCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGGTTCTATTCCTTCGACCATg          >  1:2965556/1‑151 (MQ=11)
                                                                                                                                         gcgcCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATgg         <  2:1822900/151‑1 (MQ=2)
                                                                                                                                            ccGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGaca      >  1:2410558/1‑151 (MQ=2)
                                                                                                                                             cGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGCTTCTATTCCTTCGACCATGGacac     >  2:503022/1‑151 (MQ=2)
                                                                                                                                              gACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACAcc    >  2:2642154/1‑151 (MQ=11)
                                                                                                                                              gACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACAcc    >  1:4315278/1‑151 (MQ=2)
                                                                                                                                              gACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACAcc    >  2:1377682/1‑151 (MQ=11)
                                                                                                                                               aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACAccg   >  1:1985377/1‑150 (MQ=2)
                                                                                                                                               aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACct   <  2:1532646/151‑1 (MQ=2)
                                                                                                                                               aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACct   <  1:5771364/151‑1 (MQ=2)
                                                                                                                                               aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACct   <  1:1718050/151‑1 (MQ=11)
                                                                                                                                               aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACct   <  1:1717807/151‑1 (MQ=2)
                                                                                                                                                cGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACctg  >  1:5768689/1‑151 (MQ=2)
                                                                                                                                                cGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACctg  <  2:580292/151‑1 (MQ=2)
                                                                                                                                                      |                                                                                                                                                
ACCACTGGTCGAGCACCTGGGCGAGTCGGCGAAGATCATCGACATGATCGTCGATCAGTCGCCTGAAGGCCTGGAAGTGCTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTCGACCATGGACACCTG  >  CP009974/2738285‑2738579

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