Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,570,722	A→C	E3A (GAA→GCA)	unknown →	UPF0225 protein YchJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,570,722	0	A	C	100.0%	185.1 / NA	66	E3A (GAA→GCA)	unknown	UPF0225 protein YchJ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (30/36);  total (30/36)

GGTCGGGCTGGTATGTGCGGGTGGTACGCGAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAGAGTTAAAAATGACCGAACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTAAAAAATGTGCCTTATATTGTGGAAACCACCGTGCCGAACCAACAAAAGCTGCTGAACGTGCAGGCAATTCAAACCTGGGCGGAAAATACGCAATGGCTTGG  >  1/1570480‑1570956                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           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)                    ggtggtACGCGAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTc                                                                                                                                                                                                                  <  1:369668/250‑1 (MQ=255)                        gtACGCGAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTc                                                                                                                                                                                                              >  1:281581/1‑250 (MQ=255)                             cgAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTCGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATaa                                                                                                                                                                                                        <  1:769297/251‑1 (MQ=255)                                 ggggAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGAt                                                                                                                                                                                                                                  >  1:144642/1‑221 (MQ=255)                                 ggggAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAacac                                                                                                                                                                                                     <  1:81103/250‑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:662888/251‑1 (MQ=255)                                        ttCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATgggg                                                                                                                                                                                             <  1:408838/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)                                                                          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:419626/251‑1 (MQ=255)                                                                                                     tttGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATgcg                                                                                                                                >  1:74985/1‑251 (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)                                                                                                                                                                     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:335949/207‑1 (MQ=255)                                                                                                                                                                                                ttGCATTCTCAGCTTACCAAGTTACAACAAAAACAAAGTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTAAAAAATGTGCCTTATATTGTgg                                                                                 <  1:497018/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)                                                                                                                                                                                                                                     gTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTAtg                                                                                                            <  1:813274/143‑1 (MQ=255)                                                                                                                                                                                                                                     gTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTAAAAAATGTGCCTTATATTGTGGAAACCACCGTGCCGAACCAACAAAAGCTGCTGAACGTGCAGGCAATTCAAACCTGGGCGGAAAATACGCAAtggcttgg  >  1:671381/1‑249 (MQ=255)                                                                                                                                                                                                                                     gTTAAAAATGACCGCACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTACCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGt                                                                                                               >  1:575244/1‑140 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           GGTCGGGCTGGTATGTGCGGGTGGTACGCGAGGGGGAAATTCAGCGGGGCGACAAATTAATTTTGCAGGCGCGTCCGCATCCCGACTTTACCATTCGCCATTTGAACCGCTTGCTTTCGGCGCCCGACACGATTAACGAATTGGAACAGGCGTTGGCACTTGAGGTGTTGGCGCCGGCGTTTAAACGCTCGTTGCATTCTCAGCTTACCAAGTTACAACAAAAACAGAGTTAAAAATGACCGAACTTTTAGATTCCTATTGCCCTTGTCACTCCCATAACACCTATGGGGAATGCTGCGGGCGTTTTCATGCTTACGTGCAGTTTCCCGAAACTACCGAGCAGCTCATGCGCTCGCGCTATGCCGCGTATGTGTTAAAAAATGTGCCTTATATTGTGGAAACCACCGTGCCGAACCAACAAAAGCTGCTGAACGTGCAGGCAATTCAAACCTGGGCGGAAAATACGCAATGGCTTGG  >  1/1570480‑1570956

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