Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	498,116	C→T	P21L (CCT→CTT)	hemH →	ferrochelatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	498,116	0	C	T	99.0%	378.1 / NA	98	P21L (CCT→CTT)	hemH	ferrochelatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (57/40);  minor base G (0/1);  total (57/41)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACCTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAAT  >  NC_000913/497869‑498353                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              gggTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCTGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCGCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGa                                                                                                                                                                                                                                            <  2:586313/251‑1 (MQ=255)         ggAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTaa                                                                                                                                                                                                                                    <  2:41177/251‑1 (MQ=255)              aGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACACAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGc                                                                                                                                                                                                                               >  1:241395/1‑251 (MQ=255)              aGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGc                                                                                                                                                                                                                               >  2:416910/1‑251 (MQ=255)               ggCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCt                                                                                                                                                                                                                              >  2:444174/1‑251 (MQ=255)                gCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCt                                                                                                                                                                                                                              >  1:118928/1‑250 (MQ=255)                  ccGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATc                                                                                                                                                                                                                           >  2:209496/1‑251 (MQ=255)                          gCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGGTGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACa                                                                                                                                                                                                                    >  1:164996/1‑250 (MQ=255)                            aaaCTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAAtt                                                                                                                                                                                                                 >  1:469229/1‑251 (MQ=255)                                ttATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGGAAAACGCTATCTGAAACAAttttt                                                                                                                                                                                                              >  2:695969/1‑250 (MQ=255)                                 tATCAGCAAGCTCAGGCCTGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTGTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTa                                                                                                                                                                                                             <  2:614787/250‑1 (MQ=255)                                      gCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGAc                                                                                                                                                                                                       <  1:377098/251‑1 (MQ=255)                                             cAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGAcgcg                                                                                                                                                                                                >  1:227193/1‑251 (MQ=255)                                                 ccGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGt                                                                                                                                                                                            >  1:367183/1‑251 (MQ=255)                                                           attcGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTCGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTc                                                                                                                                                                                   >  1:166272/1‑250 (MQ=255)                                                               ggCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAACCCTGGGTACGCCCGATGCCCCCACACGTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGgt                                                                                                                                                                              <  2:644779/251‑1 (MQ=255)                                                                  cTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGgttgt                                                                                                                                                                           >  1:417874/1‑251 (MQ=255)                                                                      tACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTAt                                                                                                                                                                        <  1:672830/250‑1 (MQ=255)                                                                        cAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTAtggg                                                                                                                                                                     >  2:209848/1‑250 (MQ=255)                                                                           aaaCCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGc                                                                                                                                                                  <  1:140101/251‑1 (MQ=255)                                                                             aCCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGccc                                                                                                                                                                >  2:611149/1‑250 (MQ=255)                                                                                       gcgcTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGTCATTGCTGcgcg                                                                                                                                                      <  1:481815/251‑1 (MQ=255)                                                                                           tATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGt                                                                                                                                                  >  1:601036/1‑251 (MQ=255)                                                                                            aTTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGGTGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTg                                                                                                                                                 <  2:164996/251‑1 (MQ=255)                                                                                               ccTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGAtt                                                                                                                                              >  2:163601/1‑251 (MQ=255)                                                                                                  ttttCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGAttt                                                                                                                                             <  2:510559/249‑1 (MQ=255)                                                                                                    ttCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGcc                                                                                                                                         >  1:100781/1‑251 (MQ=255)                                                                                                     tCTGATTTGCCCTCTCACAGCCAGTAGTTCTTCTTCCTCCCTTTTCCGCTACAATTATCACCCAGTTGAATCGATAAGAGGCGGTAATGCGTGAGAATAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCGCCACACTTGAAGCGGTAAAAGGCAATCTGAAACAATTTTTAAACGACAGACGCGGGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGGGATTTTGCCg                                                                                                                                        >  2:250025/1‑251 (MQ=255)                                                                                                        gATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCt                                                                                                                                      <  2:602104/250‑1 (MQ=255)                                                                                                         aTTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTgc                                                                                                                                    >  2:653863/1‑251 (MQ=255)                                                                                                         aTTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTgc                                                                                                                                    >  1:291903/1‑251 (MQ=255)                                                                                                            tGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCt                                                                                                                                 >  1:133843/1‑251 (MQ=255)                                                                                                            tGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCt                                                                                                                                 >  1:158287/1‑251 (MQ=255)                                                                                                             gACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTc                                                                                                                                <  1:444174/251‑1 (MQ=255)                                                                                                                    cacaGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTcgccg                                                                                                                            <  1:209848/248‑1 (MQ=255)                                                                                                                      caGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGGGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCgtgt                                                                                                                       >  1:639581/1‑251 (MQ=255)                                                                                                                       aGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGtgtg                                                                                                                      >  1:285719/1‑251 (MQ=255)                                                                                                                          aaTTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCg                                                                                                                   >  2:702670/1‑251 (MQ=255)                                                                                                                                    tcttcCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGt                                                                                                            <  1:695969/248‑1 (MQ=255)                                                                                                                                    tcttcCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTAt                                                                                                          <  1:473457/250‑1 (MQ=255)                                                                                                                                       tcCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGc                                                                                                        <  2:268972/249‑1 (MQ=255)                                                                                                                                         ntCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCctct                                                                                                    >  2:309751/2‑251 (MQ=255)                                                                                                                                         ntCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAAAGCTATCTGAAAAAATTTTTAAGCGACAGACGCGTGGTTGATAACTAACGGTTGGTATGGTGGCCATTGCTGCGAGGAGTGATTTTGCCGCGGCGCTCGCCGCGTGGGGCGAAGCTGTATGCctct                                                                                                    >  2:309740/2‑251 (MQ=255)                                                                                                                                            aCTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGGGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCtctgtc                                                                                                 >  2:409414/1‑251 (MQ=255)                                                                                                                                            aCTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTAAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGGACGCCCGATGCCCCCAAACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGGGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGGATGCCtctggc                                                                                                 >  1:26669/1‑249 (MQ=255)                                                                                                                                               tttCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCtctgtct                                                                                                <  2:429427/249‑1 (MQ=255)                                                                                                                                                  ccGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCtggatg                                                                                           <  2:469229/251‑1 (MQ=255)                                                                                                                                                   cGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCtggatgg                                                                                          <  2:440249/251‑1 (MQ=255)                                                                                                                                                       aCAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAg                                                                                       <  2:285719/250‑1 (MQ=255)                                                                                                                                                        cAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGt                                                                                     <  2:158287/251‑1 (MQ=255)                                                                                                                                                            tATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGc                                                                                  <  1:102272/250‑1 (MQ=255)                                                                                                                                                             aTCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTc                                                                                >  1:237841/1‑251 (MQ=255)                                                                                                                                                                aacaaGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTcgcc                                                                             <  2:100781/251‑1 (MQ=255)                                                                                                                                                                 acaaGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTcgccg                                                                            >  1:191540/1‑251 (MQ=255)                                                                                                                                                                   aaGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCt                                                                          >  1:27730/1‑251 (MQ=255)                                                                                                                                                                   aaGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGGATTCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAAAGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGGGTGGCGGAGCCGTATGCCTCTGTCTGGCTGGAAGGGGGCTcgccgcg                                                                          >  2:283580/1‑250 (MQ=255)                                                                                                                                                                         aaTCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCcgcg                                                                                                                          >  1:26128/1‑197 (MQ=255)                                                                                                                                                                         aaTCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCcgcg                                                                                                                          <  2:26128/197‑1 (MQ=255)                                                                                                                                                                         aaTCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGt                                                                    >  1:448352/1‑251 (MQ=255)                                                                                                                                                                          aTCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGGGGCTCGCCGCTGATGGtt                                                                   >  1:506165/1‑251 (MQ=255)                                                                                                                                                                            cGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCCGCGCGGCGTGAttt                                                                                                                                             <  2:609090/175‑1 (MQ=255)                                                                                                                                                                            cGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATACCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGAttt                                                                                                                                             >  1:609090/1‑175 (MQ=255)                                                                                                                                                                             gATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTAc                                                                >  1:392602/1‑251 (MQ=255)                                                                                                                                                                              aTAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGCAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTAAGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTAc                                                                >  1:337415/1‑250 (MQ=255)                                                                                                                                                                              aTAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACa                                                               <  1:209496/251‑1 (MQ=255)                                                                                                                                                                                aaGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAgc                                                             <  2:291903/251‑1 (MQ=255)                                                                                                                                                                                   agGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAgccgc                                                          >  1:352672/1‑251 (MQ=255)                                                                                                                                                                                    ggCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAgccgcc                                                         >  1:420086/1‑251 (MQ=255)                                                                                                                                                                                        gTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGgttgt                                                                                                                                                                           <  2:264551/133‑1 (MQ=255)                                                                                                                                                                                        gTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGgttgt                                                                                                                                                                           >  1:264551/1‑133 (MQ=255)                                                                                                                                                                                        gTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGc                                                      >  2:369090/1‑250 (MQ=255)                                                                                                                                                                                          aaTGCGTCAGACTAAAACCGGTATCCTGCTGGCAACCCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAAc                                                   >  1:186297/1‑251 (MQ=255)                                                                                                                                                                                          aaTGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAAc                                                   >  2:272884/1‑251 (MQ=255)                                                                                                                                                                                               gTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGACGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAAcagc                                                >  2:598798/1‑249 (MQ=255)                                                                                                                                                                                                 cAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAgg                                             <  2:352672/250‑1 (MQ=255)                                                                                                                                                                                                 cAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGgc                                            <  2:227193/251‑1 (MQ=255)                                                                                                                                                                                                 cAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGgc                                            <  2:367183/251‑1 (MQ=255)                                                                                                                                                                                                 cAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTAGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGgc                                            <  1:416910/251‑1 (MQ=255)                                                                                                                                                                                                  aGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGgcg                                           <  1:653863/251‑1 (MQ=255)                                                                                                                                                                                                       aaaaCCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGc                                      >  1:325183/1‑251 (MQ=255)                                                                                                                                                                                                         aaCCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGcac                                    >  2:698756/1‑251 (MQ=255)                                                                                                                                                                                                          aCCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGcac                                    >  1:678318/1‑250 (MQ=255)                                                                                                                                                                                                           ccGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGcaca                                   >  1:507063/1‑250 (MQ=255)                                                                                                                                                                                                            cGGTATCCTGCTGGCAACCCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCAGTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAAc                                 <  1:283580/251‑1 (MQ=255)                                                                                                                                                                                                             ggTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACg                                <  2:118928/251‑1 (MQ=255)                                                                                                                                                                                                               tATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTg                                                                                    >  2:192559/1‑197 (MQ=255)                                                                                                                                                                                                               tATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTg                                                                                    <  1:192559/197‑1 (MQ=255)                                                                                                                                                                                                                  ccTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTa                            <  2:186297/250‑1 (MQ=255)                                                                                                                                                                                                                      ctgGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCgg                        >  2:274932/1‑250 (MQ=255)                                                                                                                                                                                                                          cAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATg                   >  1:673550/1‑251 (MQ=255)                                                                                                                                                                                                                            aaCCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAACCAATTTTTAAGCGACAGACGCGCGGTTGATACCTCTCTGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGcc                 <  2:498273/251‑1 (MQ=255)                                                                                                                                                                                                                            aaCCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGcc                 >  2:81092/1‑251 (MQ=255)                                                                                                                                                                                                                                 gggTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGTCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCTTAg            <  1:309740/251‑1 (MQ=255)                                                                                                                                                                                                                                 gggTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAg            <  2:166272/251‑1 (MQ=255)                                                                                                                                                                                                                                 gggTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAg            <  1:309751/251‑1 (MQ=255)                                                                                                                                                                                                                                         ccGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTggg     >  2:280701/1‑250 (MQ=255)                                                                                                                                                                                                                                           gATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAAt  >  1:512014/1‑251 (MQ=255)                                                                                                                                                                                                                                           gATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAAt  >  1:288648/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              GGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACCTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAAT  >  NC_000913/497869‑498353

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