Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	943,761	T→C	Y116Y (TAT→TAC)	dmsB →	dimethyl sulfoxide reductase subunit B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	943,761	0	T	C	100.0%	102.5 / NA	30	Y116Y (TAT→TAC)	dmsB	dimethyl sulfoxide reductase subunit B
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/14);  total (16/14)

CGGCTTGTACTAAAGTCTGCCCGAGCGGTGCGATGCATAAACGTGAAGATGGTTTTGTTGTGGTCGATGAAGATGTGTGCATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACG  >  NC_000913/943631‑943882                                                                                                                                   |                                                                                                                          cGGCTTGTACTAAAGTCTGCCCGAGCGGTGCGATGCATAAACGTGAAGATGGTTTTGTTGTGGTCGATGAAGATATGGGCATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATACaacgaaact                                                                                                                  >  2:544500/1‑139 (MQ=255)   gCTTGTACTAAAGTCTGCCCGAGCGGTGCGATGCATAAACGTGAAGATGGTTTTGTTGTGGTCGATGAAGATGTGTGCATTGGCTGCCGCTACTGCCATATGGCGTGCCCGTATGGCCCACCGCAATACAacgaaacgaa                                                                                                                >  2:297940/1‑140 (MQ=255)          cTAAAGTCTGCCCGAGCGGTGCGATGCATAAACGTGAAGATGGTTTTGTTGTGGTCGATGAAGATGTGTGCATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATACAACGAAACGAAAGGCCat                                                                                                         <  1:520163/140‑1 (MQ=255)                  tGCCCGAGCGGTGCGATGCATAAACGTGAAGATGGTTTTGTTGTGGTCGATGAAGATGTGTGCATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATACAACGAAACGAAAGGCCATATGACCaa                                                                                                 >  1:632688/1‑140 (MQ=255)                  tGCCCGAGCGGTGCGATGCATAAACGTGAAGATGGTTTTGTTGTGGTCGATGAAGATGTGTGCATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATACAACG‑AACGAAAGGCCATATGACCaaa                                                                                                >  1:853527/1‑140 (MQ=255)                    cccGAGCGGTGCGATGCATAAACGTGAAGATGGTTTTGTTGTGGTCGATGAAGATGTGTGCATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATACAACGAAACGAAAGGCCATATGACCAAAt                                                                                               >  2:767949/1‑140 (MQ=255)                       gAGCGGTGCGATGCATAAACGTGAAGATGGTTTTGTTGTGGTCGATGAAGATGTGTGCATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATACAACG‑AACGAAAGGCCATATGACCAAATGCGa                                                                                           <  2:853527/140‑1 (MQ=255)                         gCGGTGCGATGCATAAACGTGAAGATGGTTTTGTTGTGGTCGATGAAGATGTGTGCATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCTATACAACGAAAGGAGAGGCCAGATGACCAAATGCGAt                                                                                          >  1:254966/1‑140 (MQ=255)                            gTGCGATGCATAAACGTTAAGACGGTTTTGTCGTGGTCGCTGAAGATGTGTGCATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATACAACGAAACGAAAGGCCATATGACCAAATGCGATGgt                                                                                       <  2:254966/140‑1 (MQ=255)                                    cATAAACGTGAAGATGGCTTTGTTGTGGTCGTTGACGATTTGTGCATTCGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATACAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGa                                                                               <  2:250250/140‑1 (MQ=255)                                      tAAACGTGAAGATGGTTTTGTTGTGGTCGATGAAGATGTGTGCATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATACAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGAcc                                                                             <  2:470564/140‑1 (MQ=255)                                       aaaCGTGAAGATGGTTTTGTTGTGGTCGATGAAGATGTGTGCATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATACAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCg                                                                            >  2:689752/1‑140 (MQ=255)                                            tGAAGATGGTTTTGTTGTGGTCGATGAAGATGTGTGCATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATACAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTg                                                                       <  1:498835/140‑1 (MQ=255)                                             gAAGATGGTTTTGTTGTGGTCGATGAAGATGTGTGCATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATACAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGc                                                                      <  1:689752/140‑1 (MQ=255)                                              aaGATGGTTTTGTTGTGGTCGATGAAGATGTGTGCATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATACAACGAGACGAAAGGTCATGTGAACAAATGCGATGGTGGTTATGACCGTGTTGCg                                                                     >  2:271122/1‑140 (MQ=255)                                              aaGATGGTTTTGTTGTGGTCGATGAAGATGTGTGCATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATACAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCg                                                                     >  2:757085/1‑140 (MQ=255)                                                       ttgttgTGGTCGATGAAGATGTGTGCATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATACAACGAAACGAAAGGCCATATGACCAGATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTaaa                                                            <  1:704357/140‑1 (MQ=255)                                                       ttgttgTGGTCGATGAAGATGTGTGCATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATACAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTaaa                                                            >  2:211909/1‑140 (MQ=255)                                                       ttgttgTGGTCGATGAAGATGTGTGCATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATACAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTaaa                                                            <  1:211909/140‑1 (MQ=255)                                                              ggTCGATGAAGATGTGTGCATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATACAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAg                                                         >  1:723834/1‑136 (MQ=255)                                                              ggTCGATGAAGATGTGTGCATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATACAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAg                                                         <  2:723834/136‑1 (MQ=255)                                                               gTCGATGAAGATGTGTGCATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATACAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGAt                                                    <  1:626856/140‑1 (MQ=255)                                                                                cATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATACAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc                                   >  1:557825/1‑140 (MQ=255)                                                                                   tGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATACAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGc                                >  2:820783/1‑140 (MQ=255)                                                                                             tACTGCCATATGGCTTGCCCGTATGCCGCACCGCAATCCAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCt                      <  1:297940/140‑1 (MQ=255)                                                                                                         gCTTGCCCGTATGGCGCTCCGCGATACAACGCAACCAAACGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGccc          <  2:270778/140‑1 (MQ=255)                                                                                                            tGCCCGTATGGCGCACCGCAATACAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGCAAAAAGCCGATCTGTGTTCAATCTTGGCCGATGAGGGCGCGGGATTTCGGCCCTAt       >  1:281989/1‑140 (MQ=255)                                                                                                              cccGTATGGCGCACCGCAATACAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATcg     >  1:168609/1‑140 (MQ=255)                                                                                                                 gTATGGCGCACCGCAATACAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCGGTGTTGAAACTTATCCACTGCGGGCGCTGGATTTCGGCCCTATcgacg  >  2:243596/1‑140 (MQ=255)                                                                                                                   aTGGCGCACCGCAATACAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGccc          <  1:594039/130‑1 (MQ=255)                                                                                                                   aTGGCGCACCGCAATACAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGccc          >  2:594039/1‑130 (MQ=255)                                                                                                                                   |                                                                                                                          CGGCTTGTACTAAAGTCTGCCCGAGCGGTGCGATGCATAAACGTGAAGATGGTTTTGTTGTGGTCGATGAAGATGTGTGCATTGGCTGCCGCTACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACG  >  NC_000913/943631‑943882

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: ‑