Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 79,040 A→G A138A (GCT→GCC leuD ← 3‑isopropylmalate dehydratase subunit LeuD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091379,0400AG96.4% 84.3 / ‑4.0 29A138A (GCT→GCCleuD3‑isopropylmalate dehydratase subunit LeuD
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (14/13);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-01

GATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAATGCCCAGGGCGCGTGCTCACGCGAAGAGCCACAGCCGAAGTTTTCTCGTGCCAGCAAAATGGAAGCGCCCTGATACTGCGGGAAGTTCAGCACGAAGT  >  NC_000913/78810‑79275
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           
gATAGCGGAAATCCCGGTCATGTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATTGCGTCGTCGTGCTGCAAGTTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAg                                                                                                                                                                                                                         <  1:433660/251‑1 (MQ=255)
       gAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGAAGGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCACATCCACGTCGAAATGGATCCCCGGATTGGCTTTCACCAGCCCAAACAGTTCGTcc                                                                                                                                                                                                                  >  2:462917/1‑251 (MQ=255)
            cccGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCAGCATGCAGTGGCGGCGGAAGGCATCGATGGTAAATCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTc                                                                                                                                                                                                             <  1:65162/251‑1 (MQ=255)
               ggTCATGTGAACGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCTGCGGAAGGCATCGATGGTAAAGCGATAGGTGTTCTCTCCCGCTTTCACCGCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTGGCTTACACCAGTGCAAACAGTTCGTCCACTTCTGc                                                                                                                                                                                                          <  1:36075/251‑1 (MQ=255)
                       gACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGGCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCtt                                                                                                                                                                                                   <  2:76069/250‑1 (MQ=255)
                        aCCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGCTACCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCCCTTCTGCATCGCTTa                                                                                                                                                                                                  >  2:11445/1‑250 (MQ=255)
                              aaatgggATTAAGTCAGAATCGCAGGTTGTTTCGCATCCTAAGCGGGAAAGGCGTCGTCGTGCTGCAAGGGAAGTCCAATACTGTCCAGACCGTTCACCATGCAGTTGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCCCCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCAc                                                                                                                                                                                           <  2:144088/247‑1 (MQ=255)
                                 ggggATTAGTTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGGCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGTTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGACGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCgg                                                                                                                                                                                        <  2:165238/251‑1 (MQ=255)
                                        aaTTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGAAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCagcagc                                                                                                                                                                                 >  2:464680/1‑251 (MQ=255)
                                        aaTTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAAGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTGGCTTTCACCAGCCCCAACAGTTCGTCCCCTTCTGCATCGCTTAATTTCACCGGCagcagc                                                                                                                                                                                 >  2:462047/1‑251 (MQ=255)
                                                                cTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGt                                                                                                                                                           >  1:80873/1‑249 (MQ=255)
                                                                 ttCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGTAATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGttg                                                                                                                                                         >  1:335621/1‑248 (MQ=255)
                                                                   cATAAGCGGCAAGGGCGTCGTCGTGCTGCAAGGTAAGCCCAATCCTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTGTCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTagaa                                                                                                                                                       <  1:307023/250‑1 (MQ=255)
                                                                          ggCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATTCTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAg                                                                                                                                               <  2:356424/251‑1 (MQ=255)
                                                                          ggCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAACTGGATCCCCGGATTGGCTTTCACCCACGCAAACAGTTCGGCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGGTGTTAAAGCTATAGCCGgagagaattac                                                                                                                                                 >  2:105699/1‑242 (MQ=255)
                                                                                   gtcgtcggGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCTGTTGCTGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTTGCTTTCACCAGCGCAAATAGTTCGTCCACTTCTTCAGCGCTTAATGTCACCGGCAGGAGCTGGTTGTTAAAGCTATTGCCGTAGAAGACGTCAGCAAAACTc                                                                                                                                       <  2:242053/250‑1 (MQ=255)
                                                                                         gtGCTGCAAGGTAAGCCCAATACTGTACAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCAACGTCGAAATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCa                                                                                                                                 >  1:117990/1‑250 (MQ=255)
                                                                                              gcAAGGTAAGCCCAATACTGTACAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGTTAAAGCGATAGGTTTTCTCTCCCGCTTTCTCCTCTTGCGCTTCCAGTTCCACGTCGAAATGGATTCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGGCCTCTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATca                                                                                                                             <  2:31588/249‑1 (MQ=255)
                                                                                                              aCTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGt                                                                                                           >  1:123111/1‑251 (MQ=255)
                                                                                                                                           ggcggAAGGCATCGCTGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTGGCTTTCACCAGCGCAAACCGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTTAAGCTATTGCCGTAGAAGATGTCCGCAAAACTCGGCGCAATCACCACTTTAAAACCATAGTCGGTCACTGCCCAGGGCGCGTGCTCAcgc                                                                              >  2:405048/1‑251 (MQ=255)
                                                                                                                                                                       ggTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATcac                                                                                                                            <  2:274141/177‑1 (MQ=255)
                                                                                                                                                                       ggTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATcac                                                                                                                            >  1:274141/1‑177 (MQ=255)
                                                                                                                                                                                 cccGCTTTCACCTCTGGCGTTTCCAGATCCACGTCGAAATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTGCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACGCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAATGCCCAGGGGGCGTCCTCACGCGAAGAGCCACAGCCGAAGTTTTCTCGTGCCAGCaa                                           <  2:184730/248‑1 (MQ=255)
                                                                                                                                                                                       ttCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAATGCCCAGGGCGCGTGCTCACGCGAAGAGCCACAGCCGAAGTTTTCTCGTGCCAGCAAAATGGAAgc                                  >  2:113367/1‑251 (MQ=255)
                                                                                                                                                                                                tGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTGGCTTTCACCAGCGCTAACAGTTCTTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAATGCCCAGGGCGCGTGCTCACGCGAAGAGCCACAGCCGAAGTTTTCTCGTGCCAGCAAAATGGAAGCGCCCTg                            <  1:212091/248‑1 (MQ=255)
                                                                                                                                                                                                 gcgcTTCCAGATCCACGTCGAAATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAATGCCCAGGGCGCGTGCTCACGCGAAGAGCCACAGCCGAAGTTTTCTATGGCCAGCAAAATAGAAGCGCCCTGTTAc                        >  1:195698/1‑251 (MQ=255)
                                                                                                                                                                                                                cGTCGAAATGGATCCACGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAATGCCCAGGGCGCGTGCTCACGCGAAGAGCCACAGCCGAAGTTTTCTCGTGCCAGCAAAATGGAAACGCCCTGATCATGCGGGAAGTTCAGc         >  2:267033/1‑251 (MQ=255)
                                                                                                                                                                                                                   cGAAATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAATGCCCAGGGCGCGTGCTCACGCGAAGAGCCACAGCCGAAGTTTTCTCGTGCCAGCAAAATGGAAGCGCCCTGATACTGCGGGAAGTTCAGCAc       <  1:26875/250‑1 (MQ=255)
                                                                                                                                                                                                                       aTGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAATGCCCAGGGCGCGTGCTCACGCGAAGAGCCACAGCCGAAGTTTTCTCGTGCCAGCAAAATGGAAGCGCCCTGATACTGCGGGAAGTTCAGCACGAAGt  <  1:464680/251‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           
GATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAATGCCCAGGGCGCGTGCTCACGCGAAGAGCCACAGCCGAAGTTTTCTCGTGCCAGCAAAATGGAAGCGCCCTGATACTGCGGGAAGTTCAGCACGAAGT  >  NC_000913/78810‑79275

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