BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,386,017	G→T	13.1%	*105L (TGA→TTA)	yjeO →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,386,017	0	G	T	13.1%	198.9 / 2.5	113	*105L (TGA→TTA)	yjeO	inner membrane protein
Reads supporting (aligned to +/- strand): ref base G (32/60); new base T (12/4); total (48/65)
Fisher's exact test for biased strand distribution p-value = 4.58e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

TGGATGATTGATGGACGCGGGATTAACAACGTCTGCGATGTGTTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAAT  >  NC_000913/4385791‑4386241

tGGGTGATTGATGGACGCGGGGTTAACAACGTCTGCGATGTGTTGATGTATCTCGAAGAAGACGATTCCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAg                                                                                                                                                                                                                                 <  1:55063/229‑1 (MQ=255)

tGGATGTTTGATGGGCGCGGGTTTATCACCTTTTGCGATGTGTTTATGTATCCCGAAGAAGACGATACCCTGGATGTTGGTGTGATAAGGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAg                                                                                                                                                                                                                                 <  2:343984/229‑1 (MQ=255)

 tggtgagtggtGGACGCGGGATTAACACCGTCTGCGTTGTGTTGATGTATCTCGAAGAAGACGATTCCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑c                                                                                                                                                                                                                               <  2:580421/226‑1 (MQ=255)

 ggATGATTGATGGACGCGGGATTAACAACGTCTGCGATGTGTTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTTTCGGCTTT‑CAGTTTTGTTTTTGAG‑c                                                                                                                                                                                                                               >  1:685616/1‑229 (MQ=255)

   atgatTGATGGACGCGGGATTAACAACGTCTGCGATGTGTTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGGTTTTTCTGCGCTTT‑CAGTTTTGTTTTTTAG‑cc                                                                                                                                                                                                                              >  1:35354/1‑228 (MQ=255)

        tgatGGACGCGGGATTAACAACGTCTGCGATGTGTTGATGTATCTAGAAGAAGACGATACCCGTGATGTTGGTGTGCTAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑AGGCCACCGGGCTGGCGATTTTAGGCTACTGGC‑GTGGCGGGTTTTTTTTGCGGTGT‑AGGTTTTTTTTTTATG‑CCGGGTTg                                                                                                                                                                                                                        >  2:668305/1‑229 (MQ=255)

          atGGACGCGGGATTAACAACGTCTGCGATGTGTTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTTCGCTAT‑CCGTTTTTTTTTGTAG‑CCGGGTTgg                                                                                                                                                                                                                       >  1:739221/1‑228 (MQ=255)

            ggACGCGGGATTAACAACGTCTGCGATGTGTTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGttt                                                                                                                                                                                                                    >  1:752190/1‑228 (MQ=255)

               cgcgGGATTAACAACGTCTGCGATGTGTTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTTTGCGCTAT‑CAGTTTTGTTTTTTAT‑CCGGATTGGTTCATc                                                                                                                                                                                                                 >  1:73852/1‑229 (MQ=255)

                 ggggTTTAACAATGTTTGGTATGTGGTTATGTATCCCGAAGGAGACTTTACCCGTGGTGTTGGGTTGAAAATGCCGGTCCCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCG‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCg                                                                                                                                                                                                               <  2:527480/228‑1 (MQ=255)

                 ggggATTAAGAACTTCGGGGTTGTGTTGGTTTATCTCGAAGAAGACGATACCCGGGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑CCGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCg                                                                                                                                                                                                               <  2:137768/228‑1 (MQ=255)

                  gggATTAACAACGTCTGCGATGTGTTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCCCGCTGGGTATTTTCGGCTACTGCCGGTGGC‑CGCTTTTTAT‑CGGTCTCCAGTTTTGTTTATAGG‑CCGGGTTGGTTTCTCCgg                                                                                                                                                                                                              >  1:1215/1‑229 (MQ=255)

                              gTCTGCGATGTGTTGATGTATCTCGAAGAAGCCGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGTCACCCCGCTGGACA‑TTTCTGCTACTGGCTGTGGA‑CGTTTTTTCTTCGCAAT‑CAGTTTTGTTTTTTGC‑CCCGGTTGTCTTCTCCGGCACacaaacttaa                                                                                                                                                                                                 >  2:623815/1‑228 (MQ=255)

                               tCTGCGATGTGTTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACacaaactta                                                                                                                                                                                                  >  1:548230/1‑228 (MQ=255)

                                cTGCGATGTGTTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACacaaacttaca                                                                                                                                                                                                <  2:655825/229‑1 (MQ=255)

                                       gtgtTGATTGTTTTCGAAGAAACCGATGCCCGTGTTGTTGGTTTGATAAGTACGTTACCGGTCTTTGTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTAGTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTcc                                                                                                                                                                                         <  1:297307/229‑1 (MQ=255)

                                        tgtTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGGGCTAT‑CAGTTTTTTTTGTGAG‑CCGGGTTGGTTTATCCGGGCCACAAACTTACAAACTTca                                                                                                                                                                                         >  1:352036/1‑227 (MQ=255)

                                        tgtTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGGGGC‑AGTTTTTTGTCGCCTAT‑CGGTTTTTGTTGTGAA‑CCGGGTTGGTTTCTCCTGCACcaaaaattaaaaactttc                                                                                                                                                                                         >  2:281072/1‑226 (MQ=255)

                                        tgtTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGGGGC‑AGGTTTTTGTTCGGCAT‑TCGATTTTGTTGTTAG‑CCGGATTGGTTCGTCCGGGACacaaacttaaaaaatttc                                                                                                                                                                                         >  2:321505/1‑226 (MQ=255)

                                               tttatccagagaagtCGATTCCCTGGATGTTGGTGTGAAAATTGCGATTCCGCTGTTTTTTCCATTTCTCTGGTTGGGTTTGGGGCGAAACAAACGCGGCT‑GGTTGAAGT‑ACGTCACCGCACTGGCTATTTTCGGCTGCTGGCTGTGGC‑AGTTTTTGCCGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCTTCCGGCACACAAACTTACAAACTTGCCGCCTga                                                                                                                                                                                  <  2:103255/219‑1 (MQ=255)

                                                tATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTCCTGGCTGGGGC‑AGTTTTTTCTGTGCGCT‑CCGTTTTGTTTTTGAG‑CCGGGTTGGTTTTTCCGGCACACAAACTTTCAAAATTCCCGCCCgaa                                                                                                                                                                                 >  2:722977/1‑227 (MQ=255)

                                                  tctcGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACTCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGA‑AGATTTTTCTTCCCCTT‑CACTTTT‑TTTGTGTGACCGGATTGTTTCTTCCTGCCCCCAAACTTACCAACATcccccccatagac                                                                                                                                                                              >  2:675672/1‑222 (MQ=255)

                                                  tctcGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTTGGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCTTCCGGCACACAAACTTACAAACTTCCCGCCTtacgc                                                                                                                                                                               >  1:551417/1‑227 (MQ=255)

                                                  tctcGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTTCGCTAT‑CAGTTTTGTTTTTGAG‑CCGGATTGGGTCATCCGGCACACCAACTTACAAACTTCCCGCCTgacggc                                                                                                                                                                              >  1:231636/1‑227 (MQ=255)

                                                    tcGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCAACTGGCTGTGGC‑AGTCTTTTTAGCCCTAT‑CCGTTTTTGTTGATGG‑CCGGGTTGGTTCATACTGCATAACAACTTACAAACTTCCCGCCTgccgcac                                                                                                                                                                             >  2:144503/1‑225 (MQ=255)

                                                     cGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTTCGCTAT‑CATTTTTGTTTGTGAG‑CCGGGTTGGTTACTCCGGAGCACAACCTTACAAAATTCCCCCCTGACGAACt                                                                                                                                                                            >  2:267141/1‑228 (MQ=255)

                                                       aagaagACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑CGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGGTCATCCGGCACACAAACTTAAAAACTTTCCGCCTGACGTCCttt                                                                                                                                                                          >  1:669356/1‑228 (MQ=255)

                                                            gaCGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCg                                                                                                                                                                     >  1:96916/1‑228 (MQ=255)

                                                            gaCGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGGTGC‑AGCTTTTTCTGCGCTAT‑CAGATTTGTTTTGTAT‑CCGGTTTGGTTCATCCTGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCg                                                                                                                                                                     >  1:95914/1‑228 (MQ=255)

                                                                    ccGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGGGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCAGCCTGACGACCTTTGCCCGCAGACGTc                                                                                                                                                             >  1:24659/1‑228 (MQ=255)

                                                                     cGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTc                                                                                                                                                             <  1:22223/227‑1 (MQ=255)

                                                                       tGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTTTGCGCATT‑CAGTTTTTTTTTGCAG‑CCGGATTGGTTAATCCGGCACACAAAATTCACAACTTCCCCCCTGCAGCCCTTTGCCCCCCGACGCCacag                                                                                                                                                         >  2:235675/1‑226 (MQ=255)

                                                                          tGTTGGTTTGATAATGGCGTTTACGCTCTTTTTCCCTTTTCTCTGGTTCGTTTTGTGGGGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTtc                                                                                                                                                      <  2:675485/229‑1 (MQ=255)

                                                                           gTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTTT‑CAGTTTTTTTTTTGAG‑CCCGGTTGGTTTATCCGGCCCACAAACTTAAAAATTTCCCCCCCGACTCCCTTTGCCCCCCGACCTCAGCCTtct                                                                                                                                                     >  2:392289/1‑229 (MQ=255)

                                                                                    aTAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CCGTTTTGTTTGGTAG‑CCGGGTTGGTTCATCCGGCACACCAACTTTCAAACTTTCCGCCCTGCGACCTTTTCCCCCAGACGTCAGCGTTTTCCCCctttg                                                                                                                                            >  2:537316/1‑229 (MQ=255)

                                                                                      aaggACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCgtttgtt                                                                                                                                          <  1:305706/226‑1 (MQ=255)

                                                                                       ttggcGTTACCCCTCTTTTCCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGTCT‑GGTTTATGT‑ACGCCACCGCGTTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCgtttgttt                                                                                                                                         <  2:669355/225‑1 (MQ=255)

                                                                                       gtggcGTTACCGTTTTTTCTCCCTTTTCGCGGGTTCGCTTTGTGGCGAAAAAAACTCGGCT‑GGTTTATGT‑ACGCCGCCGCGCTGGCTATTTTCGGCGACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTTTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCgtttgttt                                                                                                                                         <  1:345374/225‑1 (MQ=255)

                                                                                       aTGTCGTTACCGCTCTTTTTCCCTTTTCTCTGTTTCGCGTTGTGGCGAAAAAAACGGGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCCGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCgtttgttt                                                                                                                                         <  2:165611/229‑1 (MQ=255)

                                                                                       aTGACGTTACCGCTCCTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCgtttgttt                                                                                                                                         <  1:169378/229‑1 (MQ=255)

                                                                                              tACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTCCGCTAT‑CCGTTTTGTTTGTAAG‑CCGGATTGGTTCATCCGGCACCCAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGACAGCGTTCTCCCCGTTTGATCACCgct                                                                                                                                   >  1:146222/1‑226 (MQ=255)

                                                                                               aCCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCCCCTGACGACCTTTGCCCCCAGACGTCCACGTTCTCCCCGGTTGTTTACCGttt                                                                                                                                  >  1:728054/1‑227 (MQ=255)

                                                                                                cccctctTTTTCCCTTTTTTCTGGTTCGCTTGGGGGACAAAAAAAGGCGGGT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTGTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAg                                                                                                                                <  2:420089/226‑1 (MQ=255)

                                                                                                      tttttCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTc                                                                                                                          <  1:392289/229‑1 (MQ=255)

                                                                                                      tttttCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGCCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCtt                                                                                                                            >  1:49289/1‑227 (MQ=255)

                                                                                                          ttccTTTTCGTGGGTTTGCTTTGGCGAGAAAAGAACCCGGCT‑GTGTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGa                                                                                                                      <  2:50090/227‑1 (MQ=255)

                                                                                                          tgccTTTTCTCTGGTGCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGACCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGa                                                                                                                      <  1:675673/227‑1 (MQ=255)

                                                                                                          tCCTTTTTTGCTGTTTCCTTTTTTGGCGAAAAAAACCGGGTT‑GGTTTATGT‑ACGCCCCCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGa                                                                                                                      <  2:348744/229‑1 (MQ=255)

                                                                                                          tCCCTTTTTTCTTGTTCGGTTTGTGGGGAAAAACACGCGGGT‑GGTTTATGT‑ACTCCACCGCGCTTGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGa                                                                                                                      <  2:73852/229‑1 (MQ=255)

                                                                                                          tCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGa                                                                                                                      <  1:76813/229‑1 (MQ=255)

                                                                                                             cTTTTCTCTGGTTTGGTTTGGGGGGAAAAAAACCGGGCT‑TGGTTATGT‑ACGCCCCCGGGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGc                                                                                                                   <  2:739220/229‑1 (MQ=255)

                                                                                                             cTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGc                                                                                                                   <  2:446955/229‑1 (MQ=255)

                                                                                                             cTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCCGCT‑GGTTTATGT‑ACGCCACCGCGCTGGTTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAACCTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGc                                                                                                                   <  1:267141/229‑1 (MQ=255)

                                                                                                                   ctggtttccgcttGGGGGGAAAAAACCCGGGCT‑GGTTTAGGT‑ACGCCACCGCGCTGGCTATTTTAGGCTACTGGTTTGGGA‑AGTTGTTGCTGCGCTAG‑CAGTTTTGTTTGTGAG‑CAGGATTGGCTCAGCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGAATCTgg                                                                                                             <  2:216930/218‑1 (MQ=255)

                                                                                                                    ttGGGTCGTTTTGGGGGGAAAAAAACTGGGGT‑GGGTT‑TGTGACGGCGCCGCGCTGGCTATTTTCGGCGTCTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTGGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGa                                                                                                            <  2:35354/228‑1 (MQ=255)

                                                                                                                      ggTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTTCGCTAT‑CAGTTTTTTTTTTTTG‑CCGGGTTGGTTCATTCCGGTCACAAACTTACAAACTTTCCCCCTGTCGTCCTTTTCCCGCAGACGTCACCGTTTTCCCCGCTTTTTTTACGTTTGTGCTTTCCCAACGCaaccgggat                                                                                                          >  2:698659/1‑221 (MQ=255)

                                                                                                                        ttgtttttGGGGCGAAAACACCGCGGCT‑GTTTTATGT‑ACGCCTCCGCGCTGGCTTTTTTCGGCTCCTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAggg                                                                                                        <  1:125821/224‑1 (MQ=255)

                                                                                                                         tggCTTTGGTGCGAAAAAAACGCGGTT‑GGTTTATGT‑ACGCCAGCGGCTTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGc                                                                                                       <  2:551416/227‑1 (MQ=255)

                                                                                                                            cgttGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTA‑TTGCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTACCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGa                                                                                                   <  1:623816/227‑1 (MQ=255)

                                                                                                                            cTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACGGGCTGTGGC‑AGTTTTTTTTGCGCTCT‑CAGTTTTGTTTGTGAG‑CCGGCTTGGTTCATCCGGCCCACCACCTTACAAAATTCCCCCCCGCAGGACTTTGCACGCAGACGTCAGCGCTCTTCCCGTTTGTGTACCGTTAAGGATTTTCACGCGCCACTGGAAggggcgg                                                                                                    >  2:400078/1‑225 (MQ=255)

                                                                                                                             ttttgggCAAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCTACGCGCGGTGTATTTTCGGCTAAGGGCTTTGGC‑AATTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGa                                                                                                   <  2:231636/224‑1 (MQ=255)

                                                                                                                             tttGTGGCGAAAAAAACGCGGGT‑GGTTTATGT‑ACGTCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCGGCGCTAT‑CAGGTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGa                                                                                                   <  1:165794/229‑1 (MQ=255)

                                                                                                                                 tgGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCCCCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATgg                                                                                               <  2:548230/229‑1 (MQ=255)

                                                                                                                                 tgGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGGTACTGGCTGGGGC‑CGTTTTTTCTGCCCCAT‑CAGTTTTGCTTGTGAC‑CCGGATTGGTTCATCCGGCCCACCAACTTACCAACTTACCGCCTCACGACCTTTGTCCCCAGACGTCAGCGTTCTCCCCGCTTGTTTACCGGTTAAGCCTTCCAGACCGATCCGGAAGGGGGGGagatg                                                                                               >  2:558520/1‑225 (MQ=255)

                                                                                                                                  gtcGAAAAAATCGCGGCT‑GGTTTATGT‑TCGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTTAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATgg                                                                                               <  1:235675/226‑1 (MQ=255)

                                                                                                                                  ggtgAAAAAACCGCGGGTGGGTTTTTGT‑‑CGCCCTCGGGCTGGTTGTTTTCGGCGACGGGGTGTGGC‑AGTTTTTTCTGCGTTAT‑CCGTTTTGTTTGTGTG‑CCGGATTGGTTCATCCGGCACACAAACTTCCAACCTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATgg                                                                                               <  2:645062/225‑1 (MQ=255)

                                                                                                                                   tgaaaaacatCGGGGTT‑GTTTTATTT‑ACGCCTCCGCGCGGGCTATTTTCGGTTACTGTGTGTGGT‑AGTTTTTGCTGCACTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGGCCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCa                                                                                             <  2:294653/226‑1 (MQ=255)

                                                                                                                                   gCGAAAAAAAGGCGGCT‑GGTTTATGT‑ACGCCCCCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCa                                                                                             <  1:624951/229‑1 (MQ=255)

                                                                                                                                      aaaaaaaCGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTGCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGa                                                                                                                      <  1:604278/201‑1 (MQ=255)

                                                                                                                                                            tgttccgcACCGCCTTGGCTAGTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAag                                                                      <  1:321505/222‑1 (MQ=255)

                                                                                                                                                              ttcgcctctgtgcTGGCTATTTTCGGCTTCTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CATTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCTCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAagc                                                                     <  1:144503/217‑1 (MQ=255)

                                                                                                                                                                 gCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAagcca                                                                   >  1:417440/1‑227 (MQ=255)

                                                                                                                                                                 gCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGGTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGTTACCGTGGGCATGGAagaca                                                                   >  1:511762/1‑227 (MQ=255)

                                                                                                                                                                 gCCACCCCGGTTGCTATTGTCGGGTTCTTGTTTTTGC‑GGTTTTTTCGTCTCTAT‑CGATTTTGTTTGGTGG‑CCGGATTGGTTCCTCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTGCCAGACGCATCTGGAAGGGCGGGAATGGCATATGGATACCGTGGGCATGGAagccagc                                                                 <  2:24659/229‑1 (MQ=255)

                                                                                                                                                                      cgcgCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGGTTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCCAACGTCAGGGTTCTCCCCGTTTGCTTACCGTTTAGGCTTTCCAGACGCCACTGGAAGGGCGGGAATGGCATATCGGTACCGTGGGGAAGGAAACCAGACaaaa                                                            >  2:722105/1‑229 (MQ=255)

                                                                                                                                                                      cgcgCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCCGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGag                                                            >  1:590416/1‑228 (MQ=255)

                                                                                                                                                                        cgcTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATc                                                          >  1:557343/1‑229 (MQ=255)

                                                                                                                                                                         tgttgCTATTTTCTGCTTGTGGCTGTGGC‑AGTTTTCTCTTACCTAT‑CGGGTTTTGTTGTTAG‑CGGGAGTGGTTGATCCGGGACTCAAACTTACAAACTTCCCGGCTGACGACCTTTGCCCGCAGACGTCAGTGTTCTCCCCGTTTGTTTGCCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATTGAAGCCATCCAGAATc                                                          <  1:131066/224‑1 (MQ=255)

                                                                                                                                                                         gcTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATc                                                          <  1:129885/228‑1 (MQ=255)

                                                                                                                                                                          gttgCTATTTTCCGCTACTTGCTGGTGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTCCAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATGCAGTGGGCATGGAAGCCAGCCAGAATCAg                                                        <  1:668306/226‑1 (MQ=255)

                                                                                                                                                                               ttttttCGGTTTCTGGCTGTGGC‑ATTTTTTTCTGCGCTAT‑CCGTTTTGTTTGTGAG‑CCGGATTGGTTCCACCGGCACCCAAACTTACAAACTTCCCGACTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCtggt                                                   <  2:171777/227‑1 (MQ=255)

                                                                                                                                                                                  tttCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGAAACCGTGGGCATGGAAACCAGCCAGAATAAGCTGGTGGa                                                >  1:420712/1‑229 (MQ=255)

                                                                                                                                                                                   ttCGGCTACTGGCTGGGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGGTTGGGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGa                                                <  1:665655/228‑1 (MQ=255)

                                                                                                                                                                                        cTACTGGCTGTGGC‑AGTTGTTTCTGCGCTAG‑CAGTGTTGTTTGTGAG‑CCGGATGGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCa                                           <  1:436924/228‑1 (MQ=255)

                                                                                                                                                                                           cTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTACAGCCTGACGACCTTTGCCCGCAGAAGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCCGAATCAGCTGGTGGGTCCCTTGGc                                       >  1:396340/1‑229 (MQ=255)

                                                                                                                                                                                              ggtgtgGC‑AGTTTTTTCGGCGCTAT‑CCGTTTTGTTTGTGGG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTGCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTa                                    <  2:3737/227‑1 (MQ=255)

                                                                                                                                                                                              gCTGTGGC‑CTTTTTTTCTGCGCTCT‑CATTTTGGTTTGTCAG‑CCGGATTGGTTCAACCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAAGGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTa                                    <  1:698660/229‑1 (MQ=255)

                                                                                                                                                                                              gCTGGGGC‑AGTTTTTTCTGCGGTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGGCGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTa                                    <  1:722267/229‑1 (MQ=255)

                                                                                                                                                                                               gtttttgccgTTTTTTCTCCACTTT‑CAGTTTTGTTTGTGTG‑CCGGATTGTTTCATCCGGCACACAAACCTACAAACTTCCCGCCTGCCGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTGTGTTTACCGGTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTa                                    <  2:728053/220‑1 (MQ=255)

                                                                                                                                                                                                    gC‑AGTTTTTTCGGCGGTAT‑CAGTTTTTGTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCa                              <  1:79770/229‑1 (MQ=255)

                                                                                                                                                                                                    gC‑AGTTTTTTCCGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCa                              <  1:704401/229‑1 (MQ=255)

                                                                                                                                                                                                     tcatttttttCTGGGCTAT‑CAGGTTGTGTTGGGGG‑CCGGATGGGTTCAGCCGGCAAACAAACTTACAAACTTCCCGCCTGACGACCTGTGCCCGCAGGCGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCa                              <  2:644365/225‑1 (MQ=255)

                                                                                                                                                                                                      tatttttttCTGCGCTTT‑CAGTTTTGTTTGTGAG‑CCGGATTTGTTCATCCGGCACGCAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCAt                             <  2:166167/226‑1 (MQ=255)

                                                                                                                                                                                                         ttttttCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTgg                                                                                                             >  2:697858/1‑146 (MQ=255)

                                                                                                                                                                                                         ttttttCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTgg                                                                                                             <  1:697859/146‑1 (MQ=255)

                                                                                                                                                                                                         ttttttCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCTTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATAc                           >  1:266908/1‑228 (MQ=255)

                                                                                                                                                                                                           ttttCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCAGCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACgg                         <  1:133979/228‑1 (MQ=255)

                                                                                                                                                                                                                      aT‑CAGTTTTGTTTGGGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAACCTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGTCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCa             <  1:493501/229‑1 (MQ=255)

                                                                                                                                                                                                                       t‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGTCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAg            <  1:409963/229‑1 (MQ=255)

                                                                                                                                                                                                                                  ttgtGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGCTGGATCTCATGGCCTACCGGCATACGGTGACCCATCTCAGCCGCGTaaaa  >  1:175853/1‑228 (MQ=255)

TGGATGATTGATGGACGCGGGATTAACAACGTCTGCGATGTGTTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCT‑GGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAG‑CCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAAT  >  NC_000913/4385791‑4386241

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