Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,474,488	T→C	G530G (GGA→GGG)	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,474,488	0	T	C	100.0%	154.2 / NA	57	G530G (GGA→GGG)	unknown	tRNA cytosine(34) acetyltransferase (EC 2.3.1.193)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/31);  total (26/31)

GGGCGTTTGCAATAGGCTGTCAATAAAGGAAAATCCCGCTCCCCGCTTAATTTCATCAAGCGCCGAATAGCGGGCAAGGCAAACGCCAACGGCAGGAAGAAATCGGCAAAATCCGCCAAGCATTGCCAATCGGTTTCATTGAGCGACCAATCCGTTTCGGTTCCGACAAATTGCGTCGCCAGCGGGTGCAACGAAAGGGGCAGATTACGTTTAAATTGTTGCTCCGCCTGTTGCACCAACGCCCGTCCGCGCTCCGTCAGCGGGCACAATGCCACTACGGAATAACAGCCGCTGCCTGCATCTTTGCCTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCCAGCAGCTCGGGGGTGTAACCGAAACTGACGGAAAGGAAATCCACCGCACGTTGTTGTTTAACTC‑GTTCGATTAAACGTTGCCCCAGTCCCCGCTGCTGCCAGTTGGGTTGCACGGCGATACGGGAAATGCGCAGTGA  >  1/1474243‑1474733                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       ggacgcttGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGcccgcccgtg                                                                                                                                                                                                                                                   >  1:890464/7‑249 (MQ=255)   acgcttGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGcccgcccgtgc                                                                                                                                                                                                                                                  <  1:46631/246‑4 (MQ=255)   acgcttGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGcccgcccgtgc                                                                                                                                                                                                                                                  <  1:253026/246‑4 (MQ=255)       ttGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGt                                                                                                                                                                                                                                             >  1:475604/1‑251 (MQ=255)        tGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGtt                                                                                                                                                                                                                                            >  1:551660/1‑250 (MQ=255)         gCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGtta                                                                                                                                                                                                                                           >  1:664681/1‑249 (MQ=255)                gCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCggaca                                                                                                                                                                                                                                    >  1:410682/1‑248 (MQ=255)                         aaaGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCAc                                                                                                                                                                                                                           <  1:332777/251‑1 (MQ=255)                           aGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCAcca                                                                                                                                                                                                                         >  1:616559/1‑249 (MQ=255)                             gAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCAccac                                                                                                                                                                                                                        <  1:214156/250‑4 (MQ=255)                             gAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCAccac                                                                                                                                                                                                                        <  1:457608/250‑4 (MQ=255)                                aatcctgttCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAGCGCCACCACGGaa                                                                                                                                                                                                                    <  1:728512/243‑1 (MQ=255)                                      gttCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAg                                                                                                                                                                                                              >  1:205980/3‑251 (MQ=255)                                       ttCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAg                                                                                                                                                                                                              <  1:23021/249‑1 (MQ=255)                                       ttCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAg                                                                                                                                                                                                              <  1:858142/249‑1 (MQ=255)                                       ttCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAg                                                                                                                                                                                                              <  1:548670/249‑1 (MQ=255)                                          cccGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGc                                                                                                                                                                                                          <  1:323976/251‑1 (MQ=255)                                             gCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCAc                                                                                                                                                                                                                           >  1:752862/1‑231 (MQ=255)                                                tactttCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGctgcccg                                                                                                                                                                                                    >  1:458958/4‑249 (MQ=255)                                                  ctttCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCa                                                                                                                                                                                                  >  1:88031/2‑251 (MQ=255)                                                                                    cGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAAc                                                                                                                                                                <  1:550499/251‑1 (MQ=255)                                                                                    cGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAAc                                                                                                                                                                <  1:139986/251‑1 (MQ=255)                                                                                          cggtaGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTGGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCa                                                                                                                                                            >  1:429372/5‑249 (MQ=255)                                                                                          cggtaGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCa                                                                                                                                                            >  1:639961/5‑249 (MQ=255)                                                                                          cggtaGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCa                                                                                                                                                            >  1:672984/5‑249 (MQ=255)                                                                                            gtaGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACttt                                                                                                                                                        >  1:35151/3‑251 (MQ=255)                                                                                                gaagaaATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGc                                                                                                                                                     >  1:383121/1‑250 (MQ=255)                                                                                                                  cgccaaaccttgccggtCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCgg                                                                                                                                  >  1:41675/17‑251 (MQ=255)                                                                                                                    ccaaaccttgccggtCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCggg                                                                                                                                 <  1:748746/236‑1 (MQ=255)                                                                                                                     caaaccttgccggtCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCggggg                                                                                                                               <  1:271732/238‑1 (MQ=255)                                                                                                                     caaaccttgccggtCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCggggg                                                                                                                               <  1:38965/238‑1 (MQ=255)                                                                                                                     caaaccttgccggtCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCggggg                                                                                                                               <  1:497213/238‑1 (MQ=255)                                                                                                                     caaaccttgccggtCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCggggg                                                                                                                               <  1:766329/238‑1 (MQ=255)                                                                                                                     caaacctggccggtCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCggggg                                                                                                                               <  1:855603/238‑1 (MQ=255)                                                                                                                                  tCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAAc                                                                                                                   <  1:622352/250‑1 (MQ=255)                                                                                                                                     gcttCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCggggg                                                                                                                               >  1:571006/3‑235 (MQ=255)                                                                                                                                      cttCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACg                                                                                                              >  1:8584/2‑251 (MQ=255)                                                                                                                                        tCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACgga                                                                                                            >  1:25636/1‑251 (MQ=255)                                                                                                                                                  aCCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACGGAAAGAAAATc                                                                                                   >  1:28006/1‑250 (MQ=255)                                                                                                                                                                 ttCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACgttgttgt                                                                                   <  1:394151/251‑1 (MQ=255)                                                                                                                                                                   ccggcAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCggggg                                                                                                                               >  1:780018/5‑205 (MQ=255)                                                                                                                                                                     ggcAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTaa                                                                               >  1:422809/3‑251 (MQ=255)                                                                                                                                                                       cAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTaatc                                                                             <  1:807808/251‑3 (MQ=255)                                                                                                                                                                       cAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTaatc                                                                             <  1:140263/251‑3 (MQ=255)                                                                                                                                                                       cAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTaatc                                                                             <  1:559844/251‑3 (MQ=255)                                                                                                                                                                                       cGGGTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGt                                                             <  1:810302/250‑1 (MQ=255)                                                                                                                                                                                         ggtacaatgAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTg                                                           >  1:709743/9‑250 (MQ=255)                                                                                                                                                                                         ggTGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAAc                                                                                                                   <  1:497447/195‑1 (MQ=255)                                                                                                                                                                                           tGCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGccc                                                        <  1:476885/251‑1 (MQ=255)                                                                                                                                                                                            gCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTaa                                                                               <  1:77605/228‑1 (MQ=255)                                                                                                                                                                                            gCAATGAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGccc                                                        >  1:704049/1‑250 (MQ=255)                                                                                                                                                                                             caatgAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGt                                                             >  1:881949/5‑244 (MQ=255)                                                                                                                                                                                                 gAAAGGGGCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAgcc                                                   >  1:447361/1‑248 (MQ=255)                                                                                                                                                                                                     ggggCAAATTGCGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAgccctcg                                               >  1:465959/1‑244 (MQ=255)                                                                                                                                                                                                        gcaaattgcGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCgctg                                            <  1:819576/242‑1 (MQ=255)                                                                                                                                                                                                            attgcGTTTAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGcc                                       <  1:168332/247‑1 (MQ=255)                                                                                                                                                                                                                    tAAATTGTTGTTCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTgggc                               <  1:635730/251‑2 (MQ=255)                                                                                                                                                                                                                            tgttCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCTAGTTGGGCTGCACGGc                       <  1:562036/248‑1 (MQ=255)                                                                                                                                                                                                                            tgttCCGCCCGTTGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGc                       <  1:602337/248‑1 (MQ=255)                                                                                                                                                                                                                                       ttGCACCAACGCCCGCCCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCAc                                                                                                <  1:315291/168‑1 (MQ=255)                                                                                                                                                                                                                                                cgcccgccCGTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCaaag   <  1:289964/251‑4 (MQ=255)                                                                                                                                                                                                                                                  cccgcccgTGCTCCGTTAGCGGACACAACGCCACCACGGAATAACAGCCGCTGCCCGCATCTTTGCTTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCGAGCAGTTCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCaaaga  <  1:369461/250‑5 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       GGGCGTTTGCAATAGGCTGTCAATAAAGGAAAATCCCGCTCCCCGCTTAATTTCATCAAGCGCCGAATAGCGGGCAAGGCAAACGCCAACGGCAGGAAGAAATCGGCAAAATCCGCCAAGCATTGCCAATCGGTTTCATTGAGCGACCAATCCGTTTCGGTTCCGACAAATTGCGTCGCCAGCGGGTGCAACGAAAGGGGCAGATTACGTTTAAATTGTTGCTCCGCCTGTTGCACCAACGCCCGTCCGCGCTCCGTCAGCGGGCACAATGCCACTACGGAATAACAGCCGCTGCCTGCATCTTTGCCTTCGCTGAAATGCACCAAAATAAAACCGCACTTTTGCCAAAACGCCAGCAGCTCGGGGGTGTAACCGAAACTGACGGAAAGGAAATCCACCGCACGTTGTTGTTTAACTC‑GTTCGATTAAACGTTGCCCCAGTCCCCGCTGCTGCCAGTTGGGTTGCACGGCGATACGGGAAATGCGCAGTGA  >  1/1474243‑1474733

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