Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	255,174	(C)6→5	coding (543/1458 nt)	pepD ←	peptidase D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	255,169	0	C	.	100.0%	84.3 / NA	21	coding (548/1458 nt)	pepD	peptidase D
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (12/9);  total (12/9)

CCCAACGTGGATTTCCCCGCCGGAGTGACCGCCTTTCAGACCTTTTAAGGTTAACTTGAAGGTTTCAAAACCAGCTGGAACCGCTTCACGATCTAAATGCAGGTTGGAGGTGAAGTCGATACCCCCCGCACAACCCATGTAGATTTCACCTTCTTCTTCGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGCTTCTT  >  NC_000913/255048‑255283                                                                                                                          |                                                                                                                   cccAACGTGGATTTCCCCGCCGGAGTGACCGCCTTTCAGACCTTTTAAGGTTAACTTGAAGGTTTCAAAACCAGCTGGAACCGCTTCACGATCTAAATGCAGGTTGGAGGTGAAGTCGATAcccccgc                                                                                                              <  1:363548/128‑3 (MQ=255) cccAACGTGGATTTCCCCGCCGGAGTGACCGCCTTTCAGACCTTTTAAGGTTAACTTGAAGGTTTCAAAACCAGCTGGAACCGCTTCACGATCTAAATGCAGGTTGGAGGTGAAGTCGATAcccccgc                                                                                                              >  2:363548/1‑126 (MQ=255) cccAACGTGGATTTCCCCGCCGGAGTGACCGCCTTTCAGACCTTTTAAGGTTAACTTGAAGGTTTCAAAACCAGCTGGAACCGCTTCACGATCTAAATGCAGGTTGGAGGTGAAGTCGATA‑CCCCCGCACAAc                                                                                                        >  2:362354/1‑133 (MQ=255)      cGTGGATTTCCCCGCCGGAGTGACCGCCTTTCAGACCTTTTAAGGTTAACTTGAAGGTTTCAAAACCAGCTGGAACCGCTTCACGATCTAAATGCAGGTTGGAGGTGAAGTCGATA‑CCCCCGCACAACCCATg                                                                                                   >  1:468541/1‑133 (MQ=255)                  cgccgGAGTGACCGCCTTTCAGACCTTTTAAGGTTAACTTGAAGGTTTCAAAACCAGCTGGAACCGCTTCACGATCTAAATGCAGGTTGGAGGTGAAGTCGATA‑CCCCCGCACAACCCATGTAGATTTCACct                                                                                       >  1:627166/1‑133 (MQ=255)                               gCCTTTCAGACCTTTTAAGGTTAACTTGAAGGTTTCAAAACCAGCTGGAACCGCTTCACGATCTAAATGCAGGTTGGAGGTGAAGTCGATA‑CCCCCGCACAACCCATGTAGATTTCACCTTCTTCTTCGGAGt                                                                          <  2:468541/133‑1 (MQ=255)                                         ccTTTTAAGGTTACCTTGAAGGTTTCAAAACCAGCTGGACCCGCTTCACGATCTAAATGCAGGTTGGAGGTGAAGTCGATA‑CCCCCGCACAACCCATGTAGATTTCCCGTTCTTCTTCGGAGTCGGTGTTAAt                                                                <  2:643858/133‑1 (MQ=255)                                         ccTTTTAAGGTTAACTTGAAGGTTTCAAAACCAGCTGGAACCGCTTCACGATCTAAATGCAGGTTGGAGGTGAAGTCGATA‑CCCCCGCACAACCCATGTAGATTTCACCTTCTTCTTCGGAGTCGGTGTTAAt                                                                <  2:44142/133‑1 (MQ=255)                                                       cTTGAAGGTTTCAAAACCAGCTGGAACCGCTTCACGATCTAAATGCAGGTTGGAGGTGAAGTCGATA‑CCCCCGCACAACCCATGTAGATTTCACCTTCTTCTTCGGAGTCGGTGTTAATCAGAATATCAGCCt                                                  >  2:100153/1‑133 (MQ=255)                                                                 tCAAAACCAGCTGGAACCGCTTCACGATCTAAATGCAGGTTGGAGGTGAAGTCGATA‑CCCCCGCACAACCCATGTAGATTTCACCTTCTTCTTCGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGtt                                        >  2:548677/1‑133 (MQ=255)                                                                      aCCAGCTGGAACCGCTTCACGATCTAAATGCAGGTTGGAGGTGAAGTCGATA‑CCCCCGCACAACCCATGTAGATTTCACCTTCTTCTTCGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCt                                   <  1:377623/133‑1 (MQ=255)                                                                         aGCTGGAACCGCTTCACGATCTAAATGCAGGTTGGAGGTGAAGTCGATA‑CCCCCGCACAACCCATGTAGATTTCACCTTCTTCTTCGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTa                                >  1:611494/1‑133 (MQ=255)                                                                           cTGGAACCGCTTCACGATCTAAATGCAGGTTGGAGGTGAAGTCGATA‑CCCCCGCACAACCCATGTAGATTTCACCTTCTTCTTCGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTg                                       >  1:623555/1‑124 (MQ=255)                                                                           cTGGAACCGCTTCACGATCTAAATGCAGGTTGGAGGTGAAGTCGATA‑CCCCCGCACAACCCATGTAGATTTCACCTTCTTCTTCGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTg                                       <  2:623555/124‑1 (MQ=255)                                                                                    cTTCACGATCTAAATGCAGGTTGGAGGTGAAGTCGATA‑CCCCCGCACAACCCATGTAGATTTCACCTTCTTCTTCGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCa                     <  1:362354/133‑1 (MQ=255)                                                                                     ttCACGATCTAAATGCAGGTTGGAGGTGAAGTCGATA‑CCCCCGCACAACCCATGTAGATTTCACCTTCTTCTTCGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCAc                    >  2:557132/1‑133 (MQ=255)                                                                                      tCACGATCTAAATGCAGGTTGGAGGTGAAGTCGATA‑CCCCCGCACAACCCATGTAGATTTCACCTTCTTCTTCGGAGTCGGTGTTAATCAGAATATCAGCCTACAACCAGGTGCACTTTAAGCCTAACGCAcc                   >  1:265620/1‑133 (MQ=255)                                                                                                 aTGCAGGTTGGAGGTGAAGTCGATA‑CCCCCGCACAACCCATGTAGATTTCACCTTCTTCTTCGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCgg        <  2:627166/133‑1 (MQ=255)                                                                                                      ggTTGGAGGTGAAGTCGATA‑CCCCCGCACAACCCATGTAGATTTCACCTTCTTCTTCGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGcttct   >  2:104863/1‑133 (MQ=255)                                                                                                      ggTTGGAGGTGAAGTCGATA‑CCCCCGCACAACCCATGTAGATTTCACCTTCTTCTTCGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGc       <  2:316389/129‑1 (MQ=255)                                                                                                      ggTTGGAGGTGAAGTCGATA‑CCCCCGCACAACCCATGTAGATTTCACCTTCTTCTTCGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGc       >  1:316389/1‑129 (MQ=255)                                                                                                       gTTGGAGGTGAAGTCGATA‑CCCCCGCACAACCCATGTAGATTTCACCTTCTTCTTCGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGcttctt  <  1:557132/133‑1 (MQ=255)                                                                                                       gTTGGAGGTGAAGTCGATA‑CCCCCGCACAACCCATGTAGATTTCACCTTCTTCTTCGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGcttctt  >  1:433703/1‑133 (MQ=255)                                                                                                                          |                                                                                                                   CCCAACGTGGATTTCCCCGCCGGAGTGACCGCCTTTCAGACCTTTTAAGGTTAACTTGAAGGTTTCAAAACCAGCTGGAACCGCTTCACGATCTAAATGCAGGTTGGAGGTGAAGTCGATACCCCCCGCACAACCCATGTAGATTTCACCTTCTTCTTCGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGCTTCTT  >  NC_000913/255048‑255283

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