Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	147,117	T→C	S50S (AGT→AGC)	rpnC →	recombination‑promoting nuclease RpnC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	147,117	0	T	C	100.0%	86.2 / NA	24	S50S (AGT→AGC)	rpnC	recombination‑promoting nuclease RpnC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/14);  total (10/14)

CATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGAT  >  NC_000913/146992‑147251                                                                                                                              |                                                                                                                                       cATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGTGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGCTTCATTGAagagag                                                                                                                          >  1:173895/1‑140 (MQ=255)   tGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGCTTCATTGAag                                                                                                                              <  2:105814/134‑1 (MQ=255)   tGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGCTTCATTGAag                                                                                                                              >  1:105814/1‑134 (MQ=255)          gTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGCTTCATTGAAGAGAGCCTGAAAgg                                                                                                                 <  1:152422/140‑1 (MQ=255)          gTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGCTTCATTGAAGAGAGCCTGAAAgg                                                                                                                 <  1:390199/140‑1 (MQ=255)             tttAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGCTTCATTGAAGAGAGCCTGAAAGGaca                                                                                                              <  1:291524/140‑1 (MQ=255)                     aTTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGCTTCATTGAAGAGAGCCTGAAAGGACACAGCACgg                                                                                                      >  1:75475/1‑140 (MQ=255)                     aTTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGCTTCATTGAAGAGAGCCTGAAAGGACACAGCACgg                                                                                                      >  2:220849/1‑140 (MQ=255)                     aTTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGCTTCATTGAAGAGAGCCTGAAAGGACACAGCACgg                                                                                                      <  2:121760/140‑1 (MQ=255)                                       gaCGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGCTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGc                                                                                    <  2:87157/140‑1 (MQ=255)                                          ggCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGTGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGCTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCaaa                                                                                 <  2:173895/140‑1 (MQ=255)                                            cTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGCTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATg                                                                               >  2:290486/1‑140 (MQ=255)                                               gcgACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGCTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAg                                                                            <  1:220849/140‑1 (MQ=255)                                                                                tACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGCTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGaa                                           <  1:290486/140‑1 (MQ=255)                                                                                tACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGCTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGaa                                           <  2:75475/140‑1 (MQ=255)                                                                                   gcgAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGCTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAacac                                        <  1:237080/140‑1 (MQ=255)                                                                                    cgAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGGGCTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGTTCTATTCCGTGCAAATGCAGGGCAATACCGGGTATCTGCATTTTGTGGTTGAACAcc                                       >  2:404016/1‑140 (MQ=255)                                                                                            gtgACCTCAACACGCTTCATTTAGAGTCGGGGAGCTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagcaag                               <  1:262372/140‑1 (MQ=255)                                                                                             tgACCTCAACACGCTTCATTTAGAGTCGGGGAGCTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagcaagc                              >  1:285664/1‑140 (MQ=255)                                                                                                           ttCATTTAGAGTCGGGGAGCTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATa                        <  1:55462/132‑1 (MQ=255)                                                                                                           ttCATTTAGAGTCGGGGAGCTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATa                        >  2:55462/1‑132 (MQ=255)                                                                                                           ttCATTTAGAGTCGGGGAGCTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAAGAAAATg                >  1:404956/1‑140 (MQ=255)                                                                                                              aTTTAGAGTCGGGGAGCTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAAGAAAATGGcc             >  2:332014/1‑140 (MQ=255)                                                                                                                       cGGGGAGCTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCatg    <  2:404956/140‑1 (MQ=255)                                                                                                                        ggggAGCTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCatga   >  1:324048/1‑140 (MQ=255)                                                                                                                         gggAGCTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCatgat  >  2:219129/1‑140 (MQ=255)                                                                                                                              |                                                                                                                                       CATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGAT  >  NC_000913/146992‑147251

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