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% 196.3 / NA 56I414T (ATC→ACC) fumBanaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/29);  total (27/29)

CCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAGG  >  NC_000913/4345963‑4346217
                                                                                                                           |                                                                                                                                   
ccgTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGtt                                                                                                                      >  2:653842/1‑139 (MQ=255)
   tGGGATTGCAGCAGATCCACGTAGGGGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCtt                                                                                                                   >  2:95/1‑139 (MQ=255)
      gATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTAc                                                                                                                <  1:242582/139‑1 (MQ=255)
          gcagcaGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGc                                                                                                            <  1:1188968/139‑1 (MQ=255)
                gATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCgg                                                                                                                           >  1:1229067/1‑118 (MQ=255)
                gATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCgg                                                                                                                           <  2:1229067/118‑1 (MQ=255)
                  tCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCa                                                                                                    <  2:353247/139‑1 (MQ=255)
                        tAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCtt                                                                                              >  2:1506021/1‑139 (MQ=255)
                        tAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCtt                                                                                              <  2:57560/139‑1 (MQ=255)
                                cATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCg                                                                                                           >  1:140304/1‑118 (MQ=255)
                                cATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCg                                                                                                           <  2:140304/118‑1 (MQ=255)
                                cATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCa                                                                                                    >  1:1393929/1‑125 (MQ=255)
                                cATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCa                                                                                                    <  2:1393929/125‑1 (MQ=255)
                                  tACGGCCTGGGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATAGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGc                                                                                                            >  2:360965/1‑115 (MQ=255)
                                  tACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGc                                                                                                            <  1:360965/115‑1 (MQ=255)
                                  tACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGc                                                                                    <  1:1506021/139‑1 (MQ=255)
                                       ccTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGtgtg                                                                               >  2:901062/1‑139 (MQ=255)
                                        cTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGc                                                                                                            <  2:494488/109‑1 (MQ=255)
                                        cTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGc                                                                                                            >  1:494488/1‑109 (MQ=255)
                                         tGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCgg                                                                                                                           >  1:616768/1‑93 (MQ=255)
                                         tGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCgg                                                                                                                           <  2:616768/93‑1 (MQ=255)
                                                gtTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATAtctc                                                                      >  1:523467/1‑139 (MQ=255)
                                                    ggCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGcc                                                                  <  1:901062/139‑1 (MQ=255)
                                                    ggCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGcc                                                                  <  1:95/139‑1 (MQ=255)
                                                     gCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGccc                                                                 >  2:860947/1‑139 (MQ=255)
                                                                   aTGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCa                                                                             <  1:1342135/113‑1 (MQ=255)
                                                                   aTGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCa                                                                             >  2:1342135/1‑113 (MQ=255)
                                                                    tGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCg                                                  >  1:611237/1‑139 (MQ=255)
                                                                        tAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGa                                              >  1:12886/1‑139 (MQ=255)
                                                                         aaCCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAg                                                                                                                        <  2:628970/64‑1 (MQ=255)
                                                                         aaCCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAg                                                                                                                        >  1:628970/1‑64 (MQ=255)
                                                                            cGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGa                                          <  2:523467/139‑1 (MQ=255)
                                                                                aGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGTAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATAtct                                                                       >  1:749808/1‑106 (MQ=255)
                                                                                aGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGTAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATAtct                                                                       <  2:749808/106‑1 (MQ=255)
                                                                                aGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACaaa                                      >  1:286925/1‑139 (MQ=255)
                                                                                  gggTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTAc                                                                                                                <  2:979013/63‑1 (MQ=255)
                                                                                  gggTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTAc                                                                                                                >  1:979013/1‑63 (MQ=255)
                                                                                     ttttCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGt                                 <  2:286925/139‑1 (MQ=255)
                                                                                      tttCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAg                                                                                                                        >  2:1087577/1‑51 (MQ=255)
                                                                                      tttCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAg                                                                                                                        <  1:1087577/51‑1 (MQ=255)
                                                                                       ttCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTgg                               <  2:611237/139‑1 (MQ=255)
                                                                                                   gcgTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTAcc                                                                                                               >  2:995059/1‑47 (MQ=255)
                                                                                                   gcgTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTAcc                                                                                                               <  1:995059/47‑1 (MQ=255)
                                                                                                              cGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGAt                                                            >  2:347591/1‑87 (MQ=255)
                                                                                                              cGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGAt                                                            <  1:347591/87‑1 (MQ=255)
                                                                                                                 gTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGt                                 >  1:21827/1‑111 (MQ=255)
                                                                                                                 gTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGt                                 <  2:21827/111‑1 (MQ=255)
                                                                                                                   gATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCtt                                                                                              <  2:341408/48‑1 (MQ=255)
                                                                                                                   gATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCtt                                                                                              >  1:341408/1‑48 (MQ=255)
                                                                                                                   gATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTgg                               <  1:223108/111‑1 (MQ=255)
                                                                                                                   gATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTgg                               >  2:223108/1‑111 (MQ=255)
                                                                                                                    aTCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATaa                                                          >  1:420966/1‑83 (MQ=255)
                                                                                                                    aTCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATaa                                                          <  2:420966/83‑1 (MQ=255)
                                                                                                                    aTCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAgg  <  1:860947/139‑1 (MQ=255)
                                                                                                                       tttGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATaa                                                          <  2:803850/80‑1 (MQ=255)
                                                                                                                       tttGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATaa                                                          >  1:803850/1‑80 (MQ=255)
                                                                                                                           |                                                                                                                                   
CCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAGG  >  NC_000913/4345963‑4346217

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