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%	144.5 / NA	45	*301* (TAG→TAA)	rimK	ribosomal protein S6 modification protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/24);  total (21/24)

GGTGATGGAGGTGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAGTCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTCCTGGTCTGGATACGCTGCGTCAAT  >  NC_000913/892737‑892999                                                                                                                                     |                                                                                                                                   ggTGATGGAGGTGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCaa                                                                                                                              <  2:69783/139‑1 (MQ=255) ggTGATGGAGGTGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCaa                                                                                                                              <  1:486071/139‑1 (MQ=255)  gTGATGGAGGTGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAAt                                                                                                                             <  1:1347475/139‑1 (MQ=255)               aTGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTa                                                                                                                      <  1:1315751/133‑1 (MQ=255)               aTGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTa                                                                                                                      >  2:1315751/1‑133 (MQ=255)                tGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCaa                                                                                                                              <  2:356114/124‑1 (MQ=255)                tGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGGTAATCGCaa                                                                                                                              >  1:356114/1‑124 (MQ=255)                  cGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATTACATTACTGATCATg                                                                                                             >  2:1078693/1‑139 (MQ=255)                                           aaaaaCCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAg                                                                                    <  1:304994/139‑1 (MQ=255)                                               accaccGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCg                                                                                        <  2:570111/131‑1 (MQ=255)                                               accaccGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCg                                                                                        >  1:570111/1‑131 (MQ=255)                                               accaccGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGAC‑GAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCtgtg                                                                               >  1:1122433/1‑139 (MQ=255)                                                    cGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCgg                                                                           >  1:425057/1‑139 (MQ=255)                                                           gACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGAC‑GAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTa                                                                   <  2:1122433/139‑1 (MQ=255)                                                                cgcgGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCtt                                                                                              <  2:244392/108‑1 (MQ=255)                                                                cgcgGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCtt                                                                                              >  1:244392/1‑108 (MQ=255)                                                                cgcgGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGa                                                               >  1:94940/1‑139 (MQ=255)                                                                 gcgGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGt                                                                                       >  1:41927/1‑114 (MQ=255)                                                                 gcgGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGt                                                                                       <  2:41927/114‑1 (MQ=255)                                                                      tAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAgg                                                         <  1:460192/139‑1 (MQ=255)                                                                      tAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTATTTGCGTAAGCTGTGCCGGTCTTCTTGTCGAAAGAgg                                                         <  1:1078693/139‑1 (MQ=255)                                                                               ccGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTAc                                                   <  1:47564/136‑1 (MQ=255)                                                                               ccGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTAc                                                   >  2:47564/1‑136 (MQ=255)                                                                                   tGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTa                                            >  1:1193997/1‑139 (MQ=255)                                                                                          aaCGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAAt                                                                                                                             >  1:339438/1‑51 (MQ=255)                                                                                          aaCGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAAt                                                                                                                             <  2:339438/51‑1 (MQ=255)                                                                                          aaCGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTa                                                                                                                      >  2:103460/1‑58 (MQ=255)                                                                                          aaCGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTa                                                                                                                      <  1:103460/58‑1 (MQ=255)                                                                                          aaCGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGaa                                                              <  2:866172/114‑1 (MQ=255)                                                                                          aaCGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGaa                                                              >  1:866172/1‑114 (MQ=255)                                                                                                 cGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGc                                                                                                                                <  2:985307/41‑1 (MQ=255)                                                                                                 cGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGc                                                                                                                                >  1:985307/1‑41 (MQ=255)                                                                                                     aCGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTc                          <  2:1193997/139‑1 (MQ=255)                                                                                                           gAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTCCTGGTc                    <  2:425057/139‑1 (MQ=255)                                                                                                                ttGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCtt                                                                                              >  1:511637/1‑60 (MQ=255)                                                                                                                ttGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCtt                                                                                              <  2:511637/60‑1 (MQ=255)                                                                                                                      gAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGc                                                                                         >  1:239289/1‑59 (MQ=255)                                                                                                                      gAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGc                                                                                         <  2:239289/59‑1 (MQ=255)                                                                                                                      gAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTCCTGGTCTGGATACGCTg         >  2:564788/1‑139 (MQ=255)                                                                                                                        aaaCGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTCCTGGTCTGGATACGCTGCg       >  1:400781/1‑139 (MQ=255)                                                                                                                             ggtggtTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAgg                                                         <  1:1214658/84‑1 (MQ=255)                                                                                                                             ggtggtTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAgg                                                         >  2:1214658/1‑84 (MQ=255)                                                                                                                             ggtggtTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAAt                                              >  2:697731/1‑95 (MQ=255)                                                                                                                             ggtggtTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAAt                                              <  1:697731/95‑1 (MQ=255)                                                                                                                             ggtggtTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTCCTGGTCTGGATACGCTGCGTCAAt  <  2:94940/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                   GGTGATGGAGGTGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAGTCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTCCTGGTCTGGATACGCTGCGTCAAT  >  NC_000913/892737‑892999

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