BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	3,880,022	T→G	12.9%	V286G (GTC→GGC)	selA →	L‑seryl‑tRNA(Sec) selenium transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	3,880,022	0	T	G	12.9%	378.9 / 2.5	174	V286G (GTC→GGC)	selA	L‑seryl‑tRNA(Sec) selenium transferase
Reads supporting (aligned to +/- strand): ref base T (29/107); new base G (30/4); total (59/115)
Fisher's exact test for biased strand distribution p-value = 5.46e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 179 total aligned reads displayed.

ACGCGAATGATTATCGTCAGGCGGTGAATGAAAATACCGCACTGTTGATGAAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCCGCAGCGCAGAGGTCATTCA  >  NC_007779/3879785‑3880259

gcgcgAGTGATTATCGTCAGGCGGTGAATGAAAATACCGCGCTGTTGATGAAAGTACATACCAGTAACTACAGCATTCAGGGGTTCGCCAATGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGGCTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGCAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGt                                                                                                                                                                                                                                                     <  2:415236/237‑1 (MQ=255)

          ttATCGTCAGGCGGTGAATGAAAATACCGCACTGTTGATGAAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGGGGCCGGCT‑GGTgg                                                                                                                                                                                                                                         >  1:152092/1‑238 (MQ=255)

              cGTCAGGCGGTGAATGAAAATACCGCGCTGTTGATGAAAGGACATACCGGTAACTACAGCGTTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGttt                                                                                                                                                                                                                                     <  2:1373140/239‑1 (MQ=255)

                  aGGCGGTGAATGAAAATACCGCACTGTTGATGAAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCGGGAGTTGATTGCGGC‑GGGGGTCCGGCT‑GGTGCGTTtc‑tc                                                                                                                                                                                                                                 >  1:340750/1‑238 (MQ=255)

                         gaatgaaAATACCGCACTGTTGATGAAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGAc                                                                                                                                                                                                                         <  1:1238443/238‑1 (MQ=255)

                              aaaaTACCGCACTGTTGATGAAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGGCT‑GGTGAGGTTC‑TCCGG‑CGACACgtt‑g                                                                                                                                                                                                                  >  1:1276905/1‑239 (MQ=255)

                                    ccGCACTGTTGATGAAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCGGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGAATGCGGC‑GGGGGGTCGTCTGGGTGAGGTT‑‑ACCGG‑CGACAACGT‑GCTAggg                                                                                                                                                                                                            >  2:1232460/1‑238 (MQ=255)

                                      gCACTGTTGATGAAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGGGTCAGTCT‑GGGGAGGTTT‑TCCGG‑CGGCAAGTT‑GTTAGGC‑gg                                                                                                                                                                                                         >  1:174206/1‑239 (MQ=255)

                                                tgaAAGTACATACCAGTAACTACAGCATTCAGGTGTTCACCAAAGCATTAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCATAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGACAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTGGGC‑GGGCCgcaggc                                                                                                                                                                                                <  2:68978/238‑1 (MQ=255)

                                                tgaAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGGTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCgcaggc                                                                                                                                                                                                <  1:43478/238‑1 (MQ=255)

                                                tgaAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCgcaggc                                                                                                                                                                                                <  2:1464263/238‑1 (MQ=255)

                                                tgaAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCgcaggc                                                                                                                                                                                                <  2:1731977/238‑1 (MQ=255)

                                                tgaAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCgcaggc                                                                                                                                                                                                <  2:788735/238‑1 (MQ=255)

                                                tgaAACTACATTCCAGTAACTACGGCATTCAGGGGTTCACCAAAGCGGTAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGTCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGGTGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCgcaggc                                                                                                                                                                                                <  2:7790/238‑1 (MQ=255)

                                                    aGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGGCAGTCT‑GGTGAGGTTC‑TCCCG‑CGGCAAGAT‑GTTTTGG‑GGGGCGcaggcaggga                                                                                                                                                                                           >  1:547978/1‑237 (MQ=255)

                                                      tACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGGGGCCGGCT‑GGGGGGTTTC‑TCCCG‑GGCCAAGTT‑GGTAGGG‑GGGCCGCCGGCAGGGatt                                                                                                                                                                                         >  1:1461198/1‑239 (MQ=255)

                                                      tACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGGGGCCAGCT‑GGGGGGGTTT‑TCCGGCCGACAAGTT‑GTAAGAG‑Gggcggcggccggagtt                                                                                                                                                                                          >  1:1367377/1‑229 (MQ=255)

                                                      tACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCCGGCT‑GGTGAGTTTT‑TCCCG‑GGCCAAGTT‑GTTAGGG‑GGGCCGccggcggggaat                                                                                                                                                                                         >  1:99032/1‑235 (MQ=255)

                                                      tACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGGGGGTTTT‑TCCCG‑CGGCCAAGT‑GTTTTGC‑GGGGCGCCGGCAGGGatt                                                                                                                                                                                         >  1:1687505/1‑239 (MQ=255)

                                                      tACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGGCCGTCC‑TGTGGGGTTT‑TTCCG‑CGGCAAGGT‑GTTTGGC‑CGGCCGCAGGCAGGaaat                                                                                                                                                                                         >  1:1626405/1‑237 (MQ=255)

                                                      tACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGGCAGTCT‑GGGGGGGTTT‑TCCGG‑CGGCAAGTT‑GTTAGGC‑GGGGCGCCGGGAGGGatt                                                                                                                                                                                         >  1:1156450/1‑239 (MQ=255)

                                                      tACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCCGCAGTCT‑GGGGGGGTTT‑TCCCG‑CGGCCAGTT‑GTTTGGC‑GGGGCGccggcggggagt                                                                                                                                                                                         >  1:868525/1‑235 (MQ=255)

                                                           gccAGTAACTACAGCATTCAGGGGTTCACCAAGGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTAttgt                                                                                                                                                                                    <  2:14316/238‑1 (MQ=255)

                                                           aCCAGTAACTACAGCGTTCAGGGGTTCACCAAAGCGAGAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACGGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTAttgt                                                                                                                                                                                    <  1:862253/239‑1 (MQ=255)

                                                           aCCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTAttgt                                                                                                                                                                                    <  1:962263/239‑1 (MQ=255)

                                                           aCCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTAttgt                                                                                                                                                                                    <  2:1385868/239‑1 (MQ=255)

                                                           aCCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTAttgt                                                                                                                                                                                    <  1:937002/239‑1 (MQ=255)

                                                             cAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGGGTGGGGCT‑GGTGGGTGTCTTCCGG‑CG‑CCAGGT‑GTTTGGG‑GGGGCGCCGGCGGGGATTAttgttt                                                                                                                                                                                  >  2:1133613/1‑238 (MQ=255)

                                                             cAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCGGGCAGGAATTAttgt                                                                                                                                                                                    <  1:139824/237‑1 (MQ=255)

                                                                tAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGGGGTCGGTTG‑GGGGGGTTTT‑CCCGG‑GCACAGATT‑GTGTGGG‑GGGCGGCGGGCGGGAGTTATTGTGGGt                                                                                                                                                                                >  2:937002/1‑239 (MQ=255)

                                                                tAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGGCGGTCT‑GGGGGGTTTT‑TCCGG‑CGGCCAGTT‑GTTAGGC‑GGGGCGCAGGGAGGAATTATTGTTGGtg                                                                                                                                                                               >  1:788164/1‑238 (MQ=255)

                                                                 aaCTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGAc                                                                                                                                                                                                                         <  1:1628822/198‑1 (MQ=255)

                                                                 aaCTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTa                                                                                                                                                                               <  1:422158/238‑1 (MQ=255)

                                                                      cAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGGCCGTCT‑GGGGAGGTTT‑TCCCG‑CGGCAAGTT‑GTTTGGC‑GGGGCGCCGGGAGGGATTATTTTTTGTaaaaaaa                                                                                                                                                                         >  1:1660197/1‑238 (MQ=255)

                                                                       aGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGGCAAGCT‑GGGGGGTTTC‑CCTGG‑CCACCAGGT‑TTTTTGC‑GGGGCGCCGGGAGGGATTATTTTTTGTGAAAAaga                                                                                                                                                                        >  1:1114599/1‑239 (MQ=255)

                                                                         cATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGGGGCAGTCT‑GGGGGGTTTT‑TCCCG‑CGGCAAGTT‑GTTTGGG‑GGGGCGCCGGCAGGGATTATTTTTTGTGAAAAAaaga                                                                                                                                                                      >  1:1125411/1‑239 (MQ=255)

                                                                          aTTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGGGGGCGCGTGGGGGGGTTTG‑TCCGC‑GGAC‑AGTT‑GTTAGGC‑GGGGCGCGGGGCGGGATATTTGTTGGTGAAAAAaagaa                                                                                                                                                                     >  2:948118/1‑238 (MQ=255)

                                                                          aTTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAag                                                                                                                                                                         <  1:639844/235‑1 (MQ=255)

                                                                          aTTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGGGCGTCT‑GGTGTGGTTC‑TCCCG‑CGACAAGAT‑GTTAGGC‑GGGGCGCAGGCGGGGATTATTGGTGGGAAAAAAaagaa                                                                                                                                                                     >  1:173829/1‑238 (MQ=255)

                                                                           ttCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTAttgt                                                                                                                                                                                    <  1:1393725/223‑1 (MQ=255)

                                                                             cAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAgatgat                                                                                                                                                                  <  1:105754/239‑1 (MQ=255)

                                                                              aGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGGGGCAGGCTGGGGGGGTTTT‑TCCGC‑GGCAAAGTT‑G‑TAGGC‑GGGGCGCCGGCCGGGAATATTGTTGGTGAAAAAAAAGTGATc                                                                                                                                                                 >  1:1643534/1‑239 (MQ=255)

                                                                                gggTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGGGGTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGGATTATTGTTGGTAAAAAAGAGATGAGcg                                                                                                                                                                >  1:878576/1‑238 (MQ=255)

                                                                                   ttCACCAAAGCTATGGAT‑GAGGCGGAACTGGGGGCGCTCGGCAAAGAGTTGGTTGTTCCCGTGGTGACTGATTTGGGAGGTGGCTC‑GCTGGTCGTTCTTGGCCAGTCGGTTTTGCGGAAGGAGCCAATGCCGCAGGAGTTGATTGCGGG‑GGGCGTCAGTCT‑GGTGATTTTC‑TCCGG‑CGAAAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATcgccc                                                                                                                                                             <  2:984916/238‑1 (MQ=255)

                                                                                   ttCACAAAAGCG‑TTGATGGAGGCGGAACTGGTGGCGCCCGGCAAATGGTTGGATTTTCCCGTATTACCTGATTAAGGCAGT‑GCTCTGCTGTTCGATCTTGGCCAGTACGGTTTGCCGAAAGAGCCATTGCCGCAGGATTTGTTTGCGGG‑GGGCTTCAGTTT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGCC‑GGGCCCCAGGCAGGAATTATTGTTGTTAAAAAAGAGATGCTcgccc                                                                                                                                                             <  2:858438/238‑1 (MQ=255)

                                                                                     cACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGt                                                                                                                                                                                                                                                     >  1:1378999/1‑153 (MQ=255)

                                                                                     cACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGGCT‑GGGGGGGTTT‑TCCCG‑CGGCAAGTT‑GTTTGGG‑GGGGCGCCGGGAGGGATTATTTTTTGGGAAAAAAAGAAGATcgccccc                                                                                                                                                           >  1:1078078/1‑236 (MQ=255)

                                                                                     cACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTAACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGGGTCAGGCT‑GGGGAGGTTT‑TTCCG‑CGGCCAGGT‑GTTAGGC‑GGGGCGCCGGCAGGGAATTTTTTTTGTAAAAAAAAAAAGAAcgccccc                                                                                                                                                           >  1:1545019/1‑236 (MQ=255)

                                                                                          aaGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCa                                                                                                                                                                                                                                                   <  2:1441186/150‑1 (MQ=255)

                                                                                          aaGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCa                                                                                                                                                                                                                                                   >  1:1441186/1‑150 (MQ=255)

                                                                                              gATAGAT‑GAAGGGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGcc                                                                                                                                                 <  1:1359226/239‑1 (MQ=255)

                                                                                              gATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTAttgt                                                                                                                                                                                    <  1:940909/204‑1 (MQ=255)

                                                                                              gATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGcc                                                                                                                                                 <  1:982298/239‑1 (MQ=255)

                                                                                              gATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTCCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGcc                                                                                                                                                 <  1:191619/239‑1 (MQ=255)

                                                                                               aTAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGcc                                                                                                                                                 <  1:834916/238‑1 (MQ=255)

                                                                                               aTAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTGGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGcc                                                                                                                                                 <  1:1253798/238‑1 (MQ=255)

                                                                                               aTAGAT‑GAAGCGGAACTGGTGGCGCTCGGAAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGcc                                                                                                                                                 <  1:710487/238‑1 (MQ=255)

                                                                                                tAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCAc                                                                                                                                               <  1:705071/239‑1 (MQ=255)

                                                                                                   aT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGAc                                                                                                                                                                                                                         <  1:120527/164‑1 (MQ=255)

                                                                                                   aT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATcgccc                                                                                                                                                             <  1:1787712/222‑1 (MQ=255)

                                                                                                       aggCGGAACGGGGGGGGCTCGGAAAAGGGCTGGTGTTTCCCGTAGGGATTGATTTGGGCAGGGGCTC‑GCGGGTGGATCTTAGCCAGTCGGGTTTGCCAAAGGAGCTAATGCCGGAGGGGTTGTTTGCGGG‑GGGCGTCAGTCT‑GGTGAGTTTC‑CCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGAAGGAATTATTGTTGGAAAAAAAGAGATGACCGCCCGCCTGCAAAGCCACCCGCt                                                                                                                                          <  2:1085120/236‑1 (MQ=255)

                                                                                                           ggAGCTGGGGGCGCTGGGCAAGGTGCGGGTTGTTCCCGTGGTGACTGTTGTAGGCGGTGGCTC‑GCGGGTCGTTTTTAGCCAGTACGGTTTGCAGAAGGAGCCAGTGCCGCGGGATTTGATGGGGGC‑GGGCGTCAGTCG‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTGGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAg                                                                                                                                      <  2:1249168/238‑1 (MQ=255)

                                                                                                             aaCTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCg                                                                                                                                    <  1:19535/238‑1 (MQ=255)

                                                                                                                       gcTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGGCCGTCT‑GGGGGGTTTC‑TCCGG‑CGGCAAGTT‑GTTAGGC‑GGGGCGCCGGCAGGGATTATTGTTTGGGAAAAAAAGATGATCGCCCGCCCGCCAAGACACCCGCCGCAGCGGGCATTAcgc                                                                                                                          >  1:576537/1‑238 (MQ=255)

                                                                                                                           ggaaaaGGGCGGGATGTTCCGGTGGTAGCTGATTTGGGAGGTGGCGG‑GCGGGTGGATTTTGGCAAGTGCGGTTTGCCGAAAGAGCCAATGCCGCAGGGTTTGATGGGGGG‑GGGCGTCAGTCT‑GGTGATTTTC‑CCGGG‑CGAAAAGTT‑GTTAGGG‑GGGCCGCAGGAAGGAATTATTGTTGGAAAAAAAGAGATGACCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑gg                                                                                                                     <  2:1687505/235‑1 (MQ=255)

                                                                                                                                         gtttCCGGTGGGGACTTGTTTGGGCAGGGGCTG‑GCGGGCGGTTCTGAGCCAGTACGGTTTGCGAAAAGAGCCAATGCCGCGGGAGTTGATTGCGGG‑GGGCGTCGGTCT‑GGTGGTTTTC‑TCCGG‑CGAAAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTTTTGTTGGTAAAAAAGAGATGACCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTcg                                                                                                      <  2:1660197/237‑1 (MQ=255)

                                                                                                                                            tCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGaa                                                                                                           >  1:1041206/1‑230 (MQ=255)

                                                                                                                                              ccGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCgcggcg                                                                                                 <  1:399595/239‑1 (MQ=255)

                                                                                                                                                gTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCt                                                                                               <  1:978725/239‑1 (MQ=255)

                                                                                                                                                  aGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTAGGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCTgg                                                                                             <  1:1446469/239‑1 (MQ=255)

                                                                                                                                                          aTTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTAttgt                                                                                                                                                                                    <  1:352018/145‑1 (MQ=255)

                                                                                                                                                           tttAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCg                                                                                                                                                                                                     <  1:356557/127‑1 (MQ=255)

                                                                                                                                                                  aggggCTC‑GCGGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTc                                                                              <  1:596392/235‑1 (MQ=255)

                                                                                                                                                                    gggCTC‑GCTGGTCGATCTTGGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑gg                                                                                                                     <  1:1153362/196‑1 (MQ=255)

                                                                                                                                                                         gctgtGGTCATTCTTAGCGGGTACGGTTGGCGGAAAGAGCCAATGCCGCGGGGGTTGTTGGCGGG‑GGGCGTCAGTCT‑GGTGAGTTTC‑CCCGG‑CGAAAGGTT‑GTTAGGG‑GGCCCGCAGGAAGGAATTTTTGTTGGTAAAAAAGAGATGACCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCtttattt                                                                       <  2:1807705/235‑1 (MQ=255)

                                                                                                                                                                         c‑GCTGGTCGATCTTAGCCAGTACGGTTGGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCtttattta                                                                      <  1:1018136/239‑1 (MQ=255)

                                                                                                                                                                                    tCTTAGCCGGTACGGTTTGCCGAAAGCGCCAATGCCGCAGGAGTGGTTTGCGGG‑GGGCGTCAGTCT‑GGTGGGTTTC‑TCCGG‑CGACAGGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTGGGTAAAAAAGAGATGACCGCCCGCCt                                                                                                                                                         <  1:1818996/147‑1 (MQ=255)

                                                                                                                                                                                      ttAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGct                                                          <  1:1251244/239‑1 (MQ=255)

                                                                                                                                                                                              tACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGGGGCACAAT‑GGGGGGTTTT‑TCCGG‑CGCCAAGAT‑GTTAGGCGGGGCCCCAGG‑AGGAAGTTTTTTTGGGTAAAAAAAGAAGGACGCCCGCCTGCACAGCCCACCGCCGACACGTGCATTTACCGCGGGACAAAAGAAACCCCCCGG‑GCCGGAAGCCACGTTGCTGCTTTCTTTTCACCCCTAACCTCTgcaggt                                                  >  2:731769/1‑234 (MQ=255)

                                                                                                                                                                                              tACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGGGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCCGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGAACCCCGCGGCGCTGGGAGCCACGTTGCGTCttt                                                                           >  1:225211/1‑214 (MQ=255)

                                                                                                                                                                                              tACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGGGGTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGGATTATTTTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCCCCGAAGCGTGCCTTACGCGC‑GGGTAAAAAAACCCCCGCGGCGCTGGAAGCCACGTTGCGTCTTTTTTTaca                                                                    >  1:1206618/1‑221 (MQ=255)

                                                                                                                                                                                              tACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGGATTATTGTTGGTAAAAAAAAGATGATCGCCCGCCTGCCAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCCGGGAGCCCCGTTGCGTCTTTTTTTACACCCTGAAGCTCTGAGTGa                                                  >  1:775805/1‑239 (MQ=255)

                                                                                                                                                                                              tACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGGTAAAATGACCCTCGCGGCGCTGGGAGCCACGGTGCGGCTTTTTTTACACCCTGAAGCTCTGAGGGa                                                  >  1:770544/1‑239 (MQ=255)

                                                                                                                                                                                              tACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGGGGGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCCGGCAGGGATTATTGTTGGTAAAAAAGAGATGAACGCCCGCCTGCAAAACCACCCGCTGAAGCGTGCATTAc                                                                                                                            >  1:1104664/1‑166 (MQ=255)

                                                                                                                                                                                                      gCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGGGCGGCT‑GGTGGGGTTT‑TCCCG‑CGACCAGTT‑GTTAGGC‑GGGGCGCGGGCCGGGCGTATTTTTGGTGAAAAAAAGGTGATCGGCCGCCCGCCAAGGCACCCCCCGAAGGGGGGATTACGCGC‑GGGGAAAATAAACCCCGCGGCGCCGGGAGCCCCGGTGCGTCTTTATTTATACCCCGAAGCTCCGACTGAAAAATTaa                                          >  2:362473/1‑238 (MQ=255)

                                                                                                                                                                                                          aaaGGGCAAGTGCCGAGGGAGTGTGTTGCGGC‑GGGGGTAGTTTT‑GGTAGTTTTC‑TCCGG‑CGA‑AAGTTGGTTGGGG‑GGGCCGCAGGCAGGATTTTTTTTTGGAAAAAAAGAGAGGACCCCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGa                                       <  2:269501/238‑1 (MQ=255)

                                                                                                                                                                                                              agCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTc                                                                              <  1:233899/195‑1 (MQ=255)

                                                                                                                                                                                                                 cAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTgcg                                <  1:257273/238‑1 (MQ=255)

                                                                                                                                                                                                                       cgcgGGAGTGGATGGGGGG‑GGGCGTCAGTCT‑GGTGTTTTTC‑TCCGG‑CGAAAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTTTTGGAAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTgc                                 <  1:188085/231‑1 (MQ=255)

                                                                                                                                                                                                                             aGTTGATTGCGGC‑GGGCGTCGGTCG‑GGTGAGTTTC‑GCCGG‑CGAAAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTGGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCcg                    <  1:1283142/238‑1 (MQ=255)

                                                                                                                                                                                                                             aGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGAAAAAAAGAGATGACCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCcg                    <  1:1374420/238‑1 (MQ=255)

                                                                                                                                                                                                                             aGTTGATGGCGGC‑GGGCGTCAGTCG‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCcg                    <  1:246323/238‑1 (MQ=255)

                                                                                                                                                                                                                              gTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCcgc                   <  1:594627/238‑1 (MQ=255)

                                                                                                                                                                                                                                   ttGGGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTTTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCcgcagcgca             <  1:362473/239‑1 (MQ=255)

                                                                                                                                                                                                                                   ttGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCcgcagcgca             <  1:262204/239‑1 (MQ=255)

                                                                                                                                                                                                                                   ttGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGAAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCcgcagcgc              <  1:130877/238‑1 (MQ=255)

                                                                                                                                                                                                                                   ttGCGGC‑GGGCGCCAGATT‑GGTGAGTTTC‑TCCGG‑AGACAAGTT‑GTTAGGC‑GCGCCGGAGGCAGAATTTATGGTTGAAAAAAAAGAGATGCCCACCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCcgcagcgca             <  1:448743/239‑1 (MQ=255)

                                                                                                                                                                                                                                   gtGCGGC‑GGGCGTCAGTCT‑GGGTGTTTTC‑TCGGG‑AGACAAGTT‑GTTGGGC‑GGGCCGCAGGCAGGATTTATTGTTGGTAAAAAAGAGATGCTCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCTTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCcgcagcgc              <  2:224339/237‑1 (MQ=255)

                                                                                                                                                                                                                                        gc‑gggcgTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCCGCAGCGCAGAgg         <  1:1756796/238‑1 (MQ=255)

                                                                                                                                                                                                                                         c‑gggcgGCAGTGT‑TGTGTTTTTC‑CCGGG‑AGACAAGTT‑GTTGGGC‑GGGCCGCAGGCAGGAATTATTGTTGATAAAAAAGAGATGACCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCCGCAGCGCAGAGGt        <  2:42328/238‑1 (MQ=255)

ACGCGAATGATTATCGTCAGGCGGTGAATGAAAATACCGCACTGTTGATGAAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGAT‑GAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTC‑GCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGC‑GGGCGTCAGTCT‑GGTGAGTTTC‑TCCGG‑CGACAAGTT‑GTTAGGC‑GGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGC‑GGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCCGCAGCGCAGAGGTCATTCA  >  NC_007779/3879785‑3880259

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