Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,261,710	C→A	S14R (AGC→AGA)	dusA →	tRNA‑dihydrouridine synthase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,261,710	0	C	A	100.0%	93.2 / NA	28	S14R (AGC→AGA)	dusA	tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (16/12);  total (16/12)

GTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGCGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGT  >  NC_000913/4261573‑4261844                                                                                                                                          |                                                                                                                                       gttATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAgagc                                                                                                                                      >  1:100321/1‑140 (MQ=255)   tATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCaa                                                                                                                                    >  2:300849/1‑140 (MQ=255)    aTGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAAt                                                                                                                                   >  2:30566/1‑140 (MQ=255)               ccccGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGaaaaaa                                                                                                                        >  2:233862/1‑140 (MQ=255)                  cGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTg                                                                                                                     >  1:187861/1‑140 (MQ=255)                    cTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGAc                                                                                                                   <  1:302526/140‑1 (MQ=255)                       ttGAGCCTCCTGAATAGTGATGCTGAATAACTTAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGtt                                                                                                                <  2:298979/140‑1 (MQ=255)                               cctGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAg                                                                                                        >  2:85076/1‑140 (MQ=255)                               cctGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAg                                                                                                        >  2:310211/1‑140 (MQ=255)                               cctGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAg                                                                                                        >  1:38835/1‑140 (MQ=255)                                   aaTAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGt                                                                                                    <  2:63662/140‑1 (MQ=255)                                             cTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTgg                                                                                                     <  2:92394/129‑1 (MQ=255)                                             cTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTgg                                                                                                     >  1:92394/1‑129 (MQ=255)                                             cTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGtt                                                                                               >  2:262010/1‑135 (MQ=255)                                             cTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGtt                                                                                               <  1:262010/135‑1 (MQ=255)                                                  taacataacCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCa                                                                                     >  1:73227/1‑140 (MQ=255)                                                     cataacCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCa                                                                                  >  1:64627/1‑140 (MQ=255)                                                          ccccATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCt                                                                             >  1:189256/1‑140 (MQ=255)                                                               tGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACt                                                                        >  2:22858/1‑140 (MQ=255)                                                                  tatatCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTgga                                                                     <  1:298425/140‑1 (MQ=255)                                                                                       aCATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGTCATTGCCGCTAtt                                                >  1:33431/1‑140 (MQ=255)                                                                                           ttGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCtt                                            >  2:90215/1‑140 (MQ=255)                                                                                                               tGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATa                        <  2:187861/140‑1 (MQ=255)                                                                                                               tGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCAGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATa                        <  1:22802/140‑1 (MQ=255)                                                                                                                 aaaTGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACg                      <  1:300849/140‑1 (MQ=255)                                                                                                                 aaaTGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACg                      >  2:55165/1‑140 (MQ=255)                                                                                                                    tGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTg                   <  2:64627/140‑1 (MQ=255)                                                                                                                    tGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCCCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTg                   <  1:85076/140‑1 (MQ=255)                                                                                                                                     accaGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGt  <  2:73227/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                       GTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGCGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGT  >  NC_000913/4261573‑4261844

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