Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,850,605 C→T T590I (ACC→ATC)  sppA → protease IV, a signal peptide peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,850,6050CT100.0% 39.2 / NA 15T590I (ACC→ATC) sppAprotease IV, a signal peptide peptidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (7/8);  total (7/8)

GGGCCACGTCTGGACCGGTCAGGATGCAAAAGCTAACGGGCTGGTCGATAGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCG  >  NC_000913/1850360‑1850839
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          
gggCCACGTCTGGACGGGTCAGGATGCAAAAGCTAACGGGCTGGTCGATAGTCGCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTGGATGAACCGACCTTCGTCGACAAAGTGATGGTCAACATGTCTGGTGCTGTCCGGGCAATATTTCCAGATTGGTTCCAGGCCATTTTACCTGCACCGCTGGCCTCGGTAGCCTCTATCGTTa                                                                                                                                                                                                                                       <  1:222015/251‑1 (MQ=255)
                                             cGATAGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAACAGAGCTGGCACAAGTGAAACAGTCGCATCTGGAATATTAAGTTGATGAACCCACCTTCTTCGACCAAGTGATGGAAAACCTGTCTGGTTCTGTCCGGGCAATGATGCCAGATGCGTTCAAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTATCGTTAAAAGTGAAAGAGACAAGCTGGCCGCGTTTAATGACCCACAAAAcc                                                                                                                                                                                          >  1:113120/1‑251 (MQ=255)
                                                                                                                     ggAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTATCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGa                                                                                                                                                                        >  1:138265/1‑197 (MQ=255)
                                                                                                                                         cGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTATCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTcgacg                                                                                              <  1:91843/251‑1 (MQ=255)
                                                                                                                                            ccTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTATCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGAtcgtcg                                                                                                                      >  1:107335/1‑224 (MQ=255)
                                                                                                                                                            naTGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTATCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCt                                                                                                                                                >  1:69536/2‑182 (MQ=255)
                                                                                                                                                               ggACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGAAGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTATCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCGGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGc                                                                        <  1:172051/251‑1 (MQ=255)
                                                                                                                                                                      aTGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTATCGTTAaaanngaaaggnnc                                                                                                                                                                                                                         >  2:4935/1‑95 (MQ=255)
                                                                                                                                                                      aTGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTATCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTgcg                                                                                                                                                                  <  1:4935/154‑1 (MQ=255)
                                                                                                                                                                             gTTCTGTCCGGGCAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTATCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTgg                                                                                                                                    <  1:19571/177‑1 (MQ=255)
                                                                                                                                                                                                        cGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTATCGTTAAAAGTGAAAGTGACAAGCTGTCCGCGTTTAATGACCCACAAAACCGTTAGTCGTTTGGCCTGACCTGCGCCAACATGCGTTAAGTCGTGTACTGAGTTGCCGCCAGAGCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTcc                                                             <  1:251676/221‑1 (MQ=255)
                                                                                                                                                                                                        cGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTATCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGa                                                                                                                                                                                                     >  1:10133/1‑85 (MQ=255)
                                                                                                                                                                                                        cGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTATCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGa                                                                                                                                                                                                     <  2:10133/85‑1 (MQ=255)
                                                                                                                                                                                                                      ttACCTGCACCGCTGGCCTCGGTAGCCTCTATCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTcn                                                                                                 <  1:67830/171‑2 (MQ=255)
                                                                                                                                                                                                                                      ccTCGGTAGCCTCTATCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAACACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTATTGTCCTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCCTATTTCCTTATCCGCCTCGTATCAGGGGCGCCCCCCCTGCCTCACGTATATACTTTTGCTCTTTCg  >  1:228023/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          
GGGCCACGTCTGGACCGGTCAGGATGCAAAAGCTAACGGGCTGGTCGATAGTCTCGGGGATTTCGATGATGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCG  >  NC_000913/1850360‑1850839

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