Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	786,855	A→C	16.3%	G247G (GGT→GGG)	gpmA ←	phosphoglyceromutase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	786,855	0	A	C	16.3%	197.7 / 2.1	83	G247G (GGT→GGG)	gpmA	phosphoglyceromutase 1
Reads supporting (aligned to +/- strand):  ref base A (55/7);  new base C (2/17);  total (57/26)
Fisher's exact test for biased strand distribution p-value = 3.70e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCA‑CGTAAAGAGTTAC‑CGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGA  >  NC_000913/786626‑787089                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                gAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGGCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGtt                                                                                                                                                                                                                                         >  1:425277/1‑238 (MQ=255)          gATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑ccgc                                                                                                                                                                                                                             >  1:327193/1‑238 (MQ=255)           aTGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAACCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATAACGTTTACTTCGCTTTACCCTGGTTT‑TC‑AA‑ccccc                                                                                                                                                                                                                            >  2:701005/1‑238 (MQ=255)           aTGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGAAACAAAAAAGCCGACCCACTTCCAGTCGGCTTTCTCATTTTAAACG‑AAAGACGTTTACTTTCTTTTAACCCGGTTT‑GC‑CA‑ccccc                                                                                                                                                                                                                            >  2:787160/1‑238 (MQ=255)           aTGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCCCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCCTTTTAAACG‑AATGACGGTACTTTCCCTTTACCCTGGTTT‑GC‑AA‑ccgcc                                                                                                                                                                                                                            >  2:424756/1‑238 (MQ=255)            tGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGCCTCACTTGCAGCCGGCTTTCTCATTTTAAACC‑AATAACGTTTACTTCCCTTTACCCTGGTTT‑GC‑aaacccccc                                                                                                                                                                                                                           >  2:635937/1‑232 (MQ=255)               tctGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGACCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGACGGCTTTCTCATTTTAGACG‑AATGCCGTTTACTTCGCTTTACCCTGGTTT‑GA‑AA‑CCGCCcct                                                                                                                                                                                                                         >  1:493022/1‑237 (MQ=255)                      cGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTc                                                                                                                                                                                                                   >  1:114096/1‑236 (MQ=255)                      cGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTc                                                                                                                                                                                                                   >  1:851639/1‑236 (MQ=255)                           aCTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGAAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCCGGTTT‑GC‑AA‑CCGCCCCGGATTTCGCTCCgg                                                                                                                                                                                                            >  1:19972/1‑237 (MQ=255)                           aCTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTCCGGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTCCCCCGGTTT‑GA‑AC‑CCGCCGCGCCTTTCCCtcccat                                                                                                                                                                                                           >  1:812618/1‑234 (MQ=255)                             tGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGAtc                                                                                                                                                                                                          >  1:733696/1‑238 (MQ=255)                                 aaaTGGATTAAAGGTGCTTCGGGGGTAAGTTTTATGTTCCGGAGGGGCCGGCCAAGTGCGTTAGCCGAGGAACCCCCACATCCTTAATAGGGCGCGCGCCGTCTCATCAGTCATTTTGCTCCTTTTTTAGCAAAAAAAAAGCCGCCATTCTTTCAGTCGCTTTTTCCATTTTAACCG‑AAGGGCGTTTCTTTCTTTTTCCCTGGGTTC‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGt                                                                                                                                                                                                      <  2:27201/237‑1 (MQ=255)                                                 cgTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCTCCAGCATTACCCAGataa                                                                                                                                                                                     >  1:1004113/1‑239 (MQ=255)                                                   tCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCCGCATTACCCAGataata                                                                                                                                                                                   >  1:93904/1‑239 (MQ=255)                                                   tCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGA‑AAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCCTTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGAAAAAAg                                                                                                                                                                                  >  2:633741/1‑239 (MQ=255)                                                                   aTTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGat                                                                                                                                                                                       >  1:624894/1‑219 (MQ=255)                                                                       tCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATTTCGTCAGCATTACCCAGAAAATAGCGTTTCAGGGGTTTGAAAt                                                                                                                                                               >  2:926427/1‑239 (MQ=255)                                                                        cGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTCCTTCGCTTTCCCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCCCTGCGAGCTCGTCAGCATTACCCCGCTAATAGCGTTTCACCGGTTTGAAAtt                                                                                                                                                              >  1:248872/1‑239 (MQ=255)                                                                        cGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCCGATAAAACCGTTTCAGCGGTTTGAAAt                                                                                                                                                               >  2:216101/1‑238 (MQ=255)                                                                         ggATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTTGCTGCGATCTCGTAAGCCTTACCCAGATAAAAGCGGTTCGGCGGTTTGAAAtt                                                                                                                                                              >  2:998275/1‑238 (MQ=255)                                                                         ggATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGACAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTACACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCTCCAGCATTACACAGAAAATAACGTTTTAGCGGGTTGAAAAtc                                                                                                                                                             >  2:427178/1‑239 (MQ=255)                                                                                         tGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATac                                                                                                                                              >  1:773213/1‑238 (MQ=255)                                                                                                tCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGACGCTGCTTTCGCTGCGATTTCGTCAGCATTACCCAGATAATAGCGTTTCACCGGTTTGAACTTATCGTCAAAATCaacaacaacagg                                                                                                                                      >  2:858018/1‑228 (MQ=255)                                                                                                           cccacgggcaataTAAGCCGCGCAGCGTCGCATAGGGAAATGTCCTCAATTGTAAGCAAAAAAAAAGCCGACTCCCTGGCAGTCGGTTTTCTCATTTAAAACG‑AATGACTTTTCCTTCCTTTTACCCTGGTTT‑GA‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCg                                                                                                                           <  1:701005/227‑1 (MQ=255)                                                                                                             cacaGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCTTCGAACTCAAACACCAGCGGCACGCCAGTCggg                                                                                                                         >  2:927163/1‑239 (MQ=255)                                                                                                                  gCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCACCGGCACGCCAGTCGGGATAtt                                                                                                                    >  1:837238/1‑239 (MQ=255)                                                                                                                       aTAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGctgtcccttttacacatctccg                                                                                                                >  1:956982/1‑217 (MQ=255)                                                                                                                       aTAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATCTTAAGc                                                                                                                >  1:546541/1‑238 (MQ=255)                                                                                                                       aTAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCt                                                                                                               >  1:594681/1‑239 (MQ=255)                                                                                                                       aTAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAAGCGGGAAATTAAGCt                                                                                                               >  2:58828/1‑239 (MQ=255)                                                                                                                       aTAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGATGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATACTCGTGGAAATCATACACCAGCGGAACGCCAGGCGGGGAATTaaacc                                                                                                               >  2:462893/1‑236 (MQ=255)                                                                                                                            tCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCCCGCCAGTCGGGCTATTAAGCTCAAGa                                                                                                          >  1:539666/1‑239 (MQ=255)                                                                                                                            tCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAAATCATAAACCAACGGCACGCCAGCCGGGGAATTAAACTaaaaa                                                                                                          >  2:136194/1‑237 (MQ=255)                                                                                                                            tCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTAATTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGTTCTCGTCAGCATTACCCAGATAATCGCGTCTCAGCGGTTTGGAATTCTCGCTGAACCCCTACACCAGCGGCCCGCCAGGCGGGTAATTACGCTTAAGa                                                                                                          >  1:1022808/1‑239 (MQ=255)                                                                                                                                 agcggTCGCAGCAGGCATTGGGCTCCATTTTTAACAAAAAAAAGGCCGCCTCCTTTGCAGTCGGTTTTCTCATTTTAAACG‑AATGACTTTTACTTCGTTTTCCCCTGTTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtctc                                                                                                     <  1:952310/235‑1 (MQ=255)                                                                                                                                         caacaggaattgcgcacaatggatAGAAACAAAAACGCCGTCCCACTTGCAGTGGGCTTTCTCATTTAAAACG‑AATGACTTTTACTTCGTTTTCCCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCTGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGc                                                                                             <  1:216101/216‑1 (MQ=255)                                                                                                                                         cATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGCTCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAAACCATACCCCAACGTTCCCCCAGTCGGGAAATTAATCTCAaaaaattttttctccc                                                                                             >  2:536832/1‑227 (MQ=255)                                                                                                                                              ggCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATg                                                                                        >  1:55918/1‑239 (MQ=255)                                                                                                                                              ggCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCCCATg                                                                                        >  1:97599/1‑239 (MQ=255)                                                                                                                                              ggCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACCCCAGCCGGGAAATTAAGCCCAAAAATCTCTTCTTCGccacct                                                                                        >  2:952310/1‑234 (MQ=255)                                                                                                                                                    gtgCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCa                                                                                                                               >  1:719371/1‑194 (MQ=255)                                                                                                                                                     tttttatttGTTACAAACAAAAGACCCGCCTCTTTGGCTGCTGGTTTTCTTGTTTAAACCG‑ATTGCCTTTACCTGCGCTTTACCGTGTTTC‑GA‑AC‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCaa                                                                                 <  2:773213/232‑1 (MQ=255)                                                                                                                                                     tgCTCTTTTGTTAGCAAAAAAAAAGCCGACTCCCTTGCAGTCGGTTTTCTATTTTTAAACG‑ATTGACGTTTACTTCGCTTTCCCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCaa                                                                                 <  1:376469/239‑1 (MQ=255)                                                                                                                                                     tcccccATTTTTAAAAAAAAAAACGCCCTCTCACTTGAGGTCGGTTTTCTCATTTTAAACG‑AATGACGTTTACTTCGTTTCCCCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCaa                                                                                 <  1:635937/235‑1 (MQ=255)                                                                                                                                                     tccccaaTTGTAAGCAACAAAAAAGCCGACCCACTTGCAGTGGTCTTTCTATTTTTAACCG‑AATGTCGTTTACTTCGCTTTCCCCTGTTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCaa                                                                                 <  1:58828/233‑1 (MQ=255)                                                                                                                                                       cTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCCTGTTATCAAGa                                                                               >  1:673580/1‑239 (MQ=255)                                                                                                                                                         ccATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCGCGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGata                                                                             >  1:570906/1‑239 (MQ=255)                                                                                                                                                          taTTGTTAAAAAAAAAAAACCCGACTCCTTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTCCCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGAt                                                                                                                                                                                                           <  2:94037/111‑1 (MQ=255)                                                                                                                                                               ttAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGTAATTTAAc                                                                       >  1:388842/1‑239 (MQ=255)                                                                                                                                                               ttAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGAAATTTCAc                                                                       >  1:63009/1‑239 (MQ=255)                                                                                                                                                               ttAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGGCGGGAATTTAAGCTCCAAAATATCATATTCTCCCAAGTTATCAAAATATTTCa                                                                        >  2:381406/1‑238 (MQ=255)                                                                                                                                                                 aGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCa                                                                     >  1:693493/1‑239 (MQ=255)                                                                                                                                                                    aacaaAAAGGCCGACCCCCTGCCCGTCGTTTTCCTTATTTTAAACG‑AATGCCGTTCACTTCGTTTTACCTTGGTTT‑CC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAgcg                                                                  <  1:136194/239‑1 (MQ=255)                                                                                                                                                                    aacaaAAAGCCCGTCTCACTTGCAGTCGTCTTTCCTATTTAAAACG‑AATGCTGTTTCCTTCG‑TTTACCCCGGGTT‑GCAAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAgcg                                                                  <  1:633741/239‑1 (MQ=255)                                                                                                                                                                    aacaaAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCa                                                                                                                                                                                                 >  1:94037/1‑112 (MQ=255)                                                                                                                                                                    aacaaAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAgc                                                                   >  1:241113/1‑238 (MQ=255)                                                                                                                                                                    aaaaaaaaaGCCTCCTCGTTTGCGGTCGTTTTTCTAATTTAAACGG‑AATGCCGTTTCCTTCGCTTTCCCTTGGTTG‑CA‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAgcg                                                                  <  1:926427/236‑1 (MQ=255)                                                                                                                                                                    aaaaaaaaaCCCGTCTCATTTGGGGTGGTCTTTCTCTTTTAAAACA‑AGTCACGTTTCCCTCGTTTTCCCGGTGTTG‑AC‑AC‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGGTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATAGTAAGCTCAAGAATCGCTTCTTCGCTCATGTTATCAAGATATTTCACCAgcg                                                                  <  2:248872/236‑1 (MQ=255)                                                                                                                                                                      caaAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGTTAATAGCGGTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGTCCGCCAGTCGGGATATTAACCTCAAGAACCTCTTCTTCGCCCATGTTATCAAAAAATTTCACCaccaca                                                                >  2:164013/1‑234 (MQ=255)                                                                                                                                                                      caaAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGTTCAAGAATCCCTTCTTCCCCCATGTTATCAAAAAATTTCACACGCGCa                                                                >  2:473415/1‑239 (MQ=255)                                                                                                                                                                          acgcccgcctcctgttcatgctgCTTTCTCATTTAAAACG‑AACGACGTTTACTCCGTTTTCCCGTGTTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCA‑CGTa                                                           <  1:26326/217‑1 (MQ=255)                                                                                                                                                                          aaagaaacagcccTTGCAGTCGGTTTTCTAATTTAAAACG‑AATGACTTTTACTCCGCTTTCCCCGGTTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCA‑CGTa                                                           <  1:682951/227‑1 (MQ=255)                                                                                                                                                                           acgccgtcttgtttGCTGTCGGCTTTCTATTTAAAAACA‑AATGTCTTTTACGCCGCTTTCCCGTGTTGT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCA‑CGTa                                                           <  1:864679/226‑1 (MQ=255)                                                                                                                                                                                    cACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGGTATTAAGCTCAAGAATCTTTTCTCCGCTCATGTTACCAAGTTATTTAACCA‑CGCACCGTaaaaggtaacc                                                 >  2:864679/1‑231 (MQ=255)                                                                                                                                                                                       ttgtatgccttttttttttttaAAACG‑ATTGGTGTTT‑CTTTCCTTTAACCTGGTTTCGC‑CA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCA‑CGTAAAGAGTTAC‑Cgtg                                             <  2:93904/221‑1 (MQ=255)                                                                                                                                                                                         gcagtggtttttCTC‑TTTTAAACGAAATGCTTTTTCCTTCCTTTTACCTGGGTTG‑AC‑AT‑CCGCCGCTGCTTTCGGTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCA‑CGTAAAGAGTTAC‑CGtgtg                                           <  2:55918/230‑1 (MQ=255)                                                                                                                                                                                             tgggtttcctgtttttAAGCG‑AATGCTTTTTCCTTCGCTTTCCCCTGTTTT‑GA‑AG‑CCGCCGTTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCA‑CGTAAAGAGTTAC‑CGTGTGCAg                                        <  1:424756/226‑1 (MQ=255)                                                                                                                                                                                             ccGGTTTTCCCATTTAAAACG‑ATTTGCTTTTACTTCCCTTTGCCCTGGTTG‑AC‑AT‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCA‑CGTAAAGAGTTAC‑CGTGTGCAg                                        <  1:906491/237‑1 (MQ=255)                                                                                                                                                                                               ggCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCCTTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGAACGCCAGTCGGGATATTACGATCCACAATCCCTTCTTCGCCCAGTTTACAACAATTTTTCCCCAGAGCA‑CGTAAAAATTTCC‑CGTGTGATGCga                                     >  2:239844/1‑239 (MQ=255)                                                                                                                                                                                               ggCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCCCATGTTATCAAGATATTTCACCACGCCA‑GTTAAAGAGTAAC‑Cggggccgcgac                                     >  2:622322/1‑229 (MQ=255)                                                                                                                                                                                               ggCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCCCTCATGTTATCAAGAAATTTTCCCAGCGCA‑CGTAAAAAGTTcc‑ccctgtctcttt                                     >  2:682951/1‑228 (MQ=255)                                                                                                                                                                                                       cATTTTAAACG‑AATGGCGTTTTCTTCGCTTTTCCCTGGTTT‑TC‑AACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCA‑CGTAAAGAGTTAC‑CGTGTGCAGCGATGATCAc                              <  1:756897/239‑1 (MQ=255)                                                                                                                                                                                                              aaCG‑AGTGCGGTTTCCCTCTCTGTCCCCGGGTTT‑GA‑CA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCA‑CGTAAAGAGTTAC‑CGTGTGCAGCGATGATCACGCGCTCAc                      <  1:622322/239‑1 (MQ=255)                                                                                                                                                                                                               cgg‑AATGATTTTTACTCCTGTTCACCCTGTTCT‑GA‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCA‑CGTAAAGAGTTAC‑CGTGTGCAGCGATGATCACGCGCTCAcc                     <  1:1042731/237‑1 (MQ=255)                                                                                                                                                                                                                     tGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑CA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCCGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCGGCGGCACGCCAGTCGGGGTATTAAGCTCACGAATCTCTTCTTCGATCATTTTATCAAGGTATTTACCCAGCGCC‑CGTAAAGATTTACGCATGCGC‑GCGATGATCACCGCCTCACCGctc                 >  2:233391/1‑238 (MQ=255)                                                                                                                                                                                                                          ttgtcttctttttCCCCTTTTTA‑GC‑AC‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAATAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCA‑CGTAAAGAGTTAC‑CGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATa           <  2:1022808/229‑1 (MQ=255)                                                                                                                                                                                                                            ttcttgcttttCCCTTGGTTT‑CC‑AC‑CCCGCGCTGCTTTCGCTGGGAGCTCGTCAGCATTACCCAGAGAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATCCCCCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCA‑CGTAAAGAGTTAC‑CGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATAcg         <  2:820913/231‑1 (MQ=255)                                                                                                                                                                                                                            tccTTCGTTTTCCCCTGTTTT‑CA‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCA‑CGTAAAGAGTTAC‑CGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATAcg         <  1:1011857/237‑1 (MQ=255)                                                                                                                                                                                                                                ttgtttcaCCCTGGTTT‑GA‑AA‑CCGCCGCCGCTTTCGCTGCTATCTCGGCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCA‑CGTAAAGAGTTAC‑CGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGc     <  1:473415/232‑1 (MQ=255)                                                                                                                                                                                                                                 cGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCA‑CGTAAAGGGTTAC‑CGTGTGCAGCGATGATCACGCGCCAACCGCTCTTCATACCCggga    >  2:117684/1‑237 (MQ=255)                                                                                                                                                                                                                                 cGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCA‑CGTAAAGAGTTAC‑CGTGTGCCGCGGTGATCAAGCGCTCCCCGCCTTCCATAcccgccc    >  2:20389/1‑236 (MQ=255)                                                                                                                                                                                                                                 cGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCA‑CGTAAAGAGTTAC‑CGGGTGCAACGATGATCAAGCCCCCACCGCCCTTCCAACCCGGCa    >  2:412844/1‑239 (MQ=255)                                                                                                                                                                                                                                   cTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCA‑CGTAAAGAGTTAC‑CGTGTGCAGCGATGATCACCCGCCCACCGCTCTTCAAACGAGGcgaa  >  2:1011857/1‑236 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                GAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACG‑AATGACGTTTACTTCGCTTTACCCTGGTTT‑GC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCA‑CGTAAAGAGTTAC‑CGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGA  >  NC_000913/786626‑787089

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: ‑