Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,767,014 C→T G558S (GGT→AGT)  ydiJ ← putative FAD‑linked oxidoreductase YdiJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,767,0140CT87.5% 50.7 / 1.9 24G558S (GGT→AGT) ydiJputative FAD‑linked oxidoreductase YdiJ
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (8/13);  total (10/14)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-02

TGGATCAACGCCGCGATCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACCGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGGTCAAATGCCGCTTTCACTTTGCGCAGTTCTGCAAAAAGTTCCTCACCGAAAAACGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGTCATCAGAGATTTGCTTCATCAAA  >  NC_000913/1766802‑1767258
                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
tGGATCAACGCCGCGATCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCAATTTCATCATCGGCGCATCg                                                                                                                                                                                                                   >  1:197983/1‑250 (MQ=255)
  gATCAACGCCGCGATCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCgaga                                                                                                                                                                                                                <  1:270756/251‑1 (MQ=255)
       aCGCCGCGATCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTc                                                                                                                                                                                                           <  1:105877/251‑1 (MQ=255)
               aTCCGCCAACAACCGAAGCCATCCACGCACCAGGTTTGCGCGCCCTTTCGGGGAATGAAATCGGTGGTGGGTGATCTGCATCGACGGACACATCGGACTCCGGGCAGCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACGGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCTCCGCGTCCACTTTCATCATCGGCGTATCGAGACCTTCTGGCgg                                                                                                                                                                                                     <  2:131313/249‑1 (MQ=255)
                    ccAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGATCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGGACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGTCCTTCTGGCGGGCATATc                                                                                                                                                                                              >  1:105159/1‑251 (MQ=255)
                       acaaACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTc                                                                                                                                                                                           <  1:65945/251‑1 (MQ=255)
                                gCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGAGCGAATGTACTGCGCTTCACCGCGTcc                                                                                                                                                                                                                                        <  1:234955/197‑1 (MQ=255)
                                                     cgctcccTTGCGGTGAATGAATCCGGTTCGGGGGTATCTTCATCGACGGACACATCGGACTCCGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTCCCGCAATGGGTATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGGTCAAAt                                                                                                                                                            <  2:13715/247‑1 (MQ=255)
                                                                    aatgaatCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCCTCATCGGCGCATCGAGAACTTCTGGCGGGCAAATCTTCCCTGGGTTTAGTCGGTTAAGCGGGTCAAATGCCGCATTCCCtttt                                                                                                                                             >  1:38480/1‑250 (MQ=255)
                                                                     atgactCCGGTTAGGGGTAATCTTCATCGACGGACGCATCGGTCTACGGGATTCAAAGTTGAAGCATAAACCGTGGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGTGATCTGCCGACTGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGTGCAAATCTCCCCTGGGTGGAGTCGGTTATGCGGTTCAAATGCCGATTTCACTTTgc                                                                                                                                            <  2:140783/251‑1 (MQ=255)
                                                                          tCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATATTCCCTGGGGTGAGTCGGTTATGCGGGTCAAATGCCGCTTTCACTTTGCGCAGt                                                                                                                                       >  1:187571/1‑251 (MQ=255)
                                                                            cGGTTCTGGGTGATCTTCATCGACGGACCCATCGGACTACGGGCATCAAAGTTGAAGCATAACCCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTgg                                                                                                                                                                                      <  2:33325/203‑1 (MQ=255)
                                                                            cGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTgg                                                                                                                                                                                      >  1:33325/1‑203 (MQ=255)
                                                                               ttctggggtagtctcctCGACGGACACCTCGGACTCCGGGCATCAACGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCGGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCCCCGCGTCCTCTTCCATCAGCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGGTCAAATGCCGCTTTCACTTTGCGCAGTTCTg                                                                                                                                   <  2:57825/234‑1 (MQ=255)
                                                                                                            cgtaCTACGGTCACCATAGTCGAAGCCTAAACCGTGGCCGTTACACTCCACCGCAACGCGCCACTGCCGGCTTACCGCAAGGGGGATCTGCCGAGCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGTTCAAATGCCGCTTTCGCTTTGCGCAGTTCTGCAAAAAGTTCCTCACCGAAAAACGCCggc                                                                                                     <  2:82819/247‑2 (MQ=255)
                                                                                                               aCTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATATTCCCTGGGTTGGGTCGGTTATGCGGGTACAAAGCCGCTTTCACTTTGGGCGATTCTGACAAAAATTCCTCAcacaaaaccgacagtgtgt                                                                                                 >  2:216363/1‑233 (MQ=255)
                                                                                                                          cAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGGTCAAATGCCGCTTTCACTTTGCGCAGTTCTGCAAAAAGTTCCTCACCGAAAAACGCCGGGCTGTATTCAGcgcg                                                                                       <  1:165100/250‑1 (MQ=255)
                                                                                                                            aaGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACCGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTc                                                                                                                                                                                                           >  1:69940/1‑134 (MQ=255)
                                                                                                                            aaGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACCGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTc                                                                                                                                                                                                           <  2:69940/134‑1 (MQ=255)
                                                                                                                                         aaCCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGGTGAGTCGGTTATGCGGGTCAAATTCCGCTTTCCCTTTGCGCAGTTCTGCAAAAAGTTCCTCCCAGAAAAACGCCGTGCTTTCTTCAGCGCGAAAACCTTTGCCGTGc                                                                       >  1:2096/1‑251 (MQ=255)
                                                                                                                                                                             gctgGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGGTCAAATGCCGCTTTCACTTTGCGCAGTTCTGCAAAAAGTTCCTCACCGAAAAACGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAgc                                   <  1:108423/251‑1 (MQ=255)
                                                                                                                                                                                          aaTGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGGTCAAATGCCGCTTTCACTTTGCGCAGTTCTGCAAAAAGTTCCTCACCGAAAAACGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTACTTCGCAGTCAGCGCCACCACGtcat                      <  1:216363/251‑1 (MQ=255)
                                                                                                                                                                                               ggATCTGCCGATCGAATGTACCGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGGTCAAATGCCGCTTTCACTTTGCGCAGTTCTGCAAAAAGTTCCTCACCGAAAAACGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGTCATCag                   >  1:64134/1‑249 (MQ=255)
                                                                                                                                                                                                               tGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGGTCAAATGCCGCTTTCACTTTGCGCAGTTCTGCAAAAAGTTCCTCACCGAAAAACGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGTCATCAGAGATTTGCTTCATCaaa  >  1:203218/1‑250 (MQ=255)
                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
TGGATCAACGCCGCGATCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACCGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGGTCAAATGCCGCTTTCACTTTGCGCAGTTCTGCAAAAAGTTCCTCACCGAAAAACGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGTCATCAGAGATTTGCTTCATCAAA  >  NC_000913/1766802‑1767258

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