Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,811,671	A→G	S2P (TCC→CCC)	rpmB ←	50S ribosomal subunit protein L28

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,811,671	0	A	G	100.0%	105.4 / NA	31	S2P (TCC→CCC)	rpmB	50S ribosomal subunit protein L28
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/17);  total (14/17)

CTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTCAGTGCGTGGGAACGGTTGTTACCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGACATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGGTATTGGCGCCTCGTCAGGCTTTACAGCCCGGTCATCGCAGTTCTATGTGAACTCTCGATTGCCAGGCCCAAATGCCAAACC  >  NC_000913/3811546‑3811798                                                                                                                              |                                                                                                                                ctcGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTCAGTGCGTGGGAACGGTTGTTACCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCt                                                                                                                          >  1:416391/1‑133 (MQ=255)         cAACCCAGAAACGGAGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTCAGTGCGTGGGAACGGTTGTTACCGGTCACCGGACGGTTGCCAGTAACTGGGCCGACTCGGGTCATGATTATTCTCCa                                                                                                                  >  2:228521/1‑133 (MQ=255)                    cGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTCAGTGCGTGGGAACGGTTGTTACCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTGTCCAAAAATCAAAtt                                                                                                       >  1:160057/1‑133 (MQ=255)                      gTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTCAGTGCGTGGGAACGGTTGTTACCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTGTATGCTGCAAAAAGCAAATTAg                                                                                                     >  1:149513/1‑133 (MQ=255)                      gTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTCAGTGCGTGGGAACGGTTGTTACCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGAGGCATGTCTATTCTCCAAAAAACAAATTAg                                                                                                     >  1:225119/1‑133 (MQ=255)                         agagTGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTCAGTGCGTGGGAACGGTTGTTACCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTc                                                                                                  <  2:32062/133‑1 (MQ=255)                             tGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTCAGTGCGTGGGAACGGTTGTTACCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGc                                                                                              >  1:427407/1‑133 (MQ=255)                              gCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTCAGTGCGTGGGAACGGTTGTTACCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAgat                                                                                             <  1:286936/133‑3 (MQ=255)                                   ttCGGCAGGAAACGGCGTTTAGTCGCGTTCAGTGCGTGGGAACGGTTGTTACCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTa                                                                                        <  2:70212/133‑1 (MQ=255)                                   ttCGGCAGGAAACGGCGTTTAGTCGCGTTCAGTGCGTGGGAACGGTTGTTACCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTa                                                                                        >  2:260696/1‑133 (MQ=255)                                     cGGCAGGAAACGGCGTTTAGTCGCGTTCAGTGCGTGGGAACGGTTGTTACCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATg                                                                                      <  1:41734/133‑1 (MQ=255)                                                ggCGTTTAGTAGCGTTCAGTGCGTTGGAGCGGGTGTTACCGCTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGGTATTGgcg                                                                           <  2:225119/133‑1 (MQ=255)                                                      tAGTCGCGTCCAGTGGGTGGGAACGGGTGTTAGCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGGTATTGGCGCCTCGt                                                                     <  2:226948/133‑1 (MQ=255)                                                       aGTCGCGTTCGGTGCTTGGGAACGGTTGTTACGGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGGTATTGGCGCCTCGTc                                                                    <  2:177322/133‑1 (MQ=255)                                                       aGTCGCGTTCAGTGCGTGGGAACGGTTGTTACCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCCTATGGGGTATTGGCGCCTCGTc                                                                    >  1:226948/1‑133 (MQ=255)                                                        gTCGCGTTCAGTGCGTGGGAACGGTTGTTACCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGGTATTGGCGCCTCGTCa                                                                   <  2:427407/133‑1 (MQ=255)                                                                     cGTGGGAACGGTTGTTACCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGGTATTGGCGCCTCGTCAGGCTTTACAGccc                                                      <  1:260696/133‑1 (MQ=255)                                                                      gTGGGAACGGTTGTTACCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGGTATTGGCGCCTCGTCAGGCTTTACAGCCCg                                                     <  2:261917/133‑1 (MQ=255)                                                                       tGGGAACGGTTGTTACCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGGTATTGGCGCCTCGTCAGGCTTTACAGCCCgg                                                    >  1:451754/1‑133 (MQ=255)                                                                        gggAACGGTTGTTACCGGCCACGGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGGTATTGGCGCCTCGTCAGGCTTTACAGCCCGGt                                                   <  1:182306/133‑1 (MQ=255)                                                                         ggAACGGTTGTTACCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGGTATTGGCGCCTCGTCAGGCTTTACAGCCCGGTc                                                  >  1:378973/1‑133 (MQ=255)                                                                             cGGTTGTTACCGGTCGCCGGACGCGTGCCAGTTACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGGTATTGGCGCCTCGTCAGGCTTTACAGCCCGGTCATCg                                              <  2:160057/133‑1 (MQ=255)                                                                             cGGTTGTTACCGGTCCCCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCGTGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGGTATTGGCGCCTAGTCAGGCTTTACAGCCCGGTCATCg                                              <  1:282031/133‑1 (MQ=255)                                                                             cGGTTGTTACCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGGTATTGGCGCCTCGTCAGGCTTTACAGCCCGGTCATCg                                              >  1:101712/1‑133 (MQ=255)                                                                             cGGTTGTTACCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGGTATTGGCGCCTCGTCAGGCTTTACAGCCCGGTCATCg                                              <  2:255842/133‑1 (MQ=255)                                                                             cGGTTGTTACCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGGTATTGGCGCCTCGTCAGGCTTTACAGCCCGGTCATCg                                              <  1:513934/133‑1 (MQ=255)                                                                               gttgttACCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGGTATTGGCGCCTCGTCAGGCTTTACAGCCCGGTCATCGCa                                            >  1:367397/1‑133 (MQ=255)                                                                                          tCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGTTATTGGCGCCTCGCCAGGCTTTACAGCCCGGTCATCGCAGTTCTATGTGa                                 >  1:486821/1‑133 (MQ=255)                                                                                          tCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGGTATTGGCGCCTCGTCAGGCTTTACAGCCCGGTCATCGCAGTTCTATGTGa                                 >  2:363489/1‑133 (MQ=255)                                                                                                 aCGCTTGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGTTATTGGCGCCTCGCCAGGCTTTACAGCCCGGTCATCGCAGTTCTATGTGAACTCTCg                          <  2:486821/133‑1 (MQ=255)                                                                                                     ttGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGGTATTGGCGCCTCGTCAGGCTTTACAGCCCGGTCATCGCAGTTCTATGTGAACTCTCGATTg                      <  2:369092/133‑1 (MQ=255)                                                                                                     ttGCCAGTAACTTGGCAGACTCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGGTATTGGCGCCTCGTCAGGCTTTACAGCCCGGTCATCGCAGTTCTATGTGAACTCTCGATTg                      >  1:399062/1‑133 (MQ=255)                                                                                                                         tCGGGGCATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGGTATTGGCGCCTCGTCAGGCTTTACAGCCCGGTCATCGCAGTTCTATGTGAACTCTCGATTGCCAGGCCCAAATGCCAAAcc  >  1:488630/1‑133 (MQ=255)                                                                                                                              |                                                                                                                                CTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTCAGTGCGTGGGAACGGTTGTTACCGGTCACCGGACGCTTGCCAGTAACTTGGCAGACTCGGGACATGTCTATTCTCCAAAAATCAAATTAGCTCGAGCTTCGTATGGGGTATTGGCGCCTCGTCAGGCTTTACAGCCCGGTCATCGCAGTTCTATGTGAACTCTCGATTGCCAGGCCCAAATGCCAAACC  >  NC_000913/3811546‑3811798

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