Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,849,310	Δ1 bp	coding (474/1857 nt)	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,307	0	G	.	100.0%	180.6 / NA	43	coding (471/1857 nt)	sppA	protease IV, a signal peptide peptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (22/21);  total (22/21)

CAATATCACCGGTATTGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAG  >  NC_000913/1849169‑1849421                                                                                                                                           |                                                                                                                   cAATATCACCGGTATTGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAgg                                                                                                                   <  2:857139/140‑1 (MQ=255)  aaTATCACCGGTATTGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAggg                                                                                                                  <  2:594719/140‑1 (MQ=255)  aaTATCACCGGTATTGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAggg                                                                                                                  <  1:10804/140‑1 (MQ=255)            gTATTGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAgg                                                                                                                   >  1:196298/1‑129 (MQ=255)            gTATTGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAgg                                                                                                                   <  2:196298/129‑1 (MQ=255)             tATTGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTAtct                                                                                                      >  2:766057/1‑140 (MQ=255)                tgtgATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGc                                                                                                   >  2:653620/1‑140 (MQ=255)                       gATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTc                                                                                            >  1:1070072/1‑140 (MQ=255)                       gATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTc                                                                                            >  2:582177/1‑140 (MQ=255)                                    tCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAAttt                                                                               >  2:545238/1‑140 (MQ=255)                                                aCCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGc                                                                   <  1:582177/140‑1 (MQ=255)                                                    aCCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAgg                                                               >  2:176361/1‑140 (MQ=255)                                                    aCCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAgg                                                               >  1:972192/1‑140 (MQ=255)                                                               aGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATc                                                    >  1:529490/1‑140 (MQ=255)                                                                  aCATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAgg                                                               <  2:939980/126‑1 (MQ=255)                                                                  aCATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAgg                                                               >  1:939980/1‑126 (MQ=255)                                                                   cATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCa                                                >  1:977191/1‑140 (MQ=255)                                                                    aTCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCAc                                               <  2:972192/140‑1 (MQ=255)                                                                      cGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCGATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACgg                                             >  2:951476/1‑140 (MQ=255)                                                                      cGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTg                                                       >  1:339993/1‑130 (MQ=255)                                                                      cGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTg                                                       <  2:339993/130‑1 (MQ=255)                                                                        gCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCa                                                >  1:757902/1‑135 (MQ=255)                                                                        gCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATAATATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCa                                                <  2:757902/135‑1 (MQ=255)                                                                         cAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCtt                                          >  2:402272/1‑140 (MQ=255)                                                                          aaaGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCttt                                         <  1:545238/140‑1 (MQ=255)                                                                            aGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGc                                       <  1:1036610/140‑1 (MQ=255)                                                                             gCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGcc                                      >  1:73198/1‑140 (MQ=255)                                                                              ctctGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGcca                                     >  1:312561/1‑140 (MQ=255)                                                                              ctctGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGcca                                     >  1:138818/1‑140 (MQ=255)                                                                                 tGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCt                                           >  2:700754/1‑131 (MQ=255)                                                                                 tGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCt                                           <  1:700754/131‑1 (MQ=255)                                                                                 tGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGccacca                                  <  1:818262/140‑1 (MQ=255)                                                                                 tGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGccacca                                  <  1:738695/140‑1 (MQ=255)                                                                                 tGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGccacca                                  <  2:529490/140‑1 (MQ=255)                                                                                 tGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGccacca                                  <  1:176361/140‑1 (MQ=255)                                                                                      gagTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGACCTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGACGTGGTTGATCTGCACGGCTTTGCCACCAACGGt                             <  2:463372/140‑1 (MQ=255)                                                                                      gagTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGTCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGt                             <  2:977191/140‑1 (MQ=255)                                                                                      gagTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGt                             <  2:312561/140‑1 (MQ=255)                                                                                      gagTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGt                             <  2:138818/140‑1 (MQ=255)                                                                                      gagTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGt                             <  2:1070072/140‑1 (MQ=255)                                                                                      gagTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGt                             <  1:447880/140‑1 (MQ=255)                                                                                       agTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGTc                            <  2:588242/140‑1 (MQ=255)                                                                                        gTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGt                             <  2:782945/138‑1 (MQ=255)                                                                                        gTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGt                             >  1:782945/1‑138 (MQ=255)                                                                                         tttCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGTCTg                          <  1:402272/140‑1 (MQ=255)                                                                                             gTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGTCTGtact                      >  2:187655/1‑140 (MQ=255)                                                                                                      ggAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGTCTGTACTACAAATCGt             >  1:281696/1‑140 (MQ=255)                                                                                                                 tATGCCGTTGGCGAGAACTACAGCCA‑GGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAg  >  2:618435/1‑140 (MQ=255)                                                                                                                                           |                                                                                                                   CAATATCACCGGTATTGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGTGGTTGATCTGCACGGCTTTGCCACCAACGGTCTGTACTACAAATCGTTGCTGGATAAG  >  NC_000913/1849169‑1849421

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: ‑