Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,426,760 G→A I221I (ATC→ATT yjfZ ← DUF2686 domain‑containing protein YjfZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,426,7600GA90.9% 60.9 / ‑3.1 22I221I (ATC→ATTyjfZDUF2686 domain‑containing protein YjfZ
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (13/7);  total (13/9)
Fisher's exact test for biased strand distribution p-value = 1.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGCAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGT  >  NC_000913/4426514‑4426993
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         
gCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGGCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTg                                                                                                                                                                                                                                       <  2:141252/251‑1 (MQ=255)
             gcaaGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAAATTGCGGCTGAACTTg                                                                                                                                                                                                                           >  1:4569/3‑250 (MQ=255)
                                      gAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAAATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCaa                                                                                                                                                                                                  >  1:29241/1‑250 (MQ=255)
                                              cGTTATTCAACGAGCCGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAAATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATgggg                                                                                                                                                                                         >  1:281111/1‑251 (MQ=255)
                                                     cAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAAATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATc                                                                                                                                                                                   >  1:284314/1‑250 (MQ=255)
                                                                    aggCCGAACATTGTGCGGCGTTTTGCTTTCTCGTTTTCGTCATGTATGAGTTTCGCAAGTAGATCCAGTCGCCATTCTCTTGTGTCCATAAGTCAGTATGGTTTGCTTGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTTCGGCACGCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTGCAAATCATGGGGTAAATCAGAATACGCGTTTAGc                                                                                                                                                                   <  2:29129/249‑1 (MQ=255)
                                                                                ggTCTGCGTTTTGCTTTCTATTTTCCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGGCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAAATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAatg                                                                                                                                                       <  2:284314/251‑1 (MQ=255)
                                                                                 gTCGGCGTTTTTCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAAATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCa                                                                                                                                                                                              >  2:160119/1‑211 (MQ=255)
                                                                                 gTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAAATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCa                                                                                                                                                                                              <  1:160119/211‑1 (MQ=255)
                                                                                      cGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAAATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTc                                                                                                                                                                                                    <  1:166755/200‑1 (MQ=255)
                                                                                         tttGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAAATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAATTCATGGTGAAAATCAGAATACGCGTTTAGCAATGCTGCAATGATGAATCCg                                                                                                                                              >  2:294195/1‑251 (MQ=255)
                                                                                                             aGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAAATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGt                                                                                                                          <  1:64861/251‑1 (MQ=255)
                                                                                                                         cGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAAATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAg                                                                                                              <  1:158891/251‑1 (MQ=255)
                                                                                                                                    ccattcgccattcgCTTGTGGCCATAAGTCAGGATGGTTTGCTTGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAAATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCTAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGatat                                                                                                   <  2:304972/251‑1 (MQ=255)
                                                                                                                                           ccATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAAATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATaa                                                                                                 >  1:85229/1‑246 (MQ=255)
                                                                                                                                            cATTCTCTTGTGTACATACGGCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAAATTGCGGCTGAAATTGATAAGGCACGCAGACCAATTTTCAAATCATGGGGTAAATCAGCATACGCGTTTAGCAATGCAGCAATGATGAATCCGGGGGCATTAGAAACACCTGTAACAGAGCGAAGTTGGTCGaactaaggatcg                                                                                           >  2:225908/1‑240 (MQ=255)
                                                                                                                                                          cATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAAATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGGGCGAAGTTGGTCGATATAAAGATCGATTGTTggacgttc                                                                             >  1:67679/1‑245 (MQ=255)
                                                                                                                                                                          ttGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAAATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAAAGCCTGTAACAGTGCGAAGTTGGTCGATATCAAGTACGATTGTTGGTAGTTCAGATGACCAGCAAtt                                                              >  2:106751/1‑250 (MQ=255)
                                                                                                                                                                             cTGGCGTTAGATACATACGAGCGCAACTGCTGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAAATTGCGGCTGACCTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGCTGAATCCGAGGGAATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGAACGATTGTTGGAAGTTCAGATGCCCAGCAATTAGGc                                                          >  1:142527/1‑251 (MQ=255)
                                                                                                                                                                                                                     cTGCGGCACTAGCTCGAAACCCAAACCCGAGTAAATTGCGGCTGAACTTGATACGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGCATACGCGATTAGCAATACAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCAGTTGTTGGAAGTTCAGATGACCAGCAATTAGGCATTTGGGCAGGCAGGGTTACAAAGTTACCAATCGTTTGGc                  >  2:157545/1‑251 (MQ=255)
                                                                                                                                                                                                                         ggCACTCGCTCGAAACCCAAACCCGAGTAAATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGCAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGt              >  1:228521/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aaacccaaacccGAGTAAATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGCAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAgg  <  2:85229/251‑2 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         
GCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGCAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGT  >  NC_000913/4426514‑4426993

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