BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,760,383	A→C	22.5%	V289G (GTG→GGG)	selA ←	selenocysteine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,760,383	0	A	C	22.5%	185.9 / 18.3	109	V289G (GTG→GGG)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand): ref base A (52/31); new base C (2/24); total (54/55)
Fisher's exact test for biased strand distribution p-value = 6.51e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

GCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGACG  >  NC_000913/3760146‑3760615

gCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGGTCATCTCTTTTTTTCCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACCACTTGTCGCCGGAGAAAATCa                                                                                                                                                                                                                                           >  1:409892/1‑238 (MQ=255)

gCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGGTCACCTCTTTTTTACCACCAATAAT‑TCCTGCCCGCGGCCCGCC‑CAACAACCTTTCGCCCGAGAAACTCa                                                                                                                                                                                                                                           >  2:1185022/1‑238 (MQ=255)

gCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCa                                                                                                                                                                                                                                           <  2:1042729/238‑1 (MQ=255)

gCCTGGATTTGAATGACCTCAGCGCTGCGGGTAATCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGCGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTTCCTACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGGAAATCa                                                                                                                                                                                                                                           >  1:245143/1‑238 (MQ=255)

           aaTGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCCCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTTTCGCCGGAGAAACTAACCCGACTGAcg                                                                                                                                                                                                                                >  1:255387/1‑238 (MQ=255)

                ccTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGGTCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCCGCCCGCC‑TAACAACTTGTCCCCGGAGAAACTCCCCAGGCTGCcgccccc                                                                                                                                                                                                                           >  1:618015/1‑236 (MQ=255)

                       gcggcgGGTAGGCAGGCGCAGGGTCGGTATTTTTTCACTCAGAGCTTAAGGGTTAAAATAAA‑GACG‑AAACGTGGCTTCCGGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑CCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAACCTCCCCAGACTGACGCCCGCCGCAATc                                                                                                                                                                                                                    <  2:1242238/235‑1 (MQ=255)

                           cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACt                                                                                                                                                                                                                >  1:1044864/1‑238 (MQ=255)

                              gTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTTTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCt                                                                                                                                                                                                             >  1:81781/1‑238 (MQ=255)

                                aaGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGCTCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAAATCACCAGACTGACGCCCGCCGCACTCCACCCCTGCg                                                                                                                                                                                                          >  1:822147/1‑239 (MQ=255)

                                      gcgcAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGGCCGCC‑TAACAACTTTTCGCCCGAGAAAATCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCGttt                                                                                                                                                                                                    >  2:843438/1‑238 (MQ=255)

                                       cgcAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTTGGCGGAGAAACTCCCCAGCACGACGCCCCCCGCAATCAACTCCTTCCGCactt                                                                                                                                                                                                    >  2:1224506/1‑235 (MQ=255)

                                              tCGGTATTTTTTCGCACAGCGTTTCAGGGTTTAAAAAAA‑CACG‑CACCGTGGCTTCCAGCGCCGGGGGGGTCTTTTTATCCGCGCGTAATGCACGTTTCAGCGGGTGGCTTG‑GCGGGCGGGCGATCATCTCTTTTTTCCCAAAAATAAT‑TCCTGCCGGCGGCCCGCC‑AAACAACTTGTCGCGGGAGAACCCCACCAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCttt                                                                                                                                                                                            <  2:1038648/238‑1 (MQ=255)

                                                ggTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGGTCATTTCTTTTTTACCAACAATAAT‑TTCTGGCTGCGGCCCGCC‑TAACCACTTGTCGCCGGAGAAACTCACCAAACTGACCGCCGCCGCCATCAACTCCTGCGGCATTGGC‑TCCTTcc                                                                                                                                                                                          >  1:309969/1‑237 (MQ=255)

                                                ggTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTTCCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACCTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCCCCATCCACTCCTCTGGCATTGGC‑TCTTTCg                                                                                                                                                                                          >  2:411757/1‑238 (MQ=255)

                                                  tatttttcCCCTCAGAGCTTCAGGTTTTAAAAAAA‑GACG‑CAACGTGGCTTCCGGCGCGGCGAGGGTCTTTTTACCCGCGGGAACTGCCCGCTCCGGCGGGTGGCTTG‑GCGGGGGGGGGTTCATCTCTTTTTTACAAAAATTAAT‑TCCTCCCTGGGGCCCCCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGc                                                                                                                                                                                        <  2:285138/235‑1 (MQ=255)

                                                        ttCACTCAGAGCTTCAGGGTGTAAATAAA‑GACG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCCGCCTTGGC‑CCTTTCGGCAAACCg                                                                                                                                                                                  >  2:68396/1‑238 (MQ=255)

                                                         tccctaagagCTTCAGGGTTTAAATAAA‑GACG‑CACCGGGGTTTCCAGCGCGGCGAGGGTCTTTTTATCCGCGCGTAATGCCCCCTTCGGGGGTGGGCTTT‑GCGGGGGGGCGATCATCTCTTTTTTCCCAAAAATAAT‑CCCTCCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAAAACCTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTa                                                                                                                                                                                <  2:589023/233‑1 (MQ=255)

                                                         tcacccagggctgaggggTGAAAAAAAA‑GCCG‑CAACGGGGCTCCCAGCGCCGGGGGGGTCATTTTCCCCGCGCGTAATCCACGCTGCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑CCCTCCCGGCGCCCCCCC‑TAACAACTTGTCGCCGGAGAACCTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTa                                                                                                                                                                                <  2:1044864/224‑1 (MQ=255)

                                                         tcaccaagagCTTCAGGGTTTAAAAAAA‑GACG‑GAACGGGGCTTCCAGCGCGGGGAGGGTCGTTTTATCCGCGCGTAAGGCGCGCTTCAGCGGGTGCCGTG‑GCAGGCGGGCGCTCATCGCTTTTTGCCCAACAATATT‑TCCGGCCTGCGCCCCGCC‑AAACAACTTGCCGCCGGAGAAACTCCCCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTa                                                                                                                                                                                <  2:306931/233‑1 (MQ=255)

                                                         tCACTGAGAGCTTCAGGGTATAAATAAT‑GCCG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCCTTTTACCCGCGCGTAATGCCCGCTGCAGCGGGTGGCTTT‑GGAGGGGGGCGATCTTCTCTTTTTGCCCAACAATAAT‑TCCTCCCGGCGGCCCGCC‑TAACACCGTGTCGCCGGAGAAACTCACCCGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCTTa                                                                                                                                                                                <  2:1124615/239‑1 (MQ=255)

                                                         tCACTGAGAGCTAGAGGGTGTAAAAAAG‑GCCG‑AAACGTGGCTTCCAGCGCCGCGGGGGTCTTTTTACCCGCGCGTAATGCACGTTTAGGGGGGTGGTTTG‑GCGGGCGGCCAATCATCTCTTTTTTACAAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAAAAACTTGCCGCCGGAAAACCTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTa                                                                                                                                                                                <  2:913432/239‑1 (MQ=255)

                                                         tCACGCAGAGCTTCAGGTTGTAAATAAA‑GACG‑CACCGTGGCTTCCAGCGCCGGGGGGGTCATTTTATCCGCGCGTAATGCCCGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTCCCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTa                                                                                                                                                                                <  1:10633/239‑1 (MQ=255)

                                                         caccgcaggtgttgaGGGTAAAAAAAAA‑GGCG‑CAACTTGGCCCACGGGGCCGCGAGGGTCATTTTACCCAGGGGTACTGCACCCTCCGGGGGGTGGCTTC‑GCGGGCGGGCACTCTTCTCTTGTTAACCAACAATAAT‑CCCTGCCTGCGGCCCCAC‑AAACAACTTGTCGCCGGAAAATCCCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTGCGGCAAACCGTa                                                                                                                                                                                <  2:1072101/225‑1 (MQ=255)

                                                          ccacccagggtttaGGGGGTAAAAAAACGCCG‑CAGCGGGGCTCCCGGCCGCGCGGGGGTCATTTTACCCGCGCGTAATGCACGCTTAAGGGGGGGGCTTT‑GGGGGGGGGCGATCATCTCTTTTTTACAAACAAAAAT‑CCCTGCCGGCGGCCCGCC‑TAACACCTTGTCGCCGGAGAAACGCACCAGACGGACGCCCGCCGCAGTCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTa                                                                                                                                                                                <  2:945676/226‑1 (MQ=255)

                                                                gagCTTCAGGGTGTAAATAAA‑GACG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCACGCCGCAATCAACTCCTGCGGCCTAGGC‑TCTTTCGGGCAACCGTACTGGCTa                                                                                                                                                                         >  2:859050/1‑239 (MQ=255)

                                                                         gtttaaaataaa‑GGCT‑CAACGTGCCTTCCCCCGCCGAGTGGGTTCTTTTATCCGCCCGAAAGTCACGCTTCCGGGGGTGGCTTG‑GCGGGGGGGCCTCCATCGCTTTTTTACCAACAAAATT‑TCTCGCCGGGGCCCCGCC‑AAACACCTGCTCGCGGAAAATACCCCCAAGATTGCCCCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGAc                                                                                                                                                                 <  2:220617/233‑1 (MQ=255)

                                                                          gtgtAAAAAAA‑GCCG‑CAACGTGGTTTCCAGCGCCGCGAGGGTCATTTTACCCGCGCGTAATGCACGCTTCAGGGGGGGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGAc                                                                                                                                                                 <  1:1019048/237‑1 (MQ=255)

                                                                             taaataaa‑GACG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACCATAAT‑TCCTGCCTGCGGCCCGCC‑TAACCACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTCCTGGCTAAGATCGACCACCg                                                                                                                                                            >  2:597785/1‑239 (MQ=255)

                                                                             taaataaa‑GACG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCg                                                                                                                                                            >  1:926334/1‑239 (MQ=255)

                                                                                   aa‑GACGCCAACGTGGCTCCAGGCGCGCGGGGGGTCTTTTTACCCGGGCTAAAGGCCCCCTTCAGGGGGGGGCTTG‑GCGGGGGGGCGCTCTTCTCTTTTTTACCAACAATAAT‑CCCTGCCCGCGGCCCGCC‑TAACAACTTGTCGCCGGAAAACCTCCCCAGCCGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                                       <  2:738894/239‑1 (MQ=255)

                                                                                    aagcacacaccggcgCTTCCAGCCCCGCGGGTGTATTTTCATCCCTTTGTAATGCCCGCTGCCGCGGGCGCCTTG‑GCGGGGGGGCAATCATCTCTTTTTCACCAAAAATATT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCCCCAGACTCACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                                       <  2:688884/225‑1 (MQ=255)

                                                                                    aaataggcaaggtgtcgcccAGCGCGGGGAGGGGCTTTTTACCCGAGCGAAAGGACCGCTTCAGGGGGTGGCTTG‑GCGGGCGGGCGATCATCTCTTTTTAACAAACAATAAT‑CCCGGCCGGCGCCCCGCA‑AAACACCTTGGCGCCGGAGAAACCCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                                       <  2:625703/221‑1 (MQ=255)

                                                                                    aaaccccaaggttgctttCCGGCGCTGGGAGGTTTTTTTTTCCCGCGCGTTTTGCCCGCTACGGGGGTTGGCTTT‑GCCGGGGGGGGCTCATATCTTTTTT‑CCAAAAAAAATGTCCTGCCGGCGCCCGGCC‑TAACATCTTGTCGCCGAACAAACGCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                                       <  2:757172/223‑1 (MQ=255)

                                                                                      gACG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTTTCGCCGGAGAAACTCACCAGACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑CCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCCCt                                                                                                                                                     >  1:560232/1‑238 (MQ=255)

                                                                                        cG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGAcgcccgc                                                                                                                                                                                                                           <  1:469467/167‑1 (MQ=255)

                                                                                        cG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGAcgcccgc                                                                                                                                                                                                                           >  2:469467/1‑167 (MQ=255)

                                                                                           cAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGGTCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCCGACTGACGCCCGCCCCCATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCCTACTCGCCAACATCGCCCAGCGGGCCACTGCCtc                                                                                                                                                >  2:150251/1‑238 (MQ=255)

                                                                                                       acgccgcggggggttCATTTTTCCCGCGCGAAACGCAC‑CTTCCGGGGGTGGCTTTGGGGGGCGGGCGTCCATCTTTTTTTTCCCAAAAAAATT‑CCCTCCCGGGCCCCCCCC‑TAACAACTTGTCGCCGGAGAAACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTa                                                                                                                                    <  2:1213906/225‑1 (MQ=255)

                                                                                                               cgAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTTCCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACCACTTTTCGCCGGAGCAACCCACCAGACTCACCCCCGCCCCAATCAACTCCCGCGGCATTTGC‑TCTTTCGGCAAAACGTACTGGCCAAGATCGACCAGCGAGCCACTGTCTAAATCAGACACAACGGGAACa                                                                                                                            >  2:1014068/1‑239 (MQ=255)

                                                                                                                              ccGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAAACGTACTGGCTAAAATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACCGGGAAATCCCGCTCCTTGCCg                                                                                                             >  2:1067658/1‑239 (MQ=255)

                                                                                                                                   cgTAATGCACGCTTCCGCGGGTTGCTTT‑GCAGGCGGGCTGCCATCCCTTTTTTTCCTCCCATAAT‑TCCTGCCTGCGGCCCGAC‑TTACAACCTGTTGGCTGAGAAACTCACCAGACTGACGCCCGCCGCAATCCACCCCTCCGGCATTGGCTTCTTCCGCCAACCCTTCTTGCCCAAGTTC‑ACCGGCGGACCCCCGCCCCACACAGTCCCTACTGGAACCTCCAGCGCTTTGCCGGgcg                                                                                                         >  2:76529/1‑238 (MQ=255)

                                                                                                                                          caccctacgtcGGGTGGGTTG‑CCGGGCGGGCTCACATCTTTTTTTTACAAAAAAAAAT‑CCCTGCCCGGGGCGCGCC‑AAACAGCTTCGCGCCGGAAAACCTCCCCAAACAGTCGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGt                                                                                                 <  2:952816/229‑1 (MQ=255)

                                                                                                                                              cTTCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCTAAATCAGTCACTACGGGAACATCCCGCCCTTTGCCGCGCCCCCCCCGTTcc                                                                                              >  2:1199718/1‑238 (MQ=255)

                                                                                                                                               ttCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGGACTGGCTAAGATCGACCAGCGAGCCACCGCCTCAATCAGTCACTACGGGGACATCCAGCCCTTTGCCGCGCGCCACCAGGTCCg                                                                                             >  2:320062/1‑238 (MQ=255)

                                                                                                                                                tcggggggTGGCTGG‑GCGGGCGGGCCTACATCTTTTTTTTACCAAAAATATC‑CCCTCCCTGGGCCCCCCC‑AAAAACCTTGGCGCGGGAAAACCTCCCCAGACTGCCCCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCt                                                                                           <  2:741795/234‑1 (MQ=255)

                                                                                                                                                  aGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGGACATCCAGCTCTTTGCCGCGCGCCACCCGGTCCCCTTc                                                                                         >  1:833371/1‑239 (MQ=255)

                                                                                                                                                  aGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTc                                                                                         >  1:9713/1‑239 (MQ=255)

                                                                                                                                                       gTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTa                                                                                    >  1:208963/1‑239 (MQ=255)

                                                                                                                                                         tgttttgccgggcggcctctactctcTTTTTCCCAAAAAATACT‑CCCTCCCTGCGCCCCCCC‑TAAAATCTGCCCGCGGGAAAACCTCCCCAGACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTAt                                                                                   <  2:236188/217‑1 (MQ=255)

                                                                                                                                                         ggCTTT‑GCAGGCGGGCGGTCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TCACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCAATGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTAAGGGAAACTCCGGCTCTTTGCCGAGCGCCACCCGTTCCGCTTCacatat                                                                                   >  2:565887/1‑233 (MQ=255)

                                                                                                                                                         ggCTTT‑GCAGGCGGGCGGTCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCCGACTGACGCCCGCCGCCATCAACTCCTGCGGCATTTGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCCAAATCAGTCACCACGGGAACATCCAGCCCTTTGCCGAGCGCCACCCCTTCCGCTTCctcttt                                                                                   >  2:697029/1‑236 (MQ=255)

                                                                                                                                                         ggCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATc                                                                                  >  1:1157481/1‑239 (MQ=255)

                                                                                                                                                         ggCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCCCTTCATCTATc                                                                                  >  1:676820/1‑239 (MQ=255)

                                                                                                                                                         gccTTG‑GCGGGCGGGCGATCATCTCTTTTTACCAAAAAATAAT‑CCCTCCCGGCGCCCCGCA‑TAACAACTTGGCGCCGGAAAACCCCACCAGCCGGACCCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTAt                                                                                   <  2:113393/236‑1 (MQ=255)

                                                                                                                                                              ttggaggggggcatcctccttttttttCCCAAAAATAAT‑CCCTGCCGGGCGCCCGCC‑TAAAAACTTTCCGCGGGAAAAACCCCCCAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCtt                                                                              <  2:1135923/218‑1 (MQ=255)

                                                                                                                                                                  aGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGt                                                                          >  1:22426/1‑239 (MQ=255)

                                                                                                                                                                    gcgggcgTTCATCTTTTTTTTCCAAACAATAAT‑TCCGCCCTGCGCCCCGCC‑AAACACCTTGCCGCCGGAGAAACCCACCAGACGGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                                         <  1:68396/238‑1 (MQ=255)

                                                                                                                                                                    gcgggcgATCATCTCTTTTTACCCAACACTAAT‑TCCGGCCTGCGGCCCGCC‑TAACAACTTGTCGCGGGGGAAACTCCCCAGCCGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                                         <  1:1111379/238‑1 (MQ=255)

                                                                                                                                                                      gggcgATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTTGTGaa                                                                       >  1:1218735/1‑238 (MQ=255)

                                                                                                                                                                        tgactattcttctttttACAAAAAAAAAT‑TCCGCCCGGCGCCCGGTT‑AAACACGCTGCCGCGGGAGAATCACCCCAGCATGCCCCCCGCCGCAATCACCTCCTGTGGCATTGGC‑TCTTTCGGCAACCCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGaa                                                                       <  2:890190/224‑1 (MQ=255)

                                                                                                                                                                          gATCCTTTCTTTTTTCCCAACAATAAT‑CCCTCCCTGCGGCCCGCC‑TAACACCTGGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCt                                                                  <  1:249538/239‑1 (MQ=255)

                                                                                                                                                                              ctttctTTTTTCCCACCAATAAT‑CCCGCCCCGCGGCCCGCC‑TAACACTTTGCCGCCGGAAAACCCCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAAt                                                              <  1:1197770/236‑1 (MQ=255)

                                                                                                                                                                               tctcTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCCTACTGGCTAAGATCGACCAGCCAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGCGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAAt                                                              >  1:779908/1‑238 (MQ=255)

                                                                                                                                                                                   ttttaacccAAAAATAAT‑CCCGGCCGCCGCCCCGCC‑TAAAAACTTGCCGCGGGAGAAACCCCCCAGCCGGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGt                                                         <  2:28490/232‑1 (MQ=255)

                                                                                                                                                                                        tcccAACAATATT‑ACCGGCCGGCGGCCCGCC‑TAACAACTTGCCGCGCGAAAACCTCACAAGATTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                                    <  1:703037/237‑1 (MQ=255)

                                                                                                                                                                                         accaaaaaagactccccgcccgcGGCCCGCC‑AACCACGTTGTCGCGGAAAAACCTCACCAGACGCACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                                    <  2:245143/222‑1 (MQ=255)

                                                                                                                                                                                         aCCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGtt                                                     >  1:397285/1‑237 (MQ=255)

                                                                                                                                                                                           caacaaGAAT‑CCCTGCCGGCGGCCCCCC‑GAAAAACTTGCCGCCGGAGAAGCTAACCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGAACGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTg                                                 <  1:1014068/239‑1 (MQ=255)

                                                                                                                                                                                             acaaTAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                               >  2:163346/1‑239 (MQ=255)

                                                                                                                                                                                             acaaTAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                               >  2:1069841/1‑239 (MQ=255)

                                                                                                                                                                                             acaaTAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                               >  1:1164261/1‑239 (MQ=255)

                                                                                                                                                                                               aaaaat‑TCCGCCCGGCGGCCCGCC‑TAACACCTTGCCGCCGGAGAAACTCCCCAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTAt                                             <  1:311862/236‑1 (MQ=255)

                                                                                                                                                                                                  aat‑CCCTGCCTGCGGCCCGCT‑TAACAACTTGCCGCCGGAGAAACCCACCAGACTCCCGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                          <  1:822208/239‑1 (MQ=255)

                                                                                                                                                                                                  aat‑CCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAACCTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                          <  1:1232117/239‑1 (MQ=255)

                                                                                                                                                                                                  aat‑CCCTGCCGGCGGCCCGCC‑TACCAACTTGCCGCCGGAAAACCCCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                          <  1:995083/239‑1 (MQ=255)

                                                                                                                                                                                                   actccccgccggcGCCCCCCC‑TACAAACTTGCCGCCGGAAAAACCCACCAAACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                          <  1:667128/232‑1 (MQ=255)

                                                                                                                                                                                                   aagccccgcctgcGGCCCGCC‑TAACAAATGGCGGCAGGAAAAACCCACCAGACGGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                          <  1:116406/232‑1 (MQ=255)

                                                                                                                                                                                                     tccctgcctgCGGCCCGCT‑AAACAACTTGCCGCCGGGAAAATCCCCCAGACTCACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACt                                        <  1:215929/237‑1 (MQ=255)

                                                                                                                                                                                                     tccctccctgCGGCCCGCC‑TAACAACTTGCCGCGGGAGAACCTCACCAGCCTGACCCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACt                                        <  1:1199718/237‑1 (MQ=255)

                                                                                                                                                                                                          gcctgcGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTCCTGGCATGTACTTtca                                    >  2:84349/1‑238 (MQ=255)

                                                                                                                                                                                                          gcctgcGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                                    >  1:706882/1‑238 (MQ=255)

                                                                                                                                                                                                          gcctgcGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                                    >  1:516792/1‑238 (MQ=255)

                                                                                                                                                                                                             tgcGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                                         >  2:416969/1‑198 (MQ=255)

                                                                                                                                                                                                             tgcGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                                         <  1:416969/198‑1 (MQ=255)

                                                                                                                                                                                                             tgcGGCCCGAC‑TAACAACTTGTCGCCGGAAAAACTCACCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCaa                                <  1:843438/239‑1 (MQ=255)

                                                                                                                                                                                                              ggcgCCCGCCTTAACACCTTGACGCCGGAAAAACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCaa                                <  1:697029/236‑1 (MQ=255)

                                                                                                                                                                                                               cGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCt                                                                                           >  1:1124155/1‑178 (MQ=255)

                                                                                                                                                                                                               cGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACa                              >  1:1111683/1‑239 (MQ=255)

                                                                                                                                                                                                               cGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACa                              >  1:142229/1‑239 (MQ=255)

                                                                                                                                                                                                                  accccaaaaacaacTTGCCGCCGAAAAAACCCCCCAGACACACCCCCGCCGCAATCAACTCCTGCGGCATTGTC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGt                            <  2:397285/231‑1 (MQ=255)

                                                                                                                                                                                                                   ccGCC‑TAACAACTGGACGCCGGAGAACCCCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGc                          <  1:444833/239‑1 (MQ=255)

                                                                                                                                                                                                                          aacaactgtcggcaggAAAAACCCACCACACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTAtt                    <  1:996777/225‑1 (MQ=255)

                                                                                                                                                                                                                                ttGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTttcat               >  2:1110306/1‑238 (MQ=255)

                                                                                                                                                                                                                                  gTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTttcattc             >  2:251641/1‑238 (MQ=255)

                                                                                                                                                                                                                                  gTCGCCGGAGAAACTCACCAGACTGACACCCGCCGCAATCAACTCCTGCAGCATTGGC‑TCTTTCGACAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGACACCAGTTCCGCTTCATCTATCGCTTTGGTAAACCCCTGCATGCTGTAGTTACTGGTATGTACTTTCATCAACAGCGCGGCAGTttgattc             >  1:1036853/1‑238 (MQ=255)

                                                                                                                                                                                                                                       cGGAAAAACTCACCAGACTAGCGCCAGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCAACGcc       <  1:199031/239‑1 (MQ=255)

                                                                                                                                                                                                                                         cagaCACTCACCAGCCTGCCCCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCt      <  2:676820/237‑1 (MQ=255)

                                                                                                                                                                                                                                          agaAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTg     >  2:1155830/1‑238 (MQ=255)

                                                                                                                                                                                                                                             aaCTCACCAGACTAACCCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGACg  <  1:163842/238‑1 (MQ=255)

GCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACG‑CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTT‑GCAGGCGGGCGATCATCTCTTTTTTACCAACAATAAT‑TCCTGCCTGCGGCCCGCC‑TAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGC‑TCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGACG  >  NC_000913/3760146‑3760615

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