BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,117,123	G→T	R34S (CGC→AGC)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,117,123	0	G	T	99.2%	428.2 / NA	129	R34S (CGC→AGC)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base G (0/0); major base T (59/69); minor base C (1/0); total (60/69)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 130 total aligned reads displayed.

GATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCA  >  NC_000913/4116848‑4117390

gATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCCGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGTTTTGCTCAAATTCGctct                                                                                                                                                                                                                                                                           >  2:281554/1‑278 (MQ=255)

    acagcttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTg                                                                                                                                                                                                                                                                          >  3:659043/5‑275 (MQ=255)

      aTCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGaca                                                                                                                                                                                                                                                                     <  4:59576/278‑1 (MQ=255)

      aTCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGaca                                                                                                                                                                                                                                                                     <  4:311109/278‑1 (MQ=255)

       tcttctAAACCGTCACGTTTTAAATGCTCGCCGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGa                                                                                                                                                                                                                                                                       <  1:205559/275‑1 (MQ=255)

          tctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGaca                                                                                                                                                                                                                                                                     >  1:503634/1‑274 (MQ=255)

          tctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCGTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCGTAAATTCTCTCTGTGacac                                                                                                                                                                                                                                                                    >  1:71821/1‑275 (MQ=255)

           agaaaCCGTCACGTTTTAAAGGCTCGCAGATTTCAGCGGTACGACGGGACTTCCAGGCAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTCTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTa                                                                                                                                                                                                                                                                <  2:227199/276‑1 (MQ=255)

            tAAACCGTCACGTTTTAAATGCTCGCAGATGTCTGCGGGACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTa                                                                                                                                                                                                                                                                <  2:366504/277‑1 (MQ=255)

             aaaCCGTGCGGTTTTAAATGCTCACAGATTTTTCCGGGACGCTGGCACGTCCAGAAGATGGGGTGTTATAGAGGTTTGCCGGTTTCGTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AAGTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAAt                                                                                                                                                                                                                                                              <  4:321557/278‑1 (MQ=255)

              aaacGTCACGGTTGAAGTTGTCGCAGATTTCAGCGATACGACGGCATTGCCAGGCAATGGCTTTATAGATAGGCTGTCCGGTTTCTATTTCCCAGACAATAGTGGTTTCAC‑GATAGTTCGTAATACCGATAGCTAC‑AATTTGATCGGAACTGATAACGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGACCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTa                                                                                                                                                                                                                                                                <  3:347591/272‑1 (MQ=255)

              aaCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCCCAAATTCGCTCTTCGCCACgtca                                                                                                                                                                                                                                                                >  3:663543/1‑272 (MQ=255)

                  gTCACGTTTTAAATGCTCGCAGATTTATGCGTGACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGaca                                                                                                                                                                                                                                                                     <  1:663543/266‑1 (MQ=255)

                         tttAAATGCGAGCAGATTTCCGCGGTACTACGGAACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTATTTTTCTCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAAGACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATc                                                                                                                                                                                                                                                  <  2:122828/278‑1 (MQ=255)

                             aaTGCTCGCAGATTTCTGCGGTCCGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGa                                                                                                                                                                                                                                              <  2:597202/278‑1 (MQ=255)

                                gttcgtaacttgatGCGGTACGACTGCCTTGCCAGACAATGGCTGTATCGATTGCTTTGCCGGTTTCTTTTTCCCCAACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGTTAGCTGC‑CATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGa                                                                                                                                                                                                                                              <  3:71821/262‑1 (MQ=255)

                                gCTCGCAGATTTCTGCGGTACGACCGCACGGCCAGAAAATGGCGTTATCGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTAAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATGTCGGCTTGCGCCAGCACTTCTACCGGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCTTGa                                                                                                                                                                                                                                              <  3:179760/275‑1 (MQ=255)

                                gCTCGCAGATTTCTGAGGTACGACGGCACTGCCAGACAATGGGGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGa                                                                                                                                                                                                                                              <  1:150717/275‑1 (MQ=255)

                                  tCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGATTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTTTTCTACCCATCCTGGTTTGGGGTAGATTTGCCCAAATTCGCTCTGCGACCACCTAATGATATTGGATTCGTg                                                                                                                                                                                                                                               >  2:482111/1‑272 (MQ=255)

                                   cGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATcc                                                                                                                                                                                                                                           >  3:229056/1‑275 (MQ=255)

                                     cAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTAc                                                                                                                                                                                                                                      <  4:293533/278‑1 (MQ=255)

                                        acttCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGAc                                                                                                                                                                                                                                   <  2:182460/276‑1 (MQ=255)

                                             tGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTTATCCATTACGACcgccc                                                                                                                                                                                                                              >  2:310862/1‑276 (MQ=255)

                                                  tACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTTCCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGCCACGCTCCTTATATTGGCATCGTTATCCATTCCGACCTCGCggca                                                                                                                                                                                                                          >  2:267527/1‑275 (MQ=255)

                                                    cggcgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGTGGCTTGTTCTACCCAACCTGGTTCTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGCGATCCATTACGACcgcgcg                                                                                                                                                                                                                             >  1:393916/4‑272 (MQ=255)

                                                    cgacgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGa                                                                                                                                                                                                                          >  2:48538/1‑275 (MQ=255)

                                                       cggaaatGCCAGACAGTGGTGTTGTTGGTAGGCTAGCCGGTTTCTGTTTCCAAAACAAGGGGGGATACAC‑GCTGGTTCGTAATACCAAGGGCTGC‑AATTTGATCGAAACTGATGACGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAAGCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGc                                                                                                                                                                                                                        <  3:177737/268‑1 (MQ=255)

                                                       cgGCAGTGCCAGACAATGGCGGTATAGATAGGCTTGGCGGTTTCTTTTTCGCGGACAATAGTGTTTTCAC‑GCTGGGTCGTAATACCGATAGCTGC‑AATTTGCTCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCTTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAg                                                                                                                                                                                                                         <  1:500405/273‑1 (MQ=255)

                                                       cgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTAGTGTTCTACCCAACCTGGGTTTGGGTAGATTTGCTACAATTCGCCTTGGGCCACGCTAATGATTTTGGCATCGTGATCCATTACGCACGCGggggcgat                                                                                                                                                                                                                       >  3:59849/1‑271 (MQ=255)

                                                            cTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCAAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGCCACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTggtggtt                                                                                                                                                                                                                >  2:174883/1‑276 (MQ=255)

                                                                      aTGGCGTTATAGATAGGCGTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATTCCGATAGCTGA‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTc                                                                                                                                                                                                        <  1:161027/275‑1 (MQ=255)

                                                                        agacatataGATAGGCGTGGCGGTTTCTTTTGCCCAGACAATAGGGGTTGCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGCTCGTAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGCGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAgcg                                                                                                                                                                                                   <  4:260680/273‑1 (MQ=255)

                                                                         gCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTGTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAg                                                                                                                                                                                                     <  1:114585/275‑1 (MQ=255)

                                                                           gtttgagcttggCTTGCCTGTTTTTTTTTCCCTGACAATAGTGTTTTCAC‑CCTGGTTAGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCaa                                                                                                                                                                                                <  4:281554/268‑1 (MQ=255)

                                                                           gTTATTGTTCGGCTTGCCGGTTTCTTTTTCCCATGTAATAGTGGTTTCAC‑GCTTGTTCGTAATACCGATAGCTGC‑AATTGTATCGGAACTGCTATCGGCTTTCGCCAGCACTTCTACCAGCTTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCaa                                                                                                                                                                                                <  4:18237/278‑1 (MQ=255)

                                                                           gTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAgcg                                                                                                                                                                                                   <  1:229056/275‑1 (MQ=255)

                                                                            ttATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCATCCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAgcgc                                                                                                                                                                                                  <  1:659043/275‑1 (MQ=255)

                                                                             tataGATAGGGTTGCCGGTGTCTTTGTCCCCGACATGAGTGGTTTCAC‑GCTTGTTAGTAATACCGGTAGCGGC‑AATTTGGTCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGGGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGGCGTGTTCTACCCAACCTGGTGTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCa                                                                                                                                                                                                 <  3:112587/275‑1 (MQ=255)

                                                                               tagatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCCCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGTTGGTGCCCTGGTCGAGCGCaag                                                                                                                                                                                               >  1:672823/1‑274 (MQ=255)

                                                                                      gCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTCGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATAGTTATCCCTTACGACCGCACGGGAGCTGGTGGTGACCGGGTCGGGCGCACAGCTATAtttt                                                                                                                                                                                     >  2:356085/1‑278 (MQ=255)

                                                                                         tGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACCCGCTCATGATATTGGCTTCGTGATCCATTACGACCGCGCGGGAGCTGGTTGTGCCCTGGTCGAGCGCAACGctatctctgtc                                                                                                                                                                                   >  2:260680/1‑270 (MQ=255)

                                                                                         tGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACATGCTAATTATATTGGCATCGTGATCCATTACGACCCCGCTGGAGCTGGTGGTGCCCTGGTCTGGCGCAACGCTAtttttt                                                                                                                                                                                     >  1:216409/1‑275 (MQ=255)

                                                                                         tGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACACAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGACCTGGTGGTGCCCTGGTCGAGGGCAAAGATATAtttt                                                                                                                                                                                     >  3:478744/1‑275 (MQ=255)

                                                                                           ccGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTGGATTTGCTTAAATTCGATCTGCGACACGCTAATGATATTGGCATAGTGATCCTTTACGACCGCGCGGGAGCTGTGGGTGCCCTGGTCGAGCGCAACGATCTAttatttttt                                                                                                                                                                                >  2:376974/1‑277 (MQ=255)

                                                                                           ccGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCACTGGTCGAGGCCAACGATCTATTTTTTTTc                                                                                                                                                                                >  4:25424/1‑278 (MQ=255)

                                                                                                     tttCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCAt                                                                                                                                                                          <  3:503634/274‑1 (MQ=255)

                                                                                                           aGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAAt                                                                                                                                                                   <  1:172750/275‑1 (MQ=255)

                                                                                                                 gtttGGTTGCACAGCTGGTTCGTAATACCCATAGCTGC‑AAGTTGTTCTGAACTGATATCGGCTTTCGCCAGCACTGCTACCAGCGTGGTGGTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGt                                                                                                                                                           <  4:496256/275‑1 (MQ=255)

                                                                                                                   gTGGTTTCCG‑GCTGGTTCG‑AATAACGGTTTGTGTAAATTTGGTGGGAACTGTAATCGGTTTTCGCCAGCACTTTTACCAGGGGGGAGCTTTGGGTGGCCCAGGTTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGt                                                                                                                                                           <  3:669624/275‑1 (MQ=255)

                                                                                                                   atGGTTTCAC‑GCTGGTTCGAAATACCGATCGCTGC‑AATTTGTTCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGt                                                                                                                                                           <  1:594770/274‑1 (MQ=255)

                                                                                                                       tttCAC‑GCTGGTTCATAATACCGATAGCTCA‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTTGGTGGCCCAAATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGAGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCGTCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTc                                                                                                                                                       <  3:393916/275‑1 (MQ=255)

                                                                                                                        ttCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGatat                                                                                                                                                                                                                                                         <  2:506594/176‑1 (MQ=255)

                                                                                                                         tCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGGCACGTAGTCATAtt                                                                                                                                                  >  2:464547/1‑278 (MQ=255)

                                                                                                                           aC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGACTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAg                                                                                                                                                                                                     >  3:60238/1‑225 (MQ=255)

                                                                                                                                  gTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGGGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCGCGTAGTCATATTTCAACg                                                                                                                                            >  2:7720/1‑276 (MQ=255)

                                                                                                                                   ttCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGTTATCGGCTTTCGCCAGCCCTTCTACCGGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCg                                                                                                                                            <  3:336617/275‑1 (MQ=255)

                                                                                                                                      gTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGaa                                                                                                                                          <  1:411489/274‑1 (MQ=255)

                                                                                                                                       tAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCttt                                                                                                                                     <  2:96814/278‑1 (MQ=255)

                                                                                                                                          tACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCtttt                                                                                                                                    >  1:504798/1‑275 (MQ=255)

                                                                                                                                                                 ggcagtGATATCGGCTTTCGCCAGGAGTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATGTCCATTGGGTCGTGTTCTACCCAGCCTTGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGGGGGGGCTGGTGGTGCCCTAGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTcc                                                                                                            <  2:540939/273‑1 (MQ=255)

                                                                                                                                                                  gAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTTCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTCCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTATTTGTCCCGTGGCCAAATTACAGCGAAGCTTTGTGTTCTGAAGGAGTTGTGGttt                                                                                                              >  1:657061/1‑275 (MQ=255)

                                                                                                                                                                     ctgatcttcgCTTTCGGCAGCACTTCGACCAGCGTGGAGCTTGGGGTGGCCCAGATTTCCATTGGGTTGGGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATGCGCTCGGCGACACGCTAATGATATTGGCATCGTGAT‑CATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTttc                                                                                                       <  4:93663/269‑1 (MQ=255)

                                                                                                                                                                             ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCACATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTATAGCGACGCTTTTTTTTCTGAAGGAGTTGTGGTTTCATTTTCTTcct                                                                                                  >  4:459955/1‑275 (MQ=255)

                                                                                                                                                                             ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCAtt                                                                                                                                                                                                                                        >  1:565679/1‑142 (MQ=255)

                                                                                                                                                                             ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTTTTCTGAAGGAGTTGTGGTTGCCTTTTCTTCcacc                                                                                                >  2:522303/1‑277 (MQ=255)

                                                                                                                                                                             ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGTGTTGTGGTTTTCTTTTCTTcccc                                                                                                 >  2:398653/1‑275 (MQ=255)

                                                                                                                                                                               cTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCcac                                                                                                 <  1:59849/275‑1 (MQ=255)

                                                                                                                                                                               cTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCacaac                                                                                              <  2:431587/278‑1 (MQ=255)

                                                                                                                                                                               cTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCacaa                                                                                               <  2:459955/277‑1 (MQ=255)

                                                                                                                                                                                  tCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCCTATTAAAGCGAATCTTTTTTTTTTTATGGGGTTGTGGTTTCTTTTTTTTCAACAacact                                                                                           >  4:77247/1‑277 (MQ=255)

                                                                                                                                                                                   cGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCacaaca                                                                                             <  1:478744/275‑1 (MQ=255)

                                                                                                                                                                                          attgttACCGGCGTTGTGGTTTGGTTGGCGCAGCGTCTCATTTGGTCGTGTTCTATCAATGCTGGGTTTGGGTTGATGTGCTCAAATTCTCTCTTCGACACGCTAATGCGTTTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTGTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGc                                                                                   <  4:472949/273‑1 (MQ=255)

                                                                                                                                                                                          aCTTCTTCCTGCGTTGGGCCTTGGGTGGCCCAGGTTTCCATTGGGTCGGTTTTTACCCGACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGAAACGATAATGAGATTGGCATCGTGATCCATTACGGCCGCGCGTGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGc                                                                                   <  4:464547/278‑1 (MQ=255)

                                                                                                                                                                                           cTTCTACAAGCTTGGGGCTTTGGGTGGGCCAGATTCCCTTTGGGTCGTGTTCTACCGAACCTGGTTTTGGGTAGAGTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGc                                                                                   <  2:463149/277‑1 (MQ=255)

                                                                                                                                                                                                   aGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCTTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGGAGTTATATTACAGCGAAGCTTATTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCTAAACACATATATAGCAAGGCAAAt                                                                          >  4:116820/1‑278 (MQ=255)

                                                                                                                                                                                                      gTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACcgc                                                                                                                                                                                                                                >  4:613623/1‑125 (MQ=255)

                                                                                                                                                                                                      gTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACATTTCTCCCAAGGCAAATgtg                                                                       >  4:310069/1‑276 (MQ=255)

                                                                                                                                                                                                      gTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAAt                                                                          >  3:346626/1‑275 (MQ=255)

                                                                                                                                                                                                       tGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACcgcg                                                                                                                                                                                                                               <  2:613623/125‑1 (MQ=255)

                                                                                                                                                                                                         gggCTTTGGGTGGCCCATGGATGAAATGAGTCGTGTGTGACCCGAGCTAGTGTTGGGTAGATTGGATCAAATTCGCTGTGCGATGAGCTAATGATATTGGCATCGTGGTCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGGACTACGGTATATTTTTTTTCAGTCATGATTACTTGTCCCGTAGTCCTATTACAGCGAAGCTTTTTGTTCTGAAAGAGTGGTGGTTTCCTTGTCTTCCACAACACAGATATCGAAAGGCAAATgg                                                                        <  3:530283/272‑1 (MQ=255)

                                                                                                                                                                                                                  gTGGCTCTGATTTCCATTTGGTCGTTTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGTTCTGCGACACTCTAATGATATTTGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCTTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGt                                                            <  4:29790/277‑1 (MQ=255)

                                                                                                                                                                                                                    ggCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCCTGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTTCTCCACACCCAAGATATCCCAAGGCAAATGGCTGCCAATCCGt                                                            >  3:189401/1‑275 (MQ=255)

                                                                                                                                                                                                                       ccattgttCCATTTGTGCGTGTTCTACCAAACCTTTTTTTTGGTATCTTTGCTCAAATTTGCTCGGTGACACGCTAATGCTATTGGCATCGTGATCCATTACGGACGCGCGGGAGCTGGTAGTGCCCTGGTCGGGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGACGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATTAGTTTGCg                                                       <  4:522303/271‑1 (MQ=255)

                                                                                                                                                                                                                       ccaagtttCCGTTTGGTCGTTTTCTACCAAACCTGGTTTTGGGTAGGTTTGCTCACATTCGCTCTGCGAAACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTAGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCgg                                                      <  4:398653/273‑1 (MQ=255)

                                                                                                                                                                                                                              tCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCa                                                <  2:77247/277‑1 (MQ=255)

                                                                                                                                                                                                                                     ggTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGc                                           >  3:251418/1‑275 (MQ=255)

                                                                                                                                                                                                                                     ggTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCAGTTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGACGGCGTTGTTGGTTCCTTTTATTCCACACCACAGATATCACAAGGCAAATGGCGCCAAAACAGTTTCCGGAAGGCAAATGCAc                                         >  4:514777/1‑277 (MQ=255)

                                                                                                                                                                                                                                           gtTCTAGCCACGCTGGTTTTGTGTTGTGTTTCCTAATTTCCCTCTTGGACACGGTAAGGTTATTGTCATCGTGATCCATTATTACCGCGCTGGGGCTGGTGGTGCCCTGGTCGAGCTCAACTATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAAt                                  <  4:174883/278‑1 (MQ=255)

                                                                                                                                                                                                                                                      ccTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTTAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCCCGCCACCGa                       >  2:406114/1‑278 (MQ=255)

                                                                                                                                                                                                                                                      ccTGGTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGACGGAGTTGTGGTTTCCTTTTCTTCCACAACTATGATATCGCAAGTCAAATGGACACCAATCAGTTTGCGTTAGGCACATACACCTACAATCGCGCCAAcaa                       >  2:405787/1‑276 (MQ=255)

                                                                                                                                                                                                                                                          gTTTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCg                                                       <  3:234864/242‑1 (MQ=255)

                                                                                                                                                                                                                                                           ttttGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGAt                      >  1:372876/1‑274 (MQ=255)

                                                                                                                                                                                                                                                           ttttGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCCCAAGTCAAATGGCGACCAATCATTTTTCGGTAGGCAAATGCACCTACAACCGCCCCAACCATa                     >  3:570674/1‑275 (MQ=255)

                                                                                                                                                                                                                                                                ggTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGATGCTTTTTGTTCTGAAGGAGTTGTGGTTTACTTTTCTTCCAAAACACAGATATCGCAAGGCAAATGGAGACCAATCAGTTTGCGATATGCAAATGCACCTACAATCGCGCCAACGATAGGACCGaa             >  2:436360/1‑278 (MQ=255)

                                                                                                                                                                                                                                                                ggTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGAAAGGGcc                >  1:680049/1‑275 (MQ=255)

                                                                                                                                                                                                                                                                    gATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCATCGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGc         >  2:229631/1‑278 (MQ=255)

                                                                                                                                                                                                                                                                    gATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACCACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCtcaaa                                   >  2:176851/1‑248 (MQ=255)

                                                                                                                                                                                                                                                                          ttCAAATTCTCTCGTCGACCCGCTAATGTTATTGGTCTCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCAcc   <  2:16252/277‑1 (MQ=255)

                                                                                                                                                                                                                                                                          cTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTTCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGACGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATTGCAACCAATCAGTTTGCGGTAGGCAAATGCAACTACAAACACCCCAACGAAAGGGCCGAAAAGCGGc      >  3:582534/1‑275 (MQ=255)

                                                                                                                                                                                                                                                                               aTTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGACAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTAAAATCGCGCCAACCATCGGGCCGAAAAGGCGCAcac  >  3:174063/1‑272 (MQ=255)

GATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAC‑GCTGGTTCGTAATACCGATAGCTGC‑AATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCA  >  NC_000913/4116848‑4117390

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