Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,871,153	T→G	N51T (AAC→ACC)	truD ←	tRNA pseudouridine(13) synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,871,153	0	T	G	100.0%	107.4 / NA	28	N51T (AAC→ACC)	truD	tRNA pseudouridine(13) synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (14/14);  total (14/14)

GCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGTTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAG  >  NC_000913/2870916‑2871386                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                          gctttccttttCAGCGCGCCTAAAGGCAGCTTGCTCTTGTGCCCCCCATACTCCAGCTCCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCTTTTCCTTGCCCGGCCCGCGAGCGCATAACCCCTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATTAATTTCAAGGAATTTCGCCAGGGCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCt                                                                                                                                                                                                                              <  1:217430/243‑1 (MQ=255)                      aaaCGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTc                                                                                                                                                                                                         >  1:129130/1‑251 (MQ=255)                                              gcATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATc                                                                                                                                                                                 >  2:42151/1‑251 (MQ=255)                                                        gCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTccaccacc                                                                                                                                                                       >  2:62524/1‑251 (MQ=255)                                                          aCCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCcaccaccac                                                                                                                                                                     >  1:29448/1‑251 (MQ=255)                                                                                                 aTTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGTTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTg                                                                                                                              >  1:122585/1‑251 (MQ=255)                                                                                                   ttCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCg                                                                                                                            >  1:236654/1‑251 (MQ=255)                                                                                                     ccTTGCCCGGCACGCGCGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGt                                                                                                                          >  2:149315/1‑251 (MQ=255)                                                                                                                aCGCGAGCGAATAACCACTGTTCCGGACCAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATAAATTTTCAGGAATTTCGCCAGTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAg                                                                                                               <  1:171436/251‑1 (MQ=255)                                                                                                                        gcATACCCACTGTTCCGTAACAGCATTTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTga                                                                                                       <  1:238407/251‑1 (MQ=255)                                                                                                                          aTAACCACTGTTCCGTAACAGCAGTTTTGTCTTTTTGCCCAGCAAAGCTGCCTTCACGGGCATGAATTTTCAGAAATTTCGCCAGTCCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTgaga                                                                                                     <  1:114078/251‑1 (MQ=255)                                                                                                                                  tGTTCCGTAACAGCATGCTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAACTCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCaa                                                                                              >  2:238009/1‑250 (MQ=255)                                                                                                                                   gttccttaacatcttttttGTCTTTTTGCCCAGCGAAGTTGTTTCCCCGGGCATGAATTTTCAGGAATTCCGCCAGTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCaa                                                                                              <  1:3728/233‑1 (MQ=255)                                                                                                                                   gTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAAc                                                                                            >  1:182265/1‑251 (MQ=255)                                                                                                                                              gCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGATTTTCGCCAGTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCAtttt                                                                                 <  1:196850/251‑1 (MQ=255)                                                                                                                                              gCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCAtttt                                                                                 >  1:130509/1‑251 (MQ=255)                                                                                                                                                      gtttttctGCCCAGCGAAGCTGACTTCCCGTGCATGAATTTTCAGGAATTTCGCCAGTGCA‑TACGCAACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTACACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGACTTGCGGTTTACCGTGGAGGGAAGTGAGATTATCAAACTCCATCATTTTGTTGCCtt                                                                         >  2:212388/8‑251 (MQ=255)                                                                                                                                                       tCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCtt                                                                         >  1:98130/1‑250 (MQ=255)                                                                                                                                                          ttttGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATGTTGTTGCCTTAAt                                                                      <  1:193082/250‑1 (MQ=255)                                                                                                                                                           tttGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATg                                                                     <  1:66434/250‑1 (MQ=255)                                                                                                                                                                           gACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCtt                                                     <  1:65143/250‑1 (MQ=255)                                                                                                                                                                                     gCATGAATTTTCAGGAATTTCGCCAGTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATc                                          >  1:72391/1‑251 (MQ=255)                                                                                                                                                                                           tttttCAGGAATTTCGCC‑GTGCAATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTc                                     <  2:126677/249‑1 (MQ=255)                                                                                                                                                                                            tttgaaGGAATTGCGCCGGTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGACTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTcc                                    <  2:169216/245‑1 (MQ=255)                                                                                                                                                                                                      tttCGCCAGTGCA‑TCCGCCACAAAACGAGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGt                         >  1:57050/1‑251 (MQ=255)                                                                                                                                                                                                              gTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCa                 <  1:106439/251‑1 (MQ=255)                                                                                                                                                                                                                            aaaaCGGGTATTGCCGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAATCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACGCAATCAGTTTGTTGTCTTAATGAGTAGCGTCCCGGTTTTTCAGGCAACCCTTTCCCCACGTCCTGTTAAGCCCAGTTTTTCCgtagg    <  2:239641/251‑2 (MQ=255)                                                                                                                                                                                                                              aaCGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCGGTTTTTCCGtagtag  <  2:72391/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                          GCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCA‑TCCGCCACAAAACGGGTATTGCAGCCGTTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAG  >  NC_000913/2870916‑2871386

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