Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,850,562	T→C	F576L (TTC→CTC)	sppA →	protease IV, a signal peptide peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,850,562	0	T	C	100.0%	92.2 / NA	27	F576L (TTC→CTC)	sppA	protease IV, a signal peptide peptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/12);  total (15/12)

TGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTT  >  NC_000913/1850429‑1850691                                                                                                                                      |                                                                                                                                  tgCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGCTCCAgg                                                                                                                             <  1:75505/140‑1 (MQ=255) tgCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGCTCCAgg                                                                                                                             <  1:507928/140‑1 (MQ=255)                                  gAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCct                                                                                            <  1:124590/140‑1 (MQ=255)                                         tGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGt                                                                                     >  2:429593/1‑140 (MQ=255)                                                  gAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAaaagtga                                                                            <  1:1004/140‑1 (MQ=255)                                                          cGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAaagtgaaagtga                                                                      >  2:568134/1‑138 (MQ=255)                                                          cGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAaagtgaaagtga                                                                      <  1:568134/138‑1 (MQ=255)                                                          cGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACa                                                                    >  2:354842/1‑140 (MQ=255)                                                             tgatgaACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGc                                                                 >  2:540422/1‑140 (MQ=255)                                                                   aCCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCcg                                                           <  1:716708/140‑1 (MQ=255)                                                                       aCCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGtt                                                       >  2:322597/1‑140 (MQ=255)                                                                            cttcGTCAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATg                                                  <  2:90205/140‑1 (MQ=255)                                                                               cGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGAcc                                               <  1:322597/140‑1 (MQ=255)                                                                                  cAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGAcc                                               >  1:721530/1‑137 (MQ=255)                                                                                  cAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGAcc                                               <  2:721530/137‑1 (MQ=255)                                                                                  cAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCcac                                            >  1:634040/1‑140 (MQ=255)                                                                                       tGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAAc                                       >  1:674417/1‑140 (MQ=255)                                                                                         aTGGACAACATGTCTGGTTCTTTCCGGGCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCg                                     >  2:517396/1‑140 (MQ=255)                                                                                                 caTGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGt                             <  1:429593/140‑1 (MQ=255)                                                                                                 caTGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGt                             <  2:634040/140‑1 (MQ=255)                                                                                                      cTGGTTCTGTCCGGGCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGc                        >  1:319829/1‑140 (MQ=255)                                                                                                               tCCGGGCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTgc               >  1:556565/1‑140 (MQ=255)                                                                                                                   ggCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCa           >  2:301260/1‑140 (MQ=255)                                                                                                                   ggCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTACCTGACCTGCGCCa           >  2:627307/1‑140 (MQ=255)                                                                                                                    gCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGAc                   >  2:31350/1‑131 (MQ=255)                                                                                                                    gCAATGTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGAc                   <  1:31350/131‑1 (MQ=255)                                                                                                                         gTTGCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGc     >  2:208842/1‑140 (MQ=255)                                                                                                                            gCCAGATGCGCTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTACCTGACCTGCGCCAACATGCGtt  <  1:627307/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                  TGCGGTCGCCAAAGCAGCAGAGCTGGCAAAAGTGAAACAGTGGCATCTGGAATATTACGTTGATGAACCGACCTTCTTCGACAAAGTGATGGACAACATGTCTGGTTCTGTCCGGGCAATGTTGCCAGATGCGTTCCAGGCCATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTT  >  NC_000913/1850429‑1850691

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑