Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,685,542 (C)6→7 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,5361.C96.3% 91.3 / ‑3.2 27coding (1053/1404 nt)fumCfumarate hydratase (fumarase C),aerobic Class II
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (13/13);  total (13/14)
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

TCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAT  >  NC_000913/1685398‑1685664
                                                                                                                                           |                                                                                                                                
tCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAgcg                                                                                                                                   <  2:1143374/139‑1 (MQ=255)
   cTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccc                                                                                                                               >  1:690952/1‑139 (MQ=255)
     cACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccccc                                                                                                                             >  2:1116858/1‑139 (MQ=255)
        gATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccccccca                                                                                                                          <  2:359250/139‑3 (MQ=255)
            cGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGtt                                                                                                                       <  2:690952/139‑1 (MQ=255)
              gTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGa                                                                                                                     <  1:248855/139‑1 (MQ=255)
                       ccACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCAcgt                                                                                                            <  1:1116858/139‑1 (MQ=255)
                            gcgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTc                                                                                                       >  2:656157/1‑139 (MQ=255)
                            gcgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTc                                                                                                       >  1:585568/1‑139 (MQ=255)
                            gcgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGcccc                                                                                                                              <  2:409825/115‑1 (MQ=255)
                            gcgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGcccc                                                                                                                              >  1:409825/1‑115 (MQ=255)
                                 gtgtTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAgcgc                                                                                                                                 >  2:963042/1‑107 (MQ=255)
                                 gtgtTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAgcgc                                                                                                                                 <  1:963042/107‑1 (MQ=255)
                                 gtgtTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCAt                                                                                                  <  1:656157/139‑1 (MQ=255)
                                 gtgtTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCAt                                                                                                  <  2:714984/139‑1 (MQ=255)
                                           aaaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagc                                                                                       <  1:895486/139‑1 (MQ=255)
                                            aaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACTGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagc                                                                                       >  2:754364/1‑139 (MQ=255)
                                            aaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTccc                                                                                                     <  1:687468/125‑1 (MQ=255)
                                            aaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTccc                                                                                                     >  2:687468/1‑125 (MQ=255)
                                            aaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccccc                                                                                                                             >  2:290975/1‑100 (MQ=255)
                                            aaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccccc                                                                                                                             <  1:290975/100‑1 (MQ=255)
                                                ttCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCagag                                                                                   >  1:661382/1‑139 (MQ=255)
                                                                  cAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                 <  2:1022554/139‑1 (MQ=255)
                                                                  cAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                 <  2:805367/139‑1 (MQ=255)
                                                                                           tgGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGc                                                                                  <  2:730803/97‑1 (MQ=255)
                                                                                           tgGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGc                                                                                  >  1:730803/1‑97 (MQ=255)
                                                                                                     ttGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCagag                                                                                   >  1:1037066/1‑86 (MQ=255)
                                                                                                     ttGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCagag                                                                                   <  2:1037066/86‑1 (MQ=255)
                                                                                                          cGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcag                                                                                        >  1:296441/1‑76 (MQ=255)
                                                                                                          cGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcag                                                                                        <  2:296441/76‑1 (MQ=255)
                                                                                                             aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCAt                                                                                                  >  1:1117478/1‑63 (MQ=255)
                                                                                                             aaGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCAt                                                                                                  <  2:1117478/63‑1 (MQ=255)
                                                                                                                  gttcagttcaAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCAtt                 >  2:391837/1‑139 (MQ=255)
                                                                                                                   ttcagttcaAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCAttt                >  1:23577/1‑139 (MQ=255)
                                                                                                                        ttcaAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCgg           >  1:830214/1‑139 (MQ=255)
                                                                                                                          caAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCCTCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGa         >  2:1056856/1‑139 (MQ=255)
                                                                                                                                 aCCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAt  <  2:830214/139‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                
TCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAT  >  NC_000913/1685398‑1685664

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: