Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,688,142 G→A A79V (GCC→GTC)  fumA ← fumarate hydratase (fumarase A), aerobic Class I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,688,1420GA88.2% 39.7 / ‑1.5 17A79V (GCC→GTC) fumAfumarate hydratase (fumarase A), aerobic Class I
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (8/7);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCT  >  NC_000913/1688015‑1688273
                                                                                                                               |                                                                                                                                   
tACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCa                                                                                                                          >  2:515282/1‑139 (MQ=255)
    cGCTGCCCTTTTTTACCAACAATTATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGACCTCCGGGTCACGCa                                                                                                                      <  1:700302/139‑1 (MQ=255)
        gCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGACCTCCGGGTCACGCAGAAt                                                                                                                  >  1:183692/1‑139 (MQ=255)
            ttttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGACCTCCGGGTCACGCAGAATGTCg                                                                                                              <  2:183692/139‑1 (MQ=255)
              ttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGACCTCCGGGTCACGCAGAATGTCGGc                                                                                                            <  2:633276/139‑1 (MQ=255)
                                          aTCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACgcagc                                                                                <  1:515282/139‑1 (MQ=255)
                                                    gTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGACCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAacgacg                                                                      >  1:19874/1‑139 (MQ=255)
                                                                              gTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGACCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAAc                                            <  1:311536/139‑1 (MQ=255)
                                                                               tCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGACCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACa                                           >  2:644049/1‑139 (MQ=255)
                                                                                       aCGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGACCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGtt                                      >  2:130379/1‑136 (MQ=255)
                                                                                       aCGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGACCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGtt                                      <  1:130379/136‑1 (MQ=255)
                                                                                                 tGCAGCGCCACATATTTATCATTTTCGCTGACCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTg                         <  1:644049/139‑1 (MQ=255)
                                                                                                 tGCAGCGCCACATATTTATCATTTTCGCTGACCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTg                         <  2:19874/139‑1 (MQ=255)
                                                                                                                   tCATTTTCGCTGACCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCt       >  1:78757/1‑139 (MQ=255)
                                                                                                                     aTTTTCGCTGACCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGc     >  1:398164/1‑139 (MQ=255)
                                                                                                                       tttCGCTGACCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGccc   >  1:199846/1‑139 (MQ=255)
                                                                                                                        ttCGCTGACCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCt  >  1:456055/1‑139 (MQ=255)
                                                                                                                               |                                                                                                                                   
TACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCT  >  NC_000913/1688015‑1688273

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: