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%	143.0 / 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)

TTGGTGATGGAGGTGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAGTCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTCCTGGTCTGGATACGCTG  >  NC_000913/892735‑892992                                                                                                                                       |                                                                                                                            ttGGTGATGGAGGTGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGc                                                                                                                         <  1:561861/139‑1 (MQ=255)      gATGGAGGTGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACAGCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATca                                                                                                                    <  2:728375/139‑1 (MQ=255)               gAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGa                                                                                                           >  2:235817/1‑139 (MQ=255)                              gCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTg                                                                                            <  2:324913/139‑1 (MQ=255)                                 ggaaggaaTAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTgcgc                                                                                         >  2:1239343/1‑139 (MQ=255)                                    aggaaTAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCttt                                                                                      <  1:235817/139‑1 (MQ=255)                                        aTAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGc                                                                                                                         <  1:687641/100‑1 (MQ=255)                                        aTAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGc                                                                                                                         >  2:687641/1‑100 (MQ=255)                                                 accaccGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCTTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTCCCTGCGCTTTTTGCGTAAGCtgt                                                                         >  2:726268/1‑139 (MQ=255)                                                  ccaccGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCtgtg                                                                        <  1:592469/139‑1 (MQ=255)                                                      cGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCtgtg                                                                        >  2:498466/1‑135 (MQ=255)                                                      cGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCtgtg                                                                        <  1:498466/135‑1 (MQ=255)                                                      cGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCgg                                                                    >  1:25975/1‑139 (MQ=255)                                                                  cgcgGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCaa                                                                                                                       >  2:736429/1‑76 (MQ=255)                                                                  cgcgGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCaa                                                                                                                       <  1:736429/76‑1 (MQ=255)                                                                        tAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAgg                                                  <  2:25975/139‑1 (MQ=255)                                                                        tAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAgg                                                  <  1:1239343/139‑1 (MQ=255)                                                                         aaaaTGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAgg                                                  >  1:800568/1‑138 (MQ=255)                                                                         aaaaTGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAgg                                                  <  2:800568/138‑1 (MQ=255)                                                                             tGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTa                                             <  1:720544/139‑1 (MQ=255)                                                                                 ccGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACAGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCtt                                                                                       >  1:381434/1‑93 (MQ=255)                                                                                 ccGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACAGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCtt                                                                                       <  2:381434/93‑1 (MQ=255)                                                                                        aTCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTa                                                                                                               >  2:65967/1‑62 (MQ=255)                                                                                        aTCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTa                                                                                                               <  1:65967/62‑1 (MQ=255)                                                                                            aaCGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCaa                                                                                                                       >  2:980647/1‑50 (MQ=255)                                                                                            aaCGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTAAATCGCaa                                                                                                                       <  1:980647/50‑1 (MQ=255)                                                                                                 cACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATc                                                                                                         <  1:1154081/59‑1 (MQ=255)                                                                                                 cACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATc                                                                                                         >  2:1154081/1‑59 (MQ=255)                                                                                                        cGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTg                                                                                                >  2:546447/1‑61 (MQ=255)                                                                                                        cGACAGAATATTGCCTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTg                                                                                                <  1:546447/61‑1 (MQ=255)                                                                                                                     ccTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAagag                                                   <  2:1252678/93‑1 (MQ=255)                                                                                                                     ccTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAagag                                                   >  1:1252678/1‑93 (MQ=255)                                                                                                                     ccTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGa                                  <  1:752238/110‑1 (MQ=255)                                                                                                                     ccTGAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGa                                  >  2:752238/1‑110 (MQ=255)                                                                                                                        gAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGt                        <  1:1244820/117‑1 (MQ=255)                                                                                                                        gAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGt                        >  2:1244820/1‑117 (MQ=255)                                                                                                                        gAAAACGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTCCTGGTCTGGATACGCTg  >  1:509793/1‑139 (MQ=255)                                                                                                                          aaaCGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTa                                                            >  2:319680/1‑79 (MQ=255)                                                                                                                          aaaCGGGTGGTTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTa                                                            <  1:319680/79‑1 (MQ=255)                                                                                                                                   gtTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCg                                                                     >  1:1470284/1‑61 (MQ=255)                                                                                                                                   gtTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCg                                                                     <  2:1470284/61‑1 (MQ=255)                                                                                                                                   gtTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCg                             <  1:1056061/101‑1 (MQ=255)                                                                                                                                   gtTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCg                             >  2:1056061/1‑101 (MQ=255)                                                                                                                                   gtTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGt                        >  2:518783/1‑106 (MQ=255)                                                                                                                                   gtTAATCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGt                        <  1:518783/106‑1 (MQ=255)                                                                                                                                       |                                                                                                                            TTGGTGATGGAGGTGAATGCGTCGCCGGGGCTGGAAGGAATAGAAAAAACCACCGGTATCGACATCGCGGGTAAAATGATCCGCTGGATCGAACGCCACGCTACGACAGAATATTGCCTGAAAACGGGTGGTTAGTCGCAATCACATTACTGATCATGGTTTTGCCTGCGCTTTTTGCGTAAGCTGTGCCGGTCTTTTTATCGAAAGAGGTTGTACAAAATTATGACATCGCTGGTCGTTCCTGGTCTGGATACGCTG  >  NC_000913/892735‑892992

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