Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,457,122	C→T	L642L (CTG→CTA)	unknown ←	DNA topoisomerase I (EC 5.99.1.2)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,457,122	0	C	T	100.0%	206.0 / NA	72	L642L (CTG→CTA)	unknown	DNA topoisomerase I (EC 5.99.1.2)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/36);  total (36/36)

GTGCATATCACCGCCGCATTTATCGCATTCAACGATTGGGCCGTCATAACCTTTAATCTTAAACTTACCTTCTTCGACTAAATAGCCGTCACAATTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACCAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTGCGA‑CCG‑CAAGTCGGGCAACAAATATCCGTCTC‑AACCAAATTATTCGGGCGCATACCGCCCTCCAATTCATCTAA‑TTCCGCTTTGGTCAATTGGGTCGAGAAGTCTTTGAAGAATTGGTTTAATTCGGTTTTCCAGTTTTTATCCCCGGAA  >  1/1456876‑1457368                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           gTGCATATCACCGCCGCATTTATCGCATTCAACGATTGGGCCGTCATAACCTTTAATCTTAAACTTACCTTCTTCGACTAAATAGCCGTCACAATTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGt                                                                                                                                                                                                                                                        >  1:265876/1‑251 (MQ=255) gTGCATATCACCGCCGCATTTATCGCATTCAACGATTGGGCCGTCATAACCTTTAATCTTAAACTTACCTTCTTCGACTAAATAGCCGTCACAATTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGt                                                                                                                                                                                                                                                        >  1:686688/1‑251 (MQ=255)    cATATCACCGCCGCATTTATCGCATTCAACGATTGGGCCGTCATAACCTTTAATCTTAAACTTACCTTCTTCGACTAAATAGCCGTCACAATTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTaa                                                                                                                                                                                                                                                     <  1:734547/251‑1 (MQ=255)       atCACCGCCGCATTTATCGCATTCAACGATTGGGCCGTCATAACCTTTAATCTTAAACTTACCTTCTTCGACTAAATAGCCGTCACAATTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAAtgg                                                                                                                                                                                                                                                  <  1:187144/251‑1 (MQ=255)                                           gTCATAACCTTTAATCTTAAACTTACCTTCTTCGACTAAATAGCCGTCACAATTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAg                                                                                                                                                                                                               >  1:500518/1‑250 (MQ=255)                                            tCATAACCTTTAATCTTAAACTTACCTTCTTCGACTAAATAGCCGTCACAATTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGcc                                                                                                                                                                                                             <  1:10527/251‑1 (MQ=255)                                            tCATAACCTTTAATCTTAAACTTACCTTCTTCGACTAAATAGCCGTCACAATTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGc                                                                                                                                                                                                              >  1:444506/1‑250 (MQ=255)                                                      tAATCTTAAACTTACCTTCTTCGACTAAATAGCCGTCACAATTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCt                                                                                                                                                                                                   <  1:817496/251‑1 (MQ=255)                                                        aTCTTAAACTTACCTTCTTCGACTAAATAGCCGTCACAATTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGc                                                                                                                                                                                                                                         >  1:175306/1‑211 (MQ=255)                                                        aTCTTAAACTTACCTTCTTCGACTAAATAGCCGTCACAATTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTaa                                                                                                                                                                                                 >  1:156929/1‑251 (MQ=255)                                                        aTCTTAAACTTACCTTCTTCGACTAAATAGCCGTCACAATTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTaa                                                                                                                                                                                                 >  1:566075/1‑251 (MQ=255)                                                           ttAAACTTACCTTCTTCGACTAAATAGCCGTCACAATTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGaa                                                                                                                                                                                              <  1:687561/251‑1 (MQ=255)                                                           ttAAACTTACCTTCTTCGACTAAATAGCCGTCACAATTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGaa                                                                                                                                                                                              <  1:443817/251‑1 (MQ=255)                                                           ttAAACTTACCTTCTTCGACTAAATAGCCGTCACAATTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGaa                                                                                                                                                                                              <  1:898894/251‑1 (MQ=255)                                                           ttAAACTTACCTTCTTCGACTAAATAGCCGTCACAATTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGaa                                                                                                                                                                                              <  1:525820/251‑1 (MQ=255)                                                            tAAACTTACCTTCTTCGACTAAATAGCCGTCACAATTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAAt                                                                                                                                                                                             <  1:40468/251‑1 (MQ=255)                                                               aCTTACCTTCTTCGACTAAATAGCCGTCACAATTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGGGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACg                                                                                                                                                                                          >  1:594888/1‑251 (MQ=255)                                                                       tcttcGACTAAATAGCCGTCACAATTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTAc                                                                                                                                                                                   <  1:101521/250‑1 (MQ=255)                                                                             aCTAAATAGCCGTCACAATTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGcc                                                                                                                                                                                                             >  1:347011/1‑218 (MQ=255)                                                                                  aTAGCCGTCACAATTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCttt                                                                                                                                                                       >  1:408298/1‑249 (MQ=255)                                                                                             aaTTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTg                                                                                                                                                            <  1:834083/251‑1 (MQ=255)                                                                                                   ggTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGGTCGCCATATTGCGA‑cc                                                                                                                                                      >  1:513536/1‑250 (MQ=255)                                                                                                              cAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGcc                                                                                                                                                                                        >  1:682417/1‑206 (MQ=255)                                                                                                                aaTATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTTTGATCGCCATATTGCGA‑CCGCCAAG‑CGGGcaa                                                                                                                                        >  1:487831/1‑250 (MQ=255)                                                                                                                   atGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTGCGA‑CCG‑CAAGTCGGGCAACaaa                                                                                                                                    <  1:55008/251‑1 (MQ=255)                                                                                                                        gTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCg                                                                                                                                                                              <  1:16030/206‑1 (MQ=255)                                                                                                                                     aTCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTGCGA‑CCG‑CAAGTCGGGCAACAAATATCCGTCTC‑AACCAAat                                                                                                                 <  1:540066/251‑1 (MQ=255)                                                                                                                                         ataataTAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTGCGA‑CCG‑CAAGTCGGGCAACAAATATc                                                                                                                                <  1:315625/233‑1 (MQ=255)                                                                                                                                         ataataTAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGGTCTTATCGCCACATTGCGA‑CCG‑CAAATCGGGCAACAAAAATCCCTCTC‑CACCCAattatt                                                                                                             >  1:143714/1‑251 (MQ=255)                                                                                                                                          taataTAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTTTGTTCGCCACATTGTGA‑CCG‑CCAGCGGGGcaaccaactgtttcttattcacatctccg                                                                                                              >  1:825376/1‑225 (MQ=255)                                                                                                                                          taataTAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCCTATTTCGA‑CCG‑CAAGACGGGCAACAAAAATCCGCCTC‑AACCCAaatattt                                                                                                            >  1:545226/1‑250 (MQ=255)                                                                                                                                          taataTAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTGCGA‑CCG‑CAAGTCGGGCAACAAATATcc                                                                                                                               <  1:22290/233‑1 (MQ=255)                                                                                                                                          taataTAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTGCGA‑CCG‑CAAGTCGGGCAACAAATATCCGTCTC‑AACCAAATTATTc                                                                                                            <  1:8171/251‑1 (MQ=255)                                                                                                                                          taataTAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTGCGA‑CCG‑CAAGTCGGGCAAAAAAAATCCGCCTC‑CACCAAaatattt                                                                                                            >  1:570621/1‑250 (MQ=255)                                                                                                                                          taataTAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGG‑TTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGcc                                                                                                                                                                                        <  1:854578/177‑1 (MQ=255)                                                                                                                                          taataTAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTACATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGcc                                                                                                                                                                                        >  1:738274/1‑178 (MQ=255)                                                                                                                                          atatatAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGGTCTGATCGCCCTCTTTCGG‑CCG‑CCAGTCGGGCAACCAatattcctctccacccaaatattt                                                                                                             >  1:368666/3‑231 (MQ=255)                                                                                                                                            atatAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGaa                                                                                                                                                                                              <  1:708896/170‑1 (MQ=255)                                                                                                                                            atatAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGcc                                                                                                                                                                                        <  1:786904/176‑1 (MQ=255)                                                                                                                                            atatAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCg                                                                                                                                                                              >  1:761678/1‑186 (MQ=255)                                                                                                                                             tataACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTGCGA‑CCG‑CAAGTCGGGCAACAAATATCCGTCTC‑AACCAAATTATTCggg                                                                                                         <  1:362027/251‑1 (MQ=255)                                                                                                                                               taACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCAt                                                                                                                                                                                                                 <  1:484787/148‑1 (MQ=255)                                                                                                                                               taACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTGCGA‑CCG‑CAAGTCGGGCAACAAATATCCGTCTC‑AACCAAATTATTCGGgcg                                                                                                       <  1:390924/251‑1 (MQ=255)                                                                                                                                                    gTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTTCGA‑CCG‑CAAGTCGGGCAACAAAAATCCGTCTC‑AACCAAAATATTCTGGCGCATa                                                                                                   >  1:58458/1‑250 (MQ=255)                                                                                                                                                     tCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTGCGA‑CCG‑CAAGTCGGGCAACAAATATCCGTCTC‑AACCAAATTATTCGGGCGCATAc                                                                                                  <  1:846750/250‑1 (MQ=255)                                                                                                                                                     tCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTGCGA‑CCG‑CAAGTCGGGCAACAAATATCCGTCTC‑AACCAAATTATTCGGGCGCATAc                                                                                                  <  1:97432/250‑1 (MQ=255)                                                                                                                                                       cATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTGCGA‑CCG‑CAAGTCGGGCAACAAATATCCGTCTC‑AACCAAATTATTCGGGCGCATACCGc                                                                                               <  1:117221/251‑1 (MQ=255)                                                                                                                                                         tCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTTATCGCCCCATTTCGG‑CCG‑CCAGGCGGGCCACAAAAATCCCTCTC‑CACCAAaatatttcggggcctaccccc                                                                                              >  1:811389/1‑235 (MQ=255)                                                                                                                                                         tCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTGCGA‑CCG‑CAAGTCGGGCAACAAATATCCGTCTC‑AACCAAat                                                                                                                 <  1:753614/231‑1 (MQ=255)                                                                                                                                                              gTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCt                                                                                                                                                                                                   >  1:33643/1‑147 (MQ=255)                                                                                                                                                              gTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGTTCGCCATATTTCGA‑CCG‑CAAGGCGGGCAACAAAAAACCGTCTC‑AACCAAaatattttggggcacacctccgccttt                                                                                        >  1:600985/1‑230 (MQ=255)                                                                                                                                                                    cATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTTTGATCGCCATATTGCGG‑CCG‑CAAGTCGGGGAACAAAAATTCGTCTC‑AACCAAATTATTCGGGGGCATACCCCCCTCCAATTCAt                                                                                   >  1:352930/1‑250 (MQ=255)                                                                                                                                                                    cATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTGCGA‑CCG‑CAAGTCGGGCAACAAATATCCGTCTC‑AACCAAATTATTCGGGCGCATACCGCCCTCCAATTCATc                                                                                  <  1:98575/251‑1 (MQ=255)                                                                                                                                                                     aTTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCCTATTTCGA‑CCG‑CAAATCGGGCAACAAAAATCCGTCTC‑CACCAAATTATTTTGGCGCATACCGCCCTCCCATTCAt                                                                                   >  1:686981/1‑249 (MQ=255)                                                                                                                                                                         cGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCCTAATGCCA‑CCC‑CAAGTCGGGCAACAAATATCCGTCTC‑CACCAAATTATTCCGGGGCATACCGCCCCCCCAATCATCTAA‑tt                                                                            >  1:175252/1‑251 (MQ=255)                                                                                                                                                                           ggCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTGCGA‑CCG‑CAAGTCGGGCAACAAATATCCGTCTC‑AACCAAATTATTCGGGCGCATACCGCCCTCCAATTCATCTAA‑TTc                                                                           <  1:560731/250‑1 (MQ=255)                                                                                                                                                                               gCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGt                                                                                                                                                                                                                                                        <  1:855015/77‑1 (MQ=255)                                                                                                                                                                               gCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTAc                                                                                                                                                                                   <  1:11017/146‑1 (MQ=255)                                                                                                                                                                                        ttCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTGCGA‑CCG‑CAAGTCGGGCAACAAATATCCGTCTC‑AACCAAATTATTCGGGCGCATACCGCCCTCCAATTCATCTAA‑TTCCGCTTTGGTCAAt                                                              <  1:766951/250‑1 (MQ=255)                                                                                                                                                                                         tCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTGCGA‑CCG‑CAAGTCGGGCAACAAATATCCGTCTC‑AACCAAATTATTCGGGCGCATACCGCCCTCCAATTCATCTAA‑TTCCGCTTTGGTCAATTg                                                            <  1:233940/251‑1 (MQ=255)                                                                                                                                                                                          ccATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTTCGA‑CCG‑CAAGTCGGGGAACAAAAATTCGTCTC‑AACCAAAATATTCGGGCGCATACCCCCCCCCCAATCATCTAATTTCCCCTTTGG‑CAATTgg                                                           >  1:245562/1‑251 (MQ=255)                                                                                                                                                                                                         tttCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTGCGA‑CCG‑CAAGTCGGGCAACAAATATCCGTCTC‑AACCAAATTATTCGGGCGCATACCGCCCTCCAATTCATCTAA‑TTCCGCTTTGGTCAATTGGGTCGAGAAGTCttt                                             <  1:760050/250‑1 (MQ=255)                                                                                                                                                                                                                      gtcCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTTTGATCGCCCTATTGCGA‑CCG‑CAAAACGGGCAAAAAAAATCCGTCTC‑AACCAAAATATTTGGGCGCCTACCCCCCCCCCATTCATCTAA‑TTCCCCCTTTGTCAATTGGGTCGGGAAGGCTTTTAAGAATTgggttt                               >  1:270638/1‑247 (MQ=255)                                                                                                                                                                                                                      gtcCAACACATTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTTATCGccccattt                                                                                                                                                            >  1:85303/1‑124 (MQ=255)                                                                                                                                                                                                                            cacaTTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGaa                                                                                                                                                                                              >  1:443854/1‑90 (MQ=255)                                                                                                                                                                                                                             acaTTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTTCGA‑CCC‑CAAGGCGGGGAACAAAAATCCGTCTC‑AACCAAAATATTCGGGCGCATACCCCCCTCCCATTCATCTAA‑TTTCGCTTTGGTCCATTTGGGCGAGAAGGCTTTGTAGAATTGGTTTTATTCGGt                        >  1:120171/1‑251 (MQ=255)                                                                                                                                                                                                                              caTTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTTATCGCCATATTTCGA‑CCG‑CAAGTCGGGCAACAAATATCCCCTTC‑AACCAAATTATTTGGGGGCATACCGCCCTCCAAATCATCTTA‑TTTCGCTTTGGTCAATTGGGTCGAGAAGTCTTTTTAGAATTGGGTTAATTCTGtt                       >  1:642079/1‑251 (MQ=255)                                                                                                                                                                                                                              caTTTAATAGTTCCGCTTCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTGCGA‑CCG‑CAAGTCGGGCAACAAATATCCGTCTC‑AACCAAATTATTCGGGCGCATAc                                                                                                  <  1:864938/177‑1 (MQ=255)                                                                                                                                                                                                                                              ttCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTGCGA‑CCG‑CAAGTCGGGCAACAAATATCCGTCTC‑AACCAAATTATTCGGGCGCATAc                                                                                                  <  1:97204/161‑1 (MQ=255)                                                                                                                                                                                                                                               tCAGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCt                                                                                                                                                                                                   >  1:856392/1‑66 (MQ=255)                                                                                                                                                                                                                                                 aGGAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTTATCGCCCTATTGCGACCCC‑CAAG‑CGGGGAACAAAAATCCCGCTCAAACCAAA‑TATTCGGGGGCCTACCCCCCCCCCAATCATCTAA‑ttt                                                                           >  1:414160/1‑179 (MQ=255)                                                                                                                                                                                                                                                   gAACTAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCCTATTTCGG‑CCC‑CCAATCGGGGCACAAAAATCCCTCTC‑CACCAAAATATTTGGGCGCCTACCCCCCCCCCATTCCTCTTA‑TTCCCCTTTTGTCAATTTGGTCGGGAAGTCTTTTTAGAATTTGTTTTATTTGGTTTTTCCGTTTTTTTTCCCggga  >  1:14590/1‑249 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           GTGCATATCACCGCCGCATTTATCGCATTCAACGATTGGGCCGTCATAACCTTTAATCTTAAACTTACCTTCTTCGACTAAATAGCCGTCACAATTCGGGTTATTGCCGCAAATATGTAGTTTTCGATGAGGATCGATAATATAACTGTCCATCGCCGTACCGCATTTCGGGCAGCGTTTACGTTCCATTAGCGCTTTGGTTTCCGACGCGTCGTCCAACACATTTAATAGTTCCGCTTCAGGAACCAGGTTAATGGTGGTTTTGCAACGTGCTTTCGGTGGTAAGGCATAGCCGGTGCAGCCTAAGAATACGCCTGTACTGGCGGTTCTGATCGCCATATTGCGA‑CCG‑CAAGTCGGGCAACAAATATCCGTCTC‑AACCAAATTATTCGGGCGCATACCGCCCTCCAATTCATCTAA‑TTCCGCTTTGGTCAATTGGGTCGAGAAGTCTTTGAAGAATTGGTTTAATTCGGTTTTCCAGTTTTTATCCCCGGAA  >  1/1456876‑1457368

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