Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,849,454	T→C	G206G (GGT→GGC)	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,849,454	0	T	C	100.0%	133.2 / NA	37	G206G (GGT→GGC)	sppA	protease IV, a signal peptide peptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/19);  total (18/19)

CCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGTACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCGCTAACCGGCAGATCC  >  NC_000913/1849324‑1849584                                                                                                                                   |                                                                                                                                   ccAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATaaa                                                                                                                           >  2:362314/1‑140 (MQ=255) ccAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATaaa                                                                                                                           >  2:709168/1‑140 (MQ=255)   aGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATc                                                                                                                         <  2:227124/140‑1 (MQ=255)                       gCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTc                                                                                                     <  1:738611/140‑1 (MQ=255)                       gCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTc                                                                                                     <  1:709168/140‑1 (MQ=255)                         tGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGt                                                                                                   <  1:524446/140‑1 (MQ=255)                         tGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGt                                                                                                   <  1:127012/140‑1 (MQ=255)                           tctcCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGtga                                                                                                 <  1:362314/140‑1 (MQ=255)                           tctcCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGtga                                                                                                 <  1:750217/140‑1 (MQ=255)                                                cTGCACGGCTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCcg                                                                            <  2:130273/140‑1 (MQ=255)                                                  gCACGGCTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCcgcg                                                                          >  2:782857/1‑140 (MQ=255)                                                    aCGGCTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGaa                                                                        >  2:220178/1‑140 (MQ=255)                                                        cTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCgg                                                                                   >  2:82968/1‑125 (MQ=255)                                                        cTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCgg                                                                                   <  1:82968/125‑1 (MQ=255)                                                        cTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCgg                                                                                   <  1:28468/125‑1 (MQ=255)                                                        cTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCgg                                                                                   >  2:28468/1‑125 (MQ=255)                                                        cTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTg                                                                    >  2:27827/1‑140 (MQ=255)                                                          ttGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGAc                                                                  <  1:353641/140‑1 (MQ=255)                                                                    cGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGa                                                        >  2:181343/1‑140 (MQ=255)                                                                        cTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTgg                                                    <  1:220178/140‑1 (MQ=255)                                                                         tGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGt                                                   >  1:701034/1‑140 (MQ=255)                                                                                cAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCtgt                                            >  2:515000/1‑140 (MQ=255)                                                                                 aaaTCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCtgtg                                           >  1:164267/1‑140 (MQ=255)                                                                                     cGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTg                                                     <  2:258525/126‑1 (MQ=255)                                                                                     cGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTg                                                     >  1:258525/1‑126 (MQ=255)                                                                                     cGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCaa                                       >  2:571669/1‑140 (MQ=255)                                                                                     cGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCaa                                       >  2:695651/1‑140 (MQ=255)                                                                                      gTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTAATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCaaa                                      <  1:245931/140‑1 (MQ=255)                                                                                       ttGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCaaaa                                     <  2:96266/140‑1 (MQ=255)                                                                                        tGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAAc                                    <  1:425095/140‑1 (MQ=255)                                                                                        tGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAAc                                    <  1:228517/140‑1 (MQ=255)                                                                                                 aGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAAt                           <  2:701034/140‑1 (MQ=255)                                                                                                  gCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATa                          >  1:760157/1‑140 (MQ=255)                                                                                                  gCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATa                          >  1:302056/1‑140 (MQ=255)                                                                                                  gCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATa                          >  2:652252/1‑140 (MQ=255)                                                                                                                     gtgtTCCGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCGCTAACCGGCa       <  2:417473/140‑1 (MQ=255)                                                                                                                          ccGCGTGGGCACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCGCTAACCGGCAGATcc  >  1:528310/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                   CCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGTACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCGCTAACCGGCAGATCC  >  NC_000913/1849324‑1849584

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