Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,548,865	C→T	14.8%	A19A (GCC→GCT)	fimH →	minor component of type 1 fimbriae

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,548,865	0	C	T	14.8%	200.8 / 4.9	89	A19A (GCC→GCT)	fimH	minor component of type 1 fimbriae
Reads supporting (aligned to +/- strand):  ref base C (44/27);  new base T (2/9);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 8.79e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATG‑GCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATA  >  NC_000913/4548629‑4549103                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     gTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGAAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAGTGCTTg                                                                                                                                                                                                                                                    >  1:539936/1‑239 (MQ=255)    gcgacacATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGAACTTCCCGTTACAGGTCAGAGCATTGTCAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGCGAAGATTGTAATGAAAC‑GAGGTATTACCCTT‑CTTGCTG‑TACTGCTGAT‑GGGCTGGTAGGCAAATAGCTGGta                                                                                                                                                                                                                                                 >  2:460082/4‑238 (MQ=255)    gCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTc                                                                                                                                                                                                                                                 >  2:970014/1‑239 (MQ=255)        cacaTTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTc                                                                                                                                                                                                                                            <  2:244178/239‑1 (MQ=255)        cacaTTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTc                                                                                                                                                                                                                                            <  2:546631/239‑1 (MQ=255)            ttGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCg                                                                                                                                                                                                                                           >  1:36689/1‑236 (MQ=255)                    tGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCCTTCAGGCAGCGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAAcc                                                                                                                                                                                                                                >  1:859484/1‑239 (MQ=255)                    tGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAAcc                                                                                                                                                                                                                                >  2:80539/1‑239 (MQ=255)                      gcgcAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTTAT‑GGGCTGGTCGGGAAATTCTTGGTC‑ATTCGCCTGTAAAAcgg                                                                                                                                                                                                                               >  2:1079408/1‑236 (MQ=255)                         cAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCaa                                                                                                                                                                                                                           >  1:933897/1‑239 (MQ=255)                         cAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCaa                                                                                                                                                                                                                           >  1:149344/1‑239 (MQ=255)                         cAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCaa                                                                                                                                                                                                                           >  1:372869/1‑239 (MQ=255)                                    gTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTa                                                                                                                                                                                                                >  2:251962/1‑239 (MQ=255)                                     ttCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTa                                                                                                                                                                                                                >  1:91216/1‑238 (MQ=255)                                      tCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTAt                                                                                                                                                                                                               >  1:1002446/1‑238 (MQ=255)                                       cAGGTGGATGATTCCTCAAAATCAGCGCACTTCCTGTTACAGGTCAGAGCATTGACATTAAATGGGGGAACCACTCAGGGAACCATACAGGCATTAATTAGCATCCCATATACCTACAGCTGACCCCGAAGAAATGATTGTAAGAAAAC‑GAGTTATTCCCCTT‑TTTGCTG‑TACTGCTGAT‑GGGCTGTTGGGTAAATGCCTGGTC‑ATTCTCCTGTAAAACCGCCAATGGTACCGCTATcc                                                                                                                                                                                                             <  1:798065/239‑1 (MQ=255)                                       aaaagggATGATTCCCAAAAATCAGCGGCATCCCCGTAACAGGTCAGGGCATGGACAATAAAGGGCGGGGCCCCTCAGGGACCCAGTCGGGCAGAGATTAGCAGCACCTATACCTACAGCTGAACCCGAAGAGAGGATGGTAAGAAAAC‑GATTTATTCCCCTG‑TTTGCTG‑TACGGCTGAG‑GGGCGGGTCGGTAAATGCCTGGCC‑AGTCGCCTGAAAAACCGCCAATGGTACCGCTATc                                                                                                                                                                                                              <  2:741439/233‑1 (MQ=255)                                          gTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTa                                                                                                                                                                                                          >  1:732501/1‑239 (MQ=255)                                               tgatTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGATGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACAGCTATCCCTATTGGc                                                                                                                                                                                                     >  2:992817/1‑239 (MQ=255)                                                    ccTCACAATCGGCGCACTTCACTTTACAGGTCAGAGGATGGCCAGTAAATGGCGGACCCACTCAGGGACCCATTCAGGCAGTGTTTAACATCACTTATCCCTACAGCTGATCCCGAAGAGAGGTTTGTAAGGAACC‑GATTTATTCCCCTT‑TTTGCTG‑TACTGTGGAG‑GGGTGGGTCGGAAAATCCTGGGTC‑TTCCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTg‑g                                                                                                                                                                                               <  2:1094006/239‑1 (MQ=255)                                                     acaaaacttagcgttattgccGTTACAGGTCAGTGTATTGACAGAAAATGGCGTAGCCACTGAGGGAACCATTAAGGCATTTATAAAAATCACTTATTCCAACTGTTGACCCTGAATAAATGATTGTAAGTAAAC‑AAGTTCTTCCCTTG‑TTTGCTG‑TCCTGCTGAT‑GGGCGGTTCGTAAAACGCGTGGTT‑TTTCGCCTGTAAATCCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑Gc                                                                                                                                                                                              <  2:399866/220‑1 (MQ=255)                                                       cacaATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCg                                                                                                                                                                                                    >  1:741439/1‑232 (MQ=255)                                                          aaccggcgcCCTCCCCGTTCAAGGTGAGAGCATTGAAAGTAAAGGGCGGCCCCACCAAGGGAACAATTCAGGAAGTATTTAGAATCCCTTATACCTCCAGTTGACCCCAAAGAGTTGTTTGTAATGAAAC‑GAGTTTTTCCCCTT‑TTTGCTG‑TACTGCTGAT‑GGGCTGTTCGGTAAATCCCTGGTC‑ATTCGCCGGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAgcgc                                                                                                                                                                                         <  2:281926/234‑1 (MQ=255)                                                            tCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCt                                                                                                                                                                                                                                                     >  2:483168/1‑179 (MQ=255)                                                            tCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCt                                                                                                                                                                                                                                                     <  1:483168/179‑1 (MQ=255)                                                                 gacctccccGTAAAAGGTAAGAGCATTAACAGAAAATGGCGGAGCCACTCAGGGAACCATGCGGGCAGTGTTAAGCACCCCTTATACCTACGGTGGACCCGGAGGGGATGATGTTAAGGAAAC‑GATTTATTCCCCTT‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGTTC‑ATCCGCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGt                                                                                                                                                                                   <  2:541949/232‑1 (MQ=255)                                                                      tCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCgtgggcggacgcaatgttaattt                                                                                                                                                                              >  2:1181117/1‑217 (MQ=255)                                                                       cccGTTACAGTTCAGAGCATTCACATAAAATGGCGGACCCACTCAGGGACCCATTCAGGAAGTGATTAGCATCCCCTATACTTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTCCCCTT‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGCC‑ATTCGCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTa                                                                                                                                                                            <  1:865375/239‑1 (MQ=255)                                                                       cccGTTACAGGTCAGAGCATGGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCGGGCAGTGATTAGCATCACCTATACCAACAGCTGAACCCGAAGAGATGATTGTAATGAACC‑GGGTTATTCCCCTT‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTa                                                                                                                                                                            <  1:26585/239‑1 (MQ=255)                                                                        ccGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GACGCACAAACGTTTATGTaa                                                                                                                                                                           >  2:484003/1‑239 (MQ=255)                                                                          gTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCCCTCAGGGAACCATTAAGGCAGTGATTAGCATCACCTATACCTACAGCTGACCCCGAAGAGATGTTTGTATTGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TGCTGCTGAT‑GGGCTGGTCGGTAAAGGCCTGGTC‑ATTCGCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGt                                                                                                                                                                                   <  1:409288/229‑1 (MQ=255)                                                                               aGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCGGCGCCAATGTTTTTTTAAAccctt                                                                                                                                                                     >  1:448692/1‑235 (MQ=255)                                                                               aGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGCG‑TCCCCGCCACAGTTTTTGttataactc                                                                                                                                                                     >  2:1156158/1‑230 (MQ=255)                                                                                 gTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGTGACGCCCAAATGTTTATTTAAAATTTgc                                                                                                                                                                    >  2:865375/1‑238 (MQ=255)                                                                                 gTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGAGCCACTGCTTTTGTAAACCTTgcg                                                                                                                                                                   >  1:621917/1‑238 (MQ=255)                                                                                  tCAGAGAATGAACAGAAAAGGCGGGAGCCATACAGGGACCCATCCAGGCAGTGATGACCGTCCCCAATACTTACAGTTGAACCAAAAGAGATGATTGTAA‑GAAACAGATTTATTCCCCTG‑TTTGTTG‑TACTGCTGAT‑GGGCGGGTCGGTAAATGCTTGGTC‑ATTGGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGcc                                                                                                                                                                 <  2:175206/239‑1 (MQ=255)                                                                                           tGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCCGCGACCATGGTTTTGTAAACCTTGCGGCCGTCGTGa‑a                                                                                                                                                       >  1:504838/1‑239 (MQ=255)                                                                                            gaaagaaaaTGGCGGAGCCACGAAGGGAACCATTAAGGGAGTGATAAGCATCACCTATACCTAAAGCTGAACCGGAAGAGATGATTGTAATGAACC‑GAGTTATTCCCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCGGGTCGGTAAATGCTTGGTC‑AGTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGa                                                                                                                                                         <  2:69328/231‑1 (MQ=255)                                                                                                gTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTT‑GTAGCACCAAAGATTTTGTAAACCTTTCGccccgccggacattgg                                                                                                                                                   >  1:822744/1‑226 (MQ=255)                                                                                                  aaaTGGCGGCGCCACCAAGGGAACCATTAAGGTAGTGATTAGCATCACCTATACCTACAGCGGAACCCGAAGAGAGTATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGAAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑At‑gt                                                                                                                                                   <  1:706400/235‑1 (MQ=255)                                                                                                        cGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCGGCGCACATGGTTATTTAAACATTGCGCCCCTCGTGAGAT‑GGGGGGaaag                                                                                                                                           >  2:73453/1‑237 (MQ=255)                                                                                                          gAGCCACTCAGGGAACAATACGGGCAGTATTTAGAATCCCCTTTCCCTACAGGGGACCCGGAAGAGTTGATTGTAATGAAAC‑GAGTTTTTCCCCTG‑TTTGCTG‑TATGGCTGAG‑GGGCGGGTGGGTAAATCCTTGGTC‑ATGCGCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAAcc                                                                                                                                        <  2:149344/238‑1 (MQ=255)                                                                                                           aGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GGACCGCCAATGTTTTTGTAAAACTTTCGCCCGTCGTGG‑AT‑GTGGGGGACAAcccg                                                                                                                                      >  1:667426/1‑237 (MQ=255)                                                                                                                        cccATTCAGGCAGTGTTTAGCACCACTTATCCCTACAGCTGAACCGGAAGAGGTGATGGTAATGAAGG‑GAGTTATTACCCTT‑TTGGCTG‑TACTGCTGGT‑GGGCGGGCCGGTAAATGCCTGGTC‑TGTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCttt                                                                                                                         <  2:699754/238‑1 (MQ=255)                                                                                                                           atgcAGGCAGTGTTTCGCATCCCCTATCCCTACAGCTGACCCCGAAGAGATGATTGTAATAACAC‑GAGTTATTCCCCTG‑TTTGCTG‑TACTGCTGGT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGa                                                                                                                      <  1:460082/236‑1 (MQ=255)                                                                                                                                  cAGTGATTAGCATCCCCTATCCCTACAGATGAACCCGAAGAGATTATTGTAATGAAAC‑GATTTATTACCTTT‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTGGGTAAATGCGTGGTC‑ATTCGCTTATAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCa                                                                                                                  <  1:1725/236‑1 (MQ=255)                                                                                                                                   aGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCCCCCAAGTTTTTTTAAAACTTGTGCCCCTCGTGT‑AT‑GTGGGGCAAAAACCGGTCGGGGGACTTTCGGCGCaca                                                                                                                >  1:388634/1‑235 (MQ=255)                                                                                                                                       aTTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTACCCCTATTGGCGGTG‑GCAGCGCCAATGTTTTTTTAAACCCTTCGCCCCCCCTGA‑AT‑GTGGGGGCAAACATGGTGGGGGGTCTTTctgcgcccaactttt                                                                                                          >  2:662288/1‑225 (MQ=255)                                                                                                                                             ataccctttcccTACAGCTGAACCCGAAAGGAGGATTGTAATGAAAC‑AGGTTTTTCCCCTG‑TTTGCTT‑AACTGCTGAT‑GGGTGGTTCGGTAACGCCGTGGGC‑CTGCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCAt                                                                                                     <  2:1002446/228‑1 (MQ=255)                                                                                                                                               cACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GGACCGCCCAATTTTATGTAAAACCTTCGCCCGTCGTGG‑AT‑GTGGGGGAAAACCTTGGGGTGGGTTTTTTGGCCCCAATCTTTTTCCATAAc                                                                                                  >  1:819565/1‑239 (MQ=255)                                                                                                                                                      aCCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCCCCCATGTTTTTGTAAACCTTGTGCCCCTCGTGG‑AA‑TTGGGGGAAAACCTGGTCGGTGATCTTTTGTCCCAAATCCTTTTCCATAACCATTAt                                                                                            >  1:852389/1‑238 (MQ=255)                                                                                                                                                          cctgtggaccccGAGGGGATGTTTGTAATGAAGC‑GAGTTATTACCTTG‑‑TTGCGGTTACTGCGGAG‑GGGCGGGCCGGTAAATGCC‑GGTCAATGGGCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑Ga                                                                                      <  2:388634/230‑1 (MQ=255)                                                                                                                                                          aCAGCTGACCCGGAAGAGTTGATTGTAATGAAAC‑GATTTATTACCCTG‑TTGGTTT‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑Ga                                                                                      <  1:73453/239‑1 (MQ=255)                                                                                                                                                          aCAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCTACGCAAATCTTTTGCCATAAAGATTATCCG‑Ga                                                                                      >  1:124190/1‑239 (MQ=255)                                                                                                                                                          aCAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGCCGTGG‑AT‑GTGGGGGAAAAACTGGTCGTGGGTCTTTCGACGCAAATCTTTTGCCATCACGATTATTCG‑Ga                                                                                      >  1:300757/1‑239 (MQ=255)                                                                                                                                                           cAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTAGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTC‑GCAGCCCCAATGTTTATGTAAACCTTGCCCCCGCCGTGG‑ATGGGGGGGCCAAAAGT‑GCCGTGGAATTTTCCACCCCAATCTTTTGCCCTAACGATTAtactgca                                                                                      >  2:326271/1‑233 (MQ=255)                                                                                                                                                            tgatgaccgcgaggggatgatTGGAATGAACA‑GATTTATTCCCTTT‑TTGGCTG‑ATCTGCGTCG‑GTGTGGGTGGGTAAGCGCCGGCTC‑TTGCGCCTGTAAGACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑Gaaa                                                                                    <  2:677940/224‑1 (MQ=255)                                                                                                                                                            aGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAGTGGTACCGCTATCCCTATTGGCGGCG‑GCGGCCCCAAAGCTTATGTTAAACCTGCACCCGTCGTGA‑AT‑TTGGGGCGAAAAATGGTCGTGGGACTTTCGCCGCCAACAtttttgccttacagtctctcggcaa                                                                                     >  2:812866/1‑218 (MQ=255)                                                                                                                                                            aGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTT‑GCAGGCCCACTATTTTTTTAAACCCTGCGCCCCCCGTGA‑AT‑TTTGGGCGCAACCCTGGCGGCGAGCTTTTGCAGCCACTCTTTTTCCACAACGATTTTGCG‑Gtaa                                                                                    >  2:27668/1‑239 (MQ=255)                                                                                                                                                            aGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GGAGCGCCCATGTTTTTGTAAACCTTGCGCCCGTCGTGA‑AT‑GGGGGGGAAAACCCGGTCGGGGGTCTTTTGGAGCAAATCTTTTTCCCCAACGATTATCCG‑Gaa                                                                                     >  1:324848/1‑238 (MQ=255)                                                                                                                                                            aGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGGG‑GCAGCGCCCATGTTTTTGTAAAACTTGCGCCCGTCGGGA‑AT‑GTGGGGCAAAAACTGGGCGGGGATCTTTTGACCCAAATCTTTTTCCATAAAGGTTTTTCG‑Gaaa                                                                                    >  1:186021/1‑239 (MQ=255)                                                                                                                                                              tggcaCCGGAAGGGATGATTGTAAGGAAGC‑GGGTTATTCCCCTG‑TTGGTTT‑TACTGCTGAT‑GGGCTGGTCGGAAAATGCCTGCTC‑TTCCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAAc                                                                                   <  2:667426/234‑1 (MQ=255)                                                                                                                                                               tGACCCCGAAGAGAGGATTGTAATGAAAC‑GATTTTTTACCCGG‑TTTGCTG‑TACTGCAGAT‑GGGTGGGTGGTAAAAGGCGTGGTC‑TTTCGCCTGTAAATCCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAAc                                                                                   <  1:970014/237‑1 (MQ=255)                                                                                                                                                               ggAACCGAAAGGGATGTTTTTAATGAAA‑‑GTGTTTTTCCCCTT‑GTTGGTG‑TTCCTGTGGTGGGGCTGGTGGGTAAATGCCTGTTC‑ATCCCCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAAc                                                                                   <  2:819565/236‑1 (MQ=255)                                                                                                                                                                 aaCCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCATCGCCAATGTTTATTTCAATCCTGCGCCCGTCCGGG‑AT‑GGGGGGCGAAACCTTGGCGGGGGGCTTTCCACCCCACTATCTTGCCATCACGATTCTCCGCGAACCCAt                                                                                >  2:609001/1‑239 (MQ=255)                                                                                                                                                                   cccGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCCATGTTTATGTAAAACTTGCGCCCGTCGTGG‑AT‑GTGGGGCAAAAACTGGGCGTGGGTCTTTTGACGCAAATCTTTTTCCCTAACGAATATCCG‑GAAAACATTAc                                                                             >  1:570886/1‑239 (MQ=255)                                                                                                                                                                      gAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTACCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAACCATTACAGa                                                                          <  1:662288/239‑1 (MQ=255)                                                                                                                                                                       gagagagagttGTAATGAAAC‑GCGTGATTCCCCTG‑TTGGCTG‑TACTGCGGAT‑GGCGTGGGCGGAAAGGGTGTGGTC‑ATGCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAACCATTACAGa                                                                          <  2:807157/229‑1 (MQ=255)                                                                                                                                                                            atgatTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTT‑GCCGCGCCAATGTTTATGTAAACCTTCCCCCCGTCGTGG‑AT‑TTGGGGCAAAAACCAGTCGTTGCACTTTCTCCCCCACTCTTTTGCCAGAACGATAATcccgcacacccacctcac                                                                          >  2:229685/1‑218 (MQ=255)                                                                                                                                                                             tgatTGTAATGAACC‑GGGTTATTCCCCTT‑TTTGCTG‑TACTGCTGAG‑GGGCTGGTCGGTAAATGCTTGGTC‑ATTCGCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAACCATTACAGACTATg                                                                     <  1:304603/237‑1 (MQ=255)                                                                                                                                                                              gatTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTT‑GCAGCCCCAATGGTTATGGTAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGGAAAAACTGGTCGTTGGTCTTTCGGCGCAAAACTTTTTTCCTAACGGTTATCCC‑GCAACACATACAGAATATGTc                                                                   >  1:984327/1‑238 (MQ=255)                                                                                                                                                                               gttGTAA‑GAAAA‑GGGTTTTTTCCCTGTTTTGCTG‑TACCGCTGTT‑GGGGGGGGCGGTAATTGCTGGGTC‑TTGCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAACCATTACAGACTATGTca                                                                  <  2:933897/237‑1 (MQ=255)                                                                                                                                                                                 tGTAATGAAGG‑GTTTTATTCCCTTT‑TTTGCTT‑AACTGCGGGG‑GTGCTGGTCGCTAATTCCCTGGTC‑TTCCCCCTGTAAACCCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAACCATTACAGACTATGTCacac                                                               <  2:570886/239‑1 (MQ=255)                                                                                                                                                                                      tGAACC‑GAGTTATTACCCTG‑TTTGCTG‑TACGGTTGAG‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCATGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAACCATTACAGACTATGTCACACTGCaa                                                          <  1:240853/239‑1 (MQ=255)                                                                                                                                                                                               tttattcCCCGG‑TTTTCGG‑TATTGCGGAT‑GGGCGGGCCGGTAAATGCGTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAACCATTACAGACTATGTCACACTGCAACGAGGCt                                                   <  1:1183709/237‑1 (MQ=255)                                                                                                                                                                                               tttattaCCCTG‑TGTGCGG‑TACTGCTGAG‑GGGCTGACAGTTAAATGCCGGGTC‑GTTCGCCTGTAAAGCCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAACCATTACAGACTATGTCACACTGCAACGAGGCt                                                   <  1:326271/237‑1 (MQ=255)                                                                                                                                                                                                            tttGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCCTGCGCCCCTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAAAGATTATCCG‑GGAAACATTACAGAATATGTCACACTGCAACGAGGCTCTGCTTATg‑gcgg                                     >  1:390692/1‑239 (MQ=255)                                                                                                                                                                                                            tgttgtggtctggctgAG‑GGGATGGTCGTTAAATGCCTGGTC‑GTTCGCCTTTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTGATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATg‑gcg                                      <  1:1156158/227‑1 (MQ=255)                                                                                                                                                                                                              tttgctttatgggtggggggTTGGTCGGTAAATTCC‑GGTCTATTCGCTTGTAAATCCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATg‑gcg                                      <  1:159873/221‑1 (MQ=255)                                                                                                                                                                                                                  tcgtaggcggAT‑GGGGTGTCCGGTAAATGCGTGTGC‑AGCCGCTTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATG‑GCGGCgtg                                 <  2:984327/229‑1 (MQ=255)                                                                                                                                                                                                                               gggCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCATGGCCCATGTTTATGTAAAACTTACGCCCGTCGGGG‑AT‑GTGGGGGCAAACATAGGTGTGGGTCTTTCGTCCCCAATCTTTTTCTATAACGATTTTTCG‑GCAACACTTTCAGACACTGGCTATATGCAACAAGGCTCGGCCTTTGTGCGGCG‑GTTATATTGttttt                    >  2:805/1‑239 (MQ=255)                                                                                                                                                                                                                                    ggTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATc                                                                                                                                                                                                              >  1:24029/1‑53 (MQ=255)                                                                                                                                                                                                                                    ggTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCCCCAATGTTTTTGTAAACCCTTCCCCCCTCGGGA‑AT‑GTGGGGGCAAACCTGGTTCTGGGTCTTTCGACGCCAACCTTTTGCCCTAACGATTATCCG‑GGAACCATTACCAGCTATTTCACCCTGCCACACAGACCGGCCTCTG‑GCGGGGGGTTAtttttttttttccgg               >  2:911247/1‑234 (MQ=255)                                                                                                                                                                                                                                      tcgtgaaaTGCTTGGTG‑ATCCGCCTCTAAAACCGCCAATGGTATCGCTAGCCCTATTGGCGGTG‑GCAGCTCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTGTCGACGCAAATCTGTTGCCATAACGATTATCCG‑GAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATG‑GCGGCGTGTTATCTAATTTTTCCGGGa              <  1:220825/233‑1 (MQ=255)                                                                                                                                                                                                                                      tcctgatcgcccTGGTC‑ATTCCCCTTAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTT‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGCCAAAGCCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATG‑GCGGCGTGTTATCTAATTTTTCCGGGa              <  1:417292/228‑1 (MQ=255)                                                                                                                                                                                                                                      tCGTTAAATTCTGGGTT‑ATTCGCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATG‑GCGGCGTGTTATCTAATTTTTCCGGGa              <  1:484003/238‑1 (MQ=255)                                                                                                                                                                                                                                      tCGGTAAATGCCTGGTC‑ATTCGCCGGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATG‑GCGGCGTGTTATCTAATTTTTCCGGGa              <  1:27668/238‑1 (MQ=255)                                                                                                                                                                                                                                      gcGGAAAATGCCTGGTC‑AGTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATG‑GCGGCGTGTTATCTAATTTTTCCGGGa              <  1:609001/237‑1 (MQ=255)                                                                                                                                                                                                                                      gaggTAAATGCCTGGTC‑ATTGGCCTTTAAACCCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATG‑GCGGCGTGTTATCTAATTTTTCCGGGa              <  1:655160/236‑1 (MQ=255)                                                                                                                                                                                                                                      ccgttAAGTGCTTGGTC‑ATTCGCCTGGAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATG‑GCGGCGTGTTATCTAATTTTTCCGGGa              <  2:324848/234‑1 (MQ=255)                                                                                                                                                                                                                                                 cTGGTC‑ATTCGCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATG‑GCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAtata  <  1:80539/239‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     GTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAAC‑GAGTTATTACCCTG‑TTTGCTG‑TACTGCTGAT‑GGGCTGGTCGGTAAATGCCTGGTC‑ATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTG‑GCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGA‑AT‑GTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCG‑GAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATG‑GCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATA  >  NC_000913/4548629‑4549103

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