Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,076,055 T→G 5.6% 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 5.6% 1247.6 / 9.7 523Q22P (CAG→CCG) fdoHformate dehydrogenase‑O, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base T (174/304);  new base G (31/6);  total (212/311)
Fisher's exact test for biased strand distribution p-value = 1.89e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 524 total aligned reads displayed.

GCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATG  >  CP009273/4075989‑4076122
                                                                  |                                                                   
      cAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGgcgc                                                              >  1:31962102/1‑68 (MQ=255)
      cAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGgcg                                                               >  1:38904019/1‑67 (MQ=255)
       aGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGcgcg                                                             >  1:27979128/1‑68 (MQ=255)
       aGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCGCCGGAGGcgcg                                                             >  1:21512186/1‑68 (MQ=255)
        gCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCgg                                                            >  1:6485763/1‑68 (MQ=255)
        gCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTATTCCTGGAAGTCCCGCGCCGGAGGCGCgg                                                            >  1:9463484/1‑68 (MQ=255)
          cGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgggg                                                          <  1:36300881/68‑1 (MQ=255)
          cGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgggg                                                          <  1:12605180/68‑1 (MQ=255)
          cGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgggg                                                          <  1:22355789/68‑1 (MQ=255)
          cGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgggg                                                          <  1:19396954/68‑1 (MQ=255)
          cGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgggg                                                          <  1:1517604/68‑1 (MQ=255)
          cGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCgggg                                                          <  1:14301814/68‑1 (MQ=255)
          cGATACAGGTGGTAACGTCGAGGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCgggg                                                          <  1:696259/68‑1 (MQ=255)
           gATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGgcgc                                                              <  1:10495332/63‑1 (MQ=255)
           gATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggggg                                                         <  1:30872409/68‑1 (MQ=255)
           gATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggggg                                                         <  1:38014727/68‑1 (MQ=255)
           gATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggggg                                                         <  1:35041748/68‑1 (MQ=255)
           gATACAGGTGGTAACGTCGAGGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggggg                                                         <  1:27110158/68‑1 (MQ=255)
           gATACAGGTGGTAACGGCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgggg                                                          <  1:19991856/67‑1 (MQ=255)
             tACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgggg                                                          <  1:21778331/65‑1 (MQ=255)
             tACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGt                                                        <  1:25513046/67‑1 (MQ=255)
               cAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGCGGCGCGGGGGTCa                                                      >  1:21153854/1‑67 (MQ=255)
                 ggtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGa                                                    >  1:33604372/1‑67 (MQ=255)
                 ggtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCAGAGGCGCGGGGGTCAGa                                                    >  1:15865596/1‑67 (MQ=255)
                  gtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCCCCTGAGGCGCGGGGGTCAGAcc                                                  >  1:17394402/1‑68 (MQ=255)
                   tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGa                                                    <  1:11098498/65‑1 (MQ=255)
                   tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGa                                                    <  1:4006448/65‑1 (MQ=255)
                   tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGAcc                                                  <  1:29673904/67‑1 (MQ=255)
                   tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGAc                                                   <  1:31807661/66‑1 (MQ=255)
                   tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCg                                                 >  1:35532945/1‑68 (MQ=255)
                   tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCg                                                 >  1:28549268/1‑68 (MQ=255)
                   tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCg                                                 >  1:10841841/1‑68 (MQ=255)
                   tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCg                                                 >  1:25584186/1‑68 (MQ=255)
                   tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGGCAGACCg                                                 >  1:7601734/1‑68 (MQ=255)
                   tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGCGGCGCGGGGGTCAGAc                                                   >  1:34132035/1‑66 (MQ=255)
                        aCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGt                                            <  1:30040643/68‑1 (MQ=255)
                        aCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGt                                            <  1:16622561/68‑1 (MQ=255)
                        aCGTCGATGAGTTGCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGt                                            <  1:17857138/68‑1 (MQ=255)
                         cGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTc                                           >  1:22042622/1‑68 (MQ=255)
                         cGTCGAAGAGTTTCGCCACTTCTTCCTGGAAGTGCCGCGGCGGAGGGCCGGGGGGCAGACCGTTatta                                           >  1:4805748/1‑65 (MQ=255)
                          gTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGGCCGTTAGTc                                           >  1:3774995/1‑67 (MQ=255)
                          gTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCGCCTGAGGCGCGGGGGTCAGACCGTTAGTc                                           >  1:17116861/1‑67 (MQ=255)
                               tGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCggc                                      >  1:25778844/1‑66 (MQ=255)
                                gAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAac                                    >  1:38172880/1‑68 (MQ=255)
                                 aGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAacg                                   >  1:37119736/1‑68 (MQ=255)
                                 aGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAac                                    >  1:16304180/1‑67 (MQ=255)
                                 aGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGGCAGCCCGTTACTCGCGGaagg                                   >  1:2316250/1‑66 (MQ=255)
                                    ttCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAacgac                                 <  1:3324470/67‑1 (MQ=255)
                                    ttCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAacgac                                 <  1:5924065/67‑1 (MQ=255)
                                    ttCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAacgac                                 <  1:21005744/67‑1 (MQ=255)
                                    ttCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAacga                                  <  1:30303061/66‑1 (MQ=255)
                                     tCGCCACTTCTTCCTGGACGTCCCGCGCCGGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAAcgacga                               <  1:5427361/68‑1 (MQ=255)
                                       gCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAAt                             <  1:29305935/68‑1 (MQ=255)
                                       gCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCCGGGGTCAGCCCGTTAGTCGCGGAACGACGAAt                             >  1:8909883/1‑68 (MQ=255)
                                       gCCACTTCTTCCTGGAAGTCCCGCCCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGAGGAAt                             >  1:7244012/1‑68 (MQ=255)
                                         cACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGa                           <  1:19667765/68‑1 (MQ=255)
                                           cttcttCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGata                         <  1:12123901/68‑1 (MQ=255)
                                           cttcttCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGata                         <  1:34562825/68‑1 (MQ=255)
                                            ttcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGatat                        <  1:11562883/68‑1 (MQ=255)
                                             tcttcCTGGACGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATc                       <  1:19879985/68‑1 (MQ=255)
                                             tcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATc                       <  1:931101/68‑1 (MQ=255)
                                             tcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATc                       <  1:36023563/68‑1 (MQ=255)
                                                 ccTGTATGTCGCGCGGCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTga                   <  1:8214252/68‑1 (MQ=255)
                                                 ccTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTga                   >  1:22647228/1‑68 (MQ=255)
                                                 ccTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTga                   <  1:25100793/68‑1 (MQ=255)
                                                 ccTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTga                   <  1:22875645/68‑1 (MQ=255)
                                                   tGGACGGCCCGCGCCTGAGGCGCGGGGGTCAGACCGGTAGTCGCGGAACGACGAATGATATCTTgaga                 <  1:1297973/68‑1 (MQ=255)
                                                     gAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGAtt               >  1:28973066/1‑68 (MQ=255)
                                                     gAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATATTGAGAtt               >  1:13474005/1‑68 (MQ=255)
                                                      aaGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGAtt               <  1:31426162/67‑1 (MQ=255)
                                                      aaGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGAtt               <  1:24384713/67‑1 (MQ=255)
                                                       aGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTg              <  1:10043551/67‑1 (MQ=255)
                                                       aGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTg              <  1:24673176/67‑1 (MQ=255)
                                                       aGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTg              <  1:15410122/67‑1 (MQ=255)
                                                       aGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGa             <  1:4233478/68‑1 (MQ=255)
                                                       aGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCCCGGTACGACGAATGATATCTTGAGATTg              <  1:8498579/67‑1 (MQ=255)
                                                           ccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAg         <  1:17557119/68‑1 (MQ=255)
                                                           ccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGGCGCGGAACGACGAATGATATCTTGAGATTGATAAg         <  1:18529231/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGTGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:35371330/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:23394972/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:23645761/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:24068268/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:7880577/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:26392582/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:10618903/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:20566547/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:5012859/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:1182808/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:32542960/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:32798523/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:19151314/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:18937010/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:35826358/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGGCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:13906369/68‑1 (MQ=255)
                                                            cgcgCCTGAGGCGCGGGGGGCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc        <  1:25993614/68‑1 (MQ=255)
                                                             gcgcCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGcc       <  1:35251874/68‑1 (MQ=255)
                                                                ccTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCcta    >  1:32353175/1‑68 (MQ=255)
                                                                ccTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCata    >  1:18237390/1‑68 (MQ=255)
                                                                 cTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCatat   >  1:33905632/1‑68 (MQ=255)
                                                                  tGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATg  <  1:36959705/68‑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 < 15 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: