BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,386,018	G→T	13.2%	*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,018	0	G	T	13.2%	203.1 / 2.9	119	*105L (TGA→TTA)	yjeO	inner membrane protein
Reads supporting (aligned to +/- strand): ref base G (34/63); new base T (13/4); total (51/68)
Fisher's exact test for biased strand distribution p-value = 2.45e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

TGGATGATTGATGGACGCGGGATTAACAACGTCTGCGATGTGTTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAAT  >  NC_000913/4385792‑4386242

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

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

 tggtgagtggtGGACGCGGGATTAACACCGTCTGCGTTGTGTTGATGTATCTCGAAGAAGACGATTCCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGc                                                                                                                                                                                                                               <  2:580421/226‑1 (MQ=255)

 ggATGATTGATGGACGCGGGATTAACAACGTCTGCGATGTGTTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTTTCGGCTTT‑CAGTTTTGTTTTTGAGc                                                                                                                                                                                                                               >  1:685616/1‑229 (MQ=255)

   atgatTGATGGACGCGGGATTAACAACGTCTGCGATGTGTTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGGTTTTTCTGCGCTTT‑CAGTTTTGTTTTTTAGcc                                                                                                                                                                                                                              >  1:35354/1‑228 (MQ=255)

        tgatGGACGCGGGATTAACAACGTCTGCGATGTGTTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGTGCCGGCtt                                                                                                                                                                                                                         >  1:644366/1‑228 (MQ=255)

        tgatGGACGCGGGATTAACAACGTCTGCGATGTGTTGATGTATCTAGAAGAAGACGATACCCGTGATGTTGGTGTGCTAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑AGGCCACCGGGCTGGCGATTTTAGGCTACTGGC‑GTGGCGGGTTTTTTTTGCGGTGT‑AGGTTTTTTTTTTATGCCGGGTTg                                                                                                                                                                                                                        >  2:668305/1‑229 (MQ=255)

          atGGACGCGGGATTAACAACGTCTGCGATGTGTTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTTCGCTAT‑CCGTTTTTTTTTGTAGCCGGGTTgg                                                                                                                                                                                                                       >  1:739221/1‑228 (MQ=255)

            ggACGCGGGATTAACAACGTCTGCGATGTGTTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGttt                                                                                                                                                                                                                    >  1:752190/1‑228 (MQ=255)

                 ggggTTTAACAATGTTTGGTATGTGGTTATGTATCCCGAAGGAGACTTTACCCGTGGTGTTGGGTTGAAAATGCCGGTCCCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCGGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCg                                                                                                                                                                                                               <  2:527480/228‑1 (MQ=255)

                 ggggATTAAGAACTTCGGGGTTGTGTTGGTTTATCTCGAAGAAGACGATACCCGGGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑CCGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCg                                                                                                                                                                                                               <  2:137768/228‑1 (MQ=255)

                  gggATTAACAACGTCTGCGATGTGTTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCCCGCTGGGTATTTTCGGCTACTGCCGGTGGC‑CGCTTTTTAT‑CGGTCTCCAGTTTTGTTTATAGGCCGGGTTGGTTTCTCCgg                                                                                                                                                                                                              >  1:1215/1‑229 (MQ=255)

                              gTCTGCGATGTGTTGATGTATCTCGAAGAAGCCGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGTCACCCCGCTGGACA‑TTTCTGCTACTGGCTGTGGA‑CGTTTTTTCTTCGCAAT‑CAGTTTTGTTTTTTGCCCCGGTTGTCTTCTCCGGCACacaaacttaa                                                                                                                                                                                                 >  2:623815/1‑228 (MQ=255)

                              gTCTGCGATGTGTTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTTTGCCGTTT‑CCGTTTTGTTTGTTAGCCGGGTTggggttttcggccccaaacctac                                                                                                                                                                                                  >  2:704400/1‑208 (MQ=255)

                               tCTGCGATGTGTTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACacaaactta                                                                                                                                                                                                  >  1:548230/1‑228 (MQ=255)

                                cTGCGATGTGTTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACacaaacttaca                                                                                                                                                                                                <  2:655825/229‑1 (MQ=255)

                                   cGATGTGTTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTCTGGC‑AGGTTTTTCTTTGCTTT‑CCGCTTTGTTTTGTACACGGATTCGTTCTTCCGGCCCccaaacttccaaa                                                                                                                                                                                              >  1:95445/1‑228 (MQ=255)

                                        tgtTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGGGCTAT‑CAGTTTTTTTTGTGAGCCGGGTTGGTTTATCCGGGCCACAAACTTACAAACTTca                                                                                                                                                                                         >  1:352036/1‑227 (MQ=255)

                                        tgtTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGGGGC‑AGTTTTTTGTCGCCTAT‑CGGTTTTTGTTGTGAACCGGGTTGGTTTCTCCTGCACcaaaaattaaaaactttc                                                                                                                                                                                         >  2:281072/1‑226 (MQ=255)

                                        tgtTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGGGGC‑AGGTTTTTGTTCGGCAT‑TCGATTTTGTTGTTAGCCGGATTGGTTCGTCCGGGACacaaacttaaaaaatttc                                                                                                                                                                                         >  2:321505/1‑226 (MQ=255)

                                          ttGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTC‑TTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTTAGCCGGATTGGGTCATCCGGCACACCAACTTACAAACTTCCCGcc                                                                                                                                                                                     >  1:713318/1‑229 (MQ=255)

                                               tttatccagagaagtCGATTCCCTGGATGTTGGTGTGAAAATTGCGATTCCGCTGTTTTTTCCATTTCTCTGGTTGGGTTTGGGGCGAAACAAACGCGGCTGGTTGAAGT‑ACGTCACCGCACTGGCTATTTTCGGCTGCTGGCTGTGGC‑AGTTTTTGCCGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCTTCCGGCACACAAACTTACAAACTTGCCGCCTga                                                                                                                                                                                  <  2:103255/219‑1 (MQ=255)

                                               gTATCTCGAAGAAGACGATACCCGTTATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTga                                                                                                                                                                                  <  1:406407/228‑1 (MQ=255)

                                                ttttcccaagaatccgttccccTTGATTTTGGTGTGATAATGGCTTTACCGCTCTTTTTCCGTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTgac                                                                                                                                                                                 <  2:746949/209‑1 (MQ=255)

                                                  tctcGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATTCGGCACccaaactgcttcttttacacatttcccag                                                                                                                                                                              >  1:592226/1‑207 (MQ=255)

                                                  tctcGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGGGGC‑AGGTTTTTTTGCGCTAT‑CAATTTTGTTTTTGAACCGGGTTTGTTCGTCCGGGTCccaaaattacaaaatttacccctggctga                                                                                                                                                                              >  2:493501/1‑216 (MQ=255)

                                                    tcGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCAACTGGCTGTGGC‑AGTCTTTTTAGCCCTAT‑CCGTTTTTGTTGATGGCCGGGTTGGTTCATACTGCATAACAACTTACAAACTTCCCGCCTgccgcac                                                                                                                                                                             >  2:144503/1‑225 (MQ=255)

                                                       aagaagACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑CGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGGTCATCCGGCACACAAACTTAAAAACTTTCCGCCTGACGTCCttt                                                                                                                                                                          >  1:669356/1‑228 (MQ=255)

                                                            gaCGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCg                                                                                                                                                                     >  1:96916/1‑228 (MQ=255)

                                                                    ccGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGGGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCAGCCTGACGACCTTTGCCCGCAGACGTc                                                                                                                                                             >  1:24659/1‑228 (MQ=255)

                                                                     cGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTc                                                                                                                                                             <  1:22223/227‑1 (MQ=255)

                                                                       tGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTTTGCGCATT‑CAGTTTTTTTTTGCAGCCGGATTGGTTAATCCGGCACACAAAATTCACAACTTCCCCCCTGCAGCCCTTTGCCCCCCGACGCCacag                                                                                                                                                         >  2:235675/1‑226 (MQ=255)

                                                                         aTGTTGGAGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTTTTAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGTCGACCTTTGCCCGCAAACGTCAGCGtt                                                                                                                                                       >  1:166167/1‑229 (MQ=255)

                                                                          tGTTGGTTTGATAATGGCGTTTACGCTCTTTTTCCCTTTTCTCTGGTTCGTTTTGTGGGGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTtc                                                                                                                                                      <  2:675485/229‑1 (MQ=255)

                                                                           gTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTTT‑CAGTTTTTTTTTTGAGCCCGGTTGGTTTATCCGGCCCACAAACTTAAAAATTTCCCCCCCGACTCCCTTTGCCCCCCGACCTCAGCCTtct                                                                                                                                                     >  2:392289/1‑229 (MQ=255)

                                                                           gTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTAATGGCTGTGGC‑AGTTATTTCTGCGCTTT‑CAGGTTTGTTTTTGTGCCGGTCTTGTTCATCCGTCTACCAAACTTACACACTTCCCCCCTGACGCCCTTTTCACGCGAACGTcccgtctct                                                                                                                                                     >  2:176829/1‑221 (MQ=255)

                                                                                    aTAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CCGTTTTGTTTGGTAGCCGGGTTGGTTCATCCGGCACACCAACTTTCAAACTTTCCGCCCTGCGACCTTTTCCCCCAGACGTCAGCGTTTTCCCCctttg                                                                                                                                            >  2:537316/1‑229 (MQ=255)

                                                                                      aaggACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCgtttgtt                                                                                                                                          <  1:305706/226‑1 (MQ=255)

                                                                                       ttggcGTTACCCCTCTTTTCCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGTCTGGTTTATGT‑ACGCCACCGCGTTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCgtttgttt                                                                                                                                         <  2:669355/225‑1 (MQ=255)

                                                                                       gtggcGTTACCGTTTTTTCTCCCTTTTCGCGGGTTCGCTTTGTGGCGAAAAAAACTCGGCTGGTTTATGT‑ACGCCGCCGCGCTGGCTATTTTCGGCGACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTTTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCgtttgttt                                                                                                                                         <  1:345374/225‑1 (MQ=255)

                                                                                       aTGTCGTTACCGCTCTTTTTCCCTTTTCTCTGTTTCGCGTTGTGGCGAAAAAAACGGGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCCGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCgtttgttt                                                                                                                                         <  2:165611/229‑1 (MQ=255)

                                                                                       aTGACGTTACCGCTCCTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCgtttgttt                                                                                                                                         <  1:169378/229‑1 (MQ=255)

                                                                                               aCCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCCCCTGACGACCTTTGCCCCCAGACGTCCACGTTCTCCCCGGTTGTTTACCGttt                                                                                                                                  >  1:728054/1‑227 (MQ=255)

                                                                                                cccctctTTTTCCCTTTTTTCTGGTTCGCTTGGGGGACAAAAAAAGGCGGGTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTGTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAg                                                                                                                                <  2:420089/226‑1 (MQ=255)

                                                                                                cccccttttttcctttttccCTGGTTCGGTTTGTGGCGAAAAAAACGCGTCTTGTTTATGT‑ACCCCTCCGCGCTGGCTATTTGCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAg                                                                                                                                <  2:599532/223‑1 (MQ=255)

                                                                                                      tttttCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTc                                                                                                                          <  1:392289/229‑1 (MQ=255)

                                                                                                          ttccTTTTCGTGGGTTTGCTTTGGCGAGAAAAGAACCCGGCTGTGTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGa                                                                                                                      <  2:50090/227‑1 (MQ=255)

                                                                                                          tgccTTTTCTCTGGTGCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGACCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGa                                                                                                                      <  1:675673/227‑1 (MQ=255)

                                                                                                          tCCCTTTTTTCTTGTTCGGTTTGTGGGGAAAAACACGCGGGTGGTTTATGT‑ACTCCACCGCGCTTGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGa                                                                                                                      <  2:73852/229‑1 (MQ=255)

                                                                                                          tCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGa                                                                                                                      <  1:76813/229‑1 (MQ=255)

                                                                                                             cTTTTCTCTGGTTTGGTTTGGGGGGAAAAAAACCGGGCTTGGTTATGT‑ACGCCCCCGGGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGc                                                                                                                   <  2:739220/229‑1 (MQ=255)

                                                                                                             cTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGc                                                                                                                   <  2:446955/229‑1 (MQ=255)

                                                                                                             cTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCCGCTGGTTTATGT‑ACGCCACCGCGCTGGTTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAACCTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGc                                                                                                                   <  1:267141/229‑1 (MQ=255)

                                                                                                                   ctggtttccgcttGGGGGGAAAAAACCCGGGCTGGTTTAGGT‑ACGCCACCGCGCTGGCTATTTTAGGCTACTGGTTTGGGA‑AGTTGTTGCTGCGCTAG‑CAGTTTTGTTTGTGAGCAGGATTGGCTCAGCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGAATCTgg                                                                                                             <  2:216930/218‑1 (MQ=255)

                                                                                                                    ttGGGTCGTTTTGGGGGGAAAAAAACTGGGGTGGGTT‑TGTGACGGCGCCGCGCTGGCTATTTTCGGCGTCTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTGGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGa                                                                                                            <  2:35354/228‑1 (MQ=255)

                                                                                                                      ggTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTTCGCTAT‑CAGTTTTTTTTTTTTGCCGGGTTGGTTCATTCCGGTCACAAACTTACAAACTTTCCCCCTGTCGTCCTTTTCCCGCAGACGTCACCGTTTTCCCCGCTTTTTTTACGTTTGTGCTTTCCCAACGCaaccgggat                                                                                                          >  2:698659/1‑221 (MQ=255)

                                                                                                                        ttgtttttGGGGCGAAAACACCGCGGCTGTTTTATGT‑ACGCCTCCGCGCTGGCTTTTTTCGGCTCCTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAggg                                                                                                        <  1:125821/224‑1 (MQ=255)

                                                                                                                         tggCTTTGGTGCGAAAAAAACGCGGTTGGTTTATGT‑ACGCCAGCGGCTTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGc                                                                                                       <  2:551416/227‑1 (MQ=255)

                                                                                                                            cgttGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTA‑TTGCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTACCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGa                                                                                                   <  1:623816/227‑1 (MQ=255)

                                                                                                                            cTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACGGGCTGTGGC‑AGTTTTTTTTGCGCTCT‑CAGTTTTGTTTGTGAGCCGGCTTGGTTCATCCGGCCCACCACCTTACAAAATTCCCCCCCGCAGGACTTTGCACGCAGACGTCAGCGCTCTTCCCGTTTGTGTACCGTTAAGGATTTTCACGCGCCACTGGAAggggcgg                                                                                                    >  2:400078/1‑225 (MQ=255)

                                                                                                                             ttttgggCAAAAAAAACGCGGCTGGTTTATGT‑ACGCCTACGCGCGGTGTATTTTCGGCTAAGGGCTTTGGC‑AATTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGa                                                                                                   <  2:231636/224‑1 (MQ=255)

                                                                                                                             tttGTGGCGAAAAAAACGCGGGTGGTTTATGT‑ACGTCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCGGCGCTAT‑CAGGTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGa                                                                                                   <  1:165794/229‑1 (MQ=255)

                                                                                                                             tttGGGGCAAAAAAAACGGGGGTTGTTTTTTT‑TCGGCCGCGCTCTGGTTTTTTTCGGCTTCTGGCTGTGGC‑ATTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGGGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGa                                                                                                   <  2:335857/229‑1 (MQ=255)

                                                                                                                                 tgGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGGTACTGGCTGGGGC‑CGTTTTTTCTGCCCCAT‑CAGTTTTGCTTGTGACCCGGATTGGTTCATCCGGCCCACCAACTTACCAACTTACCGCCTCACGACCTTTGTCCCCAGACGTCAGCGTTCTCCCCGCTTGTTTACCGGTTAAGCCTTCCAGACCGATCCGGAAGGGGGGGagatg                                                                                               >  2:558520/1‑225 (MQ=255)

                                                                                                                                  tgggAAAAAAACGGGGTTGGTTTTTGT‑TGGCCGCCGGGCTGGTTTTTTTGGGCTAGTGGCTGTGGC‑AGTTGTTTTTGCGCTAC‑CAGTTTTGTTTGTGGGCCGGATTGGTTCATCCGGCACAAAAACTTACAAACTTCCCGCCTGACGCCCTTTGCCCGCAGACGTCAGCGTTCGCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATgg                                                                                               <  2:752189/225‑1 (MQ=255)

                                                                                                                                  gtcGAAAAAATCGCGGCTGGTTTATGT‑TCGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTTAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATgg                                                                                               <  1:235675/226‑1 (MQ=255)

                                                                                                                                   tgaaaaacatCGGGGTTGTTTTATTT‑ACGCCTCCGCGCGGGCTATTTTCGGTTACTGTGTGTGGT‑AGTTTTTGCTGCACTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGGCCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCa                                                                                             <  2:294653/226‑1 (MQ=255)

                                                                                                                                   gggtaaaaaggccgggttgTTTATGT‑ATGCCCTCGCCCTCGCTTTTGTGGGCTAGTGCATGTGTC‑TGTGTTTCCTGCGCTAT‑TATGTTTGTTTGTGAGCCGGATTGGTTCATCCGGCCCACAACCTTACAAACTTCCCGCCTGCCTACCTGTTCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGTCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCa                                                                                             <  2:662823/211‑1 (MQ=255)

                                                                                                                                      aaaaaaaCGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTGCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGa                                                                                                                      <  1:604278/201‑1 (MQ=255)

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

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

                                                                                                                                                                gCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGGTTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGACGGGCGGGAATGGCACATCGATACCGTGGGCATGTGagccagc                                                                 >  1:288605/1‑229 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                                                                                                                                            ttCGGCGTTAT‑CAGTGTTGTTTGTGAGCCGTATTGGTACATCCGGCACACAAGCTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGa                      <  1:512733/229‑1 (MQ=255)

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

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

                                                                                                                                                                                                                       ctaGTTTTTTTTGGTGGCCGGATTGGTTCACCCGGCACCCAAACCTCAAAAATTCCCGCCTGACGTCCTTTTCCCTCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAg            <  2:680392/227‑1 (MQ=255)

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

TGGATGATTGATGGACGCGGGATTAACAACGTCTGCGATGTGTTGATGTATCTCGAAGAAGACGATACCCGTGATGTTGGTGTGATAATGACGTTACCGCTCTTTTTCCCTTTTCTCTGGTTCGCTTTGTGGCGAAAAAAACGCGGCTGGTTTATGT‑ACGCCACCGCGCTGGCTATTTTCGGCTACTGGCTGTGGC‑AGTTTTTTCTGCGCTAT‑CAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAAT  >  NC_000913/4385792‑4386242

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