Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,475,400	G→C	P226P (CCC→CCG)	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,400	0	G	C	100.0%	77.7 / NA	29	P226P (CCC→CCG)	unknown	tRNA cytosine(34) acetyltransferase (EC 2.3.1.193)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (10/19);  total (10/19)

GAACAATTCAAAATGCCGCAAAGTGCGGTCGGTTTTCGTCATGAATTTTAACAAGAAGCCCCTGCCCGTGCCTTCGTAGCTGTGAATGGTGGTCGTCAGCACAATACGTTTAAAGGCAGTCGTTAAATCAAACAGCGTATCCAACGGAATCATCGCCGCTTCGTCCACAAATAACCACGAATCCGCAAATTGCTGTGGATTTTCACCAACCTTTTGGCGCAGTTCATCAGGTGGCATAAAGGGGATTTCGACACCGGAAAATTCATTGAAAATCTTGACCGCACTTTTGTTCGGCGCCGTCAATAGAACGGGCTTATTTGCCGCTAATTGCTGTTTGGCGAACAAGCCCGCCAAAGCGGATTTTCCCCGCCCGCGCTTTGCCGTAACGATTAACATCTCCGCGTCAGATGTCGCCATTTGCCGCAATAAATGGGCTTGTTCCAAGGTCGGTTGGTGTTGTTCTGAAAGCGAC  >  1/1475160‑1475631                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                        ctattgacggcgctgaacaaaaGTGCGGTCGGTGTTCGCTATGAATTTTAACAAGAAGCCCCTGCCCGTACCCTCGTAACTGTGAATGGTCGTGGTCAGCACAATGCGTTTAAAGACGGCAGTTAAAGCAAACAGCGTATCCAACGGGATCATCGCCGCTTCATCCACAAATAACCATAAATCCGCAAATTGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAAACGGGATTTccg                                                                                                                                                                                                                               >  1:401924/20‑249 (MQ=255)                          cggtcggtGTTCGCTATGAATTTTAACAAGAAGCCCCTGCCCGTACCCTCGTAACTGTGAATGGTCGTGGTCAGCACAATGCGTTTAAAGACGGCAGTTAAAGCAAACAGCGTAGCCAACGGGATCATCGCCGCTTCATCCACAAATAACCATAAATCCGCAAATTGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAACCGGGATTTCCGCAGCGGCAAATTCATTGAAAAtat                                                                                                                                                                                                       <  1:879924/250‑3 (MQ=255)                            gtcggtGTTCGCTATGAATTTTAACAAGAAGCCCCTGCCCGTACCCTCGTAACTGTGAATGGTCGTGGTCAGCACAATGCGTTTAAAGACGGCAGTTAAAGCAAACAGCGTATCCAACGGGATCATCGCCGCTTCATCCACAAATAACCATAAATCCGCAAATTGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTg                                                                                                                                                                                                     <  1:196735/250‑1 (MQ=255)                            gtcggtGTTCGCTATGAATTTTAACAAGAAGCCCCTGCCCGTACCCTCGTAACTGTGAATGGTCGTGGTCAGCACAATGCGTTTAAAGACGGCAGTTAAAGCAAACAGCGTATCCAACGGGATCATCGCCGCTTCATCCACAAATAACCATAAATCCGCAAATTGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTg                                                                                                                                                                                                     <  1:408353/250‑1 (MQ=255)                            gtcggtGTTCGCTATGAATTTTAACAAGAAGCCCCTGCCCGTACCCTCGTAACTGTGAATGGTCGTGGTCAGCACAATGCGTTTAAAGACGGCAGTTAAAGCAAACAGCGTATCCAACGGGATCATCGCCGCTTCATCCACAAATAACCATAAATCCGCAAATTGCTGCGGATTTTCACCAATCTTTCGACGCAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTg                                                                                                                                                                                                     <  1:777083/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)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                        GAACAATTCAAAATGCCGCAAAGTGCGGTCGGTTTTCGTCATGAATTTTAACAAGAAGCCCCTGCCCGTGCCTTCGTAGCTGTGAATGGTGGTCGTCAGCACAATACGTTTAAAGGCAGTCGTTAAATCAAACAGCGTATCCAACGGAATCATCGCCGCTTCGTCCACAAATAACCACGAATCCGCAAATTGCTGTGGATTTTCACCAACCTTTTGGCGCAGTTCATCAGGTGGCATAAAGGGGATTTCGACACCGGAAAATTCATTGAAAATCTTGACCGCACTTTTGTTCGGCGCCGTCAATAGAACGGGCTTATTTGCCGCTAATTGCTGTTTGGCGAACAAGCCCGCCAAAGCGGATTTTCCCCGCCCGCGCTTTGCCGTAACGATTAACATCTCCGCGTCAGATGTCGCCATTTGCCGCAATAAATGGGCTTGTTCCAAGGTCGGTTGGTGTTGTTCTGAAAGCGAC  >  1/1475160‑1475631

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: ‑