Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 298,191 T→G 10.1% 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,1910TG10.1% 1348.1 / 22.0 465S32R (AGC→CGC) paoAPaoABC aldehyde oxidoreductase, 2Fe‑2S subunit
Reads supporting (aligned to +/- strand):  ref base T (73/337);  new base G (52/2);  total (126/339)
Fisher's exact test for biased strand distribution p-value = 1.76e-31
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 495 total aligned reads displayed.

CTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCT  >  CP009273/298125‑298258
                                                                  |                                                                   
cTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCGGGTGCGGCGCTc                                                                    >  1:27588311/1‑68 (MQ=255)
 tGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCa                                                                   <  1:36343001/68‑1 (MQ=255)
 tGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCa                                                                   <  1:37929585/68‑1 (MQ=255)
 tGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGGTGCGGCGCTCa                                                                   <  1:6703982/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:29973073/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:5903380/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:15518235/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:22190677/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:31940908/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:24295581/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:12783369/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:38311425/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:11884426/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:11374611/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:26320098/68‑1 (MQ=255)
  gCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGGCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:32867323/68‑1 (MQ=255)
  gCCGGAACGCTTGCAGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCAc                                                                  <  1:27069096/68‑1 (MQ=255)
    cGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACtt                                                                <  1:30999785/68‑1 (MQ=255)
    cGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACtt                                                                <  1:5541312/68‑1 (MQ=255)
    cGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACt                                                                 <  1:8236095/67‑1 (MQ=255)
     ggAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACttt                                                               <  1:15102381/68‑1 (MQ=255)
            ttGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTGGCCGCTGTTGCGGCGCTCACTTTAATCAGa                                                        <  1:37189253/68‑1 (MQ=255)
            ttGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGa                                                        <  1:22819442/68‑1 (MQ=255)
            ttGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGa                                                        <  1:18226512/68‑1 (MQ=255)
            ttGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGa                                                        <  1:34966891/68‑1 (MQ=255)
            ttGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGa                                                        <  1:7209474/68‑1 (MQ=255)
                       gTAGAATGAGTATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTATTCAGATGGCGACggg                                              <  1:4885565/67‑1 (MQ=255)
                       gTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACggg                                              <  1:13498029/67‑1 (MQ=255)
                       gTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACggg                                              <  1:20480763/67‑1 (MQ=255)
                       gTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACgg                                               <  1:6419131/66‑1 (MQ=255)
                         aGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACATTAATCAGAACGCGACGGGtga                                           >  1:26033791/1‑66 (MQ=255)
                           aaTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:10788129/68‑1 (MQ=255)
                           aaTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:31699327/68‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTTGGGCGCGCACATTAATCAGCTCGACACCGGTCaaaa                                        >  1:16413638/1‑67 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:7620038/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:37464809/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:25712043/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:12306864/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:38703551/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:7897864/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:14136437/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:133681/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:24498875/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:22565082/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:14732934/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:21210305/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:35537743/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:34494907/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:17962396/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCaaa                                         <  1:30245813/67‑1 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCGGTTGCGGCGCTCACTTTAATAAGATAGCGCCGGGTCAAAc                                        >  1:18593906/1‑68 (MQ=255)
                            aTGAGGATAAACCACGGCGGTCGCCGCGGTTGCGGCGCGCACATTAATAAGATCGCGCCGGGTGaa                                          >  1:30718620/1‑66 (MQ=255)
                             tGAGGATAAACCACGGCGGGCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACt                                       <  1:27762509/68‑1 (MQ=255)
                              gAGGATAAACCACGGCGGTCGCCGCGGTTGCGGCGCGCACTTTAATAAGCTCGCGCCCGGGCaaaaat                                      >  1:17382765/1‑65 (MQ=255)
                                 gATAAACCACGGCGGTCGCCGCGGTTGCGGGGCGCACATTAATCAGATGGCGCCGGGGGAAACATaaa                                   >  1:17612029/1‑68 (MQ=255)
                                 gATAAACCACGGCGGTCGCCGCGGTTGCGGCGCTCACTTTAAt                                                            >  1:35935855/1‑43 (MQ=255)
                                 gATAAACCACGGCGGTCGCCGCGGTTGCGGCGCTCAAGTTAATCAGCTGGCGGCGGGTGAAACGTaaa                                   >  1:21902372/1‑68 (MQ=255)
                                 gATAAACCACGGCGGTCGCCCCGGTTGCGGCGCTCACTTTAAAAAGAAAGCGGCGGGGCAAAATaaaa                                   >  1:29016397/1‑68 (MQ=255)
                                 gATAAACCACGGCGGTCCCCGCGGTTGACGCGCTGACTTTAAAAATAAAGCCACGGGTGAAACTTaaa                                   >  1:24803558/1‑68 (MQ=255)
                                 gATAAACCACGGCGGTCCCCGCGGGTGCGGCCCGCACATTAATGAGATCGCGCCGGGGCCAACTTtaa                                   >  1:29220772/1‑68 (MQ=255)
                                    aaaCCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCg                                <  1:5119478/68‑1 (MQ=255)
                                     aaCCACGGCGGTCGCCGCTGTTGGGGCGCGCACGTTTATCAGCTAGCGACGGGTCAAACTTAAAACGt                               >  1:28627032/1‑68 (MQ=255)
                                     aaCCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGt                               <  1:16659453/68‑1 (MQ=255)
                                     aaCCACGGCGGTCGCCGCGGTTGCGGGGCGCACGTTTATCAGATAGCGACGGGTGAAACTTAAATAGt                               >  1:3734975/1‑68 (MQ=255)
                                      aCCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATAAGATAGCGCCGGGGCAAACTTAAATAGt                               >  1:5324592/1‑67 (MQ=255)
                                       ccACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATAGCGACGGGTCAAACTTAAATCGTGc                             >  1:38933967/1‑68 (MQ=255)
                                         aCGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCgg                           <  1:36941723/68‑1 (MQ=255)
                                         aCGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCgg                           <  1:33108228/68‑1 (MQ=255)
                                         aCGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCgg                           <  1:4051379/68‑1 (MQ=255)
                                         aCGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCgg                           <  1:28095800/68‑1 (MQ=255)
                                          cggcggTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGc                          <  1:23193934/68‑1 (MQ=255)
                                          cggcggTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGc                          <  1:15989461/68‑1 (MQ=255)
                                              ggTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGt                      <  1:32426911/68‑1 (MQ=255)
                                                 cgccgcTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCt                   <  1:33753357/68‑1 (MQ=255)
                                                  gccgcTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCtt                  <  1:33532268/68‑1 (MQ=255)
                                                    cgcTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTc                 <  1:10424114/67‑1 (MQ=255)
                                                    cgcTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTc                 <  1:3140193/67‑1 (MQ=255)
                                                    cgcTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTc                 <  1:26914298/67‑1 (MQ=255)
                                                    cgcTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTc                 <  1:2617516/67‑1 (MQ=255)
                                                    cgcTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTc                 <  1:23627098/67‑1 (MQ=255)
                                                     gcTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTccc               >  1:12039186/1‑68 (MQ=255)
                                                          tGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAcccc          >  1:19681771/1‑68 (MQ=255)
                                                          tGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAAccc          >  1:9119251/1‑68 (MQ=255)
                                                          tGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAAccc          >  1:10534795/1‑68 (MQ=255)
                                                          tGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAAccc          >  1:19118349/1‑68 (MQ=255)
                                                          tGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAAcc           <  1:25097427/67‑1 (MQ=255)
                                                          tGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAAcc           <  1:23781968/67‑1 (MQ=255)
                                                             ggcgCTCACTTTAATCAGATCGCGCCGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGAt       >  1:11056956/1‑68 (MQ=255)
                                                             ggcgCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGa        <  1:17269690/67‑1 (MQ=255)
                                                             ggcgCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGAt       >  1:32698577/1‑68 (MQ=255)
                                                             ggcgCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGAt       <  1:20045434/68‑1 (MQ=255)
                                                             ggcgCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGAt       >  1:2525846/1‑68 (MQ=255)
                                                             ggcgCTCACTTTAATCAGATCGCGACGGGGCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGAt       <  1:8772515/68‑1 (MQ=255)
                                                               cgcTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTa     >  1:26647394/1‑68 (MQ=255)
                                                               cgcTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTa     >  1:10006772/1‑68 (MQ=255)
                                                               cgcTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTa     >  1:9606225/1‑68 (MQ=255)
                                                                  tCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATc   <  1:18994506/67‑1 (MQ=255)
                                                                  tCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATc   <  1:16943053/67‑1 (MQ=255)
                                                                  tCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATc   <  1:34189285/67‑1 (MQ=255)
                                                                  tCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCt  >  1:20746153/1‑68 (MQ=255)
                                                                  |                                                                   
CTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCT  >  CP009273/298125‑298258

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