Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,076,055 T→G 6.4% 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.4% 1264.1 / 12.1 517Q22P (CAG→CCG) fdoHformate dehydrogenase‑O, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base T (173/298);  new base G (38/2);  total (217/300)
Fisher's exact test for biased strand distribution p-value = 9.62e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 520 total aligned reads displayed.

CTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATG  >  CP009273/4075990‑4076122
                                                                 |                                                                   
 tttACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTgcg                                                                  >  1:14990967/1‑66 (MQ=255)
 tttACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAg                                                                  <  1:28219858/68‑1 (MQ=255)
    aCAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGgc                                                                >  1:22200531/1‑67 (MQ=255)
    aCAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGgc                                                                >  1:29653967/1‑67 (MQ=255)
    aCAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGgc                                                                >  1:28647571/1‑67 (MQ=255)
     cAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGgcgc                                                              >  1:17821798/1‑68 (MQ=255)
     cAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGgcgc                                                              >  1:11810865/1‑68 (MQ=255)
      aGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGcgcg                                                             >  1:3500707/1‑68 (MQ=255)
      aGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGcgcg                                                             >  1:6856032/1‑68 (MQ=255)
      aGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGcgcg                                                             >  1:4775877/1‑68 (MQ=255)
      aGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGGGGcgcg                                                             >  1:2936335/1‑68 (MQ=255)
       gCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgg                                                            >  1:4617316/1‑68 (MQ=255)
       gCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGGGCgg                                                            >  1:12313673/1‑68 (MQ=255)
       gCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCCTCCTGGAAGTCCCGCGCCTGAGGCGCgg                                                            >  1:4941326/1‑68 (MQ=255)
       gCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTATTCCTGGAAGTCCCCCGCCAGAGGCGCgg                                                            >  1:9481562/1‑68 (MQ=37)
         ctaTACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGgcgc                                                              <  1:10565105/62‑1 (MQ=255)
         cGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgggg                                                          <  1:18320826/68‑1 (MQ=255)
         cGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgggg                                                          <  1:610643/68‑1 (MQ=255)
         cGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgggg                                                          <  1:8984158/68‑1 (MQ=255)
          gATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGACGTCCCGCGCCTGAGGCGCggggg                                                         <  1:23384339/68‑1 (MQ=255)
          gATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggggg                                                         <  1:18142607/68‑1 (MQ=255)
          gATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggggg                                                         <  1:5991071/68‑1 (MQ=255)
          gATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggggg                                                         <  1:2540665/68‑1 (MQ=255)
            tACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGt                                                        <  1:30972848/67‑1 (MQ=255)
            tACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGt                                                        <  1:25790738/67‑1 (MQ=255)
            tACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCGCCTGAGGCGCGGGGGTc                                                       >  1:30091051/1‑68 (MQ=255)
             aCAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCGCCTGAGGCGCGGGGGTCa                                                      >  1:15208917/1‑68 (MQ=255)
              cAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCGCCTGAGGCGCGGGGGTCAg                                                     >  1:26412290/1‑68 (MQ=255)
                ggtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGa                                                    >  1:19916296/1‑67 (MQ=255)
                ggtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCGCCTGTGCCGCGGGGGTCata                                                    >  1:5488244/1‑65 (MQ=255)
                ggtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCGCCAGAGGCGCCGGGGTCCGa                                                    >  1:19273975/1‑67 (MQ=255)
                 gtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGAcc                                                  >  1:27005629/1‑68 (MQ=255)
                  tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGa                                                    <  1:25617279/65‑1 (MQ=255)
                  tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGa                                                    <  1:919484/65‑1 (MQ=255)
                  tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGa                                                    <  1:15844991/65‑1 (MQ=255)
                  tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGCCCg                                                 >  1:25959180/1‑68 (MQ=255)
                  tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGAc                                                   <  1:7905014/66‑1 (MQ=255)
                  tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAGACCg                                                 >  1:26220796/1‑68 (MQ=255)
                  tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCGCGGGAGGCGCGGGGGTCAGCCCg                                                 >  1:1669358/1‑68 (MQ=255)
                  tcgtAACGTCGACTAGTTTCGACACTTCTACCTGGACGGCCCGCGCCTGAGGCGCGGGGGTCAGa                                                    <  1:21959197/63‑1 (MQ=255)
                  gggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAGa                                                    <  1:11558516/64‑1 (MQ=255)
                   ggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGt                                                <  1:17367045/68‑1 (MQ=255)
                    gtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCGGACCGtt                                               >  1:10786007/1‑68 (MQ=255)
                     tAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGGCAGCCCGttt                                              >  1:2362847/1‑67 (MQ=255)
                     tAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAGCCCGTTa                                              >  1:24151786/1‑68 (MQ=255)
                       aCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAg                                             >  1:10027476/1‑67 (MQ=255)
                       aCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAg                                             >  1:17125987/1‑67 (MQ=255)
                       aCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAg                                             >  1:1969521/1‑67 (MQ=255)
                       aCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAg                                             >  1:3241706/1‑67 (MQ=255)
                       aCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAGCCCGTTAg                                             >  1:20273056/1‑67 (MQ=255)
                        cGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGGCAGACAGTTAGTc                                           >  1:13893625/1‑68 (MQ=255)
                         gTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGCGGCGCGGGGGTCAGCCCGTTAGTc                                           >  1:1337242/1‑67 (MQ=255)
                         gTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTcg                                          >  1:18666794/1‑68 (MQ=255)
                         gTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTc                                           >  1:20681608/1‑67 (MQ=255)
                         gTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTc                                           >  1:11932518/1‑67 (MQ=255)
                          tCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTcgc                                         >  1:31734904/1‑68 (MQ=255)
                          tCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTagc                                         >  1:2149835/1‑68 (MQ=255)
                               gAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAac                                    >  1:24670581/1‑68 (MQ=255)
                                   ttCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAacgac                                 <  1:4052615/67‑1 (MQ=255)
                                   ttCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAacga                                  <  1:27220948/66‑1 (MQ=255)
                                   ttCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAacga                                  <  1:15444484/66‑1 (MQ=255)
                                    tCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAAcgacga                               <  1:10348544/68‑1 (MQ=255)
                                    gcGCCACTTCTTGCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAAcgacga                               <  1:7334524/67‑1 (MQ=255)
                                      gCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAAt                             <  1:6941366/68‑1 (MQ=255)
                                      gCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAAt                             <  1:14733404/68‑1 (MQ=255)
                                       ccACTTCTTCCTGGAAGTCCCGCGCCAGAGGCGCGGGGGGCAGCCCGTTAGTCGCGGAACCACGAATg                            >  1:7559365/1‑68 (MQ=255)
                                        cACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGa                           <  1:5758200/68‑1 (MQ=255)
                                        cACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGa                           >  1:21055508/1‑68 (MQ=255)
                                        cACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGa                           <  1:16135900/68‑1 (MQ=255)
                                        cACTTCTGCCTGGAAGTCCCGCGCCGGCGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGa                           <  1:7477458/68‑1 (MQ=255)
                                           ttcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGatat                        <  1:17645481/68‑1 (MQ=255)
                                           ttcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGatat                        <  1:28357071/68‑1 (MQ=255)
                                            tcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATc                       <  1:22754387/68‑1 (MQ=255)
                                            tcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATc                       <  1:8571932/68‑1 (MQ=255)
                                            tcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATc                       <  1:2438851/68‑1 (MQ=255)
                                            tcttcCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAGACCGTTATTCGCGGAACGCCGAATGata                         >  1:18892002/1‑66 (MQ=255)
                                             cttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCt                      <  1:4267520/68‑1 (MQ=255)
                                                ccTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTga                   <  1:16386436/68‑1 (MQ=255)
                                                   ggAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGAt                >  1:29069486/1‑68 (MQ=255)
                                                   ggAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGAt                >  1:12620971/1‑68 (MQ=255)
                                                    gAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGAtt               >  1:14492324/1‑68 (MQ=255)
                                                      aGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTg              <  1:14633708/67‑1 (MQ=255)
                                                      aGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGa             <  1:14140750/68‑1 (MQ=255)
                                                      aGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGa             <  1:6345725/68‑1 (MQ=255)
                                                          ccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAg         <  1:13182166/68‑1 (MQ=255)
                                                          ccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAg         <  1:25085126/68‑1 (MQ=255)
                                                           cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:12763154/68‑1 (MQ=255)
                                                           cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:11078347/68‑1 (MQ=255)
                                                           cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:3735730/68‑1 (MQ=255)
                                                           cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:18946298/68‑1 (MQ=255)
                                                           cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:31484635/68‑1 (MQ=255)
                                                           cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:21339633/68‑1 (MQ=255)
                                                           cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:30232256/68‑1 (MQ=255)
                                                           cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:22602004/68‑1 (MQ=255)
                                                           cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:22976037/68‑1 (MQ=255)
                                                           cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:23920216/68‑1 (MQ=255)
                                                           cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:28581286/68‑1 (MQ=255)
                                                            gcgcCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGcc       <  1:27952788/68‑1 (MQ=255)
                                                                cTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCata    <  1:13634417/67‑1 (MQ=255)
                                                                 tGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATg  <  1:27554448/68‑1 (MQ=255)
                                                                 |                                                                   
CTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATG  >  CP009273/4075990‑4076122

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