BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	696,325	G→T	5.6%	intergenic (+49/+105)	ubiF → / ← glnX	2‑octaprenyl‑3‑methyl‑6‑methoxy‑1,4‑benzoquinol oxygenase/tRNA‑Gln

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	696,325	0	G	T	5.6%	375.8 / 10.4	108	intergenic (+49/+105)	ubiF/glnX	2‑octaprenyl‑3‑methyl‑6‑methoxy‑1,4‑benzoquinol oxygenase/tRNA‑Gln
Reads supporting (aligned to +/- strand): ref base G (55/47); new base T (3/3); total (58/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-02

Only 100 of 110 total aligned reads displayed.

TGGGCCAGTTATCCTGTCCTCAAGCGTTACCAGATGCGGCGCATGGCGGATAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAA  >  NC_000913/696073‑696572

tGGGCCAGTTTTCCTGTCCTCAAGCGTTACCAGATGCGGGGCATGGGGGATAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTg                                                                                                                                                                                                                                     <  4:437474/274‑1 (MQ=255)

tGGGCCAGTTATCCTGTCCTCAAGCGTTACCAGATGCGGCGCATGGCGGATAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGAGTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTg                                                                                                                                                                                                                                     <  4:444191/274‑1 (MQ=255)

             cTGTCCTCAAGCGTTACCAGATGCGGCGCATGGCGGATAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCGTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGAAGCGGCAGCGGCGCAACCGACATTTAAgcggtccctta                                                                                                                                                                                                                        >  4:400053/1‑264 (MQ=255)

                 ccTCAAGCGTTACCAGATGCGGCGCATGGCGGATAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATc                                                                                                                                                                                                                                             <  2:400054/249‑1 (MQ=255)

                     aaGCGTTACCAGATGCGGCGCATGGCGGATAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCAACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCACCCGTAGGCCGGATAAGACGCGGCAGCGTCGCATCCGCCCTTGAAGGATAAGCCGTGTCAACATc                                                                                                                                                                                                               >  2:350620/1‑275 (MQ=255)

                          ttACCAGATGCGGCGCATGGCGGATAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCAtt                                                                                                                                                                                                          <  4:262074/275‑1 (MQ=255)

                          ttACCAGATGCGGCGCATGGCGGATAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACAACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGAGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGAGGTGTCAAAATCGCAtt                                                                                                                                                                                                          >  2:510167/1‑275 (MQ=255)

                           tACCAGATGCGGCGCATGGCGGATAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTc                                                                                                                                                                                                         <  2:227680/275‑1 (MQ=255)

                            aCCAGATGCGGCGCATGGCGGATAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCCCCACTGCGTTTTATTCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGAGTGTTGACACGTCAGGCGCTGAAATATGCGG‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTc                                                                                                                                                                                                         >  1:50734/1‑274 (MQ=255)

                                       cgcATGGCGGATAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATTCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATCAGACGTGCCAACATCGCATTCGACCTtgaatga                                                                                                                                                                                             >  2:72220/1‑275 (MQ=255)

                                            ggCGGATAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGCAGCGGCAGCGTCCCCTCCCACATTGCAGGGAAAGAAGTTTCAACATCGCATTTGAACTttaatgaccac                                                                                                                                                                                         >  4:495734/1‑270 (MQ=255)

                                                gATAACTTCATTATGCAAAGCGGTATGGGTCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGaa                                                                                                                                                                                     <  4:84053/274‑1 (MQ=255)

                                                         aTTATGCAAAGCGGGTTGGTTCTGTTTTATGCCGTATTCAGCAATAATCTGCCACCCCTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAg                                                                                                                                                                           <  1:1615/275‑1 (MQ=255)

                                                           tATGCAAAGCGTGATGGATCGGTTTTATGCCGGATTTAGCAATAATCAGCCACCACTGCGTTTTATGCTTAATAGCGGGTTAATGTCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGc                                                                                                                                                                          <  4:72219/274‑1 (MQ=255)

                                                             tGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TAGGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGGCGGATAATACGCGGCAGCGTCGCATCAGACATTGCAGGATAAGACGTGTCAACATCGCATTCGACATCGAATGAACGCCGAAAAGCAAAAAGATc                                                                                                                                                                        >  2:414275/1‑275 (MQ=255)

                                                                  aGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAACTATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATCAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGAATTCGACATTGGATGAACGCAGAAAAGCAAACAGCGCGccc                                                                                                                                                                    >  2:62746/1‑272 (MQ=255)

                                                                  aGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTATAACATTGCCGGATGCGTGCAAACCGTCGGACGGATAAGACGCGGCAGGGTCGCCACCGCCACTGAAGGATAAGACGTGTACAAATCACAATTGGAAGTGAANGAAACNNGAAAAACaaaaaaaaaaaccaa                                                                                                                                                                  >  3:1615/1‑265 (MQ=255)

                                                                         tGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCgc                                                                                                                                                             >  2:102310/1‑271 (MQ=255)

                                                                                gTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTaa                                                                                                                                                   <  2:528743/275‑1 (MQ=255)

                                                                                 ttttATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAAt                                                                                                                                                  >  2:51995/1‑275 (MQ=255)

                                                                                 ttttATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTAAACATCGGATTCGACATTGAATGAACGCAGAAAAGGAAAAAGCTAGCCGAAAGCGA‑CTTTATTAAt                                                                                                                                                  >  2:134652/1‑275 (MQ=255)

                                                                                      tGCCGTAGTCAGCAATAATCTGCCACCATTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGc                                                                                                                                             <  4:510166/275‑1 (MQ=255)

                                                                                      tGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTATTGTGGc                                                                                                                                             >  1:146811/1‑275 (MQ=255)

                                                                                      tGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCAACCGACATTGAAGGATAACACGTGTCAAAACTCCATTCGACATTGAATGACGCCAGAAACACAACAAGCCCGCAG‑AACCCAGCTTTTTTAAtcctg                                                                                                                                              >  4:110602/1‑270 (MQ=255)

                                                                                      tGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTAGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTTAAACGTCAGGCGCTGAATTATGCGG‑TA‑GGGGTGTAGCCTTACAACATTGCCGGATGGTTGCCAACCGTAGGTCGGATAAGACGCGGCAGCCGTGAATCCGTCATCGAAGCATAAGACGACTCAACATCGCATTCCGCATTGAATGACCGAAGCAGAGCAAAAAGCCCGGCC‑AAGCGAGCGTTTGTAAtggt                                                                                                                                               >  4:139828/1‑271 (MQ=255)

                                                                                          ggATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATtgaat                                                                                                                                                                                               <  2:462367/222‑1 (MQ=255)

                                                                                          ggATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAg                                                                                                                                                            <  2:110603/256‑1 (MQ=255)

                                                                                           gATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGCCATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGCTgggg                                                                                                                                        >  2:241949/1‑275 (MQ=255)

                                                                                           gATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAAAATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATGCGAAATTGAAGGAAAAAGAGTGTCAACATGCCAGTCGAAATTGAATGAACGACGAAAAGCAAAAAGCCCGCCG‑AAGAGAGCTTTTTTAATGTGGGTgggg                                                                                                                                        >  4:353218/1‑275 (MQ=255)

                                                                                                       atCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCg                                                                                                                            >  1:72627/1‑275 (MQ=255)

                                                                                                                 acTGCTTTTTATTCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGa                                                                                                                   <  2:485592/274‑1 (MQ=255)

                                                                                                                   tGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACCTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGGGGCGCATCCGACATTGAAGGATAAGACGTGTCAAAATAGGATTCGACATTGGATGAAAGCAGAAAAGCAAAAAACACAACG‑AAGCGAGCTTTTTTAAAGGGACTTGGGTACGCGGATTCGAAACCCGGaaa                                                                                                                 >  4:485591/1‑273 (MQ=255)

                                                                                                                                tAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATGCGACGCTGCCGCGTcttat                                                                                                                                                                                                                                              >  7:116869/1‑133 (MQ=255)

                                                                                                                                tAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATGCGACGCTGCCGCGTcttat                                                                                                                                                                                                                                              <  8:116869/137‑5 (MQ=255)

                                                                                                                                          ttAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCTGCCGCGTctta                                                                                                                                                                                                                                               >  5:679809/1‑123 (MQ=255)

                                                                                                                                          ttAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCTGCCGCGTctta                                                                                                                                                                                                                                               <  6:679809/126‑4 (MQ=255)

                                                                                                                                          ttAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGaa                                                                                                                                                                                                                                   >  8:1068609/1‑138 (MQ=255)

                                                                                                                                            aaTGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAg                                                                                                                                                                >  2:200590/1‑202 (MQ=255)

                                                                                                                                                gcggcgGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAg                                                                                                                                                            <  4:200589/202‑1 (MQ=255)

                                                                                                                                                 cggcggAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGt                                                                                                                                                                                                                                                   >  8:829015/1‑115 (MQ=255)

                                                                                                                                                 cggcggAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGt                                                                                                                                                                                                                                                   <  7:829015/115‑1 (MQ=255)

                                                                                                                                                   gcggAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTg                                                                                 >  4:50304/1‑274 (MQ=255)

                                                                                                                                                    cggAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCa                                                                                                                                                                                                                                                       >  7:83738/1‑108 (MQ=255)

                                                                                                                                                    cggAGCGTACTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCa                                                                                                                                                                                                                                                       <  8:83738/108‑1 (MQ=255)

                                                                                                                                                       aGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGa                                                                             <  2:50305/274‑1 (MQ=255)

                                                                                                                                                                   ggtTGAAACGTGAGGCGGGGAAAAATGCGTGTG‑GGGTTGTAGCCTTACAACAGTGCCGGATGCGGACCAACCGTAGGTCGGGTAAGACGCGGCAGCGTCGCATCCGACAATGAAGGATAAGACGTGTGAACATCGCATTGGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTGTTTAATGTGGCTGGGGTACGAGGATGAGAACCTCGGAAGGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTg                                                                  <  4:414274/273‑1 (MQ=38)

                                                                                                                                                                    gtTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAAATggg                                                                <  2:129593/274‑1 (MQ=255)

                                                                                                                                                                      tGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCgg                                                                                                                                                                                                                                                         >  6:409652/1‑88 (MQ=255)

                                                                                                                                                                       gAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGc                                                                                                               <  2:305746/224‑1 (MQ=255)

                                                                                                                                                                          aCGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCAtt                                                                                                                                                                                                          <  7:211267/131‑1 (MQ=255)

                                                                                                                                                                          aCGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCAtt                                                                                                                                                                                                          >  8:211267/1‑131 (MQ=255)

                                                                                                                                                                             tCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGt                                                                                                                                                                                                                                                   <  7:1091282/87‑1 (MQ=255)

                                                                                                                                                                             tCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGt                                                                                                                                                                                                                                                   >  8:1091282/1‑87 (MQ=255)

                                                                                                                                                                               aGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGaaa                                                                                                                                                                                    <  5:316162/148‑1 (MQ=255)

                                                                                                                                                                               aGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAg                                                    >  4:140407/1‑275 (MQ=255)

                                                                                                                                                                                      gAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGaaa                                                                                                                                                                                    <  7:377896/141‑1 (MQ=255)

                                                                                                                                                                                      gAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGaaa                                                                                                                                                                                    >  8:377896/1‑141 (MQ=255)

                                                                                                                                                                                      gAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCaaaa                                                                                                                                                                             >  6:682849/1‑148 (MQ=255)

                                                                                                                                                                                         atatGCGT‑TA‑GGGTTATAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATtgaatga                                                                                                                                                                                             <  7:935467/129‑1 (MQ=255)

                                                                                                                                                                                          caaaaagcgttggggTGTGGCCTTACAACAGAGCCGGAGGAGTGCCAACCGTAGGTCGGGTAAGACGCGGCCGGGTCGCATCCGACATTGAAGTATAGGACGTGTCAACGTCGCATTGGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTGTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTc                                           <  4:434698/264‑1 (MQ=37)

                                                                                                                                                                                           atGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGaaa                                                                                                                                                                                    <  5:997459/136‑1 (MQ=255)

                                                                                                                                                                                           atGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGaaa                                                                                                                                                                                    >  6:997459/1‑136 (MQ=255)

                                                                                                                                                                                            tGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGCAGGATAAGACGTGTCAACATCGCATTCGACCTTGAATAAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGATTTTTTAATGTGGCTGGGGGACGCGGATTCGAACCTAGGAATGCCGGAATCCGAATCCGGTGCCTTACCGCTTGGCAATACCCCAATCGTGTGCACTTCCAAGGTCAGCATCttta                                       >  4:414617/1‑273 (MQ=255)

                                                                                                                                                                                            tGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCAtt                                                                                                                                                                                                          <  7:379574/113‑1 (MQ=255)

                                                                                                                                                                                            tGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCAtt                                                                                                                                                                                                          >  8:379574/1‑113 (MQ=255)

                                                                                                                                                                                            tGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCAAACCCCGGAATGCCGGAATCAGAATCCGGTGCCTTAGCCCTTGGCGATACCCCAACTGGGTGCCATTCAAAGGTaccggcctata                                       >  4:88799/1‑265 (MQ=39)

                                                                                                                                                                                                   a‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAg                                                                                                                                                                <  5:682849/148‑1 (MQ=255)

                                                                                                                                                                                                         tctcgCCCTACAACAGTGCCGGATGCGTGACAACCATAGGTCGGATAAGACGCGGCAGCGTCGCATACGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCt                            <  4:134651/271‑1 (MQ=255)

                                                                                                                                                                                                         tGTAGCAGTACAACATTGCCGGATGGGTACCAGCAGTAGGTAGGAGAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCt                            <  1:11186/275‑1 (MQ=255)

                                                                                                                                                                                                          gTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTg                           <  2:261229/275‑1 (MQ=255)

                                                                                                                                                                                                           tAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGGCATTACAATGTAAGCGTCTTGAATAAATTTGCTgg                          >  2:364237/1‑275 (MQ=255)

                                                                                                                                                                                                           tAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTgg                          >  4:137052/1‑275 (MQ=255)

                                                                                                                                                                                                                ttACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATc                                                                                                                                                                                                                                             <  6:1025675/60‑1 (MQ=255)

                                                                                                                                                                                                                ttACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATc                                                                                                                                                                                                                                             >  5:1025675/1‑60 (MQ=255)

                                                                                                                                                                                                                 tACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCATCGTCGCAGCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACg                    <  2:215757/275‑1 (MQ=255)

                                                                                                                                                                                                                  acaacaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGTGTGCACTTGAATGGTAAGCGTCTTGAATAAATTGGCCGGGGTGCGa                   >  2:396628/1‑275 (MQ=255)

                                                                                                                                                                                                                   caacaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACg                                                                                                                                                                                                                         <  8:831891/77‑1 (MQ=255)

                                                                                                                                                                                                                   caacaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACg                                                                                                                                                                                                                         >  7:831891/1‑77 (MQ=255)

                                                                                                                                                                                                                    aacaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTAc                                                       <  2:492471/237‑1 (MQ=255)

                                                                                                                                                                                                                     acaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCAt                                                                                                                                                                                                           >  5:296933/1‑89 (MQ=255)

                                                                                                                                                                                                                     acaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCAt                                                                                                                                                                                                           <  6:296933/89‑1 (MQ=255)

                                                                                                                                                                                                                      caTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGa                                                                                                                                                                                                       >  5:342642/1‑92 (MQ=255)

                                                                                                                                                                                                                        ttGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGaaa                                                                                                                                                                                    >  6:949284/1‑109 (MQ=255)

                                                                                                                                                                                                                        ttGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGaaa                                                                                                                                                                                    <  5:949284/109‑1 (MQ=255)

                                                                                                                                                                                                                        ttGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAAATGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGAtt              >  2:43861/1‑274 (MQ=255)

                                                                                                                                                                                                                        ttGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCAGCCGACATTGACGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTACTTTGGCTGGGGTACGAGGATTAGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCCATACCCCAACTGGGTGCACTTACACGTTAAGCGTCTTGAATCTCTTGTCGGCAGTTCGAGGATTc             >  4:518465/1‑275 (MQ=38)

                                                                                                                                                                                                                              gATGCGTGCCAACCGTAGGTCAGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGACAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGGGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCt       <  2:414618/275‑1 (MQ=255)

                                                                                                                                                                                                                                tGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCg     >  2:364263/1‑275 (MQ=255)

                                                                                                                                                                                                                                    tGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATtgaa                                                                                                                                                                                                <  7:1063328/85‑1 (MQ=255)

                                                                                                                                                                                                                                    tGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATCCCACAACTGGGTGCACTTACACGGGAAGCGTTTTGAAGAAATTGGCTGGGGTACGTGGATTCGAAcgaaggaa  >  4:402687/1‑267 (MQ=39)

                                                                                                                                                                                                                                                 tCGGATAAGACGCGGCAGCGTCGCATCCGACATTGa                                                                                                                                                                                                                                    >  7:415335/1‑36 (MQ=255)

                                                                                                                                                                                                                                                 tCGGATAAGACGCGGCAGCGTCGCATCCGACATTGa                                                                                                                                                                                                                                    <  8:415335/36‑1 (MQ=255)

                                                                                                                                                                                                                                                 tCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACAt                                                                                                                                                                                                    <  5:865586/68‑1 (MQ=255)

                                                                                                                                                                                                                                                 tCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACAt                                                                                                                                                                                                    >  6:865586/1‑68 (MQ=255)

                                                                                                                                                                                                                                                   ggATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATtga                                                                                                                                                                                                 <  5:469427/69‑1 (MQ=255)

                                                                                                                                                                                                                                                   ggATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATtga                                                                                                                                                                                                 >  6:469427/1‑69 (MQ=255)

                                                                                                                                                                                                                                                     aTAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACAt                                                                                                                                                                                                                >  5:969356/1‑52 (MQ=255)

                                                                                                                                                                                                                                                     aTAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACAt                                                                                                                                                                                                                <  6:969356/52‑1 (MQ=255)

                                                                                                                                                                                                                                                      tgcgACGCTGCCGCGTCTTATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATtga                                                                                                                                                                                                 >  5:829794/4‑66 (MQ=25)

                                                                                                                                                                                                                                                      tgcgACGCTGCCGCGTCTTATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATtga                                                                                                                                                                                                 <  6:829794/63‑1 (MQ=25)

TGGGCCAGTTATCCTGTCCTCAAGCGTTACCAGATGCGGCGCATGGCGGATAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGT‑TA‑GGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCG‑AAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAA  >  NC_000913/696073‑696572

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 8 ≤` `ATCG/ATCG` `< 15 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`