Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	892,869	G→A	*301* (TAG→TAA)	rimK →	ribosomal protein S6 modification protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	892,869	0	G	A	100.0%	103.8 / NA	34	*301* (TAG→TAA)	rimK	ribosomal protein S6 modification protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/16);  total (18/16)

TTGGTGATGGAGGTGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAGTCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTCCTGGTCTGGATACGCTG  >  NC_000913/892735‑892992                                                                                                                                       |                                                                                                                            ttGGTGATGGAGGTGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGc                                                                                                                         >  2:1159377/1‑139 (MQ=255)   ggTGATGGAGGTGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCaa                                                                                                                       <  1:567212/139‑1 (MQ=255)    gTGATGGAGGTGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAAt                                                                                                                      <  2:1264255/139‑1 (MQ=255)    gTGATGGAGGTGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAAt                                                                                                                      <  2:103791/139‑1 (MQ=255)       aTGGAGGTGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATcac                                                                                                                   >  1:1120500/1‑139 (MQ=255)         ggAGGTGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACAt                                                                                                                 <  2:1117096/139‑1 (MQ=255)                aaTGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGAt                                                                                                          >  1:1256292/1‑139 (MQ=255)                       cGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGtt                                                                                                   >  2:491098/1‑139 (MQ=255)                                                                aTCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATc                                                          >  2:1204793/1‑139 (MQ=255)                                                                             tGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTa                                             <  2:925787/139‑1 (MQ=255)                                                                              gATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTAc                                            <  1:1204793/139‑1 (MQ=255)                                                                               aTCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACa                                           <  1:298079/139‑1 (MQ=255)                                                                                 ccGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTa                                             >  1:407084/1‑135 (MQ=255)                                                                                 ccGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTa                                             <  2:407084/135‑1 (MQ=255)                                                                                 ccGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACaaa                                         >  2:22984/1‑139 (MQ=255)                                                                                 ccGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACaaa                                         >  2:1211813/1‑139 (MQ=255)                                                                                      ggATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACAtt                                                                                                                <  2:435517/63‑1 (MQ=255)                                                                                      ggATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACAtt                                                                                                                >  1:435517/1‑63 (MQ=255)                                                                                            aaCGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATc                                                                                                         >  1:1218892/1‑64 (MQ=255)                                                                                            aaCGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATc                                                                                                         <  2:1218892/64‑1 (MQ=255)                                                                                              cGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGc                            <  2:1120500/139‑1 (MQ=255)                                                                                                 cACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAAt                                                                                                                      <  2:907964/46‑1 (MQ=255)                                                                                                 cACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAAt                                                                                                                      >  1:907964/1‑46 (MQ=255)                                                                                                   cGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTc                       >  1:6679/1‑139 (MQ=255)                                                                                                       aCGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTc                   <  2:6679/139‑1 (MQ=255)                                                                                                                     ccTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCtt                                                                                       >  1:1135777/1‑57 (MQ=255)                                                                                                                     ccTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCtt                                                                                       <  2:1135777/57‑1 (MQ=255)                                                                                                                     ccTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGtt                                                >  2:467245/1‑96 (MQ=255)                                                                                                                     ccTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGtt                                                <  1:467245/96‑1 (MQ=255)                                                                                                                     ccTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTCCTGGTCTGGATACg     >  1:435842/1‑139 (MQ=255)                                                                                                                     ccTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTCCTGGTCTGGATACg     >  2:986192/1‑139 (MQ=255)                                                                                                                        gAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTCCTGGTCTGGATACGCTg  >  1:239050/1‑139 (MQ=255)                                                                                                                         aaaaCGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGaa                                                       <  1:1010517/85‑1 (MQ=255)                                                                                                                         aaaaCGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGaa                                                       >  2:1010517/1‑85 (MQ=255)                                                                                                                                       |                                                                                                                            TTGGTGATGGAGGTGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAGTCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTCCTGGTCTGGATACGCTG  >  NC_000913/892735‑892992

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