Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,368,956	C→T	13.1%	L98L (CTA→TTA)	fxsA →	suppressor of F exclusion of phage T7

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,368,956	0	C	T	13.1%	181.7 / 2.1	81	L98L (CTA→TTA)	fxsA	suppressor of F exclusion of phage T7
Reads supporting (aligned to +/- strand):  ref base C (23/46);  new base T (5/7);  total (28/53)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGG‑TTATGCCATTTTCCGGCGTTTTTCGTTTT  >  NC_000913/4368730‑4369193                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   ggTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAAACCAGCGGCGGAGAAGGTTTAAAAG‑GTTTTGCTGGATGTTTCTTGTTTTTT‑GCTTTTTTTTCCGGGGG‑TTTTTTCCC‑GCCTCTCCCGGTTTTCtttttttttt                                                                                                                                                                                                                                          >  2:356976/1‑238 (MQ=255)    tGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTgccgc                                                                                                                                                                                                                                     <  1:712138/239‑1 (MQ=255)     gCCCATTGATTGGGGGTTTTGCTGACCCTCGGGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTgccgc                                                                                                                                                                                                                                     <  1:584382/238‑1 (MQ=255)           gTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAAACCAACGGCGGAGATGGTTAAAAGG‑GGTTCGCTGGTCATTGCTGGTTTGCT‑GGTTTTATTACC‑GGGC‑TTTTTTACC‑GAATTCCTCGGGCTTTTAATTTTATTTCCGCCCGTGCa                                                                                                                                                                                                                              >  1:950581/1‑239 (MQ=255)                      aTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTGAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGGTTAAAAGA‑GGTTCGCTGATCATTGCTGGTTTGCT‑GCCTTTTTTACC‑GGGC‑TTTTTTTCC‑CACTTACTCCGGCTTTCACTTTTTTTGCCCCCGGGGGaaaaaaaaacaa                                                                                                                                                                                                                   >  2:415928/1‑231 (MQ=255)                      aTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGGTTAAAAAT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTTCCGCCGGTGCAAAAACa                                                                                                                                                                                                                        >  1:857576/1‑234 (MQ=255)                      aTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGGAAACCCAGCGGCGGGGGTGATTTAAAAG‑GTTTCGCTGGTCTATGCTTGGTTGCT‑GCCTTTTTTACC‑CGGG‑GTTTTTTTC‑GCCTTTCTCGGGTTTTTTCTTTTTTTTTCCCCCGGGCaaaaaaaactta                                                                                                                                                                                                                   >  2:374491/1‑231 (MQ=255)                      aTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAACCCAACGGCGGAGATGATTAAAAATGGTTTCCCT‑ATTACTGCCGGTTTGTT‑GCTTTTTTTACC‑GGGCGTTTTTTCCC‑GCTTTCCCCGTTTTTCT‑CTTTTTTTTCCCTCCGTGGAAAAACATCTgc                                                                                                                                                                                                                   >  2:585509/1‑238 (MQ=255)                       tttcTGTCCCTCGTGCTGGTTTTATTCACGTGAGTTTCCGGTATGTCACTGGTACGGAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑Ac                                                                                                                                                                                                                 <  2:225421/236‑1 (MQ=255)                        tgtttcccctgtggCTGTTTTTATTCCCGTTAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACa                                                                                                                                                                                                                <  2:573748/226‑1 (MQ=255)                        tGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAAACCAACGGGGGGAGTGATTAAAAAG‑GGTTTGCTGGTCATTTGTGGTTTTCTCGTTTTTTTTCCC‑GGGT‑TTTTTTACC‑GCCTCCCCCGG‑CTTTTTCTTTTTTTGCCCCCCGGGAaaaaaaacccgcaa                                                                                                                                                                                                                 >  2:876656/1‑229 (MQ=255)                               ccTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGtcttc                                                                                                                                                                                                                                                    <  1:571302/197‑1 (MQ=255)                                         gTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAAT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACCTCCTCGGGCTTCTACTTTTATTTCCGCCGGTGCAAAAACATCTG‑GCAGTG‑GAGGTGATTCCg                                                                                                                                                                                                >  1:89972/1‑237 (MQ=255)                                             ttttcccggttgtttttggTTTGTCACTGGTACTTAACCAGGGGTTTAAGAATTTCGGGCTGATGCAGTAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTg                                                                                                                                                                                          <  2:496011/222‑1 (MQ=255)                                               ttCCCGTCATTTATCGGTATGTCACTGGTACGGAACCTGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAAGGGCGGCGGGTGAAAACCCAGGGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTgcg                                                                                                                                                                                        <  2:1005252/239‑1 (MQ=255)                                                gacaggcAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATGTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTgcg                                                                                                                                                                                        <  1:243733/232‑1 (MQ=255)                                                 cACTTCAGTTATCGGTATGTCACTGGTACGTAACCGGGGCTTTAAGAATTTCGTGCTGATGCAGCCAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTGCC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTgcg                                                                                                                                                                                        <  1:158367/237‑1 (MQ=255)                                                   cGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGGTTAAAAGG‑GTTTCGCTGATCATTTCTGGGTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTTTACTTTTTTTTCCGCCGGTGCAAAAGACTCTT‑ACCGTG‑AAGttttgtcccccttttggcttt                                                                                                                                                                                    >  1:115720/1‑220 (MQ=255)                                                   cGTCAGTTAGCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCttt                                                                                                                                                                                    <  1:374491/239‑1 (MQ=255)                                                    gTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCCCATTTTCGCtttt                                                                                                                                                                                   >  1:397494/1‑239 (MQ=255)                                                    gTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAACCCAGAGGCGGAGATGATTAAAAAG‑GTTTGGCCGGTCATTGCTGGTTTGCT‑GCCTTTATTACC‑CGGC‑TTTTTTACCTGCCTCCCCCGGT‑TTCTACATTTTTTTACGCCCGTGCAAAAACCTCTGTACTG‑G‑AAGTTGATGCCCGATTTTCGCtttt                                                                                                                                                                                   >  2:124030/1‑239 (MQ=255)                                                       aGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACCGTG‑AAGTTGATGCCGCCTTTTCGCtttttt                                                                                                                                                                                 >  1:360739/1‑236 (MQ=255)                                                                     cACTGTTACGTAACCGGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTggcgg                                                                                                                                                                   <  1:653421/238‑1 (MQ=255)                                                                     cACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTggcggc                                                                                                                                                                  <  1:867587/239‑1 (MQ=255)                                                                      aCTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCCTTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑CACTTCCTCCGGCTTCTACTTTTTTTTCCGCCGGGGCAAAAACATCTG‑ACCGTGCAAGTT‑ATCCCCCATTTTTGCCTTTCTCTCCAGCCTGggggc‑g                                                                                                                                                                >  2:867587/1‑239 (MQ=255)                                                                          gTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTTCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTTATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGttt                                                                                                                                                            >  1:875630/1‑239 (MQ=255)                                                                                  ccAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑TCAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCCTGCCTGGCGGC‑GGGTTTAGCGCCg                                                                                                                                                    >  1:460613/1‑239 (MQ=255)                                                                                    aGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTTCTTTTATTTCCGCCGGTGCAAAAACATCTG‑ACAGGG‑AAGTTTATTCCCCATTTGCGCTTTTTTTCCCTGCCTGGCCGC‑GGTTTTTGTGCCggg                                                                                                                                                  >  1:986935/1‑239 (MQ=255)                                                                                      ggCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGGTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGGGCCAAAGCCTCTG‑ACAGTG‑AAGTTGATGCCGCATTTTCGCTTTTCTCCCATGCCCGGCGGC‑GGTTTTAGAGCCcggg                                                                                                                                                 >  1:214990/1‑237 (MQ=255)                                                                                          ttAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGGCGCCGGTGCCAAAGCATCTG‑ACAGGG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGGTTTTAGGCCCgggccggg                                                                                                                                             >  1:961343/1‑234 (MQ=255)                                                                                          ttAAGAATTTCGGTCTTATGCAGCAAAAAATGGCGGCGGGTTAAACCCCCGCGGGGGAGATGATTAAAAGT‑TTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGt                                                                                                                                            <  2:857576/239‑1 (MQ=255)                                                                                          ttAAAAATTTCGGGGTGGTGGAGAAAAAAATGGGGGGGGGTGAAAACCCAGGGGCGGAGATGGTTAAAAGT‑GTTTGGCTGATCATTGCTGGTTTGCT‑TCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGt                                                                                                                                            <  2:697582/239‑1 (MQ=255)                                                                                            aaGAATTTCGTGCTGATGCAACAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTgg                                                                                                                                          <  1:884267/239‑1 (MQ=255)                                                                                                 tttCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGG‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTTCC‑GACTTCCCCGGGCTTCTACTTTTATTGCCGCCGGGGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTTCGCTTTTCTCGCATGCCTGGCGGC‑GGGTTTATGGCCGGGGCCggggggggga                                                                                                                                     >  1:12919/1‑237 (MQ=255)                                                                                                     gTGCTGATGCAACAAAAAATGGCGGCGGGTGAAAACCCAGGGGGGGGGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATAc                                                                                                                                 <  2:1038391/239‑1 (MQ=255)                                                                                                          tgtggagcaAAAAATGGGGGGGGGTGAAAACCCAGGGGCGGAAATGATTAAAAGT‑TTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑g                                                                                                                           <  2:950581/234‑1 (MQ=255)                                                                                                          ggggcagcaAAAAGTGGCGGGGGTTGAAAACCCGGCGGCGGAGATGATTAAAAGT‑GTTTCGTTGTTATTTGCTGGTTTGCT‑TCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑g                                                                                                                           <  2:946790/236‑1 (MQ=255)                                                                                                           gaggagaaaaaaaTGGGGGGGGGTGAAAACCCGGCGGCGGAGATTATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGTTGGCGGTAATACTTTT‑g                                                                                                                           <  2:175878/231‑1 (MQ=255)                                                                                                           aTGCCGGAAAAAATGGCGGCGGGTGAAATCCCGGCGGGGGGGATGATTAAAAGT‑GTTTCGCTTATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTCCTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑Ga                                                                                                                          <  2:177712/239‑1 (MQ=255)                                                                                                            agcagcaAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTCCC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑Ga                                                                                                                          <  1:585509/237‑1 (MQ=255)                                                                                                                  aaaaaaTGGGGGGGGGTGAAAACCCAGCGGGGGAGATGATTAAAAGT‑GTTTCGCTTATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAg                                                                                                                   <  2:89972/239‑1 (MQ=255)                                                                                                                  aaaaaaTGGCGGCGGGTGAAAACCCAGTGGCGGATATGATTAAAAGT‑TTTTCGCTGATCATTTCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGa                                                                                                                    <  1:358057/238‑1 (MQ=255)                                                                                                                     aaaTGGCGGCGGGTGAAAACCCAGCGGGGGATATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTAc                                                                                                                <  1:474854/239‑1 (MQ=255)                                                                                                                          aggaagggTGAAAACCCAGCGGCGGAAATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCg                                                                                                           <  1:58712/234‑1 (MQ=255)                                                                                                                           gggcggGTGAAAACCCAGGGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCGCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCg                                                                                                           <  1:570542/237‑1 (MQ=255)                                                                                                                           cggcggGTGAAACCCCAGCGGGGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GTTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCg                                                                                                           <  1:481358/238‑1 (MQ=255)                                                                                                                           cggcggGTGAAAACCCGGCGGGGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTCCC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCg                                                                                                           <  2:139897/238‑1 (MQ=255)                                                                                                                           cggcggGTGAAAACCCAGGGGCGGTGATGATTAAAAGT‑TTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCg                                                                                                           <  1:356976/238‑1 (MQ=255)                                                                                                                             gcggGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGaa                                                                                                         <  1:55601/238‑1 (MQ=255)                                                                                                                                            tgggtggtgttttttAAAAGT‑TTTTCGCGTGTTATTGCTGTTTTGCT‑GGTTTTATTCCC‑GGTT‑TTTTTTACC‑GACTTCCATGGTTGTTTCCTTTTATTGCCGCCGGTGCAAAATCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGa                                                                                         <  2:121873/226‑1 (MQ=255)                                                                                                                                            ggtgcggAAATGATTAAAAGT‑GTTTCGCTGATCATTTCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GATTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGa                                                                                         <  1:736605/236‑1 (MQ=255)                                                                                                                                            gcggggGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGa                                                                                         <  1:631343/239‑1 (MQ=255)                                                                                                                                                  gagaTGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGGGGGCGGAAAACTTTT‑TATGGTGAGTACCAGCGAAAGGATGGTGAggcgggccgcct                                                                                   >  1:647410/1‑228 (MQ=255)                                                                                                                                                   gtatgatTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATGCTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCt                                                                                   <  1:187267/236‑1 (MQ=255)                                                                                                                                                    gttgtttAAAATT‑TTTTCGCTTCTCTTTGCTGTTTTGTT‑GTTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTTTTCTACTTTTATGGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCtt                                                                                  <  2:360739/233‑1 (MQ=255)                                                                                                                                                     atgatTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTATTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTg                                                                                 <  1:602176/238‑1 (MQ=255)                                                                                                                                                        atTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTggcggc                                                                                                                                                                  <  1:595198/156‑1 (MQ=255)                                                                                                                                                        atTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTggcggc                                                                                                                                                                  >  2:595198/1‑156 (MQ=255)                                                                                                                                                            aaaatgttttCCCTGATCAGTTCTGGTTTTCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCGCGGTCTTTTACTTTGATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATaa                                                                          <  1:359225/232‑1 (MQ=255)                                                                                                                                                               gt‑gtTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTATCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGGGAAAGGATGATGAGGGCGGCCGCCTTGATCTTAAAGAc                                                                      >  1:658889/1‑239 (MQ=255)                                                                                                                                                                       tctttttttttgtgttttGTT‑GTTTTTTTTCCC‑GGGT‑TTTTTTCCC‑TATTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGcc                                                               <  2:12919/224‑1 (MQ=255)                                                                                                                                                                       tctgctcttgtttgttttGCT‑GTTTTTATTCCC‑GGGT‑TTTTTTCCC‑GCCTTCCTCGGTTTTCTACTTTTATTGCCGCGGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGTTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGcc                                                               <  2:986935/224‑1 (MQ=255)                                                                                                                                                                       gcttgtCATTGTTGGTTGGCT‑GCTTTTATTACC‑GGGT‑TTTTTTACC‑GATTTCCTCGGTTTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGcc                                                               <  1:1057855/234‑1 (MQ=255)                                                                                                                                                                                tGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑Ga                                                                                                                          <  2:913792/171‑1 (MQ=255)                                                                                                                                                                                tGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑Ga                                                                                                                          >  1:913792/1‑171 (MQ=255)                                                                                                                                                                                     tttttgttttttttttttCC‑GGGT‑TTTTTTACC‑GGTTTTCCCGTTTTTTTTCTTTTATTGCCGCCGGTTCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATCCTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCg                                                  <  2:397494/229‑1 (MQ=255)                                                                                                                                                                                     gtttttttttttttATTCCC‑GGGT‑TTTTTTCCC‑GACTCTCTGGTTCTTCTCCTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGGGCTTTTCTCGCATGCCGGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCg                                                  <  2:460613/229‑1 (MQ=255)                                                                                                                                                                                     gTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTTATCCTAAAAACGATCGCCAGGGTTAAttgcga                                                 >  1:935287/1‑234 (MQ=255)                                                                                                                                                                                     gTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGa                                                 >  1:370048/1‑239 (MQ=255)                                                                                                                                                                                     gTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGGCATACTTTTGGTTGGTG‑GGACCAGCGAAAGGATGATGAGGGCGACCGCCTTGATCATAAAAAAGATCGGCAGGGGTAATGTGGa                                                 >  2:497307/1‑239 (MQ=255)                                                                                                                                                                                      tttGCT‑GTTTTTTTTACC‑GGGT‑TTTTTTACC‑TACTTCCTGGGGTTTTTCTTTTTATTGCCGCCGGTGCAAAAGCTTGTG‑ACAGGG‑AAGTTGATGCCGCATTGGCGCTTTTCGCGCATGCCTGGCGGC‑GGTTGTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGaa                                                <  2:1083201/239‑1 (MQ=255)                                                                                                                                                                                        tgct‑gcTTTTATTCCC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAGGCATCTG‑GCAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCTTGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAAc                                              <  1:124030/239‑1 (MQ=255)                                                                                                                                                                                         ttttttttttttTCCC‑GGGT‑TTTTTTACC‑GGTTTCCTGGGTTTTTTTTTTTTTTTGCCGCTGGTGCAAAAGCATCTG‑ACAGGG‑AAGTTGATGCCGGATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAAc                                              <  2:214990/233‑1 (MQ=255)                                                                                                                                                                                          tgtttttttttTCCG‑GGGT‑TTTTTTACC‑GTTTTCCCGGGTCTTTTGTTTTTATGGCCGCCGTTTCAAAAGCATCTG‑ACAGTG‑AAGTTGAGGCCGCATTTGCGCTTTGCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATTCTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACg                                             <  2:608633/234‑1 (MQ=255)                                                                                                                                                                                          gttttttttATTACC‑GGGC‑TTTTTTAAC‑GAGTTTCGCGGTCTTTTATTTTTATTGCCGCCGGTGGAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACg                                             <  1:415928/234‑1 (MQ=255)                                                                                                                                                                                          ct‑gctTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCTGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑TCTGGTGAGAACCGGCTCAAGGAGGCGGAGCGCGACCGCCTTGAGTCTAAAAGTGACCGCCAGGATTAAGGtttgaacc                                             >  2:1058586/1‑231 (MQ=255)                                                                                                                                                                                             gtttttATTACC‑GGTC‑TTTTTTACC‑AACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGGGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCg                                          <  1:836124/237‑1 (MQ=255)                                                                                                                                                                                               ttttttTCCC‑GGTT‑TTTTTTCCC‑GTTTTCCTCGGTTTTCTACTTTTATTGCCGCCGGTGCAAAGGCATTGG‑ACGGGG‑AAGTTGAGGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGa                                        <  2:1049951/239‑1 (MQ=255)                                                                                                                                                                                                     taCC‑GGGC‑TTTTTTACG‑GCATTCCTCGGTTCTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGTAACGCCGGATTAtgt                                  <  1:30769/239‑1 (MQ=255)                                                                                                                                                                                                     taCC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATGGCCAGGGTTAATGTCGAAAAGCCGGAtttttg                                  >  1:1049951/1‑235 (MQ=255)                                                                                                                                                                                                               ttttttACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGG‑TTATg                          <  1:1047114/238‑1 (MQ=255)                                                                                                                                                                                                                               cTCGGTCTTCTTTTTTTGTTGCCGCCGGTGCAAAAGCATCGG‑ACAGTG‑AAGTTGAGGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGG‑TTATGCCATTTTCCGGCGttt          <  1:497307/239‑1 (MQ=255)                                                                                                                                                                                                                                        ctACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGGGAGTACCAGCGCAAGGAGGATGAGCGGGACCCCCTTGTCCATAAAACACATGCCCAAGGTTAATGTGACAACCCCGGCTTATGGGGTTTGTCCCTTTTCCCGGGGTtttttttttt  >  2:667696/1‑239 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   GGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑GTTTCGCTGATCATTGCTGGTTTGCT‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTT‑GATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGG‑TTATGCCATTTTCCGGCGTTTTTCGTTTT  >  NC_000913/4368730‑4369193

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