BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,956,003	G→A	D225N (GAT→AAT)	ilvA →	threonine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,956,003	0	G	A	95.5%	444.0 / NA	136	D225N (GAT→AAT)	ilvA	threonine deaminase
Reads supporting (aligned to +/- strand): ref base G (4/0); major base A (47/79); minor base C (6/0); total (57/79)
Fisher's exact test for biased strand distribution p-value = 3.51e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 136 total aligned reads displayed.

CGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTCGCG  >  NC_000913/3955755‑3956249

cGATCGAACTGTCACAGCAGCAGGGGTTCACCTAGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACCACTGATGCCGCCACTCAAAGTGGTCGCCGCATAAGCGGCAG‑ACTCCCCCACCCTGAAACCAGCGCTGGac                                                                                                                                                                                                                                                        >  2:101296/1‑250 (MQ=255)

 gATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATg                                                                                                                                                                                                                                                       <  1:560044/251‑1 (MQ=255)

  aTCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCCAATCAAAGTAATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGGGCTGGATGc                                                                                                                                                                                                                                                      >  1:154942/1‑251 (MQ=255)

    cGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGGTCAAACAACTGATGCCGCAAATCAAAATGATCGCCGTAGAAGCGGAAG‑ACTCCGCCGGCCTGAAAGCAGCGCCGAATGCgg                                                                                                                                                                                                                                                    >  1:120769/1‑251 (MQ=255)

         tGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CAt                                                                                                                                                                                                                                              >  1:295767/1‑251 (MQ=255)

          gTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CAt                                                                                                                                                                                                                                              <  1:161510/250‑1 (MQ=255)

             agagcagaatggtttccccTGGGTGCCGCCGGTCGACCATCCGATGGGGATTGCCGGGCAAGGCACGCTGGCGTTTGAACAGCGCTAGCAGGACGCTCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCAGTCGGGCTGCTGGCGTGGCGGTTCGGATCAAACAACTGATGCCGCAAATTAAAGTGATCTCCGGAGAAGCGGAAG‑ACGCCGCATGCGTGAACGCAGCGCTTAATGCGGGT‑CATCCgg                                                                                                                                                                                                                                          <  2:652951/234‑1 (MQ=255)

             acaGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGCTGGCGCTTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGCTCAAACAACTGCTGCCGCAAATCAAAATGATCGCCGTAGCAGCGGGAG‑ACTACGCCTGCCTGAAAGCAGCGCTGACAGCGGGT‑CATCCgg                                                                                                                                                                                                                                          >  2:253591/1‑251 (MQ=255)

             acaGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGGGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑AATCCCCCTGCCTGAAAGCAGCGCTGCATGCGGGT‑CATCCgg                                                                                                                                                                                                                                          >  1:551121/1‑251 (MQ=255)

                  agctggggTTCTCCGGGGTGCCGCCGTTCGGCCATCCGATGGTGATTCCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATGTCGACCGCTTAGTTGGGCCAGTCGGGGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAATTCAAAGTGATTGCCGTAGAAGCGGAAG‑TCTCCGCCTGCCTGAAAGCAGCGCTGAAGGCGGGT‑CAGCCTGTTGAt                                                                                                                                                                                                                                     <  2:308524/247‑1 (MQ=255)

                  agcagGGGTTCCCCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGAt                                                                                                                                                                                                                                     <  2:250219/251‑1 (MQ=255)

                     agGGGTGCACCTGGGTGCCGCCGTTCGACCAGCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCGCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTg                                                                                                                                                                                                                                  <  2:406054/251‑1 (MQ=255)

                          ttCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCcgcg                                                                                                                                                                                                                             >  1:382130/1‑251 (MQ=255)

                            cACCTGGGTTCCGCCGTTCGACCATCCGAGGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCcgcgcg                                                                                                                                                                                                                           <  2:620344/251‑1 (MQ=255)

                                 gggTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAggg                                                                                                                                                                                                                      <  1:642640/251‑1 (MQ=255)

                                  ggTGACGCCGTTAGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGGGGCGGTGCTGATCAAACAACTGATGCCGCAAATCCAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCCGCCTGAAAGCAGCGCCGAATGCGGGT‑CATCCGGGTGGACTGCcgcgcggggggc                                                                                                                                                                                                                     >  1:523459/1‑245 (MQ=255)

                                    ggccgccGTTCGGCCAGCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCCCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGGCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGcga                                                                                                                                                                                                                   <  2:294920/250‑3 (MQ=255)

                                      ccgccgTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCGGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑Att                                                                                                                                                                                                                <  2:451608/251‑1 (MQ=255)

                                       cgccgTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGTAGGACGCACATATCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAAAAACTGATGCCCCAAATCACAGTGATCGCCGGAGAAACGGAAG‑ACTCCGCCCGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGCTCTGCCGCGCGAAGGGCT‑Attt                                                                                                                                                                                                               >  1:460626/1‑251 (MQ=255)

                                       cgccgTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑Attt                                                                                                                                                                                                               <  1:518720/251‑1 (MQ=255)

                                         ccgTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CCTCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGc                                                                                                                                                                                                             >  1:722169/1‑251 (MQ=255)

                                          cgTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGGGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGGAAGCAGCGCGGAAAGCGGGG‑CCACCGGTTGATCTGCCGCGCGGGGGGAT‑ATTTGCt                                                                                                                                                                                                            >  1:216885/1‑251 (MQ=255)

                                                          ggTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACGTCGCCAGCAGGACGCCCGTCTCGACCGCGTATTTGTGCCAGTGGTCGGCGGTGATCTGGCGGCTGGCGTGGCTGTGCGGATCAAACAACGGATGCCGCAAATCAAAGTGATCTCCGTATAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGTGCTGAATGCGGGT‑CAGCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGAGTAGGGGTaa                                                                                                                                                                                            <  2:700973/251‑1 (MQ=255)

                                                             gATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGTACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGGGCTGGACAAACACCTGATGTCGCAAATAAAAGGGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAACAGCGCTGAATGCGGGG‑CCTCCGGGTGATCTGCCGCGCGTAGGGGG‑ATTTTCTGGAGGCGTAGCGGTaaaa                                                                                                                                                                                          >  1:731757/1‑250 (MQ=255)

                                                                  ccGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATc                                                                                                                                                                                    <  1:512490/251‑1 (MQ=255)

                                                                  ccGGGCAAGGCACGCTGGCGCGGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATc                                                                                                                                                                                    <  2:80987/251‑1 (MQ=255)

                                                                    gggCAAAGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTTCTGATCAAACAACCGATGCCCCCAATCAAAGTGATCGCCGGAGAAGCGGAAG‑AATCCGCCTGCCCGCAAGCAGCGCCGGATGCCGGT‑CATCCGGTTGGGCTGCCGCGCGTAGGGCC‑ATTTGCTGCAgggggtgcgggaaaaacaatcgg                                                                                                                                                                                  >  2:307442/1‑230 (MQ=255)

                                                                         aGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACg                                                                                                                                                                             <  1:617250/251‑1 (MQ=255)

                                                                         aGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACg                                                                                                                                                                             <  1:238739/251‑1 (MQ=255)

                                                                           gCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAAATGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCCCCCGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGa                                                                                                                                                                            >  1:546121/1‑250 (MQ=255)

                                                                            cACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTTCTGATCAAACAACTTATGCCGCAAATCAAAGGGGTCGCCGTAGAAGCGGAAG‑ACTCCCCCCGCCTTAAAGCAGCGCTGCACGCCGGG‑CATCCGGGTGTACTGGCGCGCGCTTGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACACATCGGTGACGaac                                                                                                                                                                          >  2:17921/1‑250 (MQ=255)

                                                                             aCGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGGTGATCAAACCACTGCGGCCGCAAATCCAAGATAACGCCGTAGAACCGGAAGCAC‑CCGCCTGCCTGAAAGCAGCGGTGAAAGCGGGTGCACACGGT‑GATTTGGCGCGCGTTGGGGT‑ATTTTCTGGAGGCGTAGCGGGTAAACGCCTCggggaagtac                                                                                                                                                                          >  2:28570/1‑242 (MQ=255)

                                                                                 tggcgctggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTTCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTc                                                                                                                                                                     >  1:664450/1‑251 (MQ=255)

                                                                                  ggcgctggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCtt                                                                                                                                                                      <  2:86311/249‑1 (MQ=255)

                                                                                  ggcgctggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTcc                                                                                                                                                                    <  1:427206/251‑1 (MQ=255)

                                                                                   gcgctggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCg                                                                                                                                                                   <  2:524393/251‑1 (MQ=255)

                                                                                    cgctggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCACACCACTGATGCCGCCAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCCGCCTGGAAGCAGCGCTGCATGCGGGT‑CCTCCGGTTGATCTTCCGCGCGGAGGGCC‑ATTTGCTGTAGGCGGAGCGGTAAAACGCATCGGTGACGAAACCTTCCGt                                                                                                                                                                  >  1:144950/1‑251 (MQ=255)

                                                                                    cgctggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCAGCGGTGACGAAACCTTCCGt                                                                                                                                                                  <  1:252902/251‑1 (MQ=255)

                                                                                     gctggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGGGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGtt                                                                                                                                                                 <  2:352216/251‑1 (MQ=255)

                                                                                        ggaccTGCGCCAGCAGGGCGCCCAGCGCGGCCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGAGCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTAt                                                                                                                                                              <  2:303220/247‑1 (MQ=255)

                                                                                        ggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGttt                                                                                                                                                                <  2:123738/249‑1 (MQ=255)

                                                                                          aaCTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGc                                                                                                                                                            <  1:468240/251‑1 (MQ=255)

                                                                                          aaCTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGc                                                                                                                                                            <  1:84569/251‑1 (MQ=255)

                                                                                             tGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTTCTGAAGGCGTAGCGGTAAAACGCCTCGGTGACGAAACCTTCCGTTTATGcc                                                                                                                                                           >  1:518269/1‑249 (MQ=255)

                                                                                              gCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCa                                                                                                                                                          <  1:197861/249‑1 (MQ=255)

                                                                                                    gcagGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTAtc                                                                                                                                                  >  1:247111/1‑251 (MQ=255)

                                                                                                      agGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTAtctc                                                                                                                                                <  2:363150/251‑1 (MQ=255)

                                                                                                         aCGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCTGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTcgac                                                                                                                                             <  1:287139/251‑1 (MQ=255)

                                                                                                            cccATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCgacgac                                                                                                                                          <  1:506088/251‑1 (MQ=255)

                                                                                                              cATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGG‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGGAAAACGCATCGGGGACGAAACCTTCCGTTTATGCCCGGGGGTTCTCgacgacct                                                                                                                                        >  1:647332/1‑249 (MQ=255)

                                                                                                                   cGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTTGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGGTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGGAGCGGTAAAACGCATCGGGGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACctcatcc                                                                                                                                   >  1:88599/1‑250 (MQ=255)

                                                                                                                   cGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGAcatca                                                                                                                                      <  1:375498/248‑1 (MQ=255)

                                                                                                                   cGACCGCGTATGTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACAtcatca                                                                                                                                   <  1:41699/251‑1 (MQ=255)

                                                                                                                       cgcgTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGt                                                                                                                               <  1:669845/251‑1 (MQ=255)

                                                                                                                        gcgTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑TCTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTc                                                                                                                              <  2:120769/251‑1 (MQ=255)

                                                                                                                                gtgCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgcg                                                                                                                       <  2:188276/250‑1 (MQ=255)

                                                                                                                                         ggcggcggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAACCAACTGATGCCGCAAATCAAAGTGGTCGCCGTAGAAGCGGGAA‑ACTCCGCCTGCCTGAAAGCAGCGCTGAAAGCGGGT‑CATTCGGTTGGTCTGCCGCGCGTCGGGCC‑ATTTGCTGGCGGCGCAGCGGGAAAACCCATCGGTGACGAAACCTTCCGTTTATGCCCGGAGGATCCCGACGACATCATCACCGTCGATAGCGATGCGGTCtg                                                                                                             >  1:197232/1‑251 (MQ=255)

                                                                                                                                         ggcggcggcggTCTGGCTGCTGGCGGGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCAGCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACAGCATCACCGTCGATAGCGATGCGATCtg                                                                                                             <  2:382130/251‑1 (MQ=255)

                                                                                                                                             gcggcggTCTGGCTGCTGGCGTGGCGGTGCTGGTCAAACAACTGATGGCCCAAATCCAAGGGGTCGCCGGAAAAGCGGAAA‑ACTTCGCCCGCCCGGAAACAACGCTGCAATCGGGT‑CAACCGGTTGGTCTGCCGCGCGTAGGGCTCTTTTGCT‑AAGGCGGAGCGGGAAAACGCCAAGGGGAGGAAAACTTCCTGTTATGCCAGGAGGATCTCCACGGCATCATCACCGTCGaaaacaatggggtctttgcg                                                                                                         >  2:558237/1‑232 (MQ=255)

                                                                                                                                               ggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTgcggc                                                                                                       <  2:197232/251‑1 (MQ=255)

                                                                                                                                                           gctgGCGTGGCGGTGCTGATCAAACAACTGATTCCGCAAATCAAAGTGATCTCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGTGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCAGCGGTGACGAAACCTTCCGGTTATGCCATGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGAt                                                                                             <  2:715212/249‑1 (MQ=255)

                                                                                                                                                           gctgGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGAtt                                                                                            <  1:656285/250‑1 (MQ=255)

                                                                                                                                                             tggcgtggcgGTGCTGCTCACACCACTGATTCCGCAAAACAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGGAGGGCT‑ATTTGCTGACGGCGTAGCGGGAAAAAGCACCGGGGGCGCAACCTTCCGTTTATGCCAGGAGTATCCCGGCGGCATCATCACCGGCGATAGCGATGCGCTCTGTGCGGGGAGGAAGGGTTTAt                                                                                         >  2:471116/1‑251 (MQ=255)

                                                                                                                                                             tggcgtggcgGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGGTCGCCGTAGAAGCGGAAG‑ACTCCGACTGCCTGAAAGCAGCGCTGGATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGCAGAAACCTTCCGTTTATGCCAGGAGGATCTCGACGCCATCATCACCGTCGGTAGAGATGCGATATGTGCGGCGATGAAAGATTTAt                                                                                         >  1:267799/1‑251 (MQ=255)

                                                                                                                                                              ggcgtggcgGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTAtt                                                                                        <  1:456553/251‑1 (MQ=255)

                                                                                                                                                                   ggcgGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGa                                                                                     >  1:173301/1‑249 (MQ=255)

                                                                                                                                                                    gcgGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGa                                                                                  <  1:414396/251‑1 (MQ=255)

                                                                                                                                                                            gATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACGCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTgcgcg                                                                          <  2:295316/251‑1 (MQ=255)

                                                                                                                                                                                aaaCAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCggg                                                                      >  2:335944/1‑250 (MQ=255)

                                                                                                                                                                                               cAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTgg                                                       >  2:176303/1‑251 (MQ=255)

                                                                                                                                                                                               cAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTgg                                                       <  2:518269/251‑1 (MQ=255)

                                                                                                                                                                                                   ccAAAGTGAGCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGgcgcg                                                   <  2:37115/250‑1 (MQ=255)

                                                                                                                                                                                                    cAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGgcgcgc                                                  >  2:703891/1‑251 (MQ=255)

                                                                                                                                                                                                    cAAAGTGATCGCCGTAGAAGCGGAAG‑ACTACGCCTGCCTGAAAGCAGCGCAGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGAAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAAAGCATCGGTGACGACACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCCTCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGAGGCGGAACCCTCTGgcgcgc                                                  >  1:674072/1‑251 (MQ=255)

                                                                                                                                                                                                     aaaGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCgcgct                                                 <  2:19234/251‑1 (MQ=255)

                                                                                                                                                                                                        gTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCgcgctggc                                              <  1:97776/251‑1 (MQ=255)

                                                                                                                                                                                                        gTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTGTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCGCTGGCgcgctggc                                              <  1:101296/251‑1 (MQ=255)

                                                                                                                                                                                                             cGCCGTAGAAGCGGAAG‑CCTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCgcgctggcgctg                                          <  2:154942/250‑1 (MQ=255)

                                                                                                                                                                                                             cGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCgcgcgggcgctg                                          <  2:730725/250‑1 (MQ=255)

                                                                                                                                                                                                             cGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGcgctggcgctgg                                         <  2:648036/251‑1 (MQ=255)

                                                                                                                                                                                                               ccGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCgctggcgctggc                                        >  1:76916/1‑250 (MQ=255)

                                                                                                                                                                                                                cGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCgg                                      <  1:670205/251‑1 (MQ=255)

                                                                                                                                                                                                                     aaGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATg                                 >  2:233040/1‑251 (MQ=255)

                                                                                                                                                                                                                        cGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGCATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGaa                               >  2:533440/1‑250 (MQ=255)

                                                                                                                                                                                                                         ggAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGaaaa                             <  1:335944/251‑1 (MQ=255)

                                                                                                                                                                                                                         ggAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGaaa                              >  2:557928/1‑250 (MQ=255)

                                                                                                                                                                                                                          gAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGaaaa                             >  1:720553/1‑250 (MQ=255)

                                                                                                                                                                                                                             g‑ACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGACAAAat                          >  2:548192/1‑250 (MQ=255)

                                                                                                                                                                                                                                    gcctgcctGAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGc                   >  1:713879/1‑251 (MQ=255)

                                                                                                                                                                                                                                     cctgcctgAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGcc                  <  2:546121/251‑1 (MQ=255)

                                                                                                                                                                                                                                        gcctgAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGCAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGccc                 <  2:216885/249‑1 (MQ=255)

                                                                                                                                                                                                                                        gcctgAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCTCCGTCGATAGCGATGCGATCTGTGCGTCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCGGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCt                <  2:370183/250‑1 (MQ=255)

                                                                                                                                                                                                                                        gcctgAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTg               <  1:198947/251‑1 (MQ=255)

                                                                                                                                                                                                                                         cctgAAAGCAGCGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCGGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTGTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCGCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGTCCTGc              <  2:731757/251‑1 (MQ=255)

                                                                                                                                                                                                                                          ctgAAAGCAGCGCTGCATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGAGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGc              >  2:75927/1‑250 (MQ=255)

                                                                                                                                                                                                                                            aaaaGCAGCGCTGAATGCGGGT‑CATCAGGTTGATCTGCAGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACATTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGTAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGACCCCTCTGGCGCGCTGGCGCTGGCGGGAATGaacaaaaaaatcgaccttcga            >  2:503316/2‑235 (MQ=255)

                                                                                                                                                                                                                                             aaaGCAGCGCTGAATGGGGGT‑CATCCGGTTGTTCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACCTCATCACCGTCGATAGCGATGCGTTCTGTGCGGCGTTGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGGAAGGAAAAATTTATCGGCCTGcgca           >  2:43582/1‑247 (MQ=255)

                                                                                                                                                                                                                                               agcagcGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAAGCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCacaaca        <  2:551121/251‑1 (MQ=255)

                                                                                                                                                                                                                                                gcagcGCTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACAt       >  1:386648/1‑251 (MQ=255)

                                                                                                                                                                                                                                                     gcTGAATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTcgcg  <  1:7579/251‑1 (MQ=255)

CGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAG‑ACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGT‑CATCCGGTTGATCTGCCGCGCGTAGGGCT‑ATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTCGCG  >  NC_000913/3955755‑3956249

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