Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,448,508	G→A	intergenic (‑68/‑98)	prmB ← / → smrB	50S ribosomal subunit protein L3 N(5)‑glutamine methyltransferase/putative endonuclease SmrB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,448,508	0	G	A	100.0%	74.4 / NA	27	intergenic (‑68/‑98)	prmB/smrB	50S ribosomal subunit protein L3 N(5)‑glutamine methyltransferase/putative endonuclease SmrB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/16);  total (11/16)

GACCAGCGCAACATGTCCTGAATGGTTTGCAGCTCATTTACTGCTTCATCAACGAAAATTTTATCCACGTATTCCTCCAGGGCATGATCGCAATAATTTCGGCGGCTAGTTTGCCATGAAGATGACGATAAATCAGCATTCACGCGCGGTGAGTGAGGAAAAATACGTTTAAAACGATCGATTGCGCTACGAGTCGGGTAAACTGTAGGAAAATTAGAAATAGAGACAGATAAATGAAAAAGAAAACAACACTTAGCGAGGAGGACCAGGCT  >  NC_000913/2448373‑2448644                                                                                                                                        |                                                                                                                                         gACCAGCGCAACATGTCCTGAATGGTTTGCAGCTCATTTACTGCTTCATCAACGAAAATTTTATCCACGTATTCCTCCAGGGCATGATCGCAATAATTTCGGCGGCTAGTTTGCCATGAAGATGACGATAAATCAACATTCAc                                                                                                                                   >  1:566074/1‑143 (MQ=255)              tGTCCTGAATGGTTTGCAGCTCATTTACTGCTTCATCAACGAAAATTTTATCCACGTATTCCTCCAGGGCATGATCGCAATAATTTCGGCGGCTAGTTTGCCATGAAGATGACGATAAATCAACATTCACGCGCGgtgagtga                                                                                                                      <  1:228539/143‑1 (MQ=255)                  cTGAATGGTTTGCAGCTCATTTACTGCTTCATCAACGAAAATTTTATCCACGTATTCCTCCAGGGCATGATCGCAATAATTTCGTCGGCTATGTTGCCATGAAGATGACGATAAATCAACATTCACGCGCGGTGAGTGAGGaa                                                                                                                  <  2:322084/143‑1 (MQ=255)                     aaTGGTTTGCAGCTCATTTACTGCTTCATCAACGAAAATTTTATCCACGTATTCCTCCAGGGCATGATCGCAATAATTTCGGCGGCTAGTTTGCCATGAAGATGACGATAAATCAACATTCACGCGCGGTGAGTGAGGaaaaa                                                                                                               <  1:322198/143‑1 (MQ=255)                      aTGGTTTGCAGCTCATTTACTGCTTCATCAACGAAAATTTTATCCACGTATTCCTCCAGGGCATGATCGCAATAATTTCGGCGGCTAGTTTGCCATGAAGATGACGATAAATCAACATGCACGCGCGGTGAGTGAGGAAAAAt                                                                                                              <  2:243381/143‑1 (MQ=255)                        gttttGCAGCTAATTAACGTCTTCATCAAAGAAAATTGTAGCCACGTATTCCTCCAGGGCATGATCGCAAGAATTTCGGCGGCTAGTTTGCCATGAAGATGACGATAAATCAACATTCACGCGCGGTGAGTGAGGAAAAATAc                                                                                                            <  1:286315/141‑1 (MQ=255)                                     tgtACGGCTGCAGCAACGGAAAGGTTACCCACGAAATCCGCCAGGGCATGATCGCAATAATTTCGGCGGCTAGTCTGCAATGAAGATGACGATAAATCAACATTCACGCGCGGTGAGTGAGGAAAAATACGTTTAAAAcgatc                                                                                               <  1:653480/141‑1 (MQ=255)                                            cTTCATCAACGAAAATTTTATCCACGTATTCCTCCAGGGCATGATCGCAATAATTTCGGCGGCTAGTTTGCCATGAAGATGACGATAAATCAACATTCACGCGCGGTGAGTGAGGAAAAATACGTTTAAAACGATCGATTgcg                                                                                        <  1:185459/143‑1 (MQ=255)                                                 tcaACGGAAATTTTATCCGCGTACTCCTCCAGGGCATGATCGCAATAATTTCGGCGGCTAGTTTGCCATGAAGATGACGATAAATCAACATTCACGCGCGGTGAGTGAGGAAAAATACGTTTAAAACGATCGATTGCGCTACg                                                                                   <  2:108344/143‑1 (MQ=255)                                                      gAAAATTTTATCCACGTATTCCTCCAGGGCATGATCGCAATAATTTCGGCGGCTAGTTTGCCATGAAGATGACGATAAATCAACATTCACGCGCGGTGAGTGAGGAAAAAGACGGTAAAAACGATCGAGTGGGCTACGAGACg                                                                              >  2:103575/1‑143 (MQ=255)                                                       aaaaTTTTATCCACGTATTCCTCCAGGGGATGATCGCAATAATTTCGGCGGCTAGTTTGCCATGAAGATGACGATAAATCAACATTCACGCGCGGTGAGTGAGGAAAAATACGTTTAAAACGATCGATTGCGCTACGAGTCgg                                                                             <  2:500078/143‑1 (MQ=255)                                                          aTTTGATCCACGTATTCCTCCAGGGCATGATCGCAATAATTTCGGCGGCTAGTTTGCCATGAAGATGACGATAAATCAACATTCACGCGCGGTGAGTGAGGAAAAATACGTTTAAAACGATCGATTGCGCTACGAGTCGGGTa                                                                          <  1:257570/143‑1 (MQ=255)                                                                  cACGTATTCCTCCAGGGCATGATCGCAATAATTTCGGCGGCTAGTTTGCCATGAAGATGACGATAAATCAACATTCACGCGCGGTGAGTGAGGAAAAATACGTTTAAAACGATCGATTGCGCTACGAGTCGGGGAAAccggag                                                                  >  2:569251/1‑138 (MQ=255)                                                                   aCGTATTCCTCCAGGGCATGATCGCAATAATTTCGGCGGCTAGTTTGCCATGAAGATGACGATAAATCAACATTCACGCGCGGTGAGTGAGGAAAAATACGTTTAAAACGATCGATTGCGCTACGAGACGGGTAAACTGTAgg                                                                 >  1:598348/1‑143 (MQ=255)                                                                      tATTCCTCCAGGGCATGATCGCAATAATTTCGGCGGCTAGTTCGCCATGAAGATGACGATAAATCAACATTCACGCGCGGTGAGTGAGGAAAAATACGTTTAAAACGATCGATTGCGCTACGAGTCGGGTAAACTGTAGGaaa                                                              <  1:388878/143‑1 (MQ=255)                                                                        ttCCGCCAGGGCATGATCGCAATAATTTCGGCGGCTAGTTTGCCATGAAGATGACGATAAATCAACATTCACGCGCGGTGAGTGAGGAAAAATACGTTTAAAACGATCGATTGCGCTACGAGTCGGGTAAACTGTAGGAAAAt                                                            <  1:18282/143‑1 (MQ=255)                                                                         tcctccAGGGCATGATCGCAATAATTTCGGCGGCTAGTTTGCCATGAAGATGACGATAAATCAACATTCACGCGCGGTGAGTGAGGAAAAATACGTTTAAAACGATCGATTGCGCTACGAGTCGGGTAAACTGTAGGAAAAtt                                                           >  1:660374/1‑143 (MQ=255)                                                                               aGGGCATGATCGCAATAATTTCGGCGGCTAGTTTGCCATGAAGATGACGATAAATCAACATTCACGCGCGGTGAGTGAGGAAAAATACGTTTAAAACGATCGATTGCGCTCCGAGTCGGGTAAACTGTAGGAAAATTAGAAAt                                                     >  1:681077/1‑143 (MQ=255)                                                                                              taatTTCGGCGGCTAGTTTGCCATGAAGATGACGATAAATCAACATTCACGCGCGGTGAGTGAGGAAAAATACGTTTAAAACGATCGATTGAGCAATGAGGCGGGTAAGCTGTAGGAAAATTAGAAATAGAGACAGATAAATg                                      >  1:689345/1‑143 (MQ=255)                                                                                                   tCGGCGGCTAGTTTGCCAGGAAGATGACGATAAATCAACATTCACGCGCGGTGAGTGAGGAAAAATACGTTTAAAACGATCGATTGCGCTACGAGTCGGGTAAACTGTAGGAAAATTAGAAATAGAGACAGATAAATGaaaaa                                 <  1:441936/143‑1 (MQ=255)                                                                                                       cggcTAGTTTGCCATGAAGATGACGATAAATCAACATTCACGCGCGGTGAGTGAGGAAAAATACGTTTAAAACGATCGATTGCGCTACGAGTCGGGTAAACTGTAGGAACATTAGAAATAGAGAGAGAAGAATAAAAAAGaaa                             >  2:463006/1‑143 (MQ=255)                                                                                                              tttGCCATGAAGATGACGATAAATCAACATTCACGCGCGGTGAGTGAGGAAAAATACGTTTAAAACGATCGATTGCGCTACGAGTCGGGTAAACTGTAGGAAAATTAGAAATAGAGACAGATAAATGAAAAAGAAAACAacac                      >  1:679762/1‑143 (MQ=255)                                                                                                                 gCCATGAAGATGACGATAAATCAACATTCACGCGCGGTGAGTGAGGAAAAATACGTTTAAAACGATCGGTTGCGCTACGAGTCGGGTAAACTGTAGGAAAATTAGAAATAGAGACAGATAAATGAAAAAGAAAACAACACTTa                   <  1:380984/143‑1 (MQ=255)                                                                                                                  acaagaagaagaaGATAAATCAACATTCACGCGAGGAGAAAGAGGAAAAATACGATTAAAACGAACGATTGCGATACGAGTCGAGTAAAATGTAGAAAAATTAGAAATAGAGACAGATAAATGAAAAAGAAAACAACACTTAg                  <  2:298672/139‑1 (MQ=255)                                                                                                                     tGAAGATGACGATAAATCAACATTCACGCGCGGTGAGTGAGGAAAAATACGTTTAAAACGATCGATTGCGCTACGAGTCGGGTAAACTGTAGGAAAATTAGAAATAGAGACAGATAAATGAAAAAGAAAACAACACTGAGCga               >  2:687151/1‑143 (MQ=255)                                                                                                                               gATAAATCAACATTCACGCGCGGTGAGTGAGGAAAAATACGTTTAAAACGATCGATTGCGCTACGAGTCGGGTAAACTGTAGGAAAATTAGAAATAGAGACAGATAAATGAAACAGAAAACAACACTTCTCGAGGAGGAcccg     >  1:2477/1‑141 (MQ=255)                                                                                                                                  aaaTCAACATCCAGTCGCGGTGAGTGAGGAAAAATACGTTTAAAACGATCGATTGCGCTACGAGTCGGGTACACTGTAAGAACATTAGAAATAGATTCAGATAAATGAAAAAGAATACAACACTTAGCGAGGAGGACCAGGct  <  1:442029/143‑1 (MQ=255)                                                                                                                                        |                                                                                                                                         GACCAGCGCAACATGTCCTGAATGGTTTGCAGCTCATTTACTGCTTCATCAACGAAAATTTTATCCACGTATTCCTCCAGGGCATGATCGCAATAATTTCGGCGGCTAGTTTGCCATGAAGATGACGATAAATCAGCATTCACGCGCGGTGAGTGAGGAAAAATACGTTTAAAACGATCGATTGCGCTACGAGTCGGGTAAACTGTAGGAAAATTAGAAATAGAGACAGATAAATGAAAAAGAAAACAACACTTAGCGAGGAGGACCAGGCT  >  NC_000913/2448373‑2448644

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