Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,346,086 A→G I414T (ATC→ACC)  fumB ← anaerobic class I fumarate hydratase (fumarase B)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,346,0860AG100.0% 210.6 / NA 58I414T (ATC→ACC) fumBanaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/28);  total (30/28)

ATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCG  >  NC_000913/4345951‑4346214
                                                                                                                                       |                                                                                                                                
aTCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGgtgt                                                                                                                               >  1:680411/1‑139 (MQ=255)
        gccgccGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGa                                                                                                                       <  1:64739/139‑1 (MQ=255)
         ccgccgTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGaa                                                                                                                      >  1:337659/1‑139 (MQ=255)
                            gATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAAt                                                                                                   >  2:212013/1‑139 (MQ=255)
                                    tAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCtt                                                                                           >  2:560346/1‑139 (MQ=255)
                                          tCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCt                                                                                     >  2:1424058/1‑139 (MQ=255)
                                            cATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTg                                                                                   >  1:728546/1‑139 (MQ=255)
                                             aTACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTAc                                                                                                             <  2:821040/112‑1 (MQ=255)
                                             aTACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTAc                                                                                                             >  1:821040/1‑112 (MQ=255)
                                              tACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGc                                                                                 <  2:728546/139‑1 (MQ=255)
                                                 ggCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCgt                                                                               <  1:1137992/138‑1 (MQ=255)
                                                 ggCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCgt                                                                               >  2:1137992/1‑138 (MQ=255)
                                                    cTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTAc                                                                                                             >  2:1275319/1‑105 (MQ=255)
                                                    cTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTAc                                                                                                             <  1:1275319/105‑1 (MQ=255)
                                                            gtTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGa                                                                                                                       >  2:579118/1‑87 (MQ=255)
                                                            gtTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGa                                                                                                                       <  1:579118/87‑1 (MQ=255)
                                                               gggCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGc                                                                >  2:221955/1‑139 (MQ=255)
                                                                ggCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGcc                                                               <  1:212013/139‑1 (MQ=255)
                                                                    aaGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACg                                                           >  1:580514/1‑139 (MQ=255)
                                                                      gTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGc                                                                                                                          >  1:1028628/1‑74 (MQ=255)
                                                                      gTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGc                                                                                                                          <  2:1028628/74‑1 (MQ=255)
                                                                      gTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGAt                                                         >  2:1391910/1‑139 (MQ=255)
                                                                        gAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATaa                                                       <  1:1391910/139‑1 (MQ=255)
                                                                               aTGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCgg                                                                                                                        >  1:653050/1‑67 (MQ=255)
                                                                               aTGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCgg                                                                                                                        <  2:653050/67‑1 (MQ=255)
                                                                                 ggATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTg                                                                                                                                <  1:867270/57‑1 (MQ=255)
                                                                                     aaCCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCtt                                                                                           <  1:115720/90‑1 (MQ=255)
                                                                                     aaCCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCtt                                                                                           >  2:115720/1‑90 (MQ=255)
                                                                                        cGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGa                                       <  1:221955/139‑1 (MQ=255)
                                                                                            aGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTAc                                                                                                             >  1:137228/1‑65 (MQ=255)
                                                                                            aGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTAc                                                                                                             <  2:137228/65‑1 (MQ=255)
                                                                                            aGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACaaa                                   <  2:680411/139‑1 (MQ=255)
                                                                                             ggggTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAAc                                  >  1:1226170/1‑139 (MQ=255)
                                                                                               ggTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTg                                                  <  2:1048675/121‑1 (MQ=255)
                                                                                               ggTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTg                                                  >  1:1048675/1‑121 (MQ=255)
                                                                                                 ttttCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGt                              <  1:1424058/139‑1 (MQ=255)
                                                                                                  tttCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTg                                                                                                                                >  1:632969/1‑40 (MQ=255)
                                                                                                   ttCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTgg                            <  2:337659/139‑1 (MQ=255)
                                                                                                        cGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATAcc                       <  1:560346/139‑1 (MQ=255)
                                                                                                          gACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTgg                                                                                  >  2:106879/1‑78 (MQ=255)
                                                                                                          gACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTgg                                                                                  <  1:106879/78‑1 (MQ=255)
                                                                                                                cgTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGc                                                                                               <  1:922874/59‑1 (MQ=255)
                                                                                                                cgTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGc                                                                                               >  2:922874/1‑59 (MQ=255)
                                                                                                                cgTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGTcc                                                               >  1:908018/1‑91 (MQ=255)
                                                                                                                cgTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGcc                                                               <  2:908018/91‑1 (MQ=255)
                                                                                                                         tCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGc                                                                                 >  2:944726/1‑64 (MQ=255)
                                                                                                                         tCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGc                                                                                 <  1:944726/64‑1 (MQ=255)
                                                                                                                            ggTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACaaa                                   >  2:146220/1‑107 (MQ=255)
                                                                                                                            ggTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACaaa                                   <  1:146220/107‑1 (MQ=255)
                                                                                                                            ggTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTgg                            <  2:999349/114‑1 (MQ=255)
                                                                                                                            ggTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTgg                            >  1:999349/1‑114 (MQ=255)
                                                                                                                             gTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGa                                >  2:1204733/1‑109 (MQ=255)
                                                                                                                             gTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGa                                <  1:1204733/109‑1 (MQ=255)
                                                                                                                             gTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCg  >  2:614941/1‑139 (MQ=255)
                                                                                                                             gTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCg  >  1:1037731/1‑139 (MQ=255)
                                                                                                                               gATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAg                                                                                                >  2:409862/1‑43 (MQ=255)
                                                                                                                               gATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAg                                                                                                <  1:409862/43‑1 (MQ=255)
                                                                                                                               gATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTc                                                                                         <  2:1225126/50‑1 (MQ=255)
                                                                                                                               gATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTc                                                                                         >  1:1225126/1‑50 (MQ=255)
                                                                                                                                       |                                                                                                                                
ATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCG  >  NC_000913/4345951‑4346214

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