Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,391,886	C→T	W48* (TGG→TAG)	tldD ←	metalloprotease TldD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,391,886	0	C	T	100.0%	62.7 / NA	19	W48* (TGG→TAG)	tldD	metalloprotease TldD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (6/13);  total (6/13)

CGCCAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCGGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCACT  >  NC_000913/3391765‑3392010                                                                                                                          |                                                                                                                             cgcCAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCGGTCTTCTAAAACCTAGGATTCGTGATAGCTCg                                                                                                            <  1:269479/140‑1 (MQ=255)   ccAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGGTCGCTGCTGTAAGAACCCCTTTTATTACTGCGGTCTTCTAAAACCTAGGATTCGTCATAGCTCGAc                                                                                                          >  2:12933/1‑140 (MQ=255)    cagcagGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCGGTCTTCTAAAACCTAGGATTCGTGATAGCTCGACt                                                                                                         <  1:213199/140‑1 (MQ=255)                                    cGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCGGTCTTCTAAAACCTAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGt                                                                         <  1:494431/140‑1 (MQ=255)                                      gTTTTTTCACCGCTGATTGCACACACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCGGTCTTCTAAAACCTAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTc                                                                       <  2:6111/140‑1 (MQ=255)                                         ttttCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCGGTCTTCTAAAACCTAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGc                                                                    <  1:286831/140‑1 (MQ=255)                                             cACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCGGTCTTCTAAAACCTAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGt                                                                <  1:14220/140‑1 (MQ=255)                                             cACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCGGTCTTCTAAAACCTAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGt                                                                <  2:90525/140‑1 (MQ=255)                                                    gATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCGGTCTTCTAAAACCTAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGa                                                         <  1:86947/140‑1 (MQ=255)                                                                   aaCGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCGGTCTTCTAAAACCTAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGt                                          <  1:55887/140‑1 (MQ=255)                                                                       cccTGATCGATGTTGTAAGAACCATCTTTAATAATGCGGTCTTCTAAAACCTAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTg                                      <  2:43956/140‑1 (MQ=255)                                                                       cccTGATCGATGTTGTAAGAACCATCTTTAATAATGCGGTCTTCTAAAACCTAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTg                                      <  1:120333/140‑1 (MQ=255)                                                                                           aCCATCTTTAATAATGCGGTCTTCTAAAACCTAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTcgccg                  >  1:43488/1‑140 (MQ=255)                                                                                               tCTTTAATAATGCGGTCTTCTAAAACCTAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAg              >  1:377949/1‑140 (MQ=255)                                                                                                        atGCGGTCTTCTAAAACCTAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTc     <  2:377949/140‑1 (MQ=255)                                                                                                         tGCGGTCTTCTAAAACCTAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGtt      <  1:388557/138‑1 (MQ=255)                                                                                                         tGCGGTCTTCTAAAACCTAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGtt      >  2:388557/1‑138 (MQ=255)                                                                                                         tGCGGTCTTCTAAAACCTAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCa    >  2:170793/1‑140 (MQ=255)                                                                                                           cGGTCTTCTAAAACCTAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCACt  >  1:388412/1‑140 (MQ=255)                                                                                                           cGGTCTTCTAAAACCTAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCACt  >  2:403980/1‑140 (MQ=255)                                                                                                                          |                                                                                                                             CGCCAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCGGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCACT  >  NC_000913/3391765‑3392010

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