Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,643,610	C→T	G525S (GGC→AGC)	prlC ←	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,643,610	0	C	T	100.0%	102.6 / NA	32	G525S (GGC→AGC)	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/14);  total (18/14)

CGGGCGGAACTCGGCATGAAGGCGGAAATCAAACAGGCCGAACTCCAGCTGACGCAGAATAAACAGCGCCGCCTGGTAGTTCTTCGCCGCCAGCATTTTATCCAGCAACTCTTTCGGCAGCGGTTCGCCGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCGGCAGTTCGACCGCATCCCACGGCACACCGCTGATACCGGAAACACCAGCGGTTTCGATGCGGGTC  >  NC_000913/3643482‑3643746                                                                                                                                 |                                                                                                                                         cGGGCGGAACTCGGCATGAAGGCGGAAATCAAACAGGCCGAACTCCAGCTGACGCAGAATAAACAGCGCCGCCTGGTAGTTCTTCGCCGCCAGCATTTTATCCAGCAACTCTTTCGGCAGCGGTTCGCTGGTTTCATAGt                                                                                                                               <  2:575529/140‑1 (MQ=255)                       cGGAAATCAAACAGGCCGAACTCCAGCTGACGCAGAATAAACAGCGCCGCCTGGTAGTTCTTCGCCGCCAGCATTTTATCCAGCAACTCTTTCGGCAGCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAgcgc                                                                                                         >  2:412223/1‑140 (MQ=255)                                      ccGAACTCCAGCTGACGCAGAATAAACAGCGCCGCCTGGTAGTTCTTCGCCGCCAGCATTTTATCCAACAACTCTTTCGGCAGCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCa                                                                                          >  1:595798/1‑140 (MQ=255)                                       cGAACTCCAGCTGACGCCGAATCAACAGCGCTGCCCGGTCGTTCTTCGTCGCCAGCATTTTATCCAGCAACTCTTTCGGCAGCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCAc                                                                                         <  2:224554/140‑1 (MQ=255)                                       cGAACTCCAGCTGACGCAGAATAAACAGCGCCGCCTGGTAGTTCTTCGCCGCCAGCATTTTATCCAGCAACTCTTTCGGCAGCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCAc                                                                                         <  2:426164/140‑1 (MQ=255)                                                    aCGCAGAATAAACAGCGCCGCCTGGTAGTTCTTCGCCGCCAGCATTTTATCCAGCAACTCTTTCGGCAGCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATaa                                                                            <  1:108574/140‑1 (MQ=255)                                                    aCGCAGAATAAACAGCGCCGCCTGGTAGTTCTTCGCCGCCAGCATTTTATCCAGCAACTCTTTCGGCAGCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATaa                                                                            >  2:108574/1‑140 (MQ=255)                                                              aaCAGCGCCGCCTGGTAGTTCTTCGCCGCCAGCATTTTATCCAGCAACTCTTTCGGCAGCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCgg                                                                  >  2:724563/1‑140 (MQ=255)                                                                    gccgccTGGTAGTTCTTCGCCGCCAGCATTTTATCCAGCAACTCTTTCGGCAGCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCATCGCCTCCGCCTCCCAGGACCAGTTTACCATTAACTGACTCGGCAGTTc                                                            >  2:235592/1‑140 (MQ=255)                                                                    gccgccTGGTAGTTCTTCGCCGCCAGCATTTTATCCAGCAACTCTTTCGGCAGCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCGGCAGTTc                                                            >  2:751596/1‑140 (MQ=255)                                                                    gccgccTGGTAGTTCTTCGCCGCCAGCATTTTATCCAGCAACTCTTTCGGCAGCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCGGCAGTTc                                                            >  2:668160/1‑140 (MQ=255)                                                                    gccgccTGGTAGTTCTTCGCCGCCAGCATTTTATCCAGCAACTCTTTCGGCAGCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCGGCAGTTc                                                            >  2:59252/1‑140 (MQ=255)                                                                        ccTGGTAGTTCTTCGCCGCCAGCATTTTATCCAGCAACTCTTTCGGCAGCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCGGCAGTTCGAcc                                                        >  1:495403/1‑140 (MQ=255)                                                                                ttcttcGCCGCCAGCATTTTATCCAGCAACTCTTTCGGCAGCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCGGCAGTTCGACCGCATCCCa                                                <  2:425720/140‑1 (MQ=255)                                                                                ttcttcGCCGCCAGCATTTTATCCAGCAACTCTTTCGGCAGCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCGGCAGTTCGACCGCATCCCa                                                <  1:526164/140‑1 (MQ=255)                                                                                  cttcGCCGCCAGCATTTTATCCAGCAACTCTTTCGGCAGCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAGCTGTCTCGGCAGTTCGACCGCATCCCACg                                              >  1:336199/1‑140 (MQ=255)                                                                                     cgccgcCAGCATTTTATCCAGCCACCCTTCCGGCAGCCGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCGGCAGTTCGACCGCATCCCACGGca                                           <  1:235592/140‑1 (MQ=255)                                                                                     cgccgcCAGCATTTTATCCAGCAACTCTTTCGGCAGCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAGACTGACTCGGCAGTTCGACCGCATCCCACGGca                                           >  1:791424/1‑140 (MQ=255)                                                                                      gccgccAGCATTTTATCCAGCAACTCTTTCGGCAGCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCGGCAGTTCGACCGCATCCCACGGcac                                          >  1:823400/1‑140 (MQ=255)                                                                                          ccAGCATTTTATCCAGCAACTCTTTCGGCAGCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCGGCAGTTCGACCGCATCCCACGGCACACCg                                      <  1:78491/140‑1 (MQ=255)                                                                                                ttttATCCAGCAACTCTTTCGGCAGCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCGGCAGTTCGACCGCATCCCACGGCACACCGCTGATa                                <  2:443029/140‑1 (MQ=255)                                                                                                        aGCAACTCTTTCGGCAGCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTCACTCGGCAGTTCTACCACATCCCACAGCACCCCGCTGATACCGGAAac                        >  2:254832/1‑140 (MQ=255)                                                                                                         gCAACTCTTTCGGCAGCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCGGCAGTTCGACCGCATCCCACGGCACACCGCTGATACCGGAAaca                       >  2:543798/1‑140 (MQ=255)                                                                                                           aaCTCTTTCGGCAGCGGTTCGCTGGTTTCATTGTGCCCATAGATAAACCCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCGGCAGTTCGACCGCATCCCACGGCACACCGCTGATACCGGAAACAcc                     <  2:336199/140‑1 (MQ=255)                                                                                                               ctTTCGGCAGCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCGGCAGTTCGACCGCATCCCACGGCACACCGCTGATACCGGAAACACCAGCg                 >  1:122937/1‑140 (MQ=255)                                                                                                                        gCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCGGCAGTTCGACCGCATCCCACGGCACACCGCTGATACCGGAAACACCAGCGGTTTCGATg        >  2:643396/1‑140 (MQ=255)                                                                                                                        gCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCGGCAGTTCGACCGCATCCCACGGCACACCGCTGATACCGGAAACACCAGCGGTTTCGATg        >  1:138502/1‑140 (MQ=255)                                                                                                                        gCGGTTCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCGGCAGTTCGACCGCATCCCACGGCACACCGCTGATACCGGAAACACCAGCGGTTTCGATg        >  1:12658/1‑140 (MQ=255)                                                                                                                            ttCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCGGCAGTTCGACCGCATCCCACGGCACACCGCTGATACCGGAAACACCAGCGGTTTCGATGCggg    <  2:138502/140‑1 (MQ=255)                                                                                                                            ttCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCGGCAGTTCGACCGCATCCCACGGCACACCGCTGATACCGGAAACACCAGCGGTTTCGATGCggg    <  2:122937/140‑1 (MQ=255)                                                                                                                            ttCGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCGGCAGTTCGACCGCATCCCACGGCACACCGCTGATACCGGAAACACCAGCGGTTTCGATGCggg    <  2:595798/140‑1 (MQ=255)                                                                                                                              cGCTGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCGGCAGTTCGACCGCATCCCACGGCACACCGCTGATACCGGAAACACCAGCGGTTTCGATGCGGGTc  <  1:668160/140‑1 (MQ=255)                                                                                                                                 |                                                                                                                                         CGGGCGGAACTCGGCATGAAGGCGGAAATCAAACAGGCCGAACTCCAGCTGACGCAGAATAAACAGCGCCGCCTGGTAGTTCTTCGCCGCCAGCATTTTATCCAGCAACTCTTTCGGCAGCGGTTCGCCGGTTTCATAGTGACCAGAGATAAACGCCAGCGCCTCCGGCTCCCAGCACCAGTTTTCCATAAACTGACTCGGCAGTTCGACCGCATCCCACGGCACACCGCTGATACCGGAAACACCAGCGGTTTCGATGCGGGTC  >  NC_000913/3643482‑3643746

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