Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,475,417	A→C	S221A (TCC→GCC)	unknown ←	tRNA cytosine(34) acetyltransferase (EC 2.3.1.193)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,475,417	0	A	C	100.0%	77.4 / NA	29	S221A (TCC→GCC)	unknown	tRNA cytosine(34) acetyltransferase (EC 2.3.1.193)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (10/19);  total (10/19)

CGGTCGGTTTTCGTCATGAATTTTAACAAGAAGCCCCTGCCCGTGCCTTCGTAGCTGTGAATGGTGGTCGTCAGCACAATACGTTTAAAGGCAGTCGTTAAATCAAACAGCGTATCCAACGGAATCATCGCCGCTTCGTCCACAAATAACCACGAATCCGCAAATTGCTGTGGATTTTCACCAACCTTTTGGCGCAGTTCATCAGGTGGCATAAAGGGGATTTCGACACCGGAAAATTCATTGAAAATCTTGACCGCACTTTTGTTCGGCGCCGTCAATAGAACGGGCTTATTTGCCGCTAATTGCTGTTTGGCGAACAAGCCCGCCAAAGCGGATTTTCCCCGCCCGCGCTTTGCCGTAACGATTAACATCTCCGCGTCAGATGTCGCCATTTGCCGCAATAAATGGGCTTGTTCCAAGGTCGGTTGGTGTTGTTCTGAAAGCGACTTTTCAATCTTAGGGA  >  1/1475185‑1475647                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                       cggtcggtGTTCGCTATGAATTTTAACAAGAAGCCCCTGCCCGTACCCTCGTAACTGTGAATGGTCGTGGTCAGCACAATGCGTTTAAAGACGGCAGTTAAAGCAAACAGCGTAGCCAACGGGATCATCGCCGCTTCATCCACAAATAACCATAAATCCGCAAATTGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAACCGGGATTTCCGCAGCGGCAAATTCATTGAAAAtat                                                                                                                                                                                                                       <  1:879924/250‑3 (MQ=255)   gtcggtGTTCGCTATGAATTTTAACAAGAAGCCCCTGCCCGTACCCTCGTAACTGTGAATGGTCGTGGTCAGCACAATGCGTTTAAAGACGGCAGTTAAAGCAAACAGCGTATCCAACGGGATCATCGCCGCTTCATCCACAAATAACCATAAATCCGCAAATTGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTg                                                                                                                                                                                                                     <  1:196735/250‑1 (MQ=255)   gtcggtGTTCGCTATGAATTTTAACAAGAAGCCCCTGCCCGTACCCTCGTAACTGTGAATGGTCGTGGTCAGCACAATGCGTTTAAAGACGGCAGTTAAAGCAAACAGCGTATCCAACGGGATCATCGCCGCTTCATCCACAAATAACCATAAATCCGCAAATTGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTg                                                                                                                                                                                                                     <  1:777083/250‑1 (MQ=255)   gtcggtGTTCGCTATGAATTTTAACAAGAAGCCCCTGCCCGTACCCTCGTAACTGTGAATGGTCGTGGTCAGCACAATGCGTTTAAAGACGGCAGTTAAAGCAAACAGCGTATCCAACGGGATCATCGCCGCTTCATCCACAAATAACCATAAATCCGCAAATTGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTg                                                                                                                                                                                                                     <  1:408353/250‑1 (MQ=255)   gtcggtGTTCGCTATGAATTTTAACAAGAAGCCCCTGCCCGTAACCTCGTAACTGTGAATGGTCGTGGTCAGCACAATGCGTTTAAAGACGGCAGTTAAAGCAAACAGCGTATCCAACGGGATCATCGCCGCTTCATCCACAAATAACCATAAATCCGCAAATTGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTg                                                                                                                                                                                                                     <  1:872472/250‑1 (MQ=255)    tcggtgttcgctaTGAATTTTAACAAGAAGCCCCTGCCCGTACCCTCGTAACTGTGAATGGTCGTGGTCAGCACAATGCGTTTAAAGACGGCAGTTAAAGCAAACAGCGTATCCAACGGGATCATCGCCGCTTCATCCACAAATAACCATAAATCCGCAAATTGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGAc                                                                                                                                                                                                                   <  1:265432/239‑1 (MQ=255)       gtgttcgctaTGAATTTTAACAAGAAGCCCCTGCCCGTACCCTCGTAACTGTGAATGGTCGTGGTCAGCACAATGCGTTTAAAGACGGCAGTTAAAGCAAACAGCGTATCCAACGGGATCATCGCCGCTTCATCCACAAATAACCATAAATCCGCAAATTGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCg                                                                                                                                                                                                                 >  1:803534/10‑250 (MQ=255)       gtgttcgctaTGAATTTTAACAAGAAGCCCCTGCCCGTACCCTCGTAACTGTGAATGGTCGTGGTCAGCACAATGCGTTTAAAGACGGCAGTTAAAGCAAACAGCGTATCCAACGGGATCATCGCCGCTTCATCCACAAATAACCATAAATCCGCAAATTGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCg                                                                                                                                                                                                                 >  1:348820/10‑250 (MQ=255)                 tGAATTTTAACAAGAAGCCCCTGCCCGTACCCTCGTAACTGTGAATGGTCGTGGTCAGCACAATGCGTTTAAAGACGGCAGTTAAAGCAAACAGCGTATCCAACGGGATCATCGCCGCTTCATCCACAAATAACCATAAATCCGCAAATTGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTc                                                                                                                                                                                                      <  1:856567/251‑1 (MQ=255)                  gAATTTTAACAAGAAGCCCCTGCCCGTACCCTCGTAACTGTGAATGGTCGTGGTCAGCACAATGCGTTTAAAGACGGCAGTTAAAGCAAACAGCGTATCCAACGGGATCATCGCCGCTTCATCCACAAATAACCATAAATCCGCAAATTGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTca                                                                                                                                                                                                     >  1:234747/1‑250 (MQ=255)                                                                                                                            atcatcGCCGCTTCATCCACAAATAACCATAAATCCGCAAATTGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACAtctc                                                                                           <  1:694667/251‑1 (MQ=255)                                                                                                                                    cgcTTCATCCACAAATAACCATAAATCCGCAAATTGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACAtctctgcaccgg                                                                                   <  1:443875/251‑9 (MQ=255)                                                                                                                                     gcTTCATCCACAAATAACCATAAATCCGCAAATTGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACAtctctgcaccgga                                                                                  >  1:278547/1‑242 (MQ=255)                                                                                                                                            ccACAAATAACCATAAATCCGCAAATTGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGcc                                                                           <  1:68449/251‑1 (MQ=255)                                                                                                                                                 aaTAACCATAAATCCGCAAATTGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTg                                                                      >  1:834540/1‑251 (MQ=255)                                                                                                                                                               cGCAAATTGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATggg                                                        <  1:272634/251‑1 (MQ=255)                                                                                                                                                               cGCAAATTGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATggg                                                        <  1:897422/251‑1 (MQ=255)                                                                                                                                                                     ttGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGtt                                                  >  1:14588/1‑251 (MQ=255)                                                                                                                                                                           cggATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAgg                                            <  1:450363/250‑1 (MQ=255)                                                                                                                                                                                     ccaatctttcgacgcAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTgatg                                  <  1:536541/239‑4 (MQ=255)                                                                                                                                                                                      caatctttcgacgcAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGatgt                                 <  1:529655/240‑1 (MQ=255)                                                                                                                                                                                           tttcgacgcAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATgtagtt                             >  1:536029/7‑250 (MQ=255)                                                                                                                                                                                                acgcAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGaa                        <  1:832531/249‑1 (MQ=255)                                                                                                                                                                                                 cgcAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAg                      <  1:579153/251‑1 (MQ=255)                                                                                                                                                                                                     gTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGa                   >  1:860768/1‑250 (MQ=255)                                                                                                                                                                                                     gTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGa                   >  1:177452/1‑250 (MQ=255)                                                                                                                                                                                                      ttCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGAc                  <  1:77978/250‑1 (MQ=255)                                                                                                                                                                                                                 ataTAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTc                                                 <  1:418508/206‑1 (MQ=255)                                                                                                                                                                                                                     aaacgggATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGa  >  1:891812/5‑251 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                       CGGTCGGTTTTCGTCATGAATTTTAACAAGAAGCCCCTGCCCGTGCCTTCGTAGCTGTGAATGGTGGTCGTCAGCACAATACGTTTAAAGGCAGTCGTTAAATCAAACAGCGTATCCAACGGAATCATCGCCGCTTCGTCCACAAATAACCACGAATCCGCAAATTGCTGTGGATTTTCACCAACCTTTTGGCGCAGTTCATCAGGTGGCATAAAGGGGATTTCGACACCGGAAAATTCATTGAAAATCTTGACCGCACTTTTGTTCGGCGCCGTCAATAGAACGGGCTTATTTGCCGCTAATTGCTGTTTGGCGAACAAGCCCGCCAAAGCGGATTTTCCCCGCCCGCGCTTTGCCGTAACGATTAACATCTCCGCGTCAGATGTCGCCATTTGCCGCAATAAATGGGCTTGTTCCAAGGTCGGTTGGTGTTGTTCTGAAAGCGACTTTTCAATCTTAGGGA  >  1/1475185‑1475647

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