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%	127.2 / NA	36	S14R (AGC→AGA)	dusA	tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (23/13);  total (23/13)

CCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGCGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGT  >  NC_000913/4261589‑4261844                                                                                                                          |                                                                                                                                       ccGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACt                                                                                                                      >  1:5317/1‑140 (MQ=255)  cGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTg                                                                                                                     >  1:276742/1‑140 (MQ=255)    cTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGAc                                                                                                                   <  1:287381/140‑1 (MQ=255)                         gATGCTGAATAACATAACCCCATGATATAACGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGttt                                                                                              >  2:380761/1‑140 (MQ=255)                          aTGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTa                                                                                             >  1:1529/1‑140 (MQ=255)                                          ccccATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCAc                                                                                    >  1:165142/1‑133 (MQ=255)                                          ccccATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCAc                                                                                    <  2:165142/133‑1 (MQ=255)                                          ccccATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCt                                                                             >  1:266216/1‑140 (MQ=255)                                           cccATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTc                                                                            >  1:164349/1‑140 (MQ=255)                                               tGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACt                                                                        >  1:392916/1‑140 (MQ=255)                                               tGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACt                                                                        >  2:233820/1‑140 (MQ=255)                                                    tatCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTggacg                                                                   >  2:370851/1‑140 (MQ=255)                                                     atCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTggacgg                                                                  <  2:5317/140‑1 (MQ=255)                                                     atCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTggacgg                                                                  <  1:127484/140‑1 (MQ=255)                                                     atCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTggacgg                                                                  >  1:273211/1‑140 (MQ=255)                                                     atCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTg                                                                       <  2:53821/135‑1 (MQ=255)                                                     atCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTg                                                                       >  1:53821/1‑135 (MQ=255)                                                         aTAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGgacag                                                              >  2:96575/1‑140 (MQ=255)                                                         aTAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGgacag                                                              >  1:183153/1‑140 (MQ=255)                                                                 atCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGcc                                                      <  1:4836/140‑1 (MQ=255)                                                                           ttGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCtt                                            >  1:38291/1‑140 (MQ=255)                                                                                aaTGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTc                                       >  1:344559/1‑140 (MQ=255)                                                                                   gCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGc                                    >  2:171040/1‑140 (MQ=255)                                                                                       ggTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCttt                                 >  1:387619/1‑139 (MQ=255)                                                                                       ggTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCttt                                 <  2:387619/139‑1 (MQ=255)                                                                                          aaTTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCg                             <  1:96575/140‑1 (MQ=255)                                                                                          aaTTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCg                             <  2:164349/140‑1 (MQ=255)                                                                                            ttCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCa                           >  1:361385/1‑140 (MQ=255)                                                                                                 aaaTGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACg                      <  1:380761/140‑1 (MQ=255)                                                                                                   aTGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGtt                    >  1:297945/1‑140 (MQ=255)                                                                                                   aTGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGtt                    >  1:28233/1‑140 (MQ=255)                                                                                                   aTGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGtt                    >  1:334097/1‑140 (MQ=255)                                                                                                    tGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTg                   <  2:276742/140‑1 (MQ=255)                                                                                                                accaaaccaGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGaa       <  1:233820/140‑1 (MQ=255)                                                                                                                accaaaccaGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGaa       <  2:297945/140‑1 (MQ=255)                                                                                                                     accaGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGt  >  1:5331/1‑140 (MQ=255)                                                                                                                          |                                                                                                                                       CCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGCGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGT  >  NC_000913/4261589‑4261844

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