BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	23,962	C→T	Y524Y (TAC→TAT)	ileS →	isoleucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	23,962	0	C	T	98.8%	252.9 / ‑4.5	82	Y524Y (TAC→TAT)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (34/46); total (36/46)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

AGTGGATCCCGGACTGGGGCCAGGCGCGTATCGAGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGTGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGAT‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACA‑‑CATTG‑GATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAACACCGCGGCTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAGGGTAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCA  >  NC_000913/23715‑24194

aGTGGATCCCGGACTGGGGCCAGGCGTGTATCGAGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGGACCGATGTCACTGGTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGTAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCGGACCAGTAtgtg                                                                                                                                                                                                                                           <  1:325697/251‑1 (MQ=255)

  tGGATCCCGGACTGGGGCCAGGCGCGTATCGAGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGa                                                                                                                                                                                                                                          >  1:136474/1‑250 (MQ=255)

   ggATCCCGGGCTGGTGCCAGGCGCGTATCGAGTCGATGGTTGCTAACCGTCCTGTCTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGaaa                                                                                                                                                                                                                                        <  1:355536/251‑1 (MQ=255)

       cccGGACTGGGGCCAGGCGCGTATCGAGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGc                                                                                                                                                                                                                                    <  1:345252/251‑1 (MQ=255)

            aCTGGGGCCAGGCGCGTATCGAGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGat                                                                                                                                                                                                                               <  1:427020/251‑2 (MQ=255)

             cTGGGGCCAGGCGCGTATCGAGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGat                                                                                                                                                                                                                               >  2:408489/1‑249 (MQ=255)

                  gCCAGGCGCGTATCGAGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGatactctc                                                                                                                                                                                                                         <  1:349993/251‑8 (MQ=255)

                     aGGCGCGTATCGAGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGatactctcgac                                                                                                                                                                                                                      <  1:57863/251‑11 (MQ=255)

                     aGGCGCGTATCGAGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGatactctcgac                                                                                                                                                                                                                      >  1:396684/1‑241 (MQ=255)

                             aTCGAGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCACCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGatactctcgacgtatggtt                                                                                                                                                                                                              <  2:262159/251‑19 (MQ=255)

                             aTCGAGGCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGatactctcgacgtatggt                                                                                                                                                                                                               <  2:116818/250‑18 (MQ=255)

                                      aTGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCgg                                                                                                                                                                                                  >  1:253088/1‑251 (MQ=255)

                                        ggTTGCTAACCGTCCGGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTGGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCGCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑C‑TCGCGACGTATGGTTTGACTCCGGAt                                                                                                                                                                                                <  2:277396/251‑1 (MQ=255)

                                             cTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAACGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGT‑TATCGGCATCCAGGCGTGGTGGGATCTCGCTGCGAAAGAGATCCTCGGCGACCAAGCTGACCAGTAGGTGAAAGTGCCGGATA‑‑CTCTC‑GACGGATGGTTTGACGCCGGATCtccc                                                                                                                                                                                           >  2:138466/1‑248 (MQ=255)

                                               aaCCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTAccc                                                                                                                                                                                          <  1:102422/250‑1 (MQ=255)

                                               aaCCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTAcc                                                                                                                                                                                           <  1:67659/249‑1 (MQ=255)

                                                acattCCTGACTTGTGTAGTTCCCGTCAGCGAATAGGGGACGTATCGATGTCACTGATCGTGAACAAAGAAACGGAAGAGCTGCAACCGCGTCCCCACGAACTGATGGTAGAAGTGGCAAAACTCGTGGCAGTTAAC‑GGTACCCAGGCGTGGTGGGAGATCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTAAAAGTGCCGGATATTCAT‑‑‑GACGTATGGTTTGACTCCGTATCTACCCat                                                                                                                                                                                        <  2:2411/247‑2 (MQ=255)

                                                     ccTGGCTGGGGTATCGCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACtctt                                                                                                                                                                                    <  2:101891/250‑1 (MQ=255)

                                                     ccTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACtctt                                                                                                                                                                                    <  1:110932/250‑1 (MQ=255)

                                                      cTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGGGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACtcttct                                                                                                                                                                                  >  2:2811/1‑251 (MQ=255)

                                                       tGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTAAGGTTTGACTCCGGATCTACCCACTCTTCtg                                                                                                                                                                                 >  2:432619/1‑251 (MQ=255)

                                                          cTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACCGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGCAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACCCTTCtgttg                                                                                                                                                                              >  2:232415/1‑251 (MQ=255)

                                                                     cccGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTCCCCACTCTTCTGTTGTTGACGTGCGt                                                                                                                                                                   <  1:480933/251‑1 (MQ=255)

                                                                      ccGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTc                                                                                                                                                                  >  1:471822/1‑251 (MQ=255)

                                                                                       gcgtACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCAc                                                                                                                                                 >  1:384567/3‑251 (MQ=255)

                                                                                         gtACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTTACCAGTATGTGAAATTGCCGGATA‑‑CTTTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGc                                                                                                                                               <  1:317329/251‑1 (MQ=255)

                                                                                           aCCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGAT‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACA‑‑CATTG‑GATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGc                                                                                                                                               >  1:115482/1‑249 (MQ=255)

                                                                                              gATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGTTCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGACTTTGCCGGTCACGCAGCgg                                                                                                                                          >  2:132675/1‑251 (MQ=255)

                                                                                               aTGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGa                                                                                                                                         <  2:476720/251‑1 (MQ=255)

                                                                                                   cACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑t                                                                                                                                     <  1:31894/250‑1 (MQ=255)

                                                                                                   cACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑Tg                                                                                                                                    <  1:158470/251‑1 (MQ=255)

                                                                                                   cACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACGCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑Tg                                                                                                                                    <  2:47776/251‑1 (MQ=255)

                                                                                                        ttCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTCCC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTACCGCAGCggcctagtatcc                                                                                                                                >  2:437855/1‑241 (MQ=255)

                                                                                                              cacaAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGAAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGCGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGAATCCGGATCTACCCACTCTTCTGTTGTTGCCGTGGGTCCGGAATTTGCCGGTCACGAAACGGGAA‑TGGATCTGGAAgg                                                                                                                         >  2:234815/1‑251 (MQ=255)

                                                                                                               acaAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGt                                                                                                                        <  1:193222/251‑1 (MQ=255)

                                                                                                                     acacGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGAc                                                                                                                  <  2:215461/251‑1 (MQ=255)

                                                                                                                     acacGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGAc                                                                                                                  <  1:247324/251‑1 (MQ=255)

                                                                                                                      cacGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGAcc                                                                                                                 <  1:197507/251‑1 (MQ=255)

                                                                                                                         ggAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACC‑TGTATCTGGAAGGTTCTGCCCagc                                                                                                              >  2:359122/1‑249 (MQ=255)

                                                                                                                            agagCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCAcc                                                                                                           <  1:299602/251‑1 (MQ=255)

                                                                                                                             gagCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCGAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACcg                                                                                                          <  2:464702/251‑1 (MQ=255)

                                                                                                                             gagCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACcg                                                                                                          <  1:233390/251‑1 (MQ=255)

                                                                                                                               gatGCATCCGCGTACCCTCGAACTGTAGGAAGAAGTGGCAAAACGCGTTTAAGTTAAC‑TGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACcgtg                                                                                                        <  1:287138/249‑3 (MQ=255)

                                                                                                                                  gCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GAAGTATGGTTTGGCTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTATGGCCAGCACcgtggtt                                                                                                     >  2:86077/1‑246 (MQ=255)

                                                                                                                                   cATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTGGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACcgtggttg                                                                                                    <  2:384567/251‑7 (MQ=255)

                                                                                                                                    aTCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGAGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTACGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGCCA‑TGTATCTGGAAGGTTCTGACCAGCACcgtggttgg                                                                                                   >  1:224604/1‑244 (MQ=255)

                                                                                                                                      ccGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGtggttggt                                                                                                  <  1:205248/250‑1 (MQ=255)

                                                                                                                                             cccTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCt                                                                                          <  1:363694/251‑1 (MQ=255)

                                                                                                                                                 cgAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTccc                                                                                      >  2:333697/2‑251 (MQ=255)

                                                                                                                                                  gAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCt                                                                                     <  2:150542/251‑1 (MQ=255)

                                                                                                                                                   aaCTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCt                                                                                     >  2:456419/1‑250 (MQ=255)

                                                                                                                                                    acggATGGAAGAAGGGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGCAGCTGACCAGTATGTGGAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGGTCTACCCACGCTTCTGTTGTTGTCGCGCGTCCGGCATTTGCCGGGAACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCCTCCCTat                                                                                   >  2:161011/4‑250 (MQ=255)

                                                                                                                                                       gATGGAAGAAGTGGCAAAACGCGTTTGAGTTAAC‑GGCCCCCAGGCGTGGTTGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCCCCGTTGTTGGTTCATGTCTTCCCTAAt                                                                                  <  2:297164/249‑1 (MQ=255)

                                                                                                                                                        aTGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGACTTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGAt                                                                               >  2:74095/1‑251 (MQ=255)

                                                                                                                                                           gaagaaGTGGTACAACGCGTTGAAGTCAAC‑GGCATCAAGGCGTGATGAGATCTCGATGCGCAAGAGATCCTCTGCGATGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCGGTTGTTGACTTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCATCACCGTGGTTGGTTCATGTCGTCCCTaagtatctc                                                                            <  2:202036/251‑8 (MQ=255)

                                                                                                                                                             agaagTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCa                                                                          <  2:77394/251‑1 (MQ=255)

                                                                                                                                                               aagTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCAcc                                                                        <  1:198477/251‑1 (MQ=255)

                                                                                                                                                                  tGGCAAAACGCGTTGAAGTTAAC‑GGCATCCATGAGTGGTGGGATCTCGATGCAAAAGAGATCCTCGGCGACGAATCTGACAAGTATGTGAACGTGACGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGG‑CACTGTATCTGGACGGTTCTGCACAGCCCCGTGGTTGGTTCATGCCTTCCCTAATAATCTCCACcgcg                                                                     >  2:268165/1‑251 (MQ=255)

                                                                                                                                                                                       aacggCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCagc                                                                                                              >  1:78396/4‑188 (MQ=255)

                                                                                                                                                                                       aacggCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCagc                                                                                                              <  2:78396/187‑3 (MQ=255)

                                                                                                                                                                                       aacggCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGt                                                  >  1:346015/4‑250 (MQ=255)

                                                                                                                                                                                           gCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGt                                             >  1:125609/1‑251 (MQ=255)

                                                                                                                                                                                            cATCCAGGCGTGGTGGGAGCCCGATGCGAAAGATATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATC‑‑CTCTC‑GACGGATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTc                                            <  2:224604/251‑1 (MQ=255)

                                                                                                                                                                                               ccAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGTCGAAGCTGACCAGTATGTGAAAGTGCCTGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCGCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAg                                          <  2:40357/250‑1 (MQ=255)

                                                                                                                                                                                                 aGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGt                                        >  1:413836/1‑250 (MQ=255)

                                                                                                                                                                                                          gggATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCAc                              <  2:276684/251‑1 (MQ=255)

                                                                                                                                                                                                           ggATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACg                             <  1:311493/251‑1 (MQ=255)

                                                                                                                                                                                                                gcGATGCGAAAGAGGGCCTCGGCGGTGAAGCTGACCAGTATGTTAAAGTGCCGGATA‑‑CTCTC‑GTCGTATTGTGTGACTCCGGATCTTCCCACTCTTCGTTTGTGGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGGTCATGGCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCTGTATCGTCAGGTACTGACCCACGGCtt                         <  2:431823/249‑1 (MQ=255)

                                                                                                                                                                                                                 cGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCttt                        <  2:242327/250‑1 (MQ=255)

                                                                                                                                                                                                                   aTGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTAcc                     <  2:95405/251‑1 (MQ=255)

                                                                                                                                                                                                                      cGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTg                  >  1:125602/1‑251 (MQ=255)

                                                                                                                                                                                                                          agagATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATg              >  1:241828/1‑251 (MQ=255)

                                                                                                                                                                                                                          agagATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGAACTGACCCACGGCTTTACCGTGGATg              >  1:205447/1‑251 (MQ=255)

                                                                                                                                                                                                                            agaTCCTCGGCGACGCAGCTGACCAGTATGTGCAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCGACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGTCA‑TGTCTCTGGAAGGTTCTGCCCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCCTTACCGTGGATGGt            >  2:18457/1‑251 (MQ=255)

                                                                                                                                                                                                                            agaTCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGt            >  2:40061/1‑251 (MQ=255)

                                                                                                                                                                                                                            agaTCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGt            >  1:25444/1‑251 (MQ=255)

                                                                                                                                                                                                                             gaTCCTCGGCGACGATGCTGACCAGTATGTGAAATTGCCGGATA‑‑CTCTC‑GCCGTAGGGTTTGACTCCGGATCTACCCACTCTGCTGTTGTTGACGTGCGTCCGGAAATTGCCGGTTACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGGCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTATCGTGGATGGTc           <  2:282799/251‑1 (MQ=255)

                                                                                                                                                                                                                               tCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTc           >  2:425456/1‑249 (MQ=255)

                                                                                                                                                                                                                                     gCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGc   <  1:408489/251‑1 (MQ=255)

                                                                                                                                                                                                                                      agacgaAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGAGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATAGTCAGGTACGGACCCACGGCTTTACCGGGGATGGTCacgggcgc   >  2:5231/2‑243 (MQ=255)

                                                                                                                                                                                                                                         cgaAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCa  >  2:306390/1‑248 (MQ=255)

                                                                                                                                                                                                                                              cTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATg                                                                                             <  2:320078/152‑1 (MQ=255)

                                                                                                                                                                                                                                              cTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATg                                                                                             >  1:320078/1‑152 (MQ=255)

AGTGGATCCCGGACTGGGGCCAGGCGCGTATCGAGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGTGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGAT‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACA‑‑CATTG‑GATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACA‑TGTATCTGGAAGGTTCTGACCAACACCGCGGCTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAGGGTAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCA  >  NC_000913/23715‑24194

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