Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,570,706	G→A	Q221Q (CAG→CAA)	unknown →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,570,706	0	G	A	100.0%	190.9 / NA	70	Q221Q (CAG→CAA)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/37);  total (33/37)

ACCGTGTATCGCTCCGGCTGGTCGGGCTGGTATGTGCGGGTGGTACGCGAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAGAGTTAAAAATGACCGAACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTAAAAAATGTGCCTTATATTGTGGAAACCACCGTGCCGAACCAACAAAAGCTGCTGAACGTGCAGGCAATTCAAACCTGGGCGGAAAA  >  1/1570461‑1570941                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            aCCGTGTATCGCTCCGGCTGGTCGGGCTGGTATGTGCGGGTGGTACGCGAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTa                                                                                                                                                                                                                                        >  1:196991/1‑251 (MQ=255)               cGGCTGGTCGGGCTGGTATGTGCGGGTGGTACGCGAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCgca                                                                                                                                                                                                                            >  1:626926/1‑247 (MQ=255)                    ggTCGGGCTGGTATGTGCGGGTGGTACGCGAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAg                                                                                                                                                                                                                     >  1:882633/1‑251 (MQ=255)                     gTCGGGCTGGTATGTGCGGGTGGTACGCGAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAg                                                                                                                                                                                                                     <  1:300614/250‑1 (MQ=255)                      tCGGGCTGGTATGTGCGGGTGGTACGCGAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGAt                                                                                                                                                                                                                   <  1:347825/251‑1 (MQ=255)                        gggCTGGTATGTGCGGGTGGTACGCGAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGAtt                                                                                                                                                                                                                  >  1:696803/1‑250 (MQ=255)                        gggCTGGTATGTGCGGGTGGTACGCGAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCCCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGAtt                                                                                                                                                                                                                  >  1:857513/1‑250 (MQ=255)                             ggTATGTGCGGGTGGTACGCGAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTAtt                                                                                                                                                                                                            >  1:709855/1‑251 (MQ=255)                             ggTATGTGCGGGTGGTACGCGAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTAtt                                                                                                                                                                                                            >  1:305589/1‑251 (MQ=255)                              gTATGTGCGGGTGGTACGCGAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGAt                                                                                                                                                                                                                   >  1:606998/1‑243 (MQ=255)                                aTGTGCGGGTGGTACGCGAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGcc                                                                                                                                                                                                         <  1:177175/251‑1 (MQ=255)                                aTGTGCGGGTGGTACGCGAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGcc                                                                                                                                                                                                         <  1:664826/251‑1 (MQ=255)                                 tgtgCGGGTGGTACGCGAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGccc                                                                                                                                                                                                        <  1:131521/251‑1 (MQ=255)                                  gtgCGGGTGGTACGCGAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATg                                                                                                                                                                                                                                  >  1:660879/1‑224 (MQ=255)                                       ggtggtACGCGAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTc                                                                                                                                                                                                   <  1:369668/250‑1 (MQ=255)                                           gtACGCGAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTc                                                                                                                                                                                               >  1:281581/1‑250 (MQ=255)                                                cgAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTCGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATaa                                                                                                                                                                                         <  1:769297/251‑1 (MQ=255)                                                    ggggAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAAt                                                                                                                                                                                                                                   <  1:275293/205‑1 (MQ=255)                                                    ggggAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGAt                                                                                                                                                                                                                   >  1:144642/1‑221 (MQ=255)                                                    ggggAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAacac                                                                                                                                                                                      <  1:81103/250‑1 (MQ=255)                                                     gggAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAAt                                                                                                                                                                                                                                   <  1:552670/204‑1 (MQ=255)                                                     gggAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACAcc                                                                                                                                                                                     <  1:599938/250‑1 (MQ=255)                                                         aaTTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACAcc                                                                                                                                                                                     >  1:805107/1‑246 (MQ=255)                                                         aaTTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATgg                                                                                                                                                                                <  1:566092/251‑1 (MQ=255)                                                           ttCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATgggg                                                                                                                                                                              <  1:408838/251‑1 (MQ=255)                                                           ttCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATgggg                                                                                                                                                                              <  1:662888/251‑1 (MQ=255)                                                                       aCAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCggg                                                                                                                                                                   <  1:209553/250‑1 (MQ=255)                                                                       aCAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGc                                                                                                                                                                  <  1:324971/251‑1 (MQ=255)                                                                        cAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCg                                                                                                                                                                 <  1:101860/251‑1 (MQ=255)                                                                           attaattTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTc                                                                                                                                                                                                   <  1:429652/214‑1 (MQ=255)                                                                                            gTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTaaaaa                                                                                                                                                                                                                                    >  1:382595/1‑164 (MQ=255)                                                                                             tCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGtt                                                                                                                                            >  1:806860/1‑251 (MQ=255)                                                                                                       gACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGATACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAAtgctg                                                                                                                                                                       <  1:454274/214‑1 (MQ=255)                                                                                                              ccATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTCTTCATGCTTACGTGCAGTTTCCCGAAACTACCGag                                                                                                                            >  1:715247/1‑250 (MQ=255)                                                                                                                  tCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAGACTACCGagcag                                                                                                                         >  1:840775/1‑249 (MQ=255)                                                                                                                     ccATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAAtgc                                                                                                                                                                         <  1:348979/198‑1 (MQ=255)                                                                                                                     ccATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGATGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAac                                                                                                                                                                                        >  1:123930/1‑183 (MQ=255)                                                                                                                        tttGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATgcg                                                                                                                 >  1:74985/1‑251 (MQ=255)                                                                                                                        tttGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATgcg                                                                                                                 <  1:419626/251‑1 (MQ=255)                                                                                                                                    gcttTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGc                                                                                                                                                        <  1:897197/200‑1 (MQ=255)                                                                                                                                     cttTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGc                                                                                                     <  1:442156/250‑1 (MQ=255)                                                                                                                                            cgcCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTAt                                                                                              <  1:185229/250‑1 (MQ=255)                                                                                                                                                agACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGtt                                                                                         >  1:395828/2‑251 (MQ=255)                                                                                                                                                     cGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATg                                                                                                                                                                                 <  1:12766/158‑1 (MQ=255)                                                                                                                                                     cGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATg                                                                                                                   <  1:240848/220‑1 (MQ=255)                                                                                                                                                     cGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTaaaaa                                                                                    >  1:181331/1‑251 (MQ=255)                                                                                                                                                       aTTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTaaaaaa                                                                                   <  1:812082/250‑1 (MQ=255)                                                                                                                                                               aTTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTAAAAAATGTGCCTta                                                                          >  1:857120/1‑251 (MQ=255)                                                                                                                                                                   gAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTAAAAAATGTGCCTTatat                                                                       <  1:630791/250‑1 (MQ=255)                                                                                                                                                                        ggCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGc                                                                                                     >  1:583755/1‑215 (MQ=255)                                                                                                                                                                        ggCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTAAAAAATGTGCCTTATATTGTGGa                                                                 >  1:644006/1‑251 (MQ=255)                                                                                                                                                                        ggCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTAAAAAATGTGCCTTATATTGTGGa                                                                 >  1:27358/1‑251 (MQ=255)                                                                                                                                                                            ttGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTAAAAAATGTGCCTTATATTGTGGAAacc                                                             <  1:641506/251‑1 (MQ=255)                                                                                                                                                                                 aCTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACAcc                                                                                                                                                                                     <  1:319695/126‑1 (MQ=255)                                                                                                                                                                                  cTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATgcg                                                                                                                 >  1:215101/1‑193 (MQ=255)                                                                                                                                                                                    tGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTaaaaa                                                                                                                                                                                                                                    >  1:166579/1‑76 (MQ=255)                                                                                                                                                                                        gtgtTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAacac                                                                                                                                                                                      <  1:8079/118‑1 (MQ=255)                                                                                                                                                                                          gtTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTAAAAAATGTGc                                                                              <  1:707455/220‑1 (MQ=255)                                                                                                                                                                                          gtTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTAAAAAATGTGCCTTATATTGTGGAAACCACCGTGCCGAACca                                               >  1:733712/1‑251 (MQ=255)                                                                                                                                                                                                   ggCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGccc                                                                                                                                                                                                        >  1:450971/1‑89 (MQ=255)                                                                                                                                                                                                      gTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTAAAAAATGTGCCTTATATTGTgg                                                                  >  1:708138/1‑220 (MQ=255)                                                                                                                                                                                                             cGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTAAAAAATGTGCCTTATATTGTGGAAACCACCGTGCCGAACCAACAAAAGCTGCTGAACGTg                            >  1:613639/1‑251 (MQ=255)                                                                                                                                                                                                                 cattGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTAAAAAATGTGCCTTATATTGTGGAAACCACCGTGCCGAACCAACAAAAGCTGCTGAACGTGCAgg                        <  1:795027/249‑1 (MQ=255)                                                                                                                                                                                                                   ttGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTAAAAAATGTGCCTTATATTGTgg                                                                  <  1:497018/207‑1 (MQ=255)                                                                                                                                                                                                                   ttGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTAAAAAATGTGCCTTATATTGTgg                                                                  <  1:335949/207‑1 (MQ=255)                                                                                                                                                                                                                   ttGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTAAAAAATGTGCCTTATATTGTGGAAACCACCGTGCCGAACCAACAAAAGCTGCTGAACGTGCAGGc                       <  1:809805/250‑1 (MQ=255)                                                                                                                                                                                                                          ctcAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTAAAAAATGTGCCTTATATTGTGGAAACCACCGTGCCGAACCAACAAAAGCTGCTGAACGTGCAGGCAATTCaaa               >  1:786378/1‑251 (MQ=255)                                                                                                                                                                                                                              gCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTaaaa                                                                                     <  1:204799/177‑1 (MQ=255)                                                                                                                                                                                                                                 tACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCt                                                                                                               >  1:732138/1‑148 (MQ=255)                                                                                                                                                                                                                                       gTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTAAAAAATGTGCCTTATATTGTGGAAACCACCGTGCCGAACCAACAAAAGCTGCTGAACGTGCAGGCAATTCAAACCTGGGCGGaaaa  >  1:571476/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            ACCGTGTATCGCTCCGGCTGGTCGGGCTGGTATGTGCGGGTGGTACGCGAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAGAGTTAAAAATGACCGAACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTAAAAAATGTGCCTTATATTGTGGAAACCACCGTGCCGAACCAACAAAAGCTGCTGAACGTGCAGGCAATTCAAACCTGGGCGGAAAA  >  1/1570461‑1570941

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