Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,685,542 (C)6→5 coding (1047/1404 nt) fumC ← fumarate hydratase (fumarase C),aerobic Class II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,685,5370C.100.0% 164.0 / NA 39coding (1052/1404 nt)fumCfumarate hydratase (fumarase C),aerobic Class II
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (20/19);  total (20/19)

GCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAG  >  NC_000913/1685400‑1685662
                                                                                                                                         |                                                                                                                             
gCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGcc                                                                                                                              <  1:611655/139‑1 (MQ=255)
                  tACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCAcg                                                                                                           <  1:1001174/139‑1 (MQ=255)
                          gcgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTccc                                                                                                   <  2:249858/139‑1 (MQ=255)
                          gcgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGCTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTccc                                                                                                   <  1:149309/139‑1 (MQ=255)
                                     gttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAc                                                                                        >  2:527693/1‑139 (MQ=255)
                                          aaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGa                                                                                         <  2:416223/133‑1 (MQ=255)
                                          aaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGa                                                                                         >  1:416223/1‑133 (MQ=255)
                                          aaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagcag                                                                                   >  1:303787/1‑139 (MQ=255)
                                          aaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagcag                                                                                   >  1:391303/1‑139 (MQ=255)
                                            cTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGAAGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCagag                                                                                 >  1:590111/1‑139 (MQ=255)
                                                       aTCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCAcg                                                                                                           >  1:230885/1‑102 (MQ=255)
                                                       aTCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCAcg                                                                                                           <  2:230885/102‑1 (MQ=255)
                                                         cTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGa                                                                                         <  2:991178/118‑1 (MQ=255)
                                                         cTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGa                                                                                         >  1:991178/1‑118 (MQ=255)
                                                         cTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGcc                                                                    <  2:391303/139‑1 (MQ=255)
                                                          tGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCt                                                                   <  2:940488/139‑1 (MQ=255)
                                                              aGCAAGCGCACCGATTGCAGGAAAGTGTGGATCACCATTGGACGGAAGACGTTAAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                               <  1:495647/139‑1 (MQ=255)
                                                                      cACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTg                                                       >  1:284253/1‑139 (MQ=255)
                                                                       aCCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTgg                                                      <  1:527693/139‑1 (MQ=255)
                                                                                            aTCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGAt                                 >  2:289057/1‑139 (MQ=255)
                                                                                                ccATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCAcg                                                                                                           <  1:742925/61‑1 (MQ=255)
                                                                                                ccATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCAcg                                                                                                           >  2:742925/1‑61 (MQ=255)
                                                                                                        cGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCgg                     >  2:470642/1‑139 (MQ=255)
                                                                                                        cGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCgg                     >  2:948635/1‑139 (MQ=255)
                                                                                                           aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGtt                                                                                                                     <  2:962731/40‑1 (MQ=255)
                                                                                                           aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGtt                                                                                                                     >  1:962731/1‑40 (MQ=255)
                                                                                                           aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                               >  2:352764/1‑94 (MQ=255)
                                                                                                           aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                               <  2:11060/94‑1 (MQ=255)
                                                                                                           aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                               >  1:781437/1‑94 (MQ=255)
                                                                                                           aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                               <  2:781437/94‑1 (MQ=255)
                                                                                                           aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                               <  1:352764/94‑1 (MQ=255)
                                                                                                           aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                               >  1:11060/1‑94 (MQ=255)
                                                                                                           aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCAc                                                >  2:449016/1‑109 (MQ=255)
                                                                                                           aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCAc                                                <  1:449016/109‑1 (MQ=255)
                                                                                                           aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTc                  >  1:824133/1‑139 (MQ=255)
                                                                                                            aGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCa                 <  2:824133/139‑1 (MQ=255)
                                                                                                                gttctgttcaAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                               >  1:587273/1‑89 (MQ=255)
                                                                                                                gttcagttcaAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                               <  2:587273/89‑1 (MQ=255)
                                                                                                                 ttcagttcaAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTGATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTc            >  1:639057/1‑139 (MQ=255)
                                                                                                                           aGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTgag  <  1:289057/139‑1 (MQ=255)
                                                                                                                                         |                                                                                                                             
GCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAG  >  NC_000913/1685400‑1685662

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