BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,423,232	C→A	5.8%	intergenic (+38/+191)	yhdZ → / ← rrfF	putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,423,232	0	C	A	5.8%	591.1 / 18.8	174	intergenic (+38/+191)	yhdZ/rrfF	putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand): ref base C (87/76); new base A (5/5); total (92/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.03e-01

Only 100 of 175 total aligned reads displayed.

AAAGAGGTGCTGGATACGATGATTGGGCTGGCGCAGTCGGGTATGACAATGTTGTGTGTAACACATGAGATGGGGTTTGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCAT  >  NC_000913/3422985‑3423504

                      ttGGGCTGGCGCAGTCGGGTATGACAATGTTGTGTGTAACACATGAGATGGGGTTTGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTGTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTc                                                                                                                                                                                                                                   >  2:525481/1‑274 (MQ=255)

                                                 tgttgtGTGTAACACATGAGATGGGGTTTGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCGCCACCCCCCGGGCCGTTTCtt                                                                                                                                                                                                       >  4:520817/1‑275 (MQ=255)

                                                  gttgtGTGTAACACATGAGATGGGGTTTGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCAACCTCCGGGCCGTTTCtt                                                                                                                                                                                                       >  2:337228/1‑274 (MQ=255)

                                                      tgtgcAGCACATGAGATGGGGTTGGCACGAACCGTCGCTGACCGGGTAATTTTGATTGATCGTGTGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTGTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCa                                                                                                                                                                                                   <  1:1432/274‑1 (MQ=255)

                                                             cacaTGAGATGGGGTTTGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATg                                                                                                                                                                                                                                                          >  2:514118/1‑212 (MQ=255)

                                                             cacaTGAGATGGGGTTTGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGAAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                                                                                       >  2:233429/1‑263 (MQ=255)

                                                                 tGAGATGGGGTTTGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGAAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCa                                                                                                                                                                                                   <  4:233429/263‑1 (MQ=255)

                                                                  gagaTGGGGTTTGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTATGGCAGTTTATGGCGGGGGGCCTTCCCCCCCCCCNNNCGGCCGTTTGTTTGCCACGttttcaatcg                                                                                                                                                                                      >  3:188685/1‑267 (MQ=255)

                                                                       ggggTTTGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTccc                                                                                                                                                                                >  1:87900/1‑275 (MQ=255)

                                                                                    actgtcgtggACCGGGTAAGTTCTATGGTGCGTGGTGAAATAAGGGAGCAATCCGGAATCGTTGCATTTTGTGCTCATCCTAAAACCGCGAGTACGAGTGCATTTTTATCGCAGGTAATCCATTAAGTTAATGTTAGTGCGAACATCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTc                                                                                                                                                                    <  4:361993/265‑1 (MQ=18)

                                                                                    aCCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCGCCACCCTCCGGGCCCTTGCTTCGCAACGTTCAAATCCCGCACCGGGCGGtttttcc                                                                                                                                                                    >  2:467745/1‑272 (MQ=255)

                                                                                              aCCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCataaactgccag                                                                                                                                                                                                              >  2:486372/1‑212 (MQ=255)

                                                                                                   gTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCGcc                                                                                                                                                                                                                          >  4:567127/1‑206 (MQ=255)

                                                                                                   gTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGcc                                                                                                                                                                                                                          <  2:567127/206‑1 (MQ=255)

                                                                                                        ggttAGGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAAGTTTTTGCGCATCCTAAAGCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCAc                                                                                                                                                <  2:377883/272‑1 (MQ=255)

                                                                                                          ttATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGc                                                                                                                                                                                                                                                                             <  8:835006/149‑1 (MQ=255)

                                                                                                             tGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCTTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCCGTTTATGGCGGGCGTCCTTCCCTCCACCCTCCGGGCCGTTGCTTTGCAACGTTCAAAG‑CCGCTCCCCGCGGGTTTGTCCTCCTCGCTCAAGTGTTCACCGaaca                                                                                                                                          >  4:593970/1‑272 (MQ=255)

                                                                                                             tGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCGAAAGGATGGCC‑tt                                                                                                                                                                                                                                                                         <  5:305452/149‑1 (MQ=255)

                                                                                                               gATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTc                                                                                                                                                                                                                                                                       <  8:418133/149‑1 (MQ=255)

                                                                                                                 tCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGc                                                                                                                                                                                                                                                                     >  8:126180/1‑149 (MQ=255)

                                                                                                                 tCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGc                                                                                                                                                                                                                                                                     <  6:763247/149‑1 (MQ=255)

                                                                                                                  cGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATAGCC‑TTTCGCt                                                                                                                                                                                                                                                                    <  8:311236/149‑1 (MQ=255)

                                                                                                                    tGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCTTTGCTTCCCAACGTTCACAT‑CCCCTCCCGGCCGCTTTGTCCTCCTCGGGAGAGTGTTCCCCGacacacaccc                                                                                                                                    >  4:377883/1‑271 (MQ=37)

                                                                                                                      aggAAATAGTGGAGCAAGCTGCACCTGATGAAGTTTTTGCGCTGCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGAt                                                                                                                                 <  1:65382/274‑1 (MQ=255)

                                                                                                                        gAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCGAAAGGATGGCC‑tt                                                                                                                                                                                                                                                                         <  6:938098/138‑1 (MQ=255)

                                                                                                                        gAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCGAAAGGATGGCC‑tt                                                                                                                                                                                                                                                                         >  5:938098/1‑138 (MQ=255)

                                                                                                                           aTAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCttgatttgat                                                                                                                                                                                                                                                           <  7:633928/149‑1 (MQ=255)

                                                                                                                           aTAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCttgatttgat                                                                                                                                                                                                                                                           <  8:362829/149‑1 (MQ=255)

                                                                                                                              gTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGGGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCGAAAGGATGGCC‑ttt                                                                                                                                                                                                                                                                        <  7:755456/133‑1 (MQ=255)

                                                                                                                              gTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCGAAAGGATGGCC‑ttt                                                                                                                                                                                                                                                                        >  8:755456/1‑133 (MQ=255)

                                                                                                                                   gCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGGGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCa                                                                                                                                                                                                   <  4:305626/197‑1 (MQ=255)

                                                                                                                                   gCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGc                                                                                                                                                                                                                                                                             >  7:106151/1‑124 (MQ=255)

                                                                                                                                    cAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTACTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTTAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACCGATAACACAAAAGAcca                                                                                                                   >  2:145361/1‑274 (MQ=255)

                                                                                                                                    cAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAg                                                                                                                                                                                                                                                  >  6:534958/1‑149 (MQ=255)

                                                                                                                                       gCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCGCCACCCTCCGGGCCGTTTCTTCGCAACGTTCAAAT‑CCTCTCCCTGCGGACTTTTCCTACCAGGGAGAGTGTTCACCACCAAACCACAGATAACAAAAAAGGCCCaa                                                                                                                 >  4:544337/1‑273 (MQ=255)

                                                                                                                                          gCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATg                                                                                                                                                                                                                                            <  7:908095/149‑1 (MQ=255)

                                                                                                                                            aCCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTGGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATggc                                                                                                                                                                                                                                          <  6:663291/149‑1 (MQ=255)

                                                                                                                                              cTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATggcgg                                                                                                                                                                                                                                        >  6:418727/1‑149 (MQ=255)

                                                                                                                                              cTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATggcgg                                                                                                                                                                                                                                        >  7:200978/1‑149 (MQ=255)

                                                                                                                                              cTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATggcgg                                                                                                                                                                                                                                        <  6:1003534/149‑1 (MQ=255)

                                                                                                                                              cTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATggcgg                                                                                                                                                                                                                                        >  5:468407/1‑149 (MQ=255)

                                                                                                                                                gatgaATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATggcgggc                                                                                                                                                                                                                                      >  7:999221/1‑149 (MQ=255)

                                                                                                                                                gatgaATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATggcgggc                                                                                                                                                                                                                                      >  8:591272/1‑149 (MQ=255)

                                                                                                                                                    aaTTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGc                                                                                                                                                                                                                                                                             <  5:132601/107‑1 (MQ=255)

                                                                                                                                                                   cTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGt                                                                                                                                                                                                                                                         <  8:203006/111‑1 (MQ=255)

                                                                                                                                                                       aTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCC‑CCGGGcc                                                                                                                                                                                                              <  5:418727/149‑1 (MQ=255)

                                                                                                                                                                           gagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGtt                                                                                                                                                                                                           <  7:162410/149‑1 (MQ=255)

                                                                                                                                                                           gagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGtt                                                                                                                                                                                                           <  8:440588/149‑1 (MQ=255)

                                                                                                                                                                            agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATaaa                                                                                                                              >  4:194432/1‑224 (MQ=255)

                                                                                                                                                                            agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg                                                                            >  2:598184/1‑274 (MQ=255)

                                                                                                                                                                            agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCAACTGAGCCTTTTGTTTTATTTGATTTCTgg                                                                            >  2:273027/1‑274 (MQ=255)

                                                                                                                                                                            agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTGTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTg                                                                                                                                                                                                          >  7:55738/1‑149 (MQ=255)

                                                                                                                                                                            agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑CTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg                                                                            >  2:420096/1‑274 (MQ=255)

                                                                                                                                                                             gCGTACTAGGGGATTTTTATCGCAGGTAATCCATTAATGGAAGGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCa                                                                          <  2:628548/275‑1 (MQ=255)

                                                                                                                                                                             gCGTACGAGGGCATTTTTATCGCAGGTAATGCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGc                                                                                                                                                                                                         <  5:716192/149‑1 (MQ=255)

                                                                                                                                                                             gCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGACAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGc                                                                                                                                                                                                         <  8:562208/149‑1 (MQ=255)

                                                                                                                                                                             gCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGc                                                                                                                                                                                                         <  6:135587/149‑1 (MQ=255)

                                                                                                                                                                                tACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGgaga                                                                                                                                                        <  1:126340/194‑1 (MQ=255)

                                                                                                                                                                                   ggtggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTCGTTCGAAAAGCAAAAAGGCCCTCCTGTCGTGTGGCC‑TTTCCCTTGATTTGATGTCTGTCAGTTTATGGCGGGCGTCCTGCCCTCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTccc                                                                    <  4:598184/272‑1 (MQ=255)

                                                                                                                                                                                     gggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAc                                                                                                                                                                                                 <  8:839344/149‑1 (MQ=255)

                                                                                                                                                                                       gCATTTTTATTGCAGGGAATACATTAATTGAATGTTAGGTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTAct                                                                <  1:22400/275‑1 (MQ=255)

                                                                                                                                                                                       gCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGt                                                                                                                                                                                               >  6:197846/1‑149 (MQ=255)

                                                                                                                                                                                         aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTc                                                                                                                                                                                             >  7:772629/1‑149 (MQ=255)

                                                                                                                                                                                         aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTc                                                                                                                                                                                             >  7:377376/1‑149 (MQ=255)

                                                                                                                                                                                          ttttttTCGGAGGTAATCGATTAATTNNTGTTTAGTTCGAAAAGCGAAAAGTCCCTCCNTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTNGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTAct                                                                <  3:87900/272‑1 (MQ=17)

                                                                                                                                                                                          tttttATCGCAGGTAATCCATTAATTGCATGTTAGTTCGAAAAGCAAAAAGCCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGCAGTCCTGCCCGCCACCCTCCTGGCCGGCTCTTCCCAACGTTTAAAT‑TCCCTCCCGTCCGTGTCTTCCTTCCCGGGGGAGTGTTCACCGACCAACAACCGCTAACACAAAAGCCCCCGTCTTCCGCTTTCACCTTCTGTTTTATTCTATTTCTGGCCGTTCCctcctcc                                                              >  4:300299/1‑273 (MQ=255)

                                                                                                                                                                                                   cAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTc                                                                                                                                                                                  >  5:5755/1‑149 (MQ=255)

                                                                                                                                                                                                      gTAATACATTAATTAAATGTTAGTGCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGgaga                                                  <  1:165294/274‑1 (MQ=255)

                                                                                                                                                                                                          tCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCcggcgg                                                                                                                                                                           >  8:272142/1‑149 (MQ=255)

                                                                                                                                                                                                           ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGa                                                                                                                                                                          >  5:272214/1‑149 (MQ=255)

                                                                                                                                                                                                             aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGAtt                                                                                                                                                                        >  7:699050/1‑149 (MQ=255)

                                                                                                                                                                                                             aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGAtt                                                                                                                                                                        >  6:627153/1‑149 (MQ=255)

                                                                                                                                                                                                             aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGAtt                                                                                                                                                                        >  8:662272/1‑149 (MQ=255)

                                                                                                                                                                                                              ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGAttt                                                                                                                                                                       <  5:1102537/149‑1 (MQ=255)

                                                                                                                                                                                                              ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAacaaca                                                                                                                                    >  2:288687/1‑184 (MQ=255)

                                                                                                                                                                                                              ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCacacc                                         >  1:69112/1‑274 (MQ=255)

                                                                                                                                                                                                                  ttGAATGTTAGTTCGAAAAGCAAAAATGCCATCCTTACGGATGGCC‑TTTCGCTTGATTTGATATCTGGCAGTTTATGGCGGACGGCCTTCCCGCCACCTTTCGGGCCGTTGCTTCTCAACGTTCAAAT‑CCGCTCCCCCCCGCTTTGTcc                                                                                                                                                                   >  8:929686/1‑149 (MQ=18)

                                                                                                                                                                                                                  ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTcc                                                                                                                                                                   >  6:933647/1‑149 (MQ=255)

                                                                                                                                                                                                                  ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTcc                                                                                                                                                                   <  5:1015747/149‑1 (MQ=255)

                                                                                                                                                                                                                  ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATa                                                                                                                                <  4:288687/184‑1 (MQ=255)

                                                                                                                                                                                                                  ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGgaga                                                  <  4:467755/262‑1 (MQ=255)

                                                                                                                                                                                                                  ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGACAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTcc                                                                                                                                                                   >  7:215526/1‑149 (MQ=255)

                                                                                                                                                                                                                         ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCgg                                                                                                                                                            >  6:824060/1‑149 (MQ=255)

                                                                                                                                                                                                                           aGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGga                                                                                                                                                          <  6:5755/149‑1 (MQ=255)

                                                                                                                                                                                                                             ttCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCATCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGgaga                                                                                                                                                        <  5:513278/149‑1 (MQ=255)

                                                                                                                                                                                                                              tagAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGagag                                                                                                                                                       <  6:651389/147‑1 (MQ=255)

                                                                                                                                                                                                                                gAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAgtg                                                                                                                                                     >  5:146841/1‑149 (MQ=255)

                                                                                                                                                                                                                                 aaaaGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTGCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCCCCCGGCTGATTTGTCCTTCTCCGGCGAGTGTTCANCGACNAACAACAGATAAAACa                                                                                                                           >  4:671337/1‑174 (MQ=17)

                                                                                                                                                                                                                                  aaaGCAAAAAGGCCATCCTGTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGATATTGtt                                                                                                                                                   <  8:105140/149‑1 (MQ=255)

                                                                                                                                                                                                                                   aaGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTc                                                                                                                                                  >  7:832040/1‑149 (MQ=255)

                                                                                                                                                                                                                                    aagcgaaaGGCCATCCGAAAGGATGGCCTTTTTGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAgtg                                                                                                                                                     <  6:488484/141‑1 (MQ=11)

                                                                                                                                                                                                                                       aaaaaGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTTACCg                                                                                                                                              >  5:726573/1‑149 (MQ=255)

                                                                                                                                                                                                                                        aaaaGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGa                                                                                                                                             >  7:992165/1‑149 (MQ=255)

                                                                                                                                                                                                                                           aGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaa                                                                                                                                          >  8:515901/1‑149 (MQ=255)

                                                                                                                                                                                                                                             gCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaac                                                                                                                                        <  8:992165/149‑1 (MQ=255)

                                                                                                                                                                                                                                              ccATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaa                                                                                                                                         <  8:543978/147‑1 (MQ=255)

                                                                                                                                                                                                                                              ccATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaa                                                                                                                                         >  7:543978/1‑147 (MQ=255)

                                                                                                                                                                                                                                               cATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaacaa                                                                                                                                      >  7:179256/1‑149 (MQ=255)

                                                                                                                                                                                                                                                aTCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAaacaac                                                                                                                                     >  6:443767/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       cGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCCCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCCTGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGTCAt  >  2:377037/1‑273 (MQ=255)

AAAGAGGTGCTGGATACGATGATTGGGCTGGCGCAGTCGGGTATGACAATGTTGTGTGTAACACATGAGATGGGGTTTGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAT‑CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCAT  >  NC_000913/3422985‑3423504

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