Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 606,499 G→A V295V (GTC→GTT)  ybdK ← putative glutamate‑‑cysteine ligase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913606,4990GA96.0% 78.8 / ‑4.0 25V295V (GTC→GTT) ybdKputative glutamate‑‑cysteine ligase 2
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (14/10);  total (14/11)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

GGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCC  >  NC_000913/606252‑606710
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                   
ggCAATCCCCGGCTTAGTCCCCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACg                                                                                                                                                                                                                  <  2:157105/251‑1 (MQ=255)
    aTTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATAACGCCtt                                                                                                                                                                                                              >  1:80854/1‑251 (MQ=255)
      tCCCGGCTTAGTCACCGGCCCAGATCTCCCAATGCTTTTTCACCAGCCCAATCAGCGATCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATAACGCCTTc                                                                                                                                                                                                             <  2:247276/250‑1 (MQ=255)
                                 cacaATGCTTTTTCACCAGCCCAATCAGCTAGCAGCCATCGGAGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTTTGCGGATCGGTGATAACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACgg                                                                                                                                                                                 <  1:46787/251‑1 (MQ=255)
                                 cacaATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATAACGCCTTCAAGCCCATAGCGACAGGCCTGGa                                                                                                                                                                                      <  1:292590/246‑1 (MQ=255)
                                          ttttCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATAACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAAtt                                                                                                                                                                          <  1:197274/249‑1 (MQ=255)
                                              ccccAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATAACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATAc                                                                                                                                                                     <  1:184083/248‑1 (MQ=255)
                                                              gAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATAACGCCTTCAAGCCCATAGCGACAGGcc                                                                                                                                                                                          <  1:136053/213‑1 (MQ=255)
                                                                gccgccATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATAACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCt                                                                                                                                                  >  1:251048/1‑251 (MQ=255)
                                                                        cGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATAACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGttt                                                                                                                                           >  1:7921/1‑250 (MQ=255)
                                                                                   tCGCGCATTAGCTGAGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATAACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTc                                                                                                                               >  1:85483/1‑251 (MQ=255)
                                                                                      cgcATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATAACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATAcag                                                                                                                                                                   <  1:293747/212‑1 (MQ=255)
                                                                                      cgcATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATAACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATAcag                                                                                                                                                                   >  2:293747/1‑212 (MQ=255)
                                                                                       gcATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATAACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTc                                                                                                                           <  2:80854/251‑1 (MQ=255)
                                                                                                       attcaGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATAACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTACTGATGTTTAAACGGGCGTTCCGTAAGTAACCAGTGGGCgg                                                                                                           >  1:284393/1‑251 (MQ=255)
                                                                                                                gCTGACGACCTGGCGATGAAGGGCCTCAATCGCGCTCGCTGCACCAATTTTATGTGCGGCAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTACAGTGTGCGGATCGGTGATAACGCCTTCAAGCCCCTAGCGCCAGGCCTGGAACCGGTTACATTTATACAGCAAGCAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGCAACCAGAGGGCTGAAGCCTGac                                                                                                  >  2:28663/1‑250 (MQ=255)
                                                                                                                 gcaacgacCTGGCGATGCAGGGCCCCAAGCGCGACCGATGCACCAATTTTAGGTGCGGAAGGGGCGATGTTTGCCAGCAATCGCCAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATAACGCCTTCAAGCCCATAGCGACTGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTGTTCCTTATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAAt                                                                                                 <  2:9609/248‑1 (MQ=255)
                                                                                                                        ccTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATAACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATccc                                                                                          >  1:129814/1‑251 (MQ=255)
                                                                                                                         cTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATAACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCg                                                                                         >  1:82457/1‑251 (MQ=255)
                                                                                                                                     agCCTAAATCGCGCTCGATGCACCACTTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATAACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGAGTTCCGTCAGTAACCAGTGGGCGGAAGCCTGAATTAATCCCGCCATATTTACTg                                                                             >  2:70152/2‑251 (MQ=255)
                                                                                                                                                           ggggtgccctcccgatcgggggCGATTTTTTCCAGCAATCGCAATGTATCTTCCGTTGGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATAACGCCTTCAAGCCCATAGCGACAGGCCGGGAAACGGTTAAATTTATACAGCAGGTAATCTTTTCCCTGATGTTTAAACGGGCGTTCCGTCAGTACCCAGTGGGCGGTAGCCTGAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGt                                                       <  1:217660/233‑1 (MQ=255)
                                                                                                                                                            cAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATAACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGt                                                       >  1:217368/1‑250 (MQ=255)
                                                                                                                                                                             ggggCGATTTTTTCCAGCAATCGCAAGGAATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATAACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTagatcggaagcgcgtcgg                                     >  2:217660/1‑234 (MQ=255)
                                                                                                                                                                                                       ggTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATAACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGc           >  1:284110/1‑251 (MQ=255)
                                                                                                                                                                                                                cGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATAACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATccc  >  1:146733/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                   
GGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCC  >  NC_000913/606252‑606710

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