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,4420TC96.3% 71.0 / ‑6.3 27L70P (CTA→CCA) yaaXDUF2502 domain‑containing protein YaaX
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (8/18);  minor base A (0/1);  total (8/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGC  >  NC_000913/5218‑5678
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               
aTAAAAGGAGTAACCAATAAAAAATATGCAACCTCTCGTCCCCGCATTTTCCCTGGTTCTGTTCGCTCCAAGGGCAGCCCAGGCTGCGTAATTTACGTTGTCCCCGTCAGTAAAATAACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACGGGCGCGCCCCCGGCTGGTGGAAACAACATTATGAAGGGCGAGGCACTCGCGGCCCCCCCCCCGGCCCGCCGCCACCGCCGCGcc                                                                                                                                                                                                                      <  2:146654/251‑1 (MQ=255)
     aGGAGTAACCAATTAAAAAGATGCAATCTATCGGACTCGCACTTTCCCTGGTCCTGGGCGCTCCCATGGCAGCACAGGCTGCGGAAATGACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGGTGGTGGAAACAACATTATGAAGGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCAt                                                                                                                                                                                                                 <  1:426577/251‑1 (MQ=255)
           aaCCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGaaa                                                                                                                                                                                                           <  1:120760/251‑1 (MQ=255)
                                  agcGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCATAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGGGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCAcg                                                                                                                                                                                        <  1:86102/245‑1 (MQ=255)
                                          cGCACTTTCCCTGGTTCTGGTCGTACCCATGGCAGCACAGGCGGCGGAAATTACGTTAGTCCCGGCAGTAAAATTACAGATAGGCGATCGTGATAAGCGTGGTTATTACGGGGATGGAGGTCCCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAAGGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTc                                                                                                                                                                            <  1:352354/251‑1 (MQ=255)
                                                   ccTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAAc                                                                                                                                                                   <  1:569318/251‑1 (MQ=255)
                                                                   cccATGGCAGCACAGGCGTGGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGACAGGCTATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGAc                                                                                                                                                  <  1:158798/251‑1 (MQ=255)
                                                                                 ggCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACAACACCGCTAAATAACAAAAGCCGGGTaa                                                                                                                                    >  1:435364/1‑250 (MQ=255)
                                                                                            aTTACGTAAGGCCCGCCAGGAAAATTACAGATAGGGGATCGGGATAATGGTGGCAATCCCTGGGTGGGAGGTCACTGGCGCGGCCTGGGCTGGTGGAAACAACATTATGAATGGCGGGGAAATCGCTGGCACCCACACGGCCCGCCGCCGCCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCAt                                                                                                                         <  2:86652/251‑1 (MQ=255)
                                                                                                 gTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTTCAGGCAAACATCACCGCTAAATGGCAAATGCCGGGTAACAATCCGGCATTCAgc                                                                                                                    >  1:548384/1‑251 (MQ=255)
                                                                                                 gTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAACCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAgc                                                                                                                    <  1:201722/251‑1 (MQ=255)
                                                                                                   tAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAgcg                                                                                                                   <  1:670829/250‑1 (MQ=255)
                                                                                                   tAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGAGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGGCATGGTCCCGGCAAACAAACACCCCAAATGACCAATGCCGGGTGAACAACCCgcaattacccc                                                                                                                   >  2:313034/1‑241 (MQ=255)
                                                                                                      tCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGGCATTGTCCAGGCAAACATCACCGCTAAATGACAAAAGCCGGGTAACAATTCGGCATTCCACGCCTg                                                                                                               >  1:173810/1‑251 (MQ=255)
                                                                                                       cccagTCAGTAAAAGTACATAGAGGCGATCGTGATAAGTGTGGATATTAAGGGGATGGAGGTCCGTGGCGCGACCACGGCTGTGGAAAACAACATTATGAAGGGCAAGGTAAGCCTTGCCCCCACCAGGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCGCAGGATCATCACGGCGGTCATGGTCCAGGCAAACATCGCCGCTAAAGGACAAAGGCCGGGTAAAAATCCGTCATTCAGCGCCTg                                                                                                               <  2:319119/247‑1 (MQ=255)
                                                                                                                        cAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCCGGCAAACATCACCGCTaa                                                                                                                                                       >  1:206279/1‑195 (MQ=255)
                                                                                                                        cAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGAGGGAGGTCCCGGGCGCGCCCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTaa                                                                                                                                                       <  2:206279/195‑1 (MQ=255)
                                                                                                                           aaaggcttgcGTGAAAATCGTGGGTATTGCTGGGAGGGAGGTCACTGGCGCGCCCACGGTGGGTGTAACCAACATTAGGATTGGGGGGGCAATCTCGGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACAGCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCtt                                                                                           <  2:660139/242‑1 (MQ=255)
                                                                                                                            aaGGCGATCGTGATAACCTTGGCTATTACTGGGATGGAGGTCACTGGCGCGCCCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGCCCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTa                                                                                          <  1:633388/250‑1 (MQ=255)
                                                                                                                                              ggcgctataagggggATGGGTGTCACGGGCGCGACCACGGCTGGGGAAAACACCTTTATAAAGGGCGAGGAAATCGTTGGCCCCCCCACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAAGGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAAtgc                                                                        <  2:453377/239‑3 (MQ=255)
                                                                                                                                                ggCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGACATGGGCCAGGCAAACATCACCGCCAAATTACAAAAGCCGGGGAACAATCCGGCCTTCCGCGCCTGATGCGACGCTGGCGGGTCTGATCAgggccacgcgaattatc                                                                      >  1:415822/1‑236 (MQ=255)
                                                                                                                                                            gaaggAGGTCCCTGGGGCGCCCACGGCGGGGGGAAAAAACATTAGGAATGGCGATGCAATCGCAGGCCACCAACGGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCGTATCAGGCCTACGTGAATTCTGCTATATATTGaa                                                          <  2:105024/248‑1 (MQ=255)
                                                                                                                                                               ggAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAAGGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                         <  1:607211/217‑1 (MQ=255)
                                                                                                                                                                            gcgACCACGGCTGGTGGAAACAACATGTGGAGGGGCGAGGCAATCGCTGGCCCCCACACGGACCGCGGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTa                                          <  2:80783/251‑1 (MQ=255)
                                                                                                                                                                                                                      aaTCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATACGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTTCATGCTTTTTTAGGCCGGATAAGGCGTTCACGCCGCCTCCGGCATTGACTGc  >  1:442143/1‑250 (MQ=255)
                                                                                                                                                                                                                          gCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCa                                                                                       <  1:191657/161‑1 (MQ=255)
                                                                                                                                                                                                                            tGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGGAAc                                                                                                                                   >  1:62399/1‑115 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               
ATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGC  >  NC_000913/5218‑5678

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