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

CTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAT  >  NC_000913/1685401‑1685664
                                                                                                                                        |                                                                                                                               
cTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccc                                                                                                                               >  1:1231687/1‑139 (MQ=255)
cTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccc                                                                                                                               <  1:877386/139‑1 (MQ=255)
      aTTCAGTTCAATACCCACTGCGCAGTGTTTGTTAAGACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGaccccatgt                                                                                                                         <  2:294494/139‑5 (MQ=255)
       ttCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGtt                                                                                                                       <  2:24173/139‑1 (MQ=255)
                aTACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCAc                                                                                                              <  2:1520007/139‑1 (MQ=255)
                     cACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCAcgtcgt                                                                                                         <  1:478816/139‑1 (MQ=255)
                         gcgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTccc                                                                                                     >  2:102550/1‑139 (MQ=255)
                                       aaaaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagc                                                                                       >  2:1302306/1‑139 (MQ=255)
                                             ttCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCa                                                                                 <  2:230598/139‑1 (MQ=255)
                                                   gCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTa                                                                           <  2:1191855/139‑1 (MQ=255)
                                                      aTCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCAcgtcg                                                                                                          >  2:206509/1‑105 (MQ=255)
                                                      aTCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCAcgtcg                                                                                                          <  1:206509/105‑1 (MQ=255)
                                                            cAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcac                                                                  >  1:372537/1‑139 (MQ=255)
                                                             aGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                 <  2:1231687/139‑1 (MQ=255)
                                                               cAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCACCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACt                                                               >  2:1013573/1‑139 (MQ=255)
                                                                     cACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTg                                                         >  2:1028398/1‑139 (MQ=255)
                                                                        cGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGt                                                      <  1:1028398/139‑1 (MQ=255)
                                                                         gATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCAcgtcgt                                                                                                         >  1:1112544/1‑87 (MQ=255)
                                                                         gATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCAcgtcgt                                                                                                         <  2:1112544/87‑1 (MQ=255)
                                                                         gATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGtt                                                     >  1:1112779/1‑139 (MQ=255)
                                                                            tGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccc                                                                                                                               >  2:952113/1‑63 (MQ=255)
                                                                            tGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccc                                                                                                                               <  1:952113/63‑1 (MQ=255)
                                                                            tGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTccc                                                                                                     >  2:1286847/1‑88 (MQ=255)
                                                                            tGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTccc                                                                                                     <  1:1286847/88‑1 (MQ=255)
                                                                             gCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACt                                                 <  2:1566560/139‑1 (MQ=255)
                                                                                ggAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTc                                              <  2:1365767/139‑1 (MQ=255)
                                                                                  aaaTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGAt                                                                                                                    >  2:213188/1‑67 (MQ=255)
                                                                                  aaaTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGAt                                                                                                                    <  1:213188/67‑1 (MQ=255)
                                                                                     ttGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCgg                                         <  1:80092/139‑1 (MQ=255)
                                                                                               ccATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCAcgtcgt                                                                                                         <  1:1441306/65‑1 (MQ=255)
                                                                                               ccATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCAcgtcgt                                                                                                         >  2:1441306/1‑65 (MQ=255)
                                                                                               ccATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAg                               >  2:169376/1‑139 (MQ=255)
                                                                                                     gACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTa                                                                           >  1:1488639/1‑89 (MQ=255)
                                                                                                     gACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTa                                                                           <  2:1488639/89‑1 (MQ=255)
                                                                                                       cGGAAGACGTTCAGTTCAAAGTTACCGGAAGCACCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCgg                       <  1:1013573/139‑1 (MQ=255)
                                                                                                          taGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTACCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTc                    >  2:1295395/2‑139 (MQ=255)
                                                                                                          aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCAt                                                                                                  >  2:211281/1‑61 (MQ=255)
                                                                                                          aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCAt                                                                                                  <  1:211281/61‑1 (MQ=255)
                                                                                                          aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTc                    >  1:1236981/1‑139 (MQ=255)
                                                                                                             aCGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGt                                                                             <  2:1375058/79‑1 (MQ=255)
                                                                                                             aCGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGt                                                                             >  1:1375058/1‑79 (MQ=255)
                                                                                                                    gttcaAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCa                                                                                 >  2:1250886/1‑68 (MQ=255)
                                                                                                                    gttcaAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCa                                                                                 <  1:1250886/68‑1 (MQ=255)
                                                                                                                    gttcaAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTTTTGGGTTCACTTTCCCCgg                                         <  2:331763/108‑1 (MQ=255)
                                                                                                                    gttcaAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCgg                                         >  1:331763/1‑108 (MQ=255)
                                                                                                                      tcaAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcag                                                                                        <  2:957931/59‑1 (MQ=255)
                                                                                                                      tcaAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcag                                                                                        >  1:957931/1‑59 (MQ=255)
                                                                                                                       caAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGAtt       <  1:102550/139‑1 (MQ=255)
                                                                                                                        aaaGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTc                                                                    >  2:549569/1‑77 (MQ=255)
                                                                                                                        aaaGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTc                                                                    <  1:549569/77‑1 (MQ=255)
                                                                                                                            ttACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAt  <  2:1236981/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                               
CTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAT  >  NC_000913/1685401‑1685664

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: