Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 5,442 T→C L70P (CTA→CCA)  yaaX → DUF2502 domain‑containing protein YaaX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009135,4420TC81.8% 40.4 / 0.1 22L70P (CTA→CCA) yaaXDUF2502 domain‑containing protein YaaX
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base C (12/6);  total (13/9)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

AACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTG  >  NC_000913/5207‑5677
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           
aaCGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCgccgcc                                                                                                                                                                                                                              >  1:78663/1‑251 (MQ=255)
                   aGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATa                                                                                                                                                                                                             >  1:334187/1‑249 (MQ=255)
                             gAAAAAGATGCAAGCTATCGTACTCGCACGTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGctcct                                                                                                                                                                                                  <  1:32161/250‑1 (MQ=255)
                                   gATGCAATCTATCGCACCCGCACTTTCACTGGTTCTGGTCGCCCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGGACCTACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCATCATGatc                                                                                                                                                                                           <  1:155065/251‑1 (MQ=255)
                                                    tCGCACTTTCCCTGGTTCTGGTAGCTCCCATGGCAGCACAGGCTGCGGAAATTCCGTTAGTCCCGTCAGTAAAATTACAGATCGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAAGTCACTGGCGCGACCACGGCTAGTGGAAACAACATTATGAATGGCGAGGCCAGCGCTGGCACCCACAAGGCACGCCGCCACCGCCCGGCCACCATAAGACAGATCCTCATGCTCATCAcggcgcaaacggc                                                                                                                                                                          >  2:48656/1‑243 (MQ=255)
                                                     ccgcccttCCCTGTGTTCTGGGCGCTCCCAGGGCAGCACAGGCTTCGGAAATTACGTTAGTCCCTTCAGTAAAATTACAGAGAGGCGAGCGTGATAATCGTGGCTATTACGGGGATGGAGGTCACTGGCGCGACCACGGCTGGGGAAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATgg                                                                                                                                                                           <  2:334187/244‑1 (MQ=255)
                                                      gcccttccccTGGTTCGGGTCGCTCCCATGGCATCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGAGAGGGGATCGTGATAATCGGGGCTATTACGGGGAGGGAGGTCACTGGCGCGACCCCGGCGGGGGAAACCAACATTAGGAAGGGCGAGGCAAACGCTGGCCCCTACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTcc                                                                                                                                                                        <  2:73576/244‑1 (MQ=255)
                                                      gCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCCTgggcc                                                                                                                                                                        >  1:73576/1‑248 (MQ=255)
                                                              ccTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGTCGGGCATAGTCCAGGCCAAc                                                                                                                                                                >  1:3504/1‑251 (MQ=255)
                                                                        ggtcgccccaaaggCAGCACGGGCTGGGTAAATTACGTTAGTCCCTTCAGTAAAATTACAGATAGGGTATCGTGATAAACGTGGCGATTACTGGGATGGAGGTCACTGGCGCGACCCCGGCTGGTGGAAACAACATTATGAATGGCGAGTCAATCGCTGGCCCCCACAGGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCt                                                                                                                                                       <  2:78663/239‑1 (MQ=255)
                                                                                          acaGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCCAGGGCCAGGCAAAAATAACCGCCAAAAGAAAAAAGCCgggg                                                                                                                                     >  1:44987/1‑250 (MQ=255)
                                                                                              gcttgggAAATTTCGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGTGGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACa                                                                                                                                 <  1:313329/246‑1 (MQ=255)
                                                                                                          aCGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCGGGTGGAAACACCAGTATGAATGGCGAGGCAATCGCGGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCa                                                                                                                     <  1:47793/251‑1 (MQ=255)
                                                                                                                gTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGGCCAGGCCAACATCACCCCTAAAATACAAAAGCCGGGGAACCATCCGGCATTCAGCGccc                                                                                                               >  1:178402/1‑250 (MQ=255)
                                                                                                                     gTCAGTAAAATTACAGAGAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGc                                                                                                          <  1:142080/251‑1 (MQ=255)
                                                                                                                                 aCAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCCCCAAATGACAAATGCCGGGGAACAATCCGGCATTCAACGCCCGATGCGACGCTGGcgcg                                                                                              >  1:188050/1‑251 (MQ=255)
                                                                                                                                         gCGATCGTGAGACTCGTGGCTATTCCGGGGATGGAGGTCACGGGCGCGCCCCGGGCGGGTGGAAACAACATTAGGAAGGGCGAGCAAATCGTTGGCACCCACCCGGCCCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCa                                                                                      <  2:178402/251‑1 (MQ=255)
                                                                                                                                                                           aGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCCATATATTTAATCTGc                                                    >  1:347576/1‑251 (MQ=255)
                                                                                                                                                                                                     ggAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGGCATGGGCCAGGGAAACATCACCGCTAAATAACAAAAGCCGGGGAAAAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTGATTCTGCAATATATTTAAACTGCATGCTTTTGTAGGCCGGGTAAGGCgg                          >  1:213523/1‑250 (MQ=255)
                                                                                                                                                                                                          caacaTTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCCAACATCACCGCTAAATGACAAAAGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAGTGAATTCTGCAATATATTGAATTTGCAAGATTTTGTAgggcgggtaagggggtcacg                     >  1:377747/1‑233 (MQ=255)
                                                                                                                                                                                                                 tATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACCAATGCCGGGTAACAATCCGGCATTCAGCGCCCGATGCGACGCTGGGGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATTCTTTTGTAGGccggcaaaggcggttacgccccat               >  1:1491/1‑228 (MQ=255)
                                                                                                                                                                                                                              gCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACACCACCGCTAAAAGACAAATGCCGGGTAACAATCCGGCATTCAACGCCCGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGGATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTg  >  1:254156/1‑251 (MQ=39)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           
AACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTG  >  NC_000913/5207‑5677

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