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)
       ttCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGtt                                                                                                                       <  2:24173/139‑1 (MQ=255)
                aTACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCAc                                                                                                              <  2:1520007/139‑1 (MQ=255)
                     cACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCAcgtcgt                                                                                                         <  1:478816/139‑1 (MQ=255)
                         gcgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTccc                                                                                                     >  2:102550/1‑139 (MQ=255)
                                       aaaaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagc                                                                                       >  2:1302306/1‑139 (MQ=255)
                                             ttCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCa                                                                                 <  2:230598/139‑1 (MQ=255)
                                                   gCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTa                                                                           <  2:1191855/139‑1 (MQ=255)
                                                      aTCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCAcgtcg                                                                                                          >  2:206509/1‑105 (MQ=255)
                                                      aTCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCAcgtcg                                                                                                          <  1:206509/105‑1 (MQ=255)
                                                            cAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcac                                                                  >  1:372537/1‑139 (MQ=255)
                                                             aGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                 <  2:1231687/139‑1 (MQ=255)
                                                               cAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGC‑ACCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACt                                                               >  2:1013573/1‑139 (MQ=255)
                                                                     cACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTg                                                         >  2:1028398/1‑139 (MQ=255)
                                                                        cGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGt                                                      <  1:1028398/139‑1 (MQ=255)
                                                                         gATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCAcgtcgt                                                                                                         >  1:1112544/1‑87 (MQ=255)
                                                                         gATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCAcgtcgt                                                                                                         <  2:1112544/87‑1 (MQ=255)
                                                                         gATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGtt                                                     >  1:1112779/1‑139 (MQ=255)
                                                                            tGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccc                                                                                                                               >  2:952113/1‑63 (MQ=255)
                                                                            tGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccc                                                                                                                               <  1:952113/63‑1 (MQ=255)
                                                                            tGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTccc                                                                                                     >  2:1286847/1‑88 (MQ=255)
                                                                            tGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTccc                                                                                                     <  1:1286847/88‑1 (MQ=255)
                                                                             gCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACt                                                 <  2:1566560/139‑1 (MQ=255)
                                                                                ggAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTc                                              <  2:1365767/139‑1 (MQ=255)
                                                                                  aaaTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGAt                                                                                                                    >  2:213188/1‑67 (MQ=255)
                                                                                  aaaTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGAt                                                                                                                    <  1:213188/67‑1 (MQ=255)
                                                                                     ttGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCgg                                         <  1:80092/139‑1 (MQ=255)
                                                                                               ccATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCAcgtcgt                                                                                                         <  1:1441306/65‑1 (MQ=255)
                                                                                               ccATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCAcgtcgt                                                                                                         >  2:1441306/1‑65 (MQ=255)
                                                                                               ccATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAg                               >  2:169376/1‑139 (MQ=255)
                                                                                                     gACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTa                                                                           >  1:1488639/1‑89 (MQ=255)
                                                                                                     gACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTa                                                                           <  2:1488639/89‑1 (MQ=255)
                                                                                                       cGGAAGACGTTCAGTTCAAAGTTACCGGAAGC‑ACCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCgg                       <  1:1013573/139‑1 (MQ=255)
                                                                                                          taGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTACCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTc                    >  2:1295395/2‑139 (MQ=255)
                                                                                                          aaGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCAt                                                                                                  >  2:211281/1‑61 (MQ=255)
                                                                                                          aaGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCAt                                                                                                  <  1:211281/61‑1 (MQ=255)
                                                                                                          aaGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTc                    >  1:1236981/1‑139 (MQ=255)
                                                                                                             aCGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGt                                                                             <  2:1375058/79‑1 (MQ=255)
                                                                                                             aCGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGt                                                                             >  1:1375058/1‑79 (MQ=255)
                                                                                                                    gttcaAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCa                                                                                 >  2:1250886/1‑68 (MQ=255)
                                                                                                                    gttcaAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCa                                                                                 <  1:1250886/68‑1 (MQ=255)
                                                                                                                    gttcaAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTTTTGGGTTCACTTTCCCCgg                                         <  2:331763/108‑1 (MQ=255)
                                                                                                                    gttcaAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCgg                                         >  1:331763/1‑108 (MQ=255)
                                                                                                                      tcaAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcag                                                                                        <  2:957931/59‑1 (MQ=255)
                                                                                                                      tcaAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcag                                                                                        >  1:957931/1‑59 (MQ=255)
                                                                                                                       caAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGAtt       <  1:102550/139‑1 (MQ=255)
                                                                                                                        aaaGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTc                                                                    >  2:549569/1‑77 (MQ=255)
                                                                                                                        aaaGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTc                                                                    <  1:549569/77‑1 (MQ=255)
                                                                                                                            ttACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAt  <  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: ‑