Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	147,516	C→T	R183R (CGC→CGT)	rpnC →	recombination‑promoting nuclease RpnC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	147,516	0	C	T	100.0%	108.2 / NA	31	R183R (CGC→CGT)	rpnC	recombination‑promoting nuclease RpnC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/12);  total (19/12)

TCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCGGAAGTCAGT  >  NC_000913/147440‑147589                                                                             |                                                                          tCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGTCAgc                                                                       <  1:35614/81‑1 (MQ=255) tCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGTCAgc                                                                       >  2:35646/1‑81 (MQ=255) tCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGTCAgc                                                                       >  1:35654/1‑81 (MQ=255)  cACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGTCAgcg                                                                      >  1:35658/1‑81 (MQ=255)       tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGTCAGCGCGACt                                                                 <  1:35615/81‑1 (MQ=255)       tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGTCAGCGCGACt                                                                 <  2:35616/81‑1 (MQ=255)             cGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGTCAGCGCGACTTAATGt                                                           >  1:35651/1‑81 (MQ=255)               gATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGTCAGCGCGACTTAATGTTa                                                         <  2:35617/81‑1 (MQ=255)                  gACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGTCAGCGCGACTTAATGTTATTg                                                      >  2:35665/1‑81 (MQ=255)                      aaaTCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGTCAGCGCGACTTAATGTTATTGCTTg                                                  <  2:35618/81‑1 (MQ=255)                        aTCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGTCAGCGCGACTTAATGTTATTGCTTGAg                                                >  2:35640/1‑81 (MQ=255)                          cATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGTCAGCGCGACTTAATGTTATTGCTTGAGc                                               >  1:35685/1‑80 (MQ=255)                              cAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGTCAGCGCGACTTAATGTTATTGCTTGAGCAACTg                                          <  2:35619/81‑1 (MQ=255)                               aaCATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGTCAGCGCGACTTAATGTTATTGCTTGAGCAACTgg                                         >  2:35657/1‑81 (MQ=255)                                     ggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGTCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGc                                   >  2:35703/1‑81 (MQ=255)                                     ggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGTCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGc                                   >  2:35641/1‑81 (MQ=255)                                     ggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGTCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGc                                   <  1:35620/81‑1 (MQ=255)                                                 ttCTCGAACTACTGCAAAAACATATTCGTCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAg                       <  1:35621/81‑1 (MQ=255)                                                   ctcGAACTACTGCAAAAACATATTCGTCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAggg                     >  2:35661/1‑81 (MQ=255)                                                      gAACTACTGCAAAAACATATTCGTCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGt                    <  2:35622/79‑1 (MQ=255)                                                      gAACTACTGCAAAAACATATTCGTCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTa                   >  1:35673/1‑80 (MQ=255)                                                       aaCTACTGCAAAAACATATTCGTCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTaca                 >  2:35677/1‑81 (MQ=255)                                                        actactGCAAAAACATATTCGTCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTacac                >  1:35708/1‑81 (MQ=255)                                                           actGCAAAAACATATTCGTCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAg             >  1:35714/1‑81 (MQ=255)                                                               cAAAAACATATTCGTCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCgg          >  1:35674/1‑80 (MQ=255)                                                               cAAAAACATATTCGTCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCgg          <  1:35623/80‑1 (MQ=255)                                                                aaaaaCATATTCGTCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCgg          >  1:35688/1‑79 (MQ=255)                                                                  aaaCATATTCGTCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCGGAAGt      <  1:35624/81‑1 (MQ=255)                                                                   aaCATATTCGTCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCGGAAGTc     <  2:35625/81‑1 (MQ=255)                                                                    aCATATTCGTCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCGGAAGt      >  1:35669/1‑79 (MQ=255)                                                                        atTCGTCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGATGAAGGGTACACTAGCGGAAGTCagt  >  2:35664/1‑79 (MQ=255)                                                                             |                                                                          TCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCGGAAGTCAGT  >  NC_000913/147440‑147589

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