Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,611,005 (C)5→6 coding (300/759 nt) hyfI → hydrogenase 4, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,611,0001.C95.7% 236.4 / 3.7 70coding (295/759 nt)hyfIhydrogenase 4, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (38/29);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACAC  >  NC_000913/2610902‑2611094
                                                                                                   |                                                                                                
ttAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGc                                                                                                 >  2:155578/1‑100 (MQ=255)
ttAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGc                                                                                                 >  2:2557085/1‑100 (MQ=255)
 tAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGcc                                                                                                <  2:1671199/100‑1 (MQ=255)
     gttgttTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGcccccc                                                                                          <  1:1941367/100‑2 (MQ=255)
     gttgttTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGcccccc                                                                                          >  2:2341586/1‑99 (MQ=255)
     gttgttTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGcccccc                                                                                          <  2:1909403/100‑2 (MQ=255)
       tgttTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGccccccga                                                                                        <  1:2660231/100‑4 (MQ=255)
        gttTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGccccccgat                                                                                       <  1:1056387/100‑5 (MQ=255)
              tcaCCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATaaa                                                                                  >  2:506736/1‑100 (MQ=255)
                   gcgTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAAtttn                                                                             <  1:2806581/100‑2 (MQ=255)
                   gcgTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTtg                                                                             <  2:108970/100‑1 (MQ=255)
                     gTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTtgtg                                                                           <  1:1295469/100‑1 (MQ=255)
                     gTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTtgtg                                                                           <  2:494395/100‑1 (MQ=255)
                     gTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTtgtg                                                                           <  1:2557103/100‑1 (MQ=255)
                     gTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTtgtg                                                                           <  1:2546829/100‑1 (MQ=255)
                        aCGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGttt                                                                        <  1:2800205/100‑1 (MQ=255)
                         cGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTc                                                                       <  2:2557505/100‑1 (MQ=255)
                          gCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTcc                                                                      >  2:293943/1‑100 (MQ=255)
                           ccGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCt                                                                     >  2:1002951/1‑100 (MQ=255)
                             gATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTAc                                                                   <  1:2753516/100‑1 (MQ=255)
                                atTTTGTTATNTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGCTCATAAAATTTGTGTTTCCTACGgc                                                                >  1:1090390/1‑100 (MQ=255)
                                 ttttGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGgcg                                                               >  2:624141/1‑100 (MQ=255)
                                 ttttGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGgcg                                                               <  1:2601652/100‑1 (MQ=255)
                                 ttttGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGgcg                                                               >  1:2252113/1‑100 (MQ=255)
                                  tttGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGgcgc                                                              >  1:768096/1‑100 (MQ=255)
                                  tttGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGgcgc                                                              >  1:1286188/1‑100 (MQ=255)
                                   tttttatttatTGGCGCATTGATCCGGGCGATGCTTATGCCTGCACTTCGGGCGTATGAGTTTGCCTCCCGATCATAAAATTTGTGTTTCCTACGGcgcg                                                             <  2:1006137/97‑1 (MQ=255)
                                     gTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTg                                                           >  1:785161/1‑100 (MQ=255)
                                        aTTTACTGGCGCAGTCACCCGGGCGATGGGTATGCGTGTACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCTTACGGCGCGTgggg                                                        >  2:1525717/1‑100 (MQ=255)
                                        aTTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCgg                                                        >  2:889746/1‑100 (MQ=255)
                                        aTTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCgg                                                        <  2:1077765/100‑1 (MQ=255)
                                         tttACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGgt                                                       <  1:2425141/100‑1 (MQ=255)
                                            aCTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTc                                                    >  2:1052122/1‑100 (MQ=255)
                                             cTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACATCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTcg                                                   <  1:2748540/100‑1 (MQ=255)
                                              tGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTcg                                                   >  2:840663/1‑99 (MQ=255)
                                               ggCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTcggc                                                 >  2:796759/1‑100 (MQ=255)
                                               ggCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTcggc                                                 <  2:2772501/100‑1 (MQ=255)
                                                  gcAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTcggcggc                                              <  1:1117936/100‑1 (MQ=255)
                                                   cAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTcggcggcg                                             >  1:120310/1‑100 (MQ=255)
                                                    aGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTcggcggcgg                                            <  2:2070638/100‑1 (MQ=255)
                                                       cACCCGGGCGACGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTAt                                         >  2:1516953/1‑100 (MQ=255)
                                                         cccGGGCGATGCGTATGCCTGCACTTCGGGCGTATGGGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTAttt                                       >  1:2040763/1‑100 (MQ=255)
                                                         cccGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTAttt                                       >  2:244459/1‑100 (MQ=255)
                                                           cGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTc                                     >  1:158358/1‑100 (MQ=255)
                                                            gggCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTcc                                    <  1:293944/100‑1 (MQ=255)
                                                              gCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCAc                                  >  1:2823979/1‑100 (MQ=255)
                                                              gCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCACCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCAc                                  <  1:1107488/100‑1 (MQ=255)
                                                                gATGCGAATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGa                                >  1:1105305/1‑100 (MQ=255)
                                                                 aTGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGAt                               <  2:768094/100‑1 (MQ=255)
                                                                  tGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCa                                   >  2:2231641/1‑94 (MQ=255)
                                                                   gCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGAtct                             <  1:2370860/100‑1 (MQ=255)
                                                                      tATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTa                          <  1:60028/100‑1 (MQ=255)
                                                                         gCCTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAg                       >  1:2008579/1‑100 (MQ=255)
                                                                          ccTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGc                      >  2:1396392/1‑100 (MQ=255)
                                                                          ccTGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGc                      >  1:703636/1‑100 (MQ=255)
                                                                            tccgatctgggCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGt                    <  2:1061187/92‑1 (MQ=255)
                                                                            tGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGt                    >  1:1611261/1‑100 (MQ=255)
                                                                            tGCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGt                    >  2:1331881/1‑100 (MQ=255)
                                                                             gCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTTCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTc                   >  2:870634/1‑100 (MQ=255)
                                                                             gCACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTc                   <  1:440453/100‑1 (MQ=255)
                                                                              cACTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCt                  <  2:2490094/100‑1 (MQ=255)
                                                                                cTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTgg                <  1:145759/100‑1 (MQ=255)
                                                                                cTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTgg                >  1:2612998/1‑100 (MQ=255)
                                                                                cTTCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTgg                >  2:725700/1‑100 (MQ=255)
                                                                                 ttCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTggg               <  1:1604430/100‑1 (MQ=255)
                                                                                  tCGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGAGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTgggg              >  1:865160/1‑100 (MQ=255)
                                                                                   cGGGCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGc             <  2:967189/100‑1 (MQ=255)
                                                                                    gggCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGtagatcgga            >  1:1061189/1‑92 (MQ=39)
                                                                                     ggCGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCgg           <  1:1134463/100‑1 (MQ=255)
                                                                                       cGTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTa         >  1:1690239/1‑100 (MQ=255)
                                                                                        gTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAg        >  1:903150/1‑100 (MQ=255)
                                                                                        gTATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAg        >  1:1424897/1‑100 (MQ=255)
                                                                                         tATGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGc       >  2:503325/1‑100 (MQ=255)
                                                                                          aTGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCg      <  1:885698/100‑1 (MQ=255)
                                                                                          aTGAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCg      <  1:2341601/100‑1 (MQ=255)
                                                                                            gAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGac    >  2:2073611/1‑100 (MQ=255)
                                                                                            gAGTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCATGCGTGGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGa     >  2:2595182/1‑99 (MQ=255)
                                                                                              gTCTGCCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGacac  >  1:1359549/1‑100 (MQ=255)
                                                                                                   |                                                                                                
TTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACAC  >  NC_000913/2610902‑2611094

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