Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,197	G→A	P323P (CCG→CCA)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,197	0	G	A	98.7%	255.9 / NA	80	P323P (CCG→CCA)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (34/44);  total (36/44)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAA  >  NC_000913/4618949‑4619443                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       gACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAGGCGGTGCTGGACAACGGGAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCccggg                                                                                                                                                                                                                                                      >  1:47426/1‑249 (MQ=255) gACATGAATCAGGTACTGGCCGCCAGGGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCccg                                                                                                                                                                                                                                                        <  2:275355/249‑1 (MQ=255) acaaTGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCccg                                                                                                                                                                                                                                                        <  2:47426/246‑1 (MQ=255)      gAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACt                                                                                                                                                                                                                                                 <  1:109937/251‑1 (MQ=255)      gAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACt                                                                                                                                                                                                                                                 >  1:480116/1‑251 (MQ=255)           aGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCATCCGACTTCGtt                                                                                                                                                                                                                                            <  1:577704/251‑1 (MQ=255)               aCTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTgaga                                                                                                                                                                                                                                        <  1:196144/251‑1 (MQ=255)               aCTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGAAGTTCCTGCCGGGTGAGTATCGTAACCCGCGTCTGTTTGCTGTCACGATGGCGCTGTGCGGGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGAAGCCGAAGCGCGAGAGAAATTGACAG‑CAATGCTGGACACCGGGAAAGCGGCAGAAGTCT‑TTGGGCGTATGG‑AGACGGCACAAAAAGACCCAACCGAATTCGTTgtg                                                                                                                                                                                                                                         >  2:406726/1‑248 (MQ=255)                        cAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTAcgcg                                                                                                                                                                                                                                >  1:48413/1‑250 (MQ=255)                         aGTGCAGGTAACGCGGTTGACGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCAGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTACGGCAAACTGGCGAAAGATGACGCCGAAGCGAGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGACAAAAGACT‑TTGATCGTATGGTAA‑CGGCACAAAAAGGCCCAACAGACTTCGtgtaggactagtcga                                                                                                                                                                                                                               >  2:327934/1‑236 (MQ=255)                         aGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGa                                                                                                                                                                                                                               >  1:135213/1‑250 (MQ=255)                           tGCAGGTACCGCGGTTGACGTTAGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGATGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTAT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGAATTAGTTGAGAACTACGCGAAgg                                                                                                                                                                                                                            >  2:320613/1‑250 (MQ=255)                           tGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGt                                                                                                                                                                                                                            >  1:176690/1‑251 (MQ=255)                              ccgattgtcgGTTGAAGTTCTGGAAGCGGTGCAGTTCCGGACGGGGGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATc                                                                                                                                                                                                                         <  2:601410/243‑1 (MQ=255)                                  aaCGCGGTTGAAGTTCGAGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGcc                                                                                                                                                                                                                     <  1:493174/251‑1 (MQ=255)                                   aCGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCg                                                                                                                                                                                                                    <  1:608535/251‑1 (MQ=255)                                     gcgGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATc                                                                                                                                                                                                                         >  1:601410/1‑244 (MQ=255)                                      cgGTTGAAGTTCGTGAAGCTGGGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCGGTTTGATGTCACGATGGCTCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGTGAAATTG‑CAAGCAGTGCTGGACAACGGTAAATCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACa                                                                                                                                                                                                                 <  2:675936/251‑1 (MQ=255)                                      cgGTTGAAGTTCGTGAAGCGGAGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGACTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGAAGAAGTAT‑TTGGTCGTAAGATAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAATGATCTGCCGAcc                                                                                                                                                                                                                 >  2:630675/1‑250 (MQ=255)                                                 cGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGAc                                                                                                                                                                                                      >  1:665587/1‑251 (MQ=255)                                                 cGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGAc                                                                                                                                                                                                      <  1:204290/251‑1 (MQ=255)                                                          gTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGt                                                                                                                                                                                             <  1:660192/251‑1 (MQ=255)                                                          gTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGATTATCGCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCa                                                                                                                                                                                               <  2:204937/249‑1 (MQ=255)                                                              agttccgtaCGGTTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGGTGGGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGGCTAt                                                                                                                                                                                         <  2:504690/243‑1 (MQ=255)                                                              aGTTCCTGACGGGAGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGCAGTATCTGCCGACCGCGATGCTGACGAAAGCAGTCTAt                                                                                                                                                                                         >  1:354986/1‑251 (MQ=255)                                                                     gACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGAt                                                                                                                                                                                   <  1:623405/250‑1 (MQ=255)                                                                      aCGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATAc                                                                                                                                                                                 <  1:546532/251‑1 (MQ=255)                                                                      aCGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCACCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATAc                                                                                                                                                                                 >  2:596344/1‑251 (MQ=255)                                                                                atcgataCCCGCGTCCGTTTGATGTCCCGATGGCGTTGTGCGTGGAGATGCTTATCTTCGGCAAACTGGCGAAAGATGACGCCGAAGCTCTCTCTAAAGTG‑CAAGCAGTGCTGGCCAACGGTAAAGCGGCATTAGTCT‑TTGGTCGTATGGTTT‑CGGCACAACAAGGCCCAACCGACTTCGTTGAGATCTACGCGAAGTATCTGCCGACAGCGATGCTGACGAATGCAGCCTATGCTGATACCGAAGTtttt                                                                                                                                                                       <  2:508597/245‑1 (MQ=255)                                                                                     accccGCGTCTGCTTGAAGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACGGGCGAAAGATGCCGCCGAAGCGCACGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGCACCAACCGACTTCTATGAGAAATAAGCGAAGTATCTGCCGACGGCACTGCTGACGACAGCAGTCTATGCTGCTACCGACGGTTTTGTc                                                                                                                                                                    >  2:122481/3‑249 (MQ=255)                                                                                     aacccttgTCTGTTTGATGTCACGAGGGCGCTGTGCGTGGAGATGCTTATCTCCTGCAATCTGGCGAAAGATGACGCCGAAGCGCGCGCGCAATTG‑CAAGCAGTGCTGGCCAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGTCTTCGTTTAGAACTACGCGAAGAAACTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCa                                                                                                                                                                   <  2:181289/243‑1 (MQ=255)                                                                                       cccGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTg                                                                                                                                                                <  1:262608/251‑1 (MQ=255)                                                                                        ccGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGa                                                                                                                                                               <  1:257248/251‑1 (MQ=255)                                                                                            gTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTAc                                                                                                                                                                                                                                   >  2:232524/1‑179 (MQ=255)                                                                                            gTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTAc                                                                                                                                                                                                                                   <  1:232524/179‑1 (MQ=255)                                                                                                  ttGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAGGGTTTTGTCAGTGAAATGGATAcc                                                                                                                                                     >  2:300743/1‑251 (MQ=255)                                                                                                  ttGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATAcc                                                                                                                                                     >  1:592256/1‑251 (MQ=255)                                                                                                  ttGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATAcc                                                                                                                                                     >  1:252251/1‑251 (MQ=255)                                                                                                    gATGTCACGATGGTGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGATGTATCTGCCGACAGCGATGCTGACGAAAGCAGGCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGAGACCcg                                                                                                                                                   <  2:373146/251‑1 (MQ=255)                                                                                                       gTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGCCGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGaa                                                                                                                                                                                                                              >  2:650538/1‑173 (MQ=255)                                                                                                       gTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGaa                                                                                                                                                                                                                              <  1:650538/173‑1 (MQ=255)                                                                                                                                  gccgaactcGGCAAACGTGCGAAAGATGCCGCCGAAGCGCGCGCGAAATTG‑CAAGCAGGGCTGGACAACGGTAAAGCGGCAGAAGTCG‑TTGGTCGGATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCaaggg                                                                                                                      <  2:302157/242‑4 (MQ=255)                                                                                                                                   ttaTCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCTATATTG‑CAAGCAGTGCTGGCCAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGACCTACGCGCAGTATCGGCAGACAGCGCTGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGTAGGCAGTGGTTGCAATgggg                                                                                                                    >  2:101410/3‑250 (MQ=255)                                                                                                                                     aTCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGCTGGCAGTGGTGGCAAGGggcgg                                                                                                                  >  2:259137/1‑251 (MQ=255)                                                                                                                                      tctcCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGt                                                                                                                                                                                             >  1:98937/1‑175 (MQ=255)                                                                                                                                      tctcCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCG‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGt                                                                                                                                                                                             <  2:98937/175‑1 (MQ=255)                                                                                                                                         ccGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGggcgga                                                                                                                 <  2:559928/248‑2 (MQ=255)                                                                                                                                         ccGGCAAACGGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGcggcggt                                                                                                              <  2:274604/251‑2 (MQ=255)                                                                                                                                          ctgcaaacaagCGAAGGATGACGCCGAAGCGCGCGCGAAATTA‑CAAGCAGTGCTGGACCACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGGAG‑CGGCACAAAAAGGCCCCACCGACTTCGTTGAGACCTACGCGCAGTATCTGACGACATCGATGCTGACGAAAGCAGTCTATGCTAATACCGAAAGGTTTGTCCGTGTAATGTATACCCGCGCGCTGGGGATGACAGTGGTTGCAATGGTcggcgggt                                                                                                             >  2:467374/11‑249 (MQ=255)                                                                                                                                          cGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATgg                                                                                                                      >  1:302157/1‑242 (MQ=255)                                                                                                                                             cAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGcggcggtcgtc                                                                                                          <  1:112782/251‑6 (MQ=255)                                                                                                                                                   ggCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCTCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAgg                                                                                                    <  1:521550/251‑1 (MQ=255)                                                                                                                                                      gAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCAt                                                                                                 <  1:300743/251‑1 (MQ=255)                                                                                                                                                            ggACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGaca                                                                                           <  2:176690/250‑1 (MQ=255)                                                                                                                                                                 ccGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCAtt                                                                                      <  2:171774/251‑2 (MQ=255)                                                                                                                                                                 ccGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGACC‑ATGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGCGAACTACACGGCGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTAAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGCCGTCGTCAGGCATCTGACACCAtt                                                                                      >  1:313234/1‑250 (MQ=255)                                                                                                                                                                           cgcgAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCg                                                                            <  1:6286/251‑1 (MQ=255)                                                                                                                                                                              gAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCg                                                                         >  1:84922/1‑251 (MQ=255)                                                                                                                                                                               aaaTTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCgg                                                                        <  1:74842/251‑1 (MQ=255)                                                                                                                                                                                aaTTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGc                                                                       <  1:96923/251‑1 (MQ=255)                                                                                                                                                                                       aagcagTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAACAGGCCCAACCGACTTCGTTGCGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTAt                                                                                                                                                                                         >  2:243456/6‑131 (MQ=255)                                                                                                                                                                                       aagcagTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTAt                                                                                                                                                                                         <  1:243456/126‑1 (MQ=255)                                                                                                                                                                                       aagcagTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCtttcct                                                                 >  2:190777/6‑248 (MQ=255)                                                                                                                                                                                           tgTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAACAAGACCCGACCGACTTCGTTGAGCACTACGCGAAGGATCTGCAGACAGCGATGCTGACGAAAGAAGTCTATGCTAATACCGAAGGGTTTGTCAGTGAAATGGATACCCGCGAGCTGGGGCAGGCAGTGGTTGCAATGGACGTCGGGCGCAGTCAGCCATCTCACACCAGCGGTGACATCGTCGGCTTTACTg                                                                >  2:14032/2‑248 (MQ=255)                                                                                                                                                                                           agTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGata                                                             <  1:415265/250‑1 (MQ=255)                                                                                                                                                                                                tGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGgcg                                                        >  1:484570/1‑251 (MQ=255)                                                                                                                                                                                                 ggACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGgcg                                                        <  2:135213/250‑1 (MQ=255)                                                                                                                                                                                                 ggACAACGGGAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGgc                                                         <  1:404594/249‑1 (MQ=255)                                                                                                                                                                                                    aaaCGGTAAAGCGGCAGAAGTCT‑TTGGTCGTTTGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGtt                                                    <  2:592256/250‑2 (MQ=255)                                                                                                                                                                                                      aCGGTAAGGCGGCAGAAGTCTGTTGG‑CGTATGGTAG‑CGGCACAAAAAGGCCCAGCCGACTTCGTTGAGAACTACGCGAAGTATCTGTCGACAGCGATGCTGACGAAAGCAGGCGATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGGCACCATTGATTACAGCGTCGGTTTTACTGATATGGGGCGtttg                                                  <  2:84922/251‑4 (MQ=255)                                                                                                                                                                                                      aCGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGtttg                                                  >  1:485386/1‑248 (MQ=255)                                                                                                                                                                                                            aaGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGGGGTTGCAATGGGCGGCGGACGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACGGATATGGCGCGttttggagg                                             >  2:137083/1‑242 (MQ=255)                                                                                                                                                                                                             aGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGAcc                                           >  1:675663/1‑251 (MQ=255)                                                                                                                                                                                                                               gTCGTATGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCgg                          >  1:651206/1‑250 (MQ=255)                                                                                                                                                                                                                                 cGTATGGTAG‑CGGCACAAAAAGGCCCAACCGTGTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCAGATACCGAAGGTTTTGGCAGTGAAATGGATACCCGCGCGCTGGGGATGTCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTc                         <  2:416079/250‑1 (MQ=255)                                                                                                                                                                                                                                    aTGGTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGAAACCCGCGCGCTGGGGATGGCAGGGGTTGCAATGGGCGGCGGTCGTCGTCCGGCATCTGACACCATTGTTTAACGAGTCGGCTTTACTGTTATGGCGCGTGTGGGCGTCCAGGAAGAAGGTCAGCgccagct                     >  2:487495/1‑245 (MQ=255)                                                                                                                                                                                                                                      ggTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGAAACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTgg                   >  1:250084/1‑251 (MQ=255)                                                                                                                                                                                                                                       gTAG‑CGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGGTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGc                  <  2:480116/251‑1 (MQ=255)                                                                                                                                                                                                                                             ggCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCgggtt             >  1:289460/1‑248 (MQ=255)                                                                                                                                                                                                                                              gCACAAAAAGGCCCAACCGACTTCGATTAGAACTACGCGAAGTATCTGCCGACAGCTATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCGTTACTGAGATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAt            <  2:330476/251‑1 (MQ=255)                                                                                                                                                                                                                                                 caAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCa         <  1:411152/251‑1 (MQ=255)                                                                                                                                                                                                                                                  aaaaaGGCCGAACCGGCTTCGTTGAGAGCTCCGCGAAGTAGCTGCCGACGGCTATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGAGCTGGGGAGGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCAc        <  2:133928/251‑1 (MQ=255)                                                                                                                                                                                                                                                        gCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGaaa  <  1:406726/251‑1 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       GACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCT‑TTGGTCGTATGGTAG‑CGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAA  >  NC_000913/4618949‑4619443

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