Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,543,611	C→T	L149F (CTT→TTT)	rpnA →	recombination‑promoting nuclease RpnA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,543,611	0	C	T	100.0%	87.2 / NA	25	L149F (CTT→TTT)	rpnA	recombination‑promoting nuclease RpnA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/12);  total (13/12)

CCACGGTATTGAAAGCCCCTATCCCTATTCGCTGTGTTGGCTGGATTGTTTCGCCGATCCCAAACTGGCAAGGCAGCTTTATGCCTCCGCATTTCCGCTGATTGATGTCACCGTCATGCCTGATGATGAAATCATGCAGCACCGACGTATGGC  >  NC_000913/3543535‑3543687                                                                             |                                                                             ccACGGTATTGAAAGCCCCTATCCCTATTCGCTGTGTTGGCTGGATTGTTTCGCCGATCCCAAACTGGCAAGGCAGtttt                                                                           <  1:679499/80‑1 (MQ=255)       tATTGAAAGCCCCTATCCCTATTCGCTGTGTTGGCTGGATTGTTTCGCCGATCCCAAACTGGCAAGGCAGTTTTATGCCTc                                                                    <  1:679500/81‑1 (MQ=255)           gAAAGCCCCTATCCCTATTCGCTGTGTTGGCTGGATTGTTTCGCCGATCCCAAACTGGCAAGGCAGTTTTATGCCTCcaca                                                                >  2:679536/1‑81 (MQ=255)           gAAAGCCCCTATCCCTATTCGCTGTGTTGGCTGGATTGTTTCGCCGATCCCAAACTGGCAAGGCAGTTTTATGCCTCCGCa                                                                >  2:679534/1‑81 (MQ=255)              aGCCCCTATCCCTATTCGCTGTGTTGGCTGGATTGTTTCGCCGATCCCAAACTGGCAAGGCAGTTTTATGCCTCCGCAttt                                                             >  2:679535/1‑81 (MQ=255)              aGCCCCTATCCCTATTCGCTGTGTTGGCTGGATTGTTTCGCCGATCCCAAACTGGCAAGGCAGTTTTATGCCTCCGCAttt                                                             >  1:679542/1‑81 (MQ=255)               gCCCCTATCCCTATTCGCTGTGTTGGCTGGATTGTTTCGCCGATCCCAAACTGGCAAGGCAGTTTTATGCCTCCGCATTTc                                                            <  2:679501/81‑1 (MQ=255)                   ctatccctatTCGCTGTGTTGGCTGGATTGTTTCGCCGATCCCAAACTGGCAAGGCAGTTTTATGCCTCCGCATTTCCGCt                                                        >  2:679538/1‑81 (MQ=255)                     atccctatTCGCTGTGTTGGCTGGATTGTTTCGCCGATCCCAAACTGGCAAGGCAGTTTTATGCCTCCGCATTTCCGCtga                                                      >  2:679533/1‑81 (MQ=255)                       ccctatTCGCTGTGTTGGCTGGATTGTTTCGCCGATCCCAAACTGGCAAGGCAGTTTTATGCCTCCGCATTTCCGCtgatt                                                    >  1:679540/1‑81 (MQ=255)                        cctatTCGCTGTGTTGGCTGGATTGTTTCGCCGATCCCAAACTGGCAAGGCAGTTTTATGCCTCCGCATTTCCGCtgatt                                                    <  2:679502/80‑1 (MQ=255)                         ctatTCGCTGTGTTGGCTGGATTGTTTCGCCGATCCCAAACTGGCAAGGCAGTTTTATGCCTCCGCATTTCCGCtgattga                                                  <  1:679503/81‑1 (MQ=255)                            ttCGCTGTGTTGGCTGGATTGTTTCGCCGATCCCAAACTGGCAAGGCAGTTTTATGCCTCCGCATTTCCGCtgattgat                                                 <  2:679504/79‑1 (MQ=255)                                cTGTGTTGGCTGGATTGTTTCGCCGATCCCAAACTGGCAAGGCAGTTTTATGCCTCCGCATTTCCGCTGATTGATGTCAcc                                           <  2:679505/81‑1 (MQ=255)                                       ggCTGGATTGTTTCGCCGATCCCAAACTGGCAAGGCAGTTTTATGCCTCCGCATTTCCGCTGATTGATGTCACCGTCATg                                     <  1:679506/80‑1 (MQ=255)                                          tGGATTGTTTCGCCGATCCCAAACTGGCAAGGCAGTTTTATGCCTCCGCATTTCCGCTGATTGATGTCACCGTCATGCCt                                  >  2:679551/1‑80 (MQ=255)                                                   tCGCCGATCCCAAACTGGCAAGGCAGTTTTATGCCTCCGCATTTCCGCTGATTGATGTCACCGTCATGCCTGATGATGaaa                        >  2:679548/1‑81 (MQ=255)                                                     gCCGATCCCAAACTGGCAAGGCAGTTTTATGCCTCCGCATTTCCGCTGATTGATGTCACCGTCATGCCTGATGATGAAAt                       >  1:679539/1‑80 (MQ=255)                                                     gCCGATCCCAAACTGGCAAGGCAGTTTTATGCCTCCGCATTTCCGCTGATTGATGTCACCGTCATGCCTGATGATGAAAt                       >  2:679550/1‑80 (MQ=255)                                                     gCCGATCCCAAACTGGCAAGGCAGTTTTATGCCTCCGCATTTCCGCTGATTGATGTCACCGTCATGCCTGATGATGAAATc                      <  2:679507/81‑1 (MQ=255)                                                        gATCCCAAACTGGCAAGGCAGTTTTATGCCTCCGCATTTCCGCTGATTGATGTCACCGTCATGCCTGATGATGAAATCATg                   >  2:679552/1‑81 (MQ=255)                                                          tCCCAAACTGGCAAGGCAGTTTTATGCCTCCGCATTTCCGCTGATTGATGTCACCGTCATGCCTGATGATGAAATCATg                   >  1:679543/1‑79 (MQ=255)                                                            ccAAACTGGCAAGGCAGTTTTATGCCTCCGCATTTCCGCTGATTGATGTCACCGTCATGCCTGATGATGAAATCATgcagc               <  1:679508/81‑1 (MQ=255)                                                                    gCAAGGCAGTTTTATGCCTCCGCATTTCCGCTGATTGATGTCACCGTCATGCCTGATGATGAAATCATGCAGCACCGAc         <  2:679509/79‑1 (MQ=255)                                                                         gCAGTTTTATGCCTCCGCATTTCCGCTGATTGATGTCACCGTCATGCCTGATGTTGAAATCATGCAGCACCGACGTATGgc  <  2:679510/81‑1 (MQ=255)                                                                         gCAGTTTTATGCCTCCGCATTTCCGCTGATTGATGTCACCGTCATGCCTGATGATGAAATCATGCAGCACCGACGTATGgc  <  2:679511/81‑1 (MQ=255)                                                                             |                                                                             CCACGGTATTGAAAGCCCCTATCCCTATTCGCTGTGTTGGCTGGATTGTTTCGCCGATCCCAAACTGGCAAGGCAGCTTTATGCCTCCGCATTTCCGCTGATTGATGTCACCGTCATGCCTGATGATGAAATCATGCAGCACCGACGTATGGC  >  NC_000913/3543535‑3543687

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑