Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	581,380	0	C	G	24.0%	180.0 / 81.8	100	intergenic (+1/‑154)	nohD/tfaD	putative DNA‑packaging protein NohD/putative tail fiber assembly protein TfaD
Reads supporting (aligned to +/- strand):  ref base C (31/45);  new base G (13/11);  total (44/56)
Fisher's exact test for biased strand distribution p-value = 3.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.90e-02
Rejected as consensus: Frequency below/above cutoff threshold.

AAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCT  >  NC_000913/581131‑581627                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         aaaaaTGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGGACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACa                                                                                                                                                                                                                                                         <  1:531434/251‑1 (MQ=38) aaaaaTGCCAGAGACGCCGCTGAAGTGGTGGAAACCGCATGCTGTACTTTCGTGCTGTCGCGGGTCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTGTTCCGGAACTGGAAAACCGACAGGTTGATTTCCGGAAACTGGATAGCATCAAAGGCATGAACAAAGCAGCTTCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACa                                                                                                                                                                                                                                                         <  1:76082/251‑1 (MQ=11)  aaaaTGCCAGATGCTCCTCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGGGCTGTCGCGGATCTCAGGTGAAATGGCCAGTATTCGCGACGGGATCCCCCGGTCGGTGCAGCGGCTGTGTCCGGAACTGGAAAACCGAAATGTTGATTGCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTGGCTGAGTGAATATATCGAACAGTCAGTTTAACAg                                                                                                                                                                                                                                                        <  2:482536/251‑1 (MQ=11)       gCCAGAGCGTCCGCTGAATTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATAACGAACAGTCAGGTTAACAGGCTGc                                                                                                                                                                                                                                                   <  2:174827/251‑1 (MQ=255)          aaagattCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGCTGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGcagc                                                                                                                                                                                                                                                >  2:290114/7‑248 (MQ=255)             gaCTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCAtt                                                                                                                                                                                                                                             >  1:470113/1‑251 (MQ=255)                   gCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTgcagcactctgcgc                                                                                                                                                                                                                                       >  2:117763/1‑243 (MQ=255)                      gaatggggggAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTGCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGAGACCGGGGTTGCTGAGTGAATATATCGAACAGTCATGTTAACAGGCTGCGGCATTTTGT‑Ccg                                                                                                                                                                                                                                     <  2:354344/241‑1 (MQ=255)                      gAAGTGGTGGGAACCGCATTCTGTACTTTCGTGCTGTCGCGGAATGCAGGTGAAATTGCCAGTATTCTCGGCGGGAGCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCGGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑Ccgcg                                                                                                                                                                                                                                   <  2:728834/251‑1 (MQ=255)                       aggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGATAAGGGGCTgcagcactctgcgcgcgc                                                                                                                                                                                                                                   >  2:564779/3‑239 (MQ=255)                        gggggtggAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCGGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTgcagcactctgcgcgcgcg                                                                                                                                                                                                                                  <  2:83146/248‑14 (MQ=255)                              ggAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTgcagcactctgcgcgcgcgggtct                                                                                                                                                                                                                             >  1:682108/1‑232 (MQ=255)                               gAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTgcagcactctgcgcgcgcgggtctac                                                                                                                                                                                                                           <  2:245697/251‑21 (MQ=255)                                aaaCCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCg                                                                                                                                                                                                                         <  1:424204/251‑1 (MQ=255)                                         tCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGtt                                                                                                                                                                                                                >  1:534461/1‑251 (MQ=255)                                            gTACTTTCGTTCTGTCGCGGTTCGCAGGTGAAATTGCCATTATTCTCGTCGGGTTCCCCCTGTCGGTGCAGCGGCGTTGTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGtt                                                                                                                                                                                                                <  2:136296/248‑1 (MQ=255)                                            gTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATGCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACGGGAAAACCGACAGGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTc                                                                                                                                                                                                               <  2:730507/249‑1 (MQ=255)                                             tACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAgg                                                                                                                                                                                                            >  1:669894/1‑251 (MQ=255)                                                      gCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCca                                                                                                                                                                                                   <  1:557581/251‑1 (MQ=255)                                                           cgcgGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGAc                                                                                                                                                                                              >  1:435875/1‑251 (MQ=255)                                                            gcgGATCGCAGGTGAAAGTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAAATGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATAACGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTACCAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTGCGCTCACTGTTCAGGCCGGAGACACAGAcc                                                                                                                                                                                             <  2:511180/251‑1 (MQ=255)                                                              ggATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGAccgc                                                                                                                                                                                           >  2:324072/1‑251 (MQ=255)                                                                        ggAAATTGCCAGTATTCTCGACGGGATCCACCTGTCGGGGCAGAGGCGTTTTCCGGAACTGGAAAACCGACAGGTTGATTTCCTGAAACGGGATATCAGCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTGGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGAATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATgg                                                                                                                                                                                 <  2:356224/250‑1 (MQ=21)                                                                           aagtGCCAGTACTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCg                                                                                                                                                                              <  1:507298/248‑1 (MQ=255)                                                                                ccAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTACGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTACGCTCACTGTTCAGGCCGGAGCAACAGACCGCCGTTGAATGGGCGGATGc                                                                                                                                                                         >  2:669312/1‑251 (MQ=255)                                                                                       tctcGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTa                                                                                                                                                                   <  1:324072/250‑1 (MQ=255)                                                                                         ccgccGGGAACCCCCGGTCGGTGCCGCGGCGTTTTCCGGAACGGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGTCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTTAATGGGCGGATGCTAATTACTa                                                                                                                                                                <  2:812104/247‑1 (MQ=25)                                                                                           gACGGGTTCCCCCTGTGGGTGCAGCGGCTTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCAACAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGTGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTAtc                                                                                                                                                              <  2:691629/251‑1 (MQ=255)                                                                                           gACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCt                                                                                                                                                                                                                       <  1:527854/194‑1 (MQ=255)                                                                                           gACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCt                                                                                                                                                                                                                       >  2:527854/1‑194 (MQ=255)                                                                                            aCGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGATTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTGCTAtct                                                                                                                                                             <  1:257214/251‑1 (MQ=255)                                                                                            aCGGGAACCCTCTGTCGGTTCAGCGGCGTTTTCCGGAACTGGAAAACCATCATGTTGATTTCCTGAAACGGGATATCATCAAAGCTATGAAAAAAGCAGACCCGCTGGAAGAACTGATTCCGGTGTTGCGGAGGGAATATATCGAACAGGCAGGTTAACAGGCTGCGGCATTT‑GTGCCGCGCCGGGCTTGGCTCACTGTACAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATTCTAATTACTAtct                                                                                                                                                             <  1:791656/251‑1 (MQ=21)                                                                                               ggATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTccc                                                                                                                                                          >  1:291415/1‑251 (MQ=255)                                                                                                agttcctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGC‑GCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTcccc                                                                                                                                                         >  1:446409/8‑250 (MQ=255)                                                                                                   cccccTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGaa                                                                                                                                                       >  1:346804/1‑250 (MQ=255)                                                                                                         gTCGGTTCAGCGGTGTTTTCCGGAACTGAAAAACCGACATGTTGATTTCCGGAAACGGGAGATCAACAAAGCCATGAATAAAGCACCCGCGCTGGATGAACTGATACCGTGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCTTTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACGATCCCCCGAAAGAAt                                                                                                                                                  <  2:502073/249‑1 (MQ=255)                                                                                                         gTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGC‑GCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGagt                                                                                                                                                  <  1:676012/249‑3 (MQ=255)                                                                                                          gctgTGCAGCGGATTTTTCCGGAATTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGTTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCg                                                                                                                                               <  1:179991/248‑1 (MQ=255)                                                                                                              tGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCAt                                                                                                                                            >  1:644919/1‑250 (MQ=255)                                                                                                                cAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGC‑GCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGCTAATTACTATCTCCCCAAAGagtcggcctaac                                                                                                                                         >  1:549524/1‑240 (MQ=255)                                                                                                                      cgTTTCCCGAAACTGGAAAACCGAAATGTTGATTTCCTGAAACGGGATATCATAAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCGGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAg                                                                                                                                   <  2:633425/251‑1 (MQ=255)                                                                                                                          ttCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGgc                                                                                                                                <  1:219631/250‑1 (MQ=255)                                                                                                                             gggAACTGGAAAACCGACATGTGGATTTCCTGAAACGGTATATCATCAAAGCCAGGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTg                                                                                                                            <  2:104487/250‑1 (MQ=255)                                                                                                                                 aCTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATAGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTACGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGa                                                                                                                         >  1:100136/1‑250 (MQ=255)                                                                                                                                   gggAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGC‑GCGGGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAac                                                                                                                      <  2:446409/250‑1 (MQ=255)                                                                                                                                     gAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGC‑GCGCGCGGGTCTACTGTCGCTGTACCGTCCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAacgt                                                                                                                    <  2:325113/251‑3 (MQ=255)                                                                                                                                       aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGC‑GCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAacgttg                                                                                                                  >  1:468494/1‑247 (MQ=255)                                                                                                                                       aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGC‑GCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAacgttg                                                                                                                  >  1:515253/1‑247 (MQ=255)                                                                                                                                         aCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGCTGCTACTTACTATATCCCGAACGAATCCGCATACCAGGAAGGGCGCTGGGACAAACTGcc                                                                                                                >  2:176269/1‑251 (MQ=255)                                                                                                                                         aCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGcc                                                                                                                <  1:625852/251‑1 (MQ=255)                                                                                                                                            gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGC‑GCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGGCGATAATTACTATCTCCCCAAAGAGGCGGCCTACCAGGAAGGGCGCTGGGAAacgttgccgtt                                                                                                             >  1:172398/1‑242 (MQ=255)                                                                                                                                                 gTTGCTTTCATGAAACGGGATATCATCAAAGCCATGCACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTAAGAGAATATAACGCACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CAGCGCCGGGCATCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGCTGACTGGGCGGATGCTAATTACTTACTCCCGTAAGACTCCGCATACCAGATAGGGCGCTAGGAAACGATGCCCtttacgg                                                                                                        >  2:443043/1‑247 (MQ=21)                                                                                                                                                         cctcaagaaggacatcattAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGC‑GCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCaacgcgcgat                                                                                                 <  1:279730/237‑10 (MQ=255)                                                                                                                                                               aCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGCCCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGACACACTGCCCTTTCAGCGGGCCATCATgcat                                                                                          >  1:370818/1‑248 (MQ=255)                                                                                                                                                                   gacatcattAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGC‑GCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAAAGTTGCCGTTTCAACGCGCGATCATGAATGCg                                                                                       >  1:474142/4‑250 (MQ=255)                                                                                                                                                                      aacatcaAAGCCAGGAACAAAGCAGCCGCGCTGGATGAACTGATCCCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATg                                                                                    <  1:43581/248‑1 (MQ=255)                                                                                                                                                                        cattaaaGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGC‑GCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCCGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCCGGAAGGGCGCTTGGAAACGGTGCCGTTTCAACGCGCGATCATGAATGCGATgggg                                                                                 >  1:799133/5‑250 (MQ=255)                                                                                                                                                                         atcaAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATACCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTGCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCa                                                                                <  1:536469/251‑1 (MQ=255)                                                                                                                                                                            aaaGCCATGAAAAAAGCAGCCGCGCTGGATGAACTGATACAGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGAGCCGGGCTTCGCTCAATGTTAAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGATGGGCAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCg                                                                             >  2:506494/1‑251 (MQ=255)                                                                                                                                                                              aGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGATGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGAc                                                                           >  2:709944/1‑251 (MQ=255)                                                                                                                                                                                  aTGAACAAACCTGCCGCGCTGGATGAAAGGATACGGGGGTTGCGGAGGGAATATAGCGAACAGTCAGGTTAACAGGCGGCGGCATTTTGT‑CCGCTCCGGGCTTCGCCCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTAACTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACa                                                                       <  2:59708/251‑1 (MQ=255)                                                                                                                                                                                    gAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGC‑GCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACAta                                                                     <  1:117763/251‑2 (MQ=255)                                                                                                                                                                                        aaaGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTg                                                                 >  1:400811/1‑251 (MQ=255)                                                                                                                                                                                          agcagcCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAGTCCGCATACCAgg                                                                                                                                      <  1:56956/180‑1 (MQ=255)                                                                                                                                                                                          agcagcCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAGTCCGCATACCAgg                                                                                                                                      >  2:56956/1‑180 (MQ=255)                                                                                                                                                                                           gcagcCGCGCTGGAGGAACTGATCCCGGGGTTGCGGAGGGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGC‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAgg                                                              <  2:105585/251‑1 (MQ=255)                                                                                                                                                                                             agcCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGAGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTg                                                            >  1:87592/1‑251 (MQ=255)                                                                                                                                                                                                cgcgCTGGATGAACTGATACCGGGGTTGCTGCGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAAt                                                         >  2:98048/1‑251 (MQ=255)                                                                                                                                                                                                    cTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAGTCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATgtgg                                                     >  1:35564/1‑251 (MQ=255)                                                                                                                                                                                                      ggATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGC‑GCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCCAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGGCGTTTCAACGCGCGATCATGAATGc                                                                                        >  1:585852/1‑214 (MQ=255)                                                                                                                                                                                                      ggATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGC‑GCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGc                                                                                        <  2:585852/214‑1 (MQ=255)                                                                                                                                                                                                         gtaatggttACCGGGTTTGTTGAGGGAATATATCAAACAGTCAGTTTAACGGGCTGGGGCTTTTGGT‑CCGCGTTGGGCTTCGCTCACTGTTCGGGCCGGAGCGACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCGTTCAGCGGGCCATCAGGAATGCTATGGGCAGCGACTACAGCCGTGAGGTGAAGGTGGTGaat                                                <  2:535976/243‑2 (MQ=255)                                                                                                                                                                                                           accTGATACCGGGGTTGCTGAGTGAATTTTTCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTc                                              <  1:556103/249‑1 (MQ=255)                                                                                                                                                                                                           aaCTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTGAACAGGCTGCGGCATTTTGC‑CCGCGCCGGGCTTCGCGCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTc                                              <  1:518844/251‑1 (MQ=255)                                                                                                                                                                                                            aataaTACCGGGGGTGCTGAGTGAATATATCGAACAGTCAGGTGAAGGGGCTGCAGCACTCTGC‑GCGCGCGGGTCTACTTTCGCTGTACCTACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTTCCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCt                                             <  2:223586/247‑1 (MQ=255)                                                                                                                                                                                                               gaacccGGGGTTGCTGAGTGAATATATCGAACAGGCAGGTTAACAGGCGGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTTTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGcc                                          <  2:669894/247‑1 (MQ=255)                                                                                                                                                                                                               aaTACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGC‑GCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCCGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGcc                                          >  1:141313/2‑251 (MQ=255)                                                                                                                                                                                                                  aCCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACCGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAACCACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGGGGTGAAGTCTGCCCgg                                       >  2:386884/1‑250 (MQ=255)                                                                                                                                                                                                                    cGGGGTTGCTGAGTGCATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGAGAATGTGGTGAAGTCTGCCCgtgt                                     >  1:182196/1‑251 (MQ=255)                                                                                                                                                                                                                    cGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTTCAGCCCTCTGC‑GCGCTCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGTTTGAATGGTCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgt                                     <  1:513104/251‑1 (MQ=255)                                                                                                                                                                                                                      gggtttctgaggtaagatatGGAACAGTCAGGTGAACAGGTCGGGGCATTTTGC‑CCGCGCCGGGCTTCGCTCACTGGTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACGATCGCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTTGTGAAGTCTGCTCGTGTCg                                   <  2:750131/235‑1 (MQ=255)                                                                                                                                                                                                                         gtGCTGAGGGAATAAATCGAACAGTCAGGTTAACAGGCTGCGGCAGTTGGG‑CCGCGCAGGGCTTCGCTCACTGTTCAGGCCGGAGCCAAAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGGCGGtt                                <  2:674709/250‑1 (MQ=255)                                                                                                                                                                                                                          tGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGtt                                >  1:488455/1‑250 (MQ=255)                                                                                                                                                                                                                           gCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGC‑GCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgttggctat                              <  1:732517/251‑8 (MQ=255)                                                                                                                                                                                                                           gCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTa                               <  1:622678/250‑1 (MQ=255)                                                                                                                                                                                                                            cTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTAtt                             >  1:246263/1‑251 (MQ=255)                                                                                                                                                                                                                              gAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTcc                           <  1:386884/251‑1 (MQ=255)                                                                                                                                                                                                                               ggTGAATATATCAAAACGTCAGGTTAACAGGCGGCGGCATTTGGT‑CCGCGCCGGGCTTCGCTCACTGTGCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGTCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCa                          <  2:714531/250‑1 (MQ=255)                                                                                                                                                                                                                                   aaTATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGACTGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGGGAATGGGGTGAAGTCTGCCCGTGTCGGTTTTTCCAAAAt                      >  2:684439/1‑251 (MQ=255)                                                                                                                                                                                                                                   aaTATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACGGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAt                      <  2:693317/251‑1 (MQ=255)                                                                                                                                                                                                                                     tatatCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTATGCACGTGTCGGTTATTCCACAAtgg                    >  2:361933/1‑250 (MQ=255)                                                                                                                                                                                                                                       tatCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgctg                  >  1:32065/1‑251 (MQ=255)                                                                                                                                                                                                                                       tagcgaccAGTCAGGTTAACGGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgct                   <  2:605803/243‑1 (MQ=255)                                                                                                                                                                                                                                        accGAACAGTCAGGTTAGCAGGCTGCGGCCTTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgctg                  <  2:14988/248‑1 (MQ=255)                                                                                                                                                                                                                                          cGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGctgctg               <  1:88818/251‑1 (MQ=255)                                                                                                                                                                                                                                              cAGTCAGGTGAACAGGCTGTGGCATTTGGT‑CCGCGCGGGGCTCCGTTCACTGTGCAGGCAGGAGCCTCAGACCGCCGTTGAATGGGCGGATGCGAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCAGCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGgtg           <  2:206697/251‑1 (MQ=255)                                                                                                                                                                                                                                               aGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGgtgt          >  2:410009/1‑251 (MQ=255)                                                                                                                                                                                                                                               aGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAAGGTGGTGAAGTCTGCCCGTGTCGGTTATTCCCAAATGATGCCgggggt          >  2:545887/1‑248 (MQ=255)                                                                                                                                                                                                                                                gTCAGGTTAACAGGCTGCGGCATTTTGG‑CCGCGCCGGGCTTCGCTCACTGTCCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGtt         <  1:339263/251‑1 (MQ=255)                                                                                                                                                                                                                                                       tAACCGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGAACGCCGTTGACTGGGCGGATGCTAATTACTATCTCCCGACAGAATCCGCATAACAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGAAGCAACTACATCCGTGAGGTGAATGGTGTGAAGTCTGCCCGTGTAGGTTTTTCAACAAtgcggtggggtgttgttgcct  >  2:365739/1‑233 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         AAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGT‑CCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCT  >  NC_000913/581131‑581627

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