Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 298,191 T→G 7.7% S32R (AGC→CGC)  paoA ← PaoABC aldehyde oxidoreductase, 2Fe‑2S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273298,1910TG 7.7% 1311.2 / 12.4 434S32R (AGC→CGC) paoAPaoABC aldehyde oxidoreductase, 2Fe‑2S subunit
Reads supporting (aligned to +/- strand):  ref base T (71/323);  new base G (37/2);  total (109/325)
Fisher's exact test for biased strand distribution p-value = 1.08e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 458 total aligned reads displayed.

CTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCT  >  CP009273/298125‑298258
                                                                  |                                                                   
 tGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCa                                                                   <  1:6861164/68‑1 (MQ=255)
 tGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCa                                                                   <  1:31332969/68‑1 (MQ=255)
 tGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCa                                                                   <  1:15774232/68‑1 (MQ=255)
  gCCGGTACGCTTGTCTCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:32732577/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:1732762/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:34544392/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:22892325/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:33134370/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:34991013/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:35243713/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:36901014/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:31536137/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:4389538/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:4557381/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:4820831/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:6035772/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:11110486/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:27623219/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:1460633/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:10084347/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:29933813/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:8899144/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:7140282/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:30229191/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:13889193/68‑1 (MQ=255)
    cGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACt                                                                 <  1:24709661/67‑1 (MQ=255)
       aaCGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTTCGGCGCTCCCtttta                                                             >  1:29453729/1‑66 (MQ=255)
            ttGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGa                                                        <  1:30768787/68‑1 (MQ=255)
            ttGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGa                                                        <  1:17845499/68‑1 (MQ=255)
             tGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCCGCCCGGGTGGGGGGCGCCCTTTAATcc                                                          >  1:21378262/1‑64 (MQ=255)
                cgccAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATcgc                                                    <  1:21618585/68‑1 (MQ=255)
                  ccAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCGCACTTTTATAACATagcg                                                   >  1:5339109/1‑67 (MQ=255)
                  ccAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCGGTTGCGGCGCTCACGTTAATCAGATcgcg                                                   >  1:20226531/1‑67 (MQ=255)
                  ccAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCGGTTGCGGCGCGCACATTAATGAGATggcg                                                   >  1:23402346/1‑67 (MQ=255)
                  ccAGCGTAGAATGAGGATAAACCACGGCGGTCGCCCCTGGTTCGGGGCGCACGTTAATCAGATagcg                                                   >  1:1502582/1‑67 (MQ=255)
                       gTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACggg                                              <  1:25569835/67‑1 (MQ=255)
                       gTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACggg                                              <  1:26927881/67‑1 (MQ=255)
                       gTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACggg                                              <  1:17022565/67‑1 (MQ=255)
                       gTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGt                                             <  1:1603451/68‑1 (MQ=255)
                            aTGAGGATAAACCTCGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:10820656/67‑1 (MQ=255)
                            aTGAGGATAAACCAGGGCGGTCGGCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:9044458/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaaa                                        >  1:7507198/1‑67 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:3864502/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:33499710/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:34281849/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:23633671/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:2329346/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:22590551/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:21036329/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:18864273/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:16685650/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:1548401/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:14902451/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:680518/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:11930519/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:10510846/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCGGTTGGGGCGCGCACATTTATAAGATCGCGACGGGGCAAAc                                        >  1:3712905/1‑68 (MQ=255)
                             tGAGGATAACCCACGGCGGGCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACt                                       <  1:9630656/68‑1 (MQ=255)
                             tGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACt                                       <  1:12517118/68‑1 (MQ=255)
                             tGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACt                                       <  1:23830487/68‑1 (MQ=255)
                             tGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACt                                       <  1:25348225/68‑1 (MQ=255)
                                 gATAAACCACGGCGGTCGCCCCTGTTGCGGCGCGCACTTTTATAAAATAGCGACGGGTGAAAATTaaa                                   >  1:3246222/1‑68 (MQ=255)
                                  aTAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTaa                                    <  1:24414025/66‑1 (MQ=255)
                                  aTAAACCACGGCGGTCGCCGCTGGTGCGGGGCGCACTTTTATAAGATCGCGCCGGGGGAAAATTaaaa                                  >  1:21985879/1‑67 (MQ=255)
                                       ccACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGc                             >  1:7017361/1‑68 (MQ=255)
                                         aCGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCgg                           <  1:18615858/68‑1 (MQ=255)
                                         aCGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCgg                           <  1:35723214/68‑1 (MQ=255)
                                         aCGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCgg                           <  1:26731941/68‑1 (MQ=255)
                                         aCGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCgg                           <  1:28159064/68‑1 (MQ=255)
                                              ggTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGt                      <  1:18233757/68‑1 (MQ=255)
                                              ggTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGt                      <  1:7305240/68‑1 (MQ=255)
                                                  gccgcTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCtt                  <  1:8655310/68‑1 (MQ=255)
                                                  gccgcTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCtt                  <  1:20746658/68‑1 (MQ=255)
                                                  gccgcTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCtt                  <  1:7934730/68‑1 (MQ=255)
                                                    cgcTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTc                 <  1:19160563/67‑1 (MQ=255)
                                                    cgcTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTc                 <  1:27252038/67‑1 (MQ=255)
                                                     gcTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTccc               >  1:29200982/1‑68 (MQ=255)
                                                        gTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCaa             >  1:13514864/1‑67 (MQ=255)
                                                        gTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAAc            <  1:36526033/68‑1 (MQ=255)
                                                        gTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAAc            <  1:13675272/68‑1 (MQ=255)
                                                          tGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAAccc          <  1:6402718/68‑1 (MQ=255)
                                                          tGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAAccc          >  1:11652061/1‑68 (MQ=255)
                                                          tGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAAccc          >  1:11462771/1‑68 (MQ=255)
                                                          tGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAAcc           <  1:8434496/67‑1 (MQ=255)
                                                          tGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTACCAAccc          >  1:24838680/1‑68 (MQ=255)
                                                          tGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCCGCTCGTGCTTCCCAAccc          >  1:6269488/1‑68 (MQ=255)
                                                             ggcgCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGAt       >  1:12337290/1‑68 (MQ=255)
                                                             ggcgCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGAt       <  1:3998286/68‑1 (MQ=255)
                                                             ggcgCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGAt       <  1:1973009/68‑1 (MQ=255)
                                                               cgcTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTa     >  1:16436397/1‑68 (MQ=255)
                                                               cgcTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTa     >  1:27449516/1‑68 (MQ=255)
                                                                gcTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTAt    <  1:31850402/68‑1 (MQ=255)
                                                                  tCAGTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATc   <  1:20580198/67‑1 (MQ=255)
                                                                  tCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATc   <  1:14156495/67‑1 (MQ=255)
                                                                  tCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATc   <  1:4991284/67‑1 (MQ=255)
                                                                  tCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATc   <  1:19486567/67‑1 (MQ=255)
                                                                  tCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATc   <  1:20032307/67‑1 (MQ=255)
                                                                  tCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATc   <  1:33955781/67‑1 (MQ=255)
                                                                  tCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATc   <  1:25976495/67‑1 (MQ=255)
                                                                  tCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCt  <  1:28610004/68‑1 (MQ=255)
                                                                  |                                                                   
CTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCT  >  CP009273/298125‑298258

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