Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,076,055 T→G 6.0% Q22P (CAG→CCG)  fdoH ← formate dehydrogenase‑O, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,076,0550TG 6.0% 2161.7 / 18.1 906Q22P (CAG→CCG) fdoHformate dehydrogenase‑O, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base T (308/522);  new base G (59/8);  total (375/531)
Fisher's exact test for biased strand distribution p-value = 1.92e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 912 total aligned reads displayed.

GCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATG  >  CP009273/4075989‑4076122
                                                                  |                                                                   
    tACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGgc                                                                <  1:46111927/68‑1 (MQ=255)
       aGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGcgcg                                                             >  1:11530468/1‑68 (MQ=255)
       aGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGcgcg                                                             >  1:1388380/1‑68 (MQ=255)
       aGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCCGGGGcgcg                                                             >  1:47115690/1‑68 (MQ=255)
       aGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGACTGAGGcgcg                                                             >  1:30133649/1‑68 (MQ=255)
        gCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGCGGCGCgg                                                            >  1:18892400/1‑68 (MQ=255)
        gCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgg                                                            >  1:54757028/1‑68 (MQ=255)
        gCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgg                                                            >  1:10036285/1‑68 (MQ=255)
        gCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgg                                                            >  1:55279397/1‑68 (MQ=255)
        gCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGCCGCGCgg                                                            >  1:43171986/1‑68 (MQ=37)
        gCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCgg                                                            >  1:30979717/1‑68 (MQ=255)
        gCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCgg                                                            >  1:57964699/1‑68 (MQ=255)
        gCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCCCCTGAGGCGCgg                                                            >  1:40951154/1‑68 (MQ=255)
          cGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgggg                                                          <  1:33069828/68‑1 (MQ=255)
          cGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgggg                                                          <  1:379485/68‑1 (MQ=255)
           gATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggggg                                                         <  1:17896526/68‑1 (MQ=255)
             tACAGCTGGTAACGGCGATGAGCTTCGCCACTTCTTCCTGGCAGACCCGCGCCTGAGGCGCGGGGGt                                                        <  1:55905655/67‑1 (MQ=255)
               cAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCa                                                      >  1:40141914/1‑67 (MQ=255)
               cAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAg                                                     >  1:34161939/1‑68 (MQ=255)
               cAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAg                                                     >  1:12027362/1‑68 (MQ=255)
               cAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCa                                                      >  1:48769256/1‑67 (MQ=255)
               cAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAg                                                     >  1:51898843/1‑68 (MQ=255)
                aGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAg                                                     >  1:26300251/1‑67 (MQ=255)
                 ggtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGCGGCGCGGGGGTCAGAc                                                   >  1:57151139/1‑68 (MQ=255)
                 ggtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAgc                                                    >  1:23443725/1‑66 (MQ=255)
                 ggtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGa                                                    >  1:33716295/1‑67 (MQ=255)
                 ggtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGa                                                    >  1:50308173/1‑67 (MQ=255)
                 ggtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGa                                                    >  1:18483485/1‑67 (MQ=255)
                 ggtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGa                                                    >  1:14295065/1‑67 (MQ=255)
                 ggtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGGCAGa                                                    >  1:3441187/1‑67 (MQ=255)
                   tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGa                                                    <  1:12695028/65‑1 (MQ=255)
                   tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGa                                                    <  1:19275682/65‑1 (MQ=255)
                   tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGa                                                    <  1:40544359/65‑1 (MQ=255)
                   tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGAc                                                   <  1:15099332/66‑1 (MQ=255)
                   tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGAc                                                   <  1:7246880/66‑1 (MQ=255)
                   tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGAc                                                   <  1:53685016/66‑1 (MQ=255)
                   tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGAc                                                   <  1:46537583/66‑1 (MQ=255)
                   tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCg                                                 >  1:23691194/1‑68 (MQ=255)
                   tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCg                                                 >  1:21998968/1‑68 (MQ=255)
                   tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCCCCAGAGGCGCGGGGGGCAGACCg                                                 >  1:1566104/1‑68 (MQ=255)
                      tAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTa                                              >  1:52426635/1‑68 (MQ=255)
                        aCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGGCGGAGGCGCGGGGGTCAGAACGTTAGt                                            >  1:879622/1‑68 (MQ=255)
                        aCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGCGGCGCGGGGGGCAGACCGTTAg                                             >  1:20041656/1‑67 (MQ=255)
                        aCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAg                                             >  1:32112690/1‑67 (MQ=255)
                        aCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAg                                             >  1:6411686/1‑67 (MQ=255)
                        aCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAg                                             >  1:27235704/1‑67 (MQ=255)
                        aCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGt                                            <  1:3056790/68‑1 (MQ=255)
                        aCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGGGGCGCGGGGGTCAGCCCGTTAg                                             >  1:50866346/1‑67 (MQ=255)
                        aCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCAGAGGCGCGGGGGTCAGACCGTTatt                                            >  1:2703970/1‑66 (MQ=255)
                          gTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTc                                           >  1:31423334/1‑67 (MQ=255)
                                gAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAac                                    <  1:29614322/68‑1 (MQ=255)
                                 aGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAacg                                   >  1:20778608/1‑68 (MQ=255)
                                 aGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAacg                                   >  1:56788318/1‑68 (MQ=255)
                                 aGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCCGGGGTCAGACCGTTAGTCGCGGAac                                    >  1:25574649/1‑67 (MQ=255)
                                    ttCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAacgacg                                >  1:6662129/1‑68 (MQ=255)
                                    ttCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAacgac                                 <  1:7603930/67‑1 (MQ=255)
                                     tCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAAcgacga                               <  1:45708125/68‑1 (MQ=255)
                                     tCGCCACATCTTCCTGGAAGGCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAAcgacga                               <  1:49739079/68‑1 (MQ=255)
                                      cGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGaa                              <  1:27747478/68‑1 (MQ=255)
                                       gCCACTTCTTCCTGGAAGGCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAAt                             <  1:48180304/68‑1 (MQ=255)
                                        ccACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGCCCGTTAGTCGCGGAACGACGAATg                            >  1:22601376/1‑68 (MQ=255)
                                         cACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGa                           <  1:51446805/68‑1 (MQ=255)
                                         cACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGa                           <  1:37347654/68‑1 (MQ=255)
                                         cACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGa                           <  1:3673367/68‑1 (MQ=255)
                                         cACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGa                           <  1:35856410/68‑1 (MQ=255)
                                         cACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGa                           <  1:6288956/68‑1 (MQ=255)
                                         cACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGa                           <  1:3254415/68‑1 (MQ=255)
                                         cACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGa                           >  1:16474187/1‑68 (MQ=255)
                                         cACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGa                           <  1:49064313/68‑1 (MQ=255)
                                           cttcttCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGata                         <  1:14786638/68‑1 (MQ=255)
                                            ttcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGatat                        <  1:39756823/68‑1 (MQ=255)
                                            ttcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGatat                        <  1:9308658/68‑1 (MQ=255)
                                            ttcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGatat                        <  1:3908998/68‑1 (MQ=255)
                                             tcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATc                       <  1:43801738/68‑1 (MQ=255)
                                             tcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATc                       <  1:19691299/68‑1 (MQ=255)
                                             tcttcCTGGAACTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATc                       <  1:53011174/68‑1 (MQ=255)
                                              cttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCt                      <  1:56482302/68‑1 (MQ=255)
                                                      aaGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGAtt               <  1:4225540/67‑1 (MQ=255)
                                                      aaGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTg              >  1:34861752/1‑68 (MQ=255)
                                                       aGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGa             <  1:21632803/68‑1 (MQ=255)
                                                       aGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGa             <  1:38588616/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:57611389/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:8068018/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:36237864/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:54147499/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:47596929/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:45078093/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:28447867/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:24925858/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:24413365/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:17501566/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:16146070/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:13432879/68‑1 (MQ=255)
                                                             gcgcCGGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGcc       <  1:10509608/68‑1 (MQ=255)
                                                               gcCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTga                   <  1:42747490/54‑1 (MQ=255)
                                                               gcCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCa      <  1:16987466/67‑1 (MQ=255)
                                                                ccTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCata    >  1:11142325/1‑68 (MQ=255)
                                                                 cTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCata    <  1:44578888/67‑1 (MQ=255)
                                                                 cTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCata    <  1:4162029/67‑1 (MQ=255)
                                                                 cTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCata    <  1:28856826/67‑1 (MQ=255)
                                                                  |                                                                   
GCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATG  >  CP009273/4075989‑4076122

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: