Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,082,696 G→C 6.6% Q22E (CAG→GAG)  fdoH ← formate dehydrogenase O subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,082,6960GC 6.6% 134.9 / 3.0 53Q22E (CAG→GAG) fdoHformate dehydrogenase O subunit beta
Reads supporting (aligned to +/- strand):  ref base G (22/27);  new base C (2/2);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-01

GTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCA  >  NC_000913/4082556‑4082835
                                                                                                                                            |                                                                                                                                           
gtaggcGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAg                                                                                                                                           <  1:315733/138‑1 (MQ=255)
gTTGTCGTTCGACCCAATAGTATTGGCGACGGGAACGCGAATGTCGTTCCACTCTGAACACGCCGCCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAg                                                                                                                                           <  1:273648/143‑1 (MQ=25)
gTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATAAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTgcg                                                                                                                                           >  2:95832/1‑141 (MQ=255)
       tgcacCCCAATGTTTTTCGCGACGGTATCGCGACTGTCGTTAAACGAGCCGCTGTCCGCGGCACTGGGTTTACGGCCTATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCggg                                                                                                                                    <  2:333424/141‑1 (MQ=255)
        acacCCCAATGTCATTGCCGACGGTACCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgggg                                                                                                                                   <  1:65144/143‑1 (MQ=255)
            gccAATGTTAACGCCGCCGGTACCGGGAATGACGTTCCACTCCGACCACGCTCCCTGACAGGCTTTTCAATCGACTCAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGGCCCGCGCCGGAGGCGCGGGGGTCa                                                                                                                               <  2:141118/142‑1 (MQ=255)
            cccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTcg                                                                                                                               >  1:164971/1‑142 (MQ=255)
            cccAATGTTACTGCCGACGGTATCGCGAATGGCGGTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCa                                                                                                                               <  1:360843/143‑1 (MQ=255)
                  gtgcttcccgccgGTAGCGCGAATGTCGTGCCGCTCTGAACACGCCACCTGACAGGCTTCACAGCCGAGACAGGTGGTAACGTCTATGAGCTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGt                                                                                                                         <  2:233006/136‑1 (MQ=25)
                    cttcgcggACGGTATCGCGACTGGCGTTCCACTCGGAACACGCCACCTTGCAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTa                                                                                                                       <  2:406765/136‑1 (MQ=37)
                        gCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAATCCAGCGCCTGAGGCGCGGGGGTCAGACCTTTAGTcg                                                                                                                   >  1:163451/1‑143 (MQ=255)
                         cctgcGGTACGGCGAATGTCGTTGCACTCGGAACAGGCGAGCTGACAGGCGTGACAGCTGGGAGAGGTGGTAACGTCGATGAGTTTCGCCACTTCGTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTcgc                                                                                                                  <  1:162507/139‑1 (MQ=255)
                         ccGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTcgc                                                                                                                  >  1:269828/1‑143 (MQ=255)
                            aCGGTATGGCGAATGGCGTTGCACGCTGAGCGGGCCACTTGACAGGCTTGACAGCCGATACGGGGGGTACCGTCGATGAGTGTCGCCACTGCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGa                                                                                                               <  2:13262/143‑1 (MQ=255)
                            aCGGTAGCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGa                                                                                                               <  2:362108/143‑1 (MQ=255)
                              ggTATCGCGAATGTCGTTCCTCTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAac                                                                                                             <  1:244012/143‑1 (MQ=255)
                                   cgcgAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGACTGAGGCGCGGGGGTCAGACCGTTCGTCGCGGGAcgaaga                                                                                                        >  1:29493/1‑143 (MQ=255)
                                             gTTCCACTCTGAACACGCCACCTGACAGACTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCAGCGCCTGAGGCGCGGAGGTCAGACGGTTGGTCACGGAACGACGAACGATATCtt                                                                                              >  2:105152/1‑143 (MQ=25)
                                               tCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGACCGGGGATGCGGCGGTTAGTCGCGGAACGACGACTGatacactgc                                                                                            >  1:195084/1‑137 (MQ=25)
                                                   ttctGAACTCGACACTTGACAGGCTTTACAGCCGATACAGGTGGTAGCGTCGATGAGTTTCGCCACTTCTTCCTGGCAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGAtt                                                                                        <  1:228880/142‑1 (MQ=38)
                                                   gtctGTACACGCCACCTGACCGGCTTGACGGCCCATACAGGGCGTGACGCCTATGAGTTTCGCCACTTCTGCCTGGAAGTCCCGCGCCTGAGTCGCGGGGGTCAGACCGTGAGTCGCGGAACGACGAATGATATCTTGAGAtt                                                                                        <  2:199255/142‑1 (MQ=255)
                                                     ctGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTGTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGa                                                                                      <  2:51009/143‑1 (MQ=255)
                                                      tGAACACGCCGCCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGAt                                                                                     <  1:153603/143‑1 (MQ=255)
                                                        tccacCCCACCTGGCGGGCTTTACAGCCGAGAGAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATaa                                                                                   <  2:50083/141‑1 (MQ=255)
                                                         acacGCCACCTGTCAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAg                                                                                  <  2:122817/143‑1 (MQ=255)
                                                             gCCACCTGACAGGCTTTACAGCCGATACAGGTGGTACCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGCGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGGTGCGCCat                                                                              >  1:300779/1‑143 (MQ=38)
                                                             gCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCat                                                                              >  1:299171/1‑143 (MQ=255)
                                                                 ccTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATCTCTTGTGCTTGCTAAGCCATCTgt                                                                          >  1:49652/1‑143 (MQ=39)
                                                                  cTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATgtc                                                                         >  2:440570/1‑143 (MQ=255)
                                                                       aGGCTTTACAGCCGATATAGGTGGTAGCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGTGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGtctc                                                                    <  2:179776/143‑1 (MQ=255)
                                                                        ggCTTTACGGCCGAGACAGGTGGTAACGTCGATGCGTTTCGCCACTTCTTCCTGGAAGTCCCGCGTCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTcc                                                                   <  2:193496/143‑1 (MQ=255)
                                                                                  ccGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCtt                                                         >  1:258331/1‑143 (MQ=255)
                                                                                  ccGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCGCCGTTACCCCtt                                                         >  1:366491/1‑143 (MQ=255)
                                                                                    gATACAGGTGGTAACGTCGCTGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCtttt                                                       <  1:396647/143‑1 (MQ=255)
                                                                                         aGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTGTTCcaca                                                  >  2:414512/1‑143 (MQ=255)
                                                                                          ggtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACAt                                                 >  2:363661/1‑143 (MQ=255)
                                                                                            tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCCTGAGATTGATAAGCCATATGTCAACTCTTGTACACCGTTTcaaacttc                                               >  2:350719/1‑136 (MQ=255)
                                                                                            tggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAAGCGCGGGGGTCAGACCGTTAGTCTCGGAACGACGAATGATCTCGTGGGGGCGACAAGCCATATGTCCTCTCCGTTACACCTTTTCaacacta                                               >  2:331920/1‑140 (MQ=255)
                                                                                            tcgtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCTCGGGGGTCAGCCCGGTAGTCGCGGAACGACGAACGATATCTTGAGATTGATAAGCCATATGTGCTCTCCGTTACACGTTTTCTACATTc                                               >  2:226895/3‑143 (MQ=255)
                                                                                                 aCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaag                                          >  1:2211/1‑143 (MQ=255)
                                                                                                    tCGATGAGTTTCGCCACTTCTTCCTGGAAGTCGCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaaggaa                                       <  2:191463/143‑1 (MQ=255)
                                                                                                        tcaGTTTCGTCACTTCTTCCTGGAAGTTCCGCGCCTGGGTCGCGGGGCTCAGACCGTTAGTCTCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGAc                                   <  2:163451/141‑1 (MQ=255)
                                                                                                             ttCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTaaa                              >  1:231048/1‑143 (MQ=255)
                                                                                                             gtctccACTTCTTCCGGGCAGTCCCGCTTCTCGGGCGCGGGGGTCAGACCGTTAGTCGCTGAACGACGAATGATATCTTGAGATTGATTAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTaaa                              <  2:399954/139‑1 (MQ=255)
                                                                                                                        ttcCTGGAAGTCCCTCTCCTGAGGGGTGGAGGTCATGGCGTTAGTTGCGCATCTACTAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCt                   <  2:115550/143‑1 (MQ=255)
                                                                                                                         tcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGCCTTTAGCTCCCCCGTCTg                  >  1:171288/1‑143 (MQ=255)
                                                                                                                           cTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCg                <  1:453981/143‑1 (MQ=255)
                                                                                                                            tGGAGGGCCCGCGCCTGAAGCGCGGGGGTCAGGCCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCgt               <  2:26178/143‑1 (MQ=255)
                                                                                                                             ggAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGGCGGACTTAAACTCCGGCGTCTGCgtg              >  1:42995/1‑143 (MQ=255)
                                                                                                                                 gTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCg          <  1:233016/143‑1 (MQ=255)
                                                                                                                                  ccccGCGCCTCAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGc         >  2:451199/2‑143 (MQ=255)
                                                                                                                                    ccGCGCCTGAGGCGAGGGGGTCAGACCGTTAGTAGCGGAACGACGAATGACATCCTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCAt       <  1:22298/143‑1 (MQ=255)
                                                                                                                                     cgcgCCTCAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATc      >  1:59046/1‑143 (MQ=255)
                                                                                                                                         ccTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGACGTCTCCGTTACACCTTTTCCACATTCACAGGGAAGGACTCACACTCCGGAGTCACCGCGTCCGCATCAcaa  >  2:170707/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                           
GTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCA  >  NC_000913/4082556‑4082835

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑