BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	747,902	C→G	5.7%	intergenic (‑133/+19)	abrB ← / ← ybgO	regulator of aidB expression; inner membrane protein/putative fimbrial protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	747,902	0	C	G	5.7%	437.8 / 12.5	124	intergenic (‑133/+19)	abrB/ybgO	regulator of aidB expression; inner membrane protein/putative fimbrial protein
Reads supporting (aligned to +/- strand): ref base C (52/65); new base G (3/4); total (55/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-01

Only 100 of 127 total aligned reads displayed.

AGTAATTCCGCGCATACTAAAGATGATCCCGGCGATCATCGGCCCAAGTAATAGCGCCGCCGGGAGATGTACCGCGAGGAAACCAATAGAAAGGAGAAGTGATAGCACACATAACATTCCCCAC‑TGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCATATTTATTGAAGGTATAAGGCGAGGCGTTTTCATCCAGTAAACGTAGCTTCAGACCATTACCAATCAACAATGCTTTATCTTCTTCAACCAACGTTTCTTCGGTATAAAAGGAG  >  NC_000913/747632‑748173

agTAATTCCGCGCATACTAAAGATGATCCCGGCGATCATCGGCCCAAGTAATAGCGCCGCCGGGAGATGTACCGCGAGGAAACAAATAGAAAGGAGAAGTGATAGCACACATAACATTCCCCAC‑TGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCaaa                                                                                                                                                                                                                                                                              <  2:404441/274‑1 (MQ=255)

     ttCCGCGCATACTAAAGATGATCCCGGCGATCATCGGCCCAAGTAATAGCGCCGCCGGGAGATGTACCGCGAGGAAACCAATAGAAAGGAGAAGTGATAGCACACATAACATTCCCCAC‑TGCAAAACTGGCATCCTTCGCTCCCTAATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATCACCGAAATCAAGGGAAAGGGGAAAATCAGAAATTTTCCGAAAGCGCTGCCCTTTCCGGCGGCAAAGGAcat                                                                                                                                                                                                                                                                        >  4:40249/1‑275 (MQ=255)

          tgcATACTAAAGATGATCCCGGCGATCATCGGCACAAGTAATAGCGCCGCCGGGAAATGTACCGCGAGGAAACCAATAGAAAGGAGAAGTGATAGCACACATAACATTCCCCAC‑TGCAAAACTGGCATCCTTCGCTCCATCATTATTTTCTCTTTGGTGTATCGATCCAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAAGTTTACGAAAGAGATGCCCTTTACGGCGGCAAAGGGCAACATgg                                                                                                                                                                                                                                                                   >  4:407324/2‑275 (MQ=255)

                   aagaagAGCCCCGCGGTCATCCGCCCAAGAAAGAGAGCCGCCGGGGGATGTACCCGGAGGAAACCAATAGAAAGGAGAAGTGATAGCACACATAACATTCCCCAC‑TGCAAAACTGGCAGCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAAGTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGt                                                                                                                                                                                                                                                          <  4:325774/275‑1 (MQ=255)

                                     accGGCCCAAGTAAAAGGGCCGCCGGGGGATGTACCGGGAGGAAACCAAAAGAAAGGAGAAGGGATAGCACACATAACATTCCCCAC‑TGCAAAACTGGCATCCGTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCTAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTaaa                                                                                                                                                                                                                                          <  4:515945/271‑1 (MQ=255)

                                      taggCCCAAGAAATAGGGCAGCCGGGAGATGTACCGGGAGGAAACCAATAGAAAGGAGACTTGATAGCACACATAACATTCCCCAC‑TGCAAAACTGGCATCCTTCGCCCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAg                                                                                                                                                                                                                                       <  4:109547/273‑1 (MQ=255)

                                               gTAATAGCGCCGCCGGGAGATGTACCGCGAGGAAACCAATAGAAAGGAGAAGTGATAGCACACATAACATTCCCCAC‑TGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCCATCAAGGAAAAGAGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATCGTTAATCTTAAAAACCACCACTg                                                                                                                                                                                                                               >  4:140667/1‑274 (MQ=255)

                                                        ccgccgGGAGATGTACCGCGAGGAAACCAATAGAAAGGAGAAGTGATAGCACACATAACATTCCCCAC‑TGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcac                                                                                                                                                                                                                                    >  4:419008/1‑260 (MQ=255)

                                                         cgccgGGAGATGTACCGCGAGGAAACCAATAGAAAGGAGAAGTGATAGCACACATAACATTCCCCAC‑TGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAgggg                                                                                                                                                                                                                    >  4:480372/1‑275 (MQ=255)

                                                            cgGGAGATGTACCGCGAGGAAACCAATAGAAAGGAGAAGTGATAGCACACATAACATTCCCCAC‑TGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACt                                                                                                                                                                                                                                <  2:419009/260‑1 (MQ=255)

                                                                         gcgAGGAAACCAATTGAAAGGAGAAGTGATAGTACACATAACATTCCCCAC‑TGGAAAGCTGGCATCCTTCGCTCCCTCATTATGTTCTCGTTGGTGGATCGATACAGGATAAAGAAAGATGGGCCAGGAAAAAAGAACAGCAAGCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGGGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgggg                                                                                                                                                                                                     <  4:216391/274‑1 (MQ=255)

                                                                                       aGAAAGGAGAAGTGATAGCACACATAACATTCCCCAC‑TGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGAc                                                                                                                                                                                      <  2:187409/275‑1 (MQ=255)

                                                                                             gagaAGTGATAGCACACATAACATTCCCCAC‑TGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTTTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATTCCCTTTCCGGCGGCAAAGGGCAACATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGAGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGt                                                                                                                                                                                >  2:234123/1‑275 (MQ=255)

                                                                                               aaaatcaTAACACACATTAC‑TTNNCCACANGCAAAACTGTCATCCTTCACTACCTTAATATTTTATCTATGGGGTATAGATACAGAATAAAAAAAGCTGAGCCAGAAAAAAAAAACAGCAACCAAGGAAAACGAGAAAATCCGCAATTTTCTGAAAGCGATCGCCTTACCGGCGGCAAAGGGAATCATGGTAAATCAATTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTc                                                                                                                                                                              <  3:7230/269‑1 (MQ=255)

                                                                                                 aGTGATAGCACACATAACATTCCCCAC‑TGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAg                                                                                                                                                                            <  2:480373/275‑1 (MQ=255)

                                                                                                       aGCACACATAACATTCCCCAC‑TGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAgg                                                                                                                                                                                                                                                                            >  4:148389/1‑173 (MQ=255)

                                                                                                        gCACACATAACATTCCCCAC‑TGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGtat                                                                                                                                                                      >  2:141827/1‑274 (MQ=255)

                                                                                                           cacaTAACATTCCCCAC‑TGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGcat                                                                                                                                                                                                                                                                        <  2:148390/173‑1 (MQ=255)

                                                                                                            acaTAACATTCCCCAC‑TGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGc                                                                                                                                                                                                                   >  4:368499/1‑225 (MQ=255)

                                                                                                             caTAACATTCCCCAC‑TGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAgg                                                                                                                                                                                                                                                                            >  4:401616/1‑167 (MQ=255)

                                                                                                                aaCATTCCCCAC‑TGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGc                                                                                                                                                                                                               <  2:368500/225‑1 (MQ=255)

                                                                                                                 aCATTCCCCAC‑TGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGcat                                                                                                                                                                                                                                                                        <  2:401617/167‑1 (MQ=255)

                                                                                                                       cccAC‑TGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGAc                                                                                                                                                      >  4:77599/1‑275 (MQ=255)

                                                                                                                                 aaaCTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAgg                                                                                                                                                                                                                                                                            >  8:799064/1‑148 (MQ=255)

                                                                                                                                 aaaCTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGa                                                                                                                                                                                    <  2:265894/236‑1 (MQ=255)

                                                                                                                                 aaaCTGGCATCCTTCGCTCCCTCATTATTTTCGCTTTGGTGGATCGATACAGCATAAAGAAAACTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAAGTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGAGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGAc                                                                                                                                                      <  2:508765/266‑1 (MQ=255)

                                                                                                                                    cTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAggg                                                                                                                                                                                                                     >  3:138925/1‑200 (MQ=255)

                                                                                                                                        cATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCa                                                                                                                                                                                                                 <  1:138925/200‑1 (MQ=255)

                                                                                                                                         aTCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATCGGTTTTTTCGTAACCCAACATGTCCATGGTGa                                                                                                                                     >  4:389838/1‑275 (MQ=255)

                                                                                                                                              tCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGc                                                                                                                                                                                                                                      >  2:466123/1‑173 (MQ=255)

                                                                                                                                               cGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATATAGAAAGCTGGGCCAGTAAAACATAACAGCAAACAAGGAAAAGGGGCAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCCAAGGGGCCGGCCTCGATGGGGTTGCCTGGATTGGACGAAAGTTCAGCCGTATAGGTTTTTTGGTAGACCAACATGTCCCGGGTTAACTCga                                                                                                                               >  4:51190/1‑274 (MQ=255)

                                                                                                                                                  tCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcacca                                                                                                                                                                                                                                  <  4:466122/173‑1 (MQ=255)

                                                                                                                                                       cATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTAt                                                                                                                                                                                                                            >  4:317482/1‑174 (MQ=255)

                                                                                                                                                           attTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGaa                                                                                                                                                                                                                        <  2:317483/174‑1 (MQ=255)

                                                                                                                                                                  ctTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGc                                                                                                                                                                                                                                      >  2:545587/1‑153 (MQ=255)

                                                                                                                                                                   tttGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAgg                                                                                                                                                                                                                                                                            >  4:77675/1‑114 (MQ=255)

                                                                                                                                                                   tttGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGcatc                                                                                                                                                                                                                                                                       >  4:569774/1‑119 (MQ=255)

                                                                                                                                                                    ttGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGc                                                                                                                                                                                                                                      >  2:547839/1‑151 (MQ=255)

                                                                                                                                                                    ttGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAAGAGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCATATTTATTGa                                                                                                          <  4:89418/275‑1 (MQ=255)

                                                                                                                                                                      ggTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcacca                                                                                                                                                                                                                                  <  4:545586/153‑1 (MQ=255)

                                                                                                                                                                       gtgtATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGcat                                                                                                                                                                                                                                                                        <  2:77676/114‑1 (MQ=255)

                                                                                                                                                                       gtgtATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATgg                                                                                                                                                                                                                                                                   <  2:569775/119‑1 (MQ=255)

                                                                                                                                                                          tATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGcatca                                                                                                                                                                                                                                                                      >  2:168153/1‑113 (MQ=255)

                                                                                                                                                                              gATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGt                                                                                                                                                                                                                                                                  <  4:168152/113‑1 (MQ=255)

                                                                                                                                                                               aTACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAg                                                                                                                                                                                                                                                                             <  6:1066682/101‑1 (MQ=255)

                                                                                                                                                                                  aaaaaTAAAAAAAGCTGGACCAAGAAAAAAAAACAGCAATCAAGGCAAAGGGGAAAATCAGCAATTTTCGTAAAGAGNTGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCATATTTATTGAAGGTATAaggcgag                                                                                            <  3:117321/271‑1 (MQ=255)

                                                                                                                                                                                    gCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGa                                                                                                                                                                                                                         <  7:799064/148‑1 (MQ=255)

                                                                                                                                                                                     cATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCaa                                                                                                                                                                                                                                                                               <  7:728137/93‑1 (MQ=255)

                                                                                                                                                                                     cATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCaa                                                                                                                                                                                                                                                                               >  8:728137/1‑93 (MQ=255)

                                                                                                                                                                                     cATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCAtattta                                                                                                              >  4:504998/1‑254 (MQ=255)

                                                                                                                                                                                     cATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTANAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGGGTAGTGAGCATATTTATTGAAGGTATAgggggaggag                                                                                         >  3:116822/1‑273 (MQ=255)

                                                                                                                                                                                         aagaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCATATTTATTGa                                                                                                          <  2:504999/254‑1 (MQ=255)

                                                                                                                                                                                            aaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATg                                                                                                                                                                                                                                                                    <  6:509413/97‑1 (MQ=255)

                                                                                                                                                                                            aaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTaaaa                                                                                                                                                                                                                                         >  6:186716/1‑124 (MQ=255)

                                                                                                                                                                                            aaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTaaaa                                                                                                                                                                                                                                         <  5:186716/124‑1 (MQ=255)

                                                                                                                                                                                              agCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTc                                                                                                                                             <  1:62746/214‑1 (MQ=255)

                                                                                                                                                                                                   ggCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCg                                                                                                                                                                         >  3:121903/1‑181 (MQ=255)

                                                                                                                                                                                                       aGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGtata                                                                                                                                                                     <  1:121903/181‑1 (MQ=255)

                                                                                                                                                                                                       aGGAAAAAATAACAGCAAGCCATGAAAAGGGGTAAAGCAGCAATTGTCGGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCATATTTATTGAAGGTATAAGGCGAGGCGTTTTCATCCAGTAAACGt                                                                       <  2:389839/275‑1 (MQ=255)

                                                                                                                                                                                                          aaaaaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTGCCGCCGGAAAGGGc                                                                                                                                                                                                                                                                          <  8:679965/77‑1 (MQ=37)

                                                                                                                                                                                                          aaaaaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTGCCGCCGGAAAGGGc                                                                                                                                                                                                                                                                          >  7:679965/1‑77 (MQ=37)

                                                                                                                                                                                                          aaaaaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAggg                                                                                                                                                                                                                                                                           <  6:735154/76‑1 (MQ=255)

                                                                                                                                                                                                          aaaaaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAggg                                                                                                                                                                                                                                                                           >  5:735154/1‑76 (MQ=255)

                                                                                                                                                                                                            aaaaaaACAGCAAGCAAGGGAAGGGGGAAAATCAGCAATTTTCAGAAAGGGATGCCCTTTCCGGCGGCAAAGGGCAGCATGGTAAATCAGTTATAGTTAAACTTAAAGAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCATATTTATTGAAGGTATAAGGCGAGGCGTTTTCAGCCAGTAAACGTAGCtt                                                                  <  1:116822/275‑1 (MQ=255)

                                                                                                                                                                                                            aaaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATAGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCATATTTATTGAAGGTATAAGGCGAGGCGTTTTCATCCAGTAAACGTAGCtt                                                                  <  2:428611/275‑1 (MQ=255)

                                                                                                                                                                                                            aaaaTAACAGCAATCAAGGAAAAGGGGAAAAACAGCAATTGTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGAATCATGGTAAACCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGGGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCATATTTATTGAAGGTATAAGGCGAGGCGTTTTCATCCAGTAAACGTAGCtt                                                                  <  2:215276/275‑1 (MQ=255)

                                                                                                                                                                                                                  aCAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgg                                                                                                                                                                                                       >  8:810613/1‑136 (MQ=255)

                                                                                                                                                                                                                  aCAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgg                                                                                                                                                                                                       <  7:810613/136‑1 (MQ=255)

                                                                                                                                                                                                                     gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATc                                                                                                                                                                                                                                                             <  6:991130/79‑1 (MQ=255)

                                                                                                                                                                                                                     gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGt                                                                                                                                                                                                                                                          >  5:306382/1‑82 (MQ=255)

                                                                                                                                                                                                                     gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGt                                                                                                                                                                                                                                                          <  6:306382/82‑1 (MQ=255)

                                                                                                                                                                                                                     gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATggg                                                                                                                                                                                                      >  6:932131/1‑134 (MQ=255)

                                                                                                                                                                                                                     gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATggg                                                                                                                                                                                                      <  5:932131/134‑1 (MQ=255)

                                                                                                                                                                                                                       aaTCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcaccac                                                                                                                                                                                                                                 >  5:461465/1‑105 (MQ=255)

                                                                                                                                                                                                                         tCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGaa                                                                                                                                                                                                                        >  8:109085/1‑112 (MQ=255)

                                                                                                                                                                                                                         tCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGaa                                                                                                                                                                                                                        <  7:109085/112‑1 (MQ=255)

                                                                                                                                                                                                                          cAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATgg                                                                                                                                                                                                                                                                   <  6:862312/68‑1 (MQ=255)

                                                                                                                                                                                                                          cAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATgg                                                                                                                                                                                                                                                                   >  5:862312/1‑68 (MQ=255)

                                                                                                                                                                                                                              gAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGa                                                                                                                                                                                                          <  6:912091/121‑1 (MQ=255)

                                                                                                                                                                                                                              gAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGa                                                                                                                                                                                                          >  5:912091/1‑121 (MQ=255)

                                                                                                                                                                                                                                aaaGGGGAAAATCAGCAATTTTCTGAAAGTGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCATATTTATTGAAGGTATAAGGCGAGGCGTTTTCATCCAGTAAACGTAGCTTCAGACCATTACCAATcaaca                                              <  2:407325/275‑1 (MQ=255)

                                                                                                                                                                                                                                  aGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCt                                                                                                                                                                                               <  6:538301/128‑1 (MQ=255)

                                                                                                                                                                                                                                  aGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCt                                                                                                                                                                                               >  5:538301/1‑128 (MQ=255)

                                                                                                                                                                                                                                       aaaaTCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAg                                                                                                                                                                                                                       >  7:569128/1‑99 (MQ=255)

                                                                                                                                                                                                                                       aaaaTCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAg                                                                                                                                                                                                                       <  8:569128/99‑1 (MQ=255)

                                                                                                                                                                                                                                                             agaTGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAg                                                                                                                                                                                                                                       >  7:149913/1‑61 (MQ=255)

                                                                                                                                                                                                                                                             agaTGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAg                                                                                                                                                                                                                                       <  8:149913/61‑1 (MQ=255)

                                                                                                                                                                                                                                                                 gCCCTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGa                                                                                                                                                                                                          <  8:836753/86‑1 (MQ=38)

                                                                                                                                                                                                                                                                 gCCCTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGa                                                                                                                                                                                                          >  7:836753/1‑86 (MQ=38)

                                                                                                                                                                                                                                                                  cccTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATggg                                                                                                                                                                                                      <  8:212358/89‑1 (MQ=39)

                                                                                                                                                                                                                                                                  cccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAACGCACCACTGTATCGAAGGGGCCAGCCTCGATggggg                                                                                                                                                                                                    >  5:984104/1‑91 (MQ=255)

                                                                                                                                                                                                                                                                  cccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATgg                                                                                                                                                                                        >  5:640374/1‑103 (MQ=255)

                                                                                                                                                                                                                                                                  cccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATgg                                                                                                                                                                                        <  6:640374/103‑1 (MQ=255)

                                                                                                                                                                                                                                                                   ccTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGtat                                                                                                                                                                      <  8:342309/120‑1 (MQ=255)

                                                                                                                                                                                                                                                                   ccTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGtat                                                                                                                                                                      >  7:342309/1‑120 (MQ=255)

                                                                                                                                                                                                                                                                   ccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATggg                                                                                                                                                                                                      <  8:63591/88‑1 (MQ=255)

                                                                                                                                                                                                                                                                   ccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATggg                                                                                                                                                                                                      >  7:63591/1‑88 (MQ=255)

                                                                                                                                                                                                                                                                   ccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGaa                                                                                                                                                                                   <  5:288351/107‑1 (MQ=255)

                                                                                                                                                                                                                                                                   ccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTAGCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGaa                                                                                                                                                                                   >  6:288351/1‑107 (MQ=255)

                                                                                                                                                                                                                                                                            cggcAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATATCACTGGTGAAACCGGCGTATTCAGCATATTTATTGAAGGTATAAGGCGAGGCGTTTTCATCCAGTAAACGTAGCTTCAGACCATTACCAAACAACAATGCTTTATCTTCTTTAACCAACGTTTTTGTGGTATAAAagggg  >  4:21591/1‑273 (MQ=255)

AGTAATTCCGCGCATACTAAAGATGATCCCGGCGATCATCGGCCCAAGTAATAGCGCCGCCGGGAGATGTACCGCGAGGAAACCAATAGAAAGGAGAAGTGATAGCACACATAACATTCCCCAC‑TGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCATATTTATTGAAGGTATAAGGCGAGGCGTTTTCATCCAGTAAACGTAGCTTCAGACCATTACCAATCAACAATGCTTTATCTTCTTCAACCAACGTTTCTTCGGTATAAAAGGAG  >  NC_000913/747632‑748173

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