Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,514,287	G→C	L8V (CTT→GTT)	fecB ←	ferric citrate ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,514,287	0	G	C	99.0%	317.4 / NA	96	L8V (CTT→GTT)	fecB	ferric citrate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (36/59);  minor base A (1/0);  total (37/59)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGT‑CATACGCCACGCGTGC  >  NC_000913/4514043‑4514523                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              tGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCAAGCGCATCGGCGACCGAGGGTTCAAGCACCACAATCCGTTGTGGCGTTTTTTCGGGTGTAAACGTGCCGTGT‑TCGTC‑CTGCACCGTG‑GCGGCAAAGGCAGGGGTGATCCCCAGCAGCAGGGCTGCAAAAACaaac                                                                                                                                                                                                                                           >  2:556956/1‑248 (MQ=255)      gcgcgTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATaa                                                                                                                                                                                                                                    <  1:49902/251‑1 (MQ=255)           ttCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCACACACAAAACGGATAAATGcc                                                                                                                                                                                                                               >  2:175692/1‑251 (MQ=255)                 cGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGGGGCTGATCACCAGCAGCAGGCCTGCAAAAACACAACGGATAAATGCCAACAt                                                                                                                                                                                                                          >  1:270894/1‑250 (MQ=255)                   gACTGCCACGGTTTCAGGTGCGCATGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAAt                                                                                                                                                                                                                       <  1:395876/251‑1 (MQ=255)                     cggcaccGGTTTCAGGTGCGCTCGCATTTCGGGCAGTATGGGTTTTGCATCGCTATCGTCGGACATGCCGGTCGGTCTGCCGCGCACGGCGGCCAGTGCAGCGGCGAACGAGAGTTCCAGCACCACAATCCGGTTTGGCGTATTGTCGAGAGTATACGTGCCGTGT‑TCGTC‑CGGAACCGTT‑GCGGCAAAGGCGTGGCTGATCCCCAGCAGCTGGCCTGCATATACAATACGGATAAAAGCCAACATAATca                                                                                                                                                                                                                     <  2:185594/245‑1 (MQ=255)                           cGGTTTCAGGTGCGCACGCACTTCGGGCAGGATTCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCCCGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTcc                                                                                                                                                                                                               <  2:270894/251‑1 (MQ=255)                             gTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAg                                                                                                                                                                                                             <  1:10152/251‑1 (MQ=255)                              cctcctccgatcacgcacCTCGGGCAGGATGCGTTGTGCATCGGTAGCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCTGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGc                                                                                                                                                                                                            <  2:508344/240‑1 (MQ=255)                                   ggTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGCTCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGGGTTCCAGCACCACAATCCGTTGTGGCGTTTTCTCGAGTGTAAACGTGCCGTGT‑TCGTA‑CTGAACCGTG‑GCGGCAAAGGCGTGGATGATCAACAGCAGCAGGCCTGCCAAAACAAACCGGCTAAATGCCAACATAATCAAATTCCAGCTaaa                                                                                                                                                                                                        >  1:653490/1‑250 (MQ=255)                                      gcgcACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGcc                                                                                                                                                                                                    <  2:194049/251‑1 (MQ=255)                                       cgcacgcaCTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAACACCTATCGGGCTGACGGCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGccc                                                                                                                                                                                                   <  2:224372/251‑1 (MQ=255)                                        gcacgcacTGCGGGCAGGATGCGTTTTGCATCGTTATCGGCGGCAATACCGTGCGGGCTGACGTCCACGGCGGCCAGTGCATCGGCGAACGAGAGTGCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCg                                                                                                                                                                                                  <  2:302021/251‑1 (MQ=255)                                           cgcacTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACCTAATCACATTCCAGc                                                                                                                                                                                                            >  1:508344/1‑238 (MQ=255)                                            gcacTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACa                                                                                                                                                                                                                           <  2:148120/222‑1 (MQ=255)                                            gcacTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACa                                                                                                                                                                                                                           >  1:148120/1‑222 (MQ=255)                                             cacTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCAGCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAg                                                                                                                                                                                             <  2:485891/251‑1 (MQ=255)                                                  cGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCAAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATaa                                                                                                                                                                                                                        <  1:100394/219‑1 (MQ=255)                                                   gggCAGGATGCTTTTTGCATCGTTATCGCCGGCAATCCCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCCCCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGc                                                                                                                                                                                       <  2:612033/251‑1 (MQ=255)                                                   gggCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCgggg                                                                                                                                                                                       >  1:295663/1‑250 (MQ=255)                                                          aTGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAaca                                                                                                                                                                                <  1:175692/251‑1 (MQ=255)                                                            gCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAacac                                                                                                                                                                               <  1:332783/250‑1 (MQ=255)                                                             cGTTTTTCAGCGTTATCGTCGGTAATACCGATCTGGCTGACGCCCACGGCGGCCATCTCATCGGCGAACGTTAGGTCCAGCCTGACAATCCGTTGGGGCGTTTGTTAGAGTGTAAACGGGCCG‑GTGTCGAC‑CTGAACCGTG‑TCGGCAAATGCGTGGCTGATCACCATCAGCCGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAATAGCCCGGCATGCCGGGCTTTAacccag                                                                                                                                                                             <  2:382053/251‑5 (MQ=255)                                                                 ttGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGGCCACGGCGGCCAGCGCATCGGCGAACGAGGGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACTTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCCCCAGCAGCAGGCCTGCAAAAACAAACCGGAAAAATGCCAACAAAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGCCGTTAACACATc                                                                                                                                                                            >  1:62499/1‑248 (MQ=255)                                                                    cATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACt                                                                                                                                                                      >  1:602850/1‑251 (MQ=255)                                                                     aTCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACtt                                                                                                                                                                     <  1:507309/251‑1 (MQ=255)                                                                          tATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAAGCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAAc                                                                                                                                                                 <  2:253915/250‑1 (MQ=255)                                                                           aTCGTCGGCAATACCGATTGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACg                                                                                                                                                                >  1:525742/1‑250 (MQ=255)                                                                               tcgGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGAcccc                                                                                                                                                           >  1:35312/1‑251 (MQ=255)                                                                                ccgtccgtgctaTCGGGCTGCCGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGTCCCCt                                                                                                                                                          <  2:194825/240‑1 (MQ=255)                                                                                      tACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGAACCCTGCatat                                                                                                                                                    >  1:602363/1‑251 (MQ=255)                                                                                      tACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGGGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCACAAAAAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGAAAGCCGGGGGTTAACACATCAGAACTTCAACGACCCCTGCatat                                                                                                                                                    >  1:28644/1‑251 (MQ=255)                                                                                       aCCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCatat                                                                                                                                                    >  1:503710/1‑250 (MQ=255)                                                                                        ccGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAACCGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCCAAGGCGTGGCTGATCACCAGCAGCCGGCCTGCAAAAACAAAACGGATAAATGCCAACCTAATCACAATCCAGCAAAAAACCCGGCCAGCCGGGCGTTAACACATCCGAAATTAAACGAACCCTgaaaattc                                                                                                                                                  >  2:371466/1‑244 (MQ=255)                                                                                        ccGATCGGGCTGACGGTCACGGCGTCCAGCGCATCGGCGAACGATAGTTCCAGCCCCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTT‑GCGGAAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGGTAACACATCAGAACTTCAACGACCCCTGCATATAc                                                                                                                                                  <  2:434112/251‑1 (MQ=255)                                                                                          gATCGGGCTGGCGCCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTGTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAg                                                                                                                                                <  2:507925/251‑1 (MQ=255)                                                                                          gATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATAc                                                                                                                                                  >  1:558310/1‑249 (MQ=255)                                                                                           agcggactGACGTCCTCGGCGGCCAGCGTATCGGCGAACGAGAGTTCCAGCACCACAATCCGGTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑GCGTC‑CTGAACCGTT‑GCGGCAAAGGCGAGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTTCATATACAgt                                                                                                                                               <  2:58930/244‑2 (MQ=255)                                                                                           aTCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCt                                                                                                                                                          >  1:194825/1‑240 (MQ=255)                                                                                           aTCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGc                                                                                                                                               <  1:241526/251‑1 (MQ=255)                                                                                                gctGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTgcg                                                                                                                                          <  1:40197/249‑1 (MQ=255)                                                                                                      gtcaaCGGAGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATACGTTGGGGCGTTTTTTCGAGAGTAAACGTGCCGTCT‑TCGGA‑CTGAACAGGG‑GCGGCAAAGGAGGGGATGATCACCAGCAACAGGACTTCAAAAACAAAAAGGATAAATGCAAAAATAATAAAATTACATATAAAAGCACCGAAAGCACGGCGATAACAAAAAAGCACTACAACGACACCAGCATATACAGCCGGCGCGGcgg                                                                                                                                    >  2:546994/5‑249 (MQ=255)                                                                                                      gTCCCCGGCGGCCAGCTCATCGGCGAACGAGAGTTCCAGCACCCCAATCCGTTGTGGCGGTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGGTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTg                                                                                                                                    <  2:192688/251‑1 (MQ=255)                                                                                                          aCGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCt                                                                                                                                <  1:549890/251‑1 (MQ=255)                                                                                                                gcCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAg                                                                                                                          <  1:467403/251‑1 (MQ=255)                                                                                                                  cagcgcagcgTCGACCGAGGGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑GCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATGCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACAGCAGAACTACAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGAt                                                                                                                        <  2:634743/251‑1 (MQ=255)                                                                                                                   aGCGCATCGGCGAACGAGAGATCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAACCGGGCGTTAACACATCAGAACTTCAACGACACCTGCATATACAGCGTGCGCGGATGGCCTGCATAGATg                                                                                                                       >  1:174919/1‑251 (MQ=255)                                                                                                                     cgcATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGcc                                                                                                                     >  1:86536/1‑251 (MQ=255)                                                                                                                      gcATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCt                                                                                                                    <  1:572793/251‑1 (MQ=255)                                                                                                                       gagcGGCGAGCGAGAGTTCCAGCTCCACAAACCGTTGTGGCGGTTTTTCGAGTGTAAACGTGCCGTGG‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACCCATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCtt                                                                                                                   <  2:24075/248‑1 (MQ=255)                                                                                                                            gcgcaCGAGAGCTCCAGCTCCACAATCCGTTGCGGCGTTTGTTCGAGGTTACACGTGCCGTGT‑TCGTC‑CTGCACCGGG‑GCGGTAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACTGATAAATGCCCACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACGCCTGCATATACATCGTGCGCGGCTGGCCTGCATATATGCCGttgttg                                                                                                              <  2:56306/247‑1 (MQ=255)                                                                                                                            gCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTTCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgtgg                                                                                                              <  2:270015/251‑3 (MQ=255)                                                                                                                              gAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgttgtt                                                                                                            <  1:68949/251‑1 (MQ=255)                                                                                                                                 cGAGAGTTCCAGCATCACAATCCGTTGTGGCGGTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCATAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAACCAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgttgttggc                                                                                                         <  2:466618/251‑3 (MQ=255)                                                                                                                                  gagaGTTCCAGCACCACAATCCGTTGGGGCGTTTTTGCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGACCCGTG‑GCTGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTgtcg                                                                                                        <  1:522813/251‑1 (MQ=255)                                                                                                                                       ttCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACTTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCCGCAGCAGGCCTGCAAAAACAACACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACCCATCAGAACTTCAACGACCCCTGCATCTACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTggcgtctt                                                                                                    >  2:327942/1‑248 (MQ=255)                                                                                                                                           aGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAAAACTTCACCGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTTTCGTCAtcag                                                                                                 >  2:131069/1‑246 (MQ=255)                                                                                                                                            gcaccatcagccGTTGTGGCGTTGTTTCGAGTTTAAACGTGCCGTGT‑TCGTC‑CCGAACCGGG‑GCGGCAAAGGCGTGGCTGATCTCCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGTCTTTGTTGTTGTCGTCATAAGAGc                                                                                              <  2:18990/241‑1 (MQ=255)                                                                                                                                            gCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGGGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGc                                                                                              <  1:286535/251‑1 (MQ=255)                                                                                                                                               ccacAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCgg                                                                                            >  1:377193/1‑250 (MQ=255)                                                                                                                                                cacaATCCGTTGTGGCGGTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTTTCGTCATAAGAGCGGAt                                                                                          <  1:456704/251‑1 (MQ=255)                                                                                                                                                cacaATCCGTGGGGGTGTGTTGTCGGGTGTAAACGTGCCGTGT‑TCGGC‑CTGAACCGTT‑GCGGCAAAGGCGGGGCTGATCTCCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATACCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTGAACTACCCCTGCATATACAGCGTGCGCGGCTTGCCTGCAGAGATGCCGTTGTTGTTGTCGTCTTAAGAGCGGAt                                                                                          <  2:77148/251‑1 (MQ=255)                                                                                                                                                      gcgtagtgGCGATTTTTAGAGTGTAAAAGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACCGATAAATGCCAACATAATCACATTACAGCTAACAGCCAGGCAGGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCACATACAGCGAGCGCGACGGGCCTGCATGGTTGCCTTTGTAGTTgtcgcactaggactggtggtagtg                                                                                    >  2:210111/6‑231 (MQ=255)                                                                                                                                                      ccGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAagta                                                                                    <  1:371466/251‑1 (MQ=255)                                                                                                                                                      ccGTTGGGGCTTTTTTTAGAGTGTAAACGTGCCGTGT‑TCGTCTCT‑AACC‑TGTGCGGCAAAGGAGTGGCTTACCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAACAGCCTGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCTCTGCATATACAGCTTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGAGGAagta                                                                                    <  2:419558/251‑1 (MQ=255)                                                                                                                                                          tgtgGCGTTTTTTCGCGTGTAAACGTGCCGTGT‑TCGTC‑CTGACCCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGAAAAATGCCAACTTCATCACATTCCAACTACAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTTTCGCCCTCAGGACGGAGGAAGTAGTcc                                                                                >  2:481407/1‑251 (MQ=255)                                                                                                                                                             ggCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTgg                                                                             <  1:63243/251‑1 (MQ=255)                                                                                                                                                                  tttttCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGa                                                                         <  1:485059/250‑1 (MQ=255)                                                                                                                                                                    tttCGGGTGTAAACGTGCCGTGT‑TCGGC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAg                                                                       <  2:212026/250‑1 (MQ=255)                                                                                                                                                                       cGAGTGTAAACGTGCCGGGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCGGACCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGGGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGGCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGt                                                                   <  2:52861/251‑1 (MQ=255)                                                                                                                                                                       cGAGGGTAAACGTGCCGGGG‑TCGCC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCATGACTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTGTGTTGTTGTCGTTATAAGCGCGGATGAAGTAGTCCTGGTCGAAGATGt                                                                   <  2:402917/251‑1 (MQ=255)                                                                                                                                                                            gtAAACGGGCCGTGT‑TCGTC‑CCTAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCATCAGCAGGCCTTCAAAAACAAAACGGATAAAGCCCAACAGAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCATACCTTCAACGTCCCCTGCCTATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTTTCGTCATAAGAGCGGATGAAGTAGTCCTTGTCGAATATGtttgtc                                                              <  2:28644/251‑4 (MQ=255)                                                                                                                                                                               aaCGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGGAGGAGTCCTGGTCGAAGATGTTTTTcac                                                            >  2:71602/1‑250 (MQ=255)                                                                                                                                                                                 cGTGCCGTGT‑TCGTC‑CTGAACCGTG‑ACGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACAcc                                                         >  1:546264/1‑251 (MQ=255)                                                                                                                                                                                      cGTGT‑TCGTC‑CTGAACCGGG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATg                                                    <  2:558238/251‑1 (MQ=255)                                                                                                                                                                                         gt‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCa                                                 <  2:602363/251‑1 (MQ=255)                                                                                                                                                                                          t‑tCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATAAAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGTTGAAGTAGTCCTGGTCGCAGATGTTTTTAACACCGAATGcc                                                  >  2:547560/1‑249 (MQ=255)                                                                                                                                                                                                  cTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTc                                            >  1:46651/1‑249 (MQ=255)                                                                                                                                                                                                      aCCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAATCCCGGCAAGCCGGGCGTTAACTCATCAGAACTTCAACGACCCCTGCATATACATCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTTGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGTTTCAGATCt                                      <  2:497464/251‑1 (MQ=255)                                                                                                                                                                                                           ttgcggcAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATCAATGACAACATAATAACATTCCAGCTAAAAGCCCGGCAACCCGAGCGATAACACATCAGACCTTCCACGAACCCGGCATATACAGCGAGCACAGCTGGCCAGAATAGATGCCTTTGTTGTTGTCTTCCTAAGAGGTGAGCAAGACGGCCTGGTCGAGGATGTTTTTAGCACCGAATGCCAGGTTCATATCttccc                                  >  2:72618/3‑247 (MQ=255)                                                                                                                                                                                                              cggcAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCt                               <  1:536251/251‑1 (MQ=255)                                                                                                                                                                                                                 cAAAGGCGTGGCCTATAACCAGCAGCAGGACTGCAAAAACAAAACGGATAAAAGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGTCCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCg                            <  2:222518/251‑1 (MQ=255)                                                                                                                                                                                                                       cGTGGCTGATTAACAGCAGCAGGCCTGCAAAAACAAAACGGATAAAGGCCAACATAATCACAGTCCAGCTAAAAGCCCGGCAAGCCGGGCGGTAACACAGTAGAACTTCAACGACCCCTGCATAGACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGATGAGGGTGTTCACACCGAATTCCAGGTTCAGATGTGCCATCTGGTGGCCaa                      <  2:448818/251‑1 (MQ=255)                                                                                                                                                                                                                        gTGGCTGCTCCCCAGCAGCAAGACCGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGATAAAAGACCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCAGATGAAGTAGTCCTGTCCGAAGATGTTTTTCACACCGCATGCCAGGTTCAGATCTGCCCTCTGCGGGCCaa                      >  2:74347/1‑250 (MQ=255)                                                                                                                                                                                                                        gTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAg                                                                                                                                                <  1:163308/128‑1 (MQ=255)                                                                                                                                                                                                                        gTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTACCGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAg                                                                                                                                                >  2:163308/1‑128 (MQ=255)                                                                                                                                                                                                                         tGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCaa                      >  1:423742/1‑249 (MQ=255)                                                                                                                                                                                                                           gCTGATCACCAGCAGCAGGACTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGT‑c                 <  1:626968/251‑1 (MQ=255)                                                                                                                                                                                                                             tGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGAAGTCCTGGTCGAAGATGTTTTTca                                                             >  1:161333/1‑206 (MQ=255)                                                                                                                                                                                                                              gATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGT‑CATa              <  1:546994/251‑1 (MQ=255)                                                                                                                                                                                                                               atccccAGCAGCTGGCCTGCAAAAACCAAGCGGATTAATGCCATCATATTCACAATCCAGCTAAAATCCCGGCATGTCGGGCTTTACCACAACAGAACTTCAACGGCTCCTACATATACAGCGTGCGCTGCTGGCCTTCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATTAAGTAGTCCTGGTCTAAGATCGTTATCCCACCGAATACCAGGTTCAGATCTTCCATCTGCGTGCCATAGT‑CATAc             <  2:653490/247‑1 (MQ=255)                                                                                                                                                                                                                                    cagcagcagGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGT‑CATACGCCAc        <  1:356322/251‑1 (MQ=255)                                                                                                                                                                                                                                        agcagGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAagta                                                                                    <  2:217897/172‑1 (MQ=255)                                                                                                                                                                                                                                        agcagGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAagta                                                                                    >  1:217897/1‑172 (MQ=255)                                                                                                                                                                                                                                        agcagGCCTGCAAAAACAAAAAGGATAAATGCCAACATAATCACATTCCAGCTACAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACACCTGCCTATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGCAGATGTTTTTAACACCGAAGGCCAGTTTCAGATCTGCCATCTGCGGGCCAAAGTACA‑ACGCCACGCGt    >  2:434698/1‑251 (MQ=255)                                                                                                                                                                                                                                          cagGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCAAGATAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTGCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGAATAGATGCCTTTGTTGTGGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGGCAGGTTCAGATCTGCCATCTGCGGGCCAAAGT‑CATACGCCAcgcgagc  >  1:176658/1‑248 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              TGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGT‑TCGTC‑CTGAACCGTG‑GCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGT‑CATACGCCACGCGTGC  >  NC_000913/4514043‑4514523

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