Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	531,186	A→G	L14P (CTG→CCG)	selU ←	tRNA 2‑selenouridine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	531,186	0	A	G	100.0%	103.7 / NA	29	L14P (CTG→CCG)	selU	tRNA 2‑selenouridine synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/18);  total (11/18)

TTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAA  >  NC_000913/531053‑531315                                                                                                                                      |                                                                                                                                  ttGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCGGCAggg                                                                                                                             <  1:840486/140‑1 (MQ=255)      gTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCGGCAGGGCACGa                                                                                                                        <  1:66408/140‑1 (MQ=255)      gTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCGGCAGGGCACGa                                                                                                                        <  1:530465/140‑1 (MQ=255)           taGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCGGCAGGGCACGATAGTc                                                                                                                   <  1:766507/140‑1 (MQ=255)           taGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCGGCAGGGCACGATAGTc                                                                                                                   <  1:755891/140‑1 (MQ=255)            aGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAAATCGGTAGGGGGGCGGACATCAAAAAAGGGGGGATCAGCAAACGGCAGGGGACGATAGTcc                                                                                                                  >  2:530465/1‑140 (MQ=255)                 gTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGGATTGCGCCGGGCTCAAACTCGATAGGGGCGCGAACATCAATAAAGGGGGTATCAGCAATCGGCAGGGGACGAAAGTCCTGTTc                                                                                                             >  2:152866/1‑140 (MQ=255)                        cggcggCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCgtgtgt                                                                                                      <  2:952709/140‑1 (MQ=255)                        cggcggCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCgtgtgt                                                                                                      <  1:949681/140‑1 (MQ=255)                           cggcgCGTTCATCGTTATTCATTAACGACAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATAGGCGTATCAGCAATCGGCAGGGGACGATAGTCCTGTTCCCTGTGtctc                                                                                                   >  2:222424/1‑140 (MQ=255)                                 gTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGc                                                                                             <  1:799436/140‑1 (MQ=255)                                 gTTCATCGTTATTCATTAACGACAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGc                                                                                             <  1:222424/140‑1 (MQ=255)                                  ttCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCa                                                                                            <  1:152866/140‑1 (MQ=255)                                  ttCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCa                                                                                            <  2:858533/140‑1 (MQ=255)                                   tCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCAt                                                                                           <  1:1220144/140‑1 (MQ=255)                                              cATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTg                                                                                <  2:82229/140‑1 (MQ=255)                                                    cGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTc                                                                          <  1:1246311/140‑1 (MQ=255)                                                    aggCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTc                                                                          <  1:786292/139‑1 (MQ=255)                                                                 gcggcgGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCGGCAGGGGACGATAGTCCTGTTCCGGGTGTCTCTCTTGCATAGGTAAAAAGGAAACTCAAATCAGCtttct                                                             >  2:786292/1‑137 (MQ=255)                                                                 gcggcgGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGCTTgcg                                                             <  1:79776/140‑1 (MQ=255)                                                                   ggcggGCATTGCGCCGTGCTCAAACTCGATAGGGGGGCGAACATCAATAATGGGCGGATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGGGTGTCTCTCTTTCATAGTTAAAAAGGAAACTCAAATCAGCTTtcgcg                                                           >  2:62234/1‑140 (MQ=255)                                                                                     cTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAAcg                                         >  1:23734/1‑140 (MQ=255)                                                                                             cGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGa                                 >  2:583703/1‑140 (MQ=255)                                                                                                      cgcgAACATCAATAATGGGCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGAtt                        >  2:1141670/1‑140 (MQ=255)                                                                                                      cgcgAACATCAATAATGGGCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGAtt                        >  2:1176084/1‑140 (MQ=255)                                                                                                        cgAACATCAATAATGGGCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGAtttt                      >  1:850274/1‑140 (MQ=255)                                                                                                                    atGGGCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTggg          <  1:62234/140‑1 (MQ=255)                                                                                                                       ggCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTaa                   <  1:1188251/128‑1 (MQ=255)                                                                                                                       ggCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTaa                   >  2:1188251/1‑128 (MQ=255)                                                                                                                       ggCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTg       >  2:391805/1‑140 (MQ=255)                                                                                                                       ggCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTg       >  1:141715/1‑140 (MQ=255)                                                                                                                       ggCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTg       >  1:444660/1‑140 (MQ=255)                                                                                                                       ggCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTg       >  2:871494/1‑140 (MQ=255)                                                                                                                       ggCGTATCAGCAATCGGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTg       >  1:613365/1‑140 (MQ=255)                                                                                                                            atcagcaatcggcaGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGaa  >  2:1242235/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                  TTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAA  >  NC_000913/531053‑531315

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