Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 82,957 C→T V191V (GTG→GTA leuA ← 2‑isopropylmalate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091382,9570CT94.4% 46.7 / NA 18V191V (GTG→GTAleuA2‑isopropylmalate synthase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (8/9);  total (9/9)
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

AACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACCTCGTCCAGCGTGCTGCGC  >  NC_000913/82715‑83197
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                
aactGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGtt                                                                                                                                                                                                                                           <  1:36608/250‑1 (MQ=255)
                      ccTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCGGATGATTc                                                                                                                                                                                                                      <  1:62466/249‑1 (MQ=255)
                                           ggACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGgt                                                                                                                                                                                                <  1:195327/250‑1 (MQ=255)
                                                               cGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAAt                                                                                                                                                                           >  2:119428/1‑251 (MQ=255)
                                                                gAACTTTGATCCGCATACTGACTTCTTCAAGGGCACAGTATACGGCACGCTCGCCGATCCCCTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGCCCGCCAGGCCTCAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTCGGCACGCGTTCATACAGGCCCTTGATGATTTCGGCGAACTCAAACGGCATGGTGTAGCCAACGGTGTTCGGAATg                                                                                                                                                                          >  2:157233/1‑251 (MQ=255)
                                                                                                gAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGCTGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGGGGTGGCACCGGCATTAATcgccgc                                                                                                                                          >  1:182371/1‑251 (MQ=255)
                                                                                                                          cccGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCAGTAATCGCGGCTTCGACCACGCGCGCCAGATCGGCaa                                                                                                                <  1:86418/251‑1 (MQ=255)
                                                                                                                                           ccACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCTCCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGc                                                                                                <  1:156407/250‑1 (MQ=255)
                                                                                                                                           ccACCTGGCGTGCAACGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGc                                                                                                >  1:52687/1‑250 (MQ=255)
                                                                                                                                                        aCCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGGGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGaa                                                                                   >  2:52071/1‑250 (MQ=255)
                                                                                                                                                                                     cGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATGCCGGCGAACTCAAACGGCATGGTTTAGCCCGCTGTGTCCGGAATGTTGATGGTGGTGGCCCCGGCATTAATCGCCGCTTCGACCTCTCGCGCCAGATCGGCAATGGGTGTGGGCGCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGgcg                                                     <  1:174201/251‑1 (MQ=255)
                                                                                                                                                                                         cGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCAc                                                                                                                                                                                       <  2:115209/117‑1 (MQ=255)
                                                                                                                                                                                         cGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCAc                                                                                                                                                                                       >  1:115209/1‑117 (MQ=255)
                                                                                                                                                                                                cctaaaTCGTCGGGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCTGTGTAATTACGGGCGCGTTTCACCat                                          <  2:250100/248‑1 (MQ=255)
                                                                                                                                                                                                           gtgagTATGGACGGAAATAATGGCTTTGGCGATGTTCAGTTCGCGTTCATACAGGCCGCTGCTGATTCCGGCTAACTCAAGCTTCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCCCCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACTGGCGCGTGTCACCATATAGATTgcgc                               <  2:116666/247‑1 (MQ=255)
                                                                                                                                                                                                            tGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGcgcg                              >  1:79925/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCgg         >  2:81067/1‑250 (MQ=255)
                                                                                                                                                                                                                                        agATGTTAGGTACGCGTTCATACAGTCCGCTGATGATTCCGGCGAACTCACACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATCTAGATAGCGCGTTCGATCTCTTCATCCAGCGTGCTgcgc  >  2:139794/2‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                
AACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACCTCGTCCAGCGTGCTGCGC  >  NC_000913/82715‑83197

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