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.95.6% 168.0 / 2.5 45coding (1052/1404 nt)fumCfumarate hydratase (fumarase C),aerobic Class II
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (20/23);  minor base A (1/1);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

CGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCA  >  NC_000913/1685399‑1685647
                                                                                                                                          |                                                                                                              
cGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAgcgc                                                                                                                >  2:899533/1‑139 (MQ=255)
                       cACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCAcgtcgt                                                                                        <  2:249519/139‑1 (MQ=255)
                           gcgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTccc                                                                                    <  2:660591/139‑1 (MQ=255)
                                      gttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCACCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGa                                                                          <  1:1109516/138‑1 (MQ=255)
                                      gttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCACCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGa                                                                          >  2:1109516/1‑138 (MQ=255)
                                           aaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagcag                                                                    >  2:853994/1‑139 (MQ=255)
                                               ttCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCt                                                                            <  2:417650/126‑1 (MQ=255)
                                               ttCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCt                                                                            >  1:417650/1‑126 (MQ=255)
                                                              cAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcac                                                 <  1:899533/139‑1 (MQ=255)
                                                               aGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                <  1:230871/139‑1 (MQ=255)
                                                               aGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                <  1:407685/139‑1 (MQ=255)
                                                                 cAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCt                                                                            >  1:1241230/1‑109 (MQ=255)
                                                                 cAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCt                                                                            <  2:1241230/109‑1 (MQ=255)
                                                                     cgcACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCAcgtcg                                                                                         >  1:185868/1‑92 (MQ=255)
                                                                     cgcACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCAcgtcg                                                                                         <  2:185868/92‑1 (MQ=255)
                                                                       cACGGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAACGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTa                                                          <  2:896836/121‑1 (MQ=255)
                                                                       cACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTa                                                          >  1:896836/1‑121 (MQ=255)
                                                                       cACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGcc                                                     <  2:1022680/126‑1 (MQ=255)
                                                                       cACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGcc                                                     >  1:1022680/1‑126 (MQ=255)
                                                                            aTTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTgg                                       <  1:165823/135‑1 (MQ=255)
                                                                            aTTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTgg                                       >  2:165823/1‑135 (MQ=255)
                                                                                  ggAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTc                             <  1:113294/139‑1 (MQ=255)
                                                                                      aTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCg                         >  2:655930/1‑139 (MQ=255)
                                                                                       ttGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCgg                        <  1:513157/139‑1 (MQ=255)
                                                                                        tgtgGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGc                       >  1:500291/1‑139 (MQ=255)
                                                                                            gATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccccca                                                                                                           <  1:256238/52‑2 (MQ=255)
                                                                                            gATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccccca                                                                                                           >  2:256238/1‑51 (MQ=255)
                                                                                             aTCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCAt                                                                                 >  2:565632/1‑76 (MQ=255)
                                                                                             aTCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCAt                                                                                 <  1:565632/76‑1 (MQ=255)
                                                                                               cACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                >  1:641765/1‑107 (MQ=255)
                                                                                               cACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                <  2:641765/107‑1 (MQ=255)
                                                                                                       gACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagca                                                                     <  2:95779/78‑1 (MQ=255)
                                                                                                       gACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagca                                                                     >  1:95779/1‑78 (MQ=255)
                                                                                                         cGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTc                             >  2:702251/1‑116 (MQ=255)
                                                                                                         cGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTc                             <  1:702251/116‑1 (MQ=255)
                                                                                                          ggAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCAcgtcg                                                                                         <  1:428124/55‑1 (MQ=255)
                                                                                                          ggAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCAcgtcg                                                                                         >  2:428124/1‑55 (MQ=255)
                                                                                                           gAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCt    <  1:655930/139‑1 (MQ=255)
                                                                                                            aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCgt                        >  1:901063/1‑117 (MQ=255)
                                                                                                            aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTc   >  2:1491730/1‑139 (MQ=255)
                                                                                                            aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTc   >  1:1383512/1‑139 (MQ=255)
                                                                                                            aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGAGTTGGGTTCACTTTCCCCgg                        <  2:901063/118‑1 (MQ=255)
                                                                                                            aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCACCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTc   >  2:1085607/1‑139 (MQ=255)
                                                                                                             aGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCa  <  2:500291/139‑1 (MQ=255)
                                                                                                                         caAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTa                                                          >  2:1125012/1‑71 (MQ=255)
                                                                                                                         caAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTa                                                          <  1:1125012/71‑1 (MQ=255)
                                                                                                                                aCCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTccc                           <  2:566220/95‑1 (MQ=255)
                                                                                                                                aCCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTccc                           >  1:566220/1‑95 (MQ=255)
                                                                                                                                          |                                                                                                              
CGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCA  >  NC_000913/1685399‑1685647

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