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%	87.9 / NA	27	N51T (AAC→ACC)	truD	tRNA pseudouridine(13) synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (18/9);  total (18/9)

CGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATG‑TTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGTTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCA  >  NC_000913/2870932‑2871395                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   ccctaaaccgaagtttgccctggggCCGCGCATCCTCCACCCCCGGGCAGCCTTCCAGTTGAAAGCCGCTCAGTGCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCCCTGTTCCGTAACAGCATT‑TTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGATTTTGCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCa                                                                                                                                                                                                                       <  2:323981/227‑1 (MQ=255)   ccTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATG‑TTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCAcc                                                                                                                                                                                                                     >  2:60840/1‑251 (MQ=255)       aaCGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATG‑TTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTcc                                                                                                                                                                                                                 >  2:280924/1‑250 (MQ=255)               ttGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGAGCATAACCACTGTTCCGTAACAGCATG‑TTTGTATTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGCTTCTAACCAGCATATGCTCACCTTCACCATc                                                                                                                                                                                                            >  2:216764/1‑248 (MQ=255)                                aTACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATG‑TTTGTCTTTTTGCCCAGCGAAGCTGACTTCCCGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCtt                                                                                                                                                                                        >  1:354053/1‑251 (MQ=255)                                          aCCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTCACAGCATG‑TTTGTATTTTTGCCCAGCGAAGCTGACTTCGCGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCTGCCGGTTTTGAGGCTTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAGGCCCAAATCTTCcaccaccac                                                                                                                                                                              >  1:399396/1‑251 (MQ=255)                                                  gCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATG‑TTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCAGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCtt                                                                                                                                                                      >  1:248516/1‑251 (MQ=255)                                                             aaaaGGCGCTAAGTGCGGGCATTTCCTGCCCCGGCACGCGAGCGAATAACCACTGTCCCTTAACAGCATG‑TTTGTCTTTTTGCCCAGCGAAGCTGACTTCCCGGGCATGAATTTTCAGGAATTCCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGc                                                                                                                                                            <  1:352587/249‑1 (MQ=255)                                                              aaaGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACAAACGTTCCGTAACAGCATG‑TTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGATAACCTTCACCATCAGGATCCAAGCCCCAATATACCACCCCCACGAAGTCTTTCggctctcctt                                                                                                                                                          >  1:248196/1‑243 (MQ=255)                                                                   cgcTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATG‑TTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTcagc                                                                                                                                                     >  2:336974/1‑251 (MQ=255)                                                                            cGGGCATTTCCTTGCTGGGCACGCGACCGCATAACCACGTTTCCGTAACAG‑ATGGTTTGTTTTTTTGCCCAGCGAAGCTGACTACGCGGGAATGATTTTTCAGGAATTTGGCCAGTACATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGt                                                                                                                                             <  1:72408/250‑1 (MQ=255)                                                                            cGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATG‑TTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTg                                                                                                                                            >  1:110942/1‑251 (MQ=255)                                                                                                    gAGCGCATAACCACTGTACAGTAACAGCATG‑TTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATAGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAACGACAAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCATTGCGGTTTACCGTGGAGGt                                                                                                                      >  2:365496/1‑249 (MQ=255)                                                                                                            aaCCACTGTTCCGTAACAGCATG‑TTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCAAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCAGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGATTACCGTGGAGGTAAGTGAGAtt                                                                                                            >  2:272305/1‑251 (MQ=255)                                                                                                              ccccctcgccgttacaaccatttttGTTTTTTTGCCGAGCGAAGCGGACTTCACGGCCAGAAATGTTCAGGATTTTCGCATGTCCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCTCCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTa                                                                                                           <  2:115929/229‑1 (MQ=255)                                                                                                                   gTTCCGTAACAGCATG‑TTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAAc                                                                                                     >  1:8313/1‑251 (MQ=255)                                                                                                                                          tttttGCCCATCGAATCTGACTTCACGGGCATGAATTTTCAGGAATTTCCCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTaa                                                                                <  1:309657/250‑1 (MQ=255)                                                                                                                                                      gAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGCGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCAc                                                                   >  1:197184/1‑251 (MQ=255)                                                                                                                                                                    gggCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCa                                                      >  2:278481/1‑250 (MQ=255)                                                                                                                                                                           aattttgaggatgttcgcacgtgactCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCGCCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGTTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTGTCAAACTCAATCATGTTGTTGCCTTAATGAGTAGCGCCACCGCTGCACAGTCAATCCCTTc                                              <  2:376860/226‑1 (MQ=255)                                                                                                                                                                            aTTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTcc                                             >  1:46014/1‑251 (MQ=255)                                                                                                                                                                                               gTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTATAACCAGCATATGCTCACCTTCACCATCAGGCTCAAAGCACAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATccc                           >  1:379710/1‑250 (MQ=255)                                                                                                                                                                                                                 gggTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTgg        >  1:152708/1‑251 (MQ=255)                                                                                                                                                                                                                   ggaTTGCAGCCGGTTTTGAGGATTATAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAgc                                                                        >  2:327233/3‑185 (MQ=255)                                                                                                                                                                                                                   gTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAgc                                                                        <  1:327233/185‑1 (MQ=255)                                                                                                                                                                                                                     aTTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCttt    <  1:387708/251‑1 (MQ=255)                                                                                                                                                                                                                        gCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCa  <  1:109813/250‑1 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   CGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATG‑TTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGTTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCA  >  NC_000913/2870932‑2871395

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