BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,100,791	G→C	5.2%	intergenic (+4/+60)	ycdZ → / ← csgG	DUF1097 family inner membrane protein/curli production assembly/transport outer membrane lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,100,791	0	G	C	5.2%	838.8 / 25.3	239	intergenic (+4/+60)	ycdZ/csgG	DUF1097 family inner membrane protein/curli production assembly/transport outer membrane lipoprotein
Reads supporting (aligned to +/- strand): ref base G (112/114); new base C (8/5); total (120/119)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.27e-01

Only 100 of 245 total aligned reads displayed.

CATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGAAAATGACCCCTGTTTCGATAGCCGACATCAGGCACAGCATAACAGGTTCGTTCGAGGTGTAACCCACTTCCCCTTCAAGCAAGCGCTGGTAGTCAATAA  >  NC_000913/1100521‑1101062

  tCTGGAAATTCTCGGTTATGTCATAACCGGTATGGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCt                                                                                                                                                                                                                                                                       <  1:269288/279‑1 (MQ=255)

                gTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTaaa                                                                                                                                                                                                                                                             >  2:242597/1‑275 (MQ=255)

                              ggTATTTTCGCGTTTTGTGGGTGTAGTCGGGCCAAAACGCTGCTGCTTTCAATTGTTCCGGGGACGTTTATAGGGTCATGCGCGACGTTTGCCGGGGAGGGTGACTGGAAACTGGTGTTACCTTCGCGGGCGCTTGGGCTGATATTTGGTTACGCAAGGAAAAACAG‑TGGCCTGTGGCTGGCGGGGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCg                                                                                                                                                                                                                                                <  4:343978/274‑1 (MQ=255)

                                    gTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGCAACTGGTGTTACCTTCGCTGGCTCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCCCGTAGTGCAAAGACCGAGCACCGTGAGCAGTAAAGCAAAAATAAAGCGTGATGGGCACTCACGCTTTCGCTTAAACCGTAAAATGCCGgactcta                                                                                                                                                                                                                                         >  4:92856/1‑270 (MQ=255)

                                     tCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATaa                                                                                                                                                                                                                                        >  2:529775/1‑275 (MQ=255)

                                         tttttTGATG‑GTCTTCGGGCCAACCGGCTGCTGCTTTCTTTGGTTCCGGGGGTGTTTGTAGGCGGCGGGGGGACGTTTGCCGGGGAGGGTGGCTGGAAACTGTTGGTACCTTGGCTGGCGCGTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGGGTAGTGCAAAGGCGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCg                                                                                                                                                                                                                                   <  4:682477/274‑1 (MQ=255)

                                                 taagcggcagggCCATACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTAt                                                                                                                                                                                                                         <  1:215055/268‑1 (MQ=255)

                                                     tATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTGGTGCACAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGTTTAAACAGTAAAATGACGGCTCATAAATCCCGCTTTTTTATCTGTc                                                                                                                                                                                                                    >  1:89987/1‑279 (MQ=255)

                                                            gCCAAACAGCTGCTGCTTTCATTTGTTCCGGGCACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGCCGCTTGGGCTGATATTTGGTTACGCCATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCACAGACGGCGCACCGTGACCAGGCAATCCACAATAAATCGTGAGGGGCACTCACGCTTTCGCTTAACCAGTAAAATGCCGGATGCTAATTTCGACTTTTTTATCTTTCAg                                                                                                                                                                                                                  >  4:120066/1‑274 (MQ=255)

                                                                 aCAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGATGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTcc                                                                                                                                                                                                            >  2:508277/1‑275 (MQ=255)

                                                                                    tttccctggACGGTTGTAGCGGCATGCGCGCCGGTTTCCGGGCTGGTGGACTTGAAACTGGTGTTACCTTTGCTGGCGCTTGGGCTATTATTTGGTGACGCAATGAAAAACAG‑TGGCCTGTGTCTGGCGGCTCGTTGTGTAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGc                                                                                                                                                                                         <  4:180573/268‑1 (MQ=255)

                                                                                    gTTCCGGGGACGTTGATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGc                                                                                                                                                                                         <  2:153769/275‑1 (MQ=255)

                                                                                     ttCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACAGATGGCGGTAt                                                                                                                                                                                    >  1:113681/1‑279 (MQ=255)

                                                                                                      ggCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGCCGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGCATg                                                                                                                                                                        >  2:669101/1‑274 (MQ=255)

                                                                                                          cATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGCCGGCGNCCCGTGAGCCGGAAAGCAACAATAAAGCGTGAGGGGCACCCAAGCTTTCGCTTAAAAAGTAAAATGCAGNATGAGAATTCCGGCTGTATTATCTGTCAGGCTTATGGTGGAACCCNCNCCTGGCGGAATTCCACCAGAACGTGATTCTGc                                                                                                                                                               >  3:83784/1‑279 (MQ=255)

                                                                                                          cATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACACTAACATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGCCATATGGCGGTATTTCACCAGTATGTCAtt                                                                                                                                                                   >  2:190064/1‑275 (MQ=255)

                                                                                                              gggggagTTTGCCGGGCGAGTGGACTGGAAATTGGGTTTACCTTCGCGTGGGGTTGGGCGGAAAGTTGGTTACCCAATAAAAAACAG‑GAGCCTGGGGCTGGGGGCTAGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCTACATATGGCGGTATTTCACCAGAATGTCATTCTg                                                                                                                                                                <  4:529771/268‑1 (MQ=255)

                                                                                                                              cAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACt                                                                                                                                                                                                                                                                             <  8:911422/149‑1 (MQ=255)

                                                                                                                               aGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTc                                                                                                                                                                                                                                                                            >  8:1646691/1‑149 (MQ=255)

                                                                                                                                    gACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGAGTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCaaa                                                                                                                                          <  2:120067/274‑1 (MQ=255)

                                                                                                                                       tGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTc                                                                                                                                                                                                                                                                    >  7:1330379/1‑149 (MQ=255)

                                                                                                                                       tGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTc                                                                                                                                                                                                                                                                    <  7:1497244/149‑1 (MQ=255)

                                                                                                                                         gAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCAGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGc                                                                                                                                                                                                                                                                  >  6:1628106/1‑149 (MQ=255)

                                                                                                                                              tGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTaaa                                                                                                                                                                                                                                                             >  8:1385774/1‑149 (MQ=255)

                                                                                                                                                gtggTACCTTCGCGGGCGCTTGGGCTGAGATTTGGTTACGCAATGAAAAACAG‑TGGCCGGTGGCTGGCGGCGCGTAGGGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCa                                                                                                                              <  2:258953/274‑1 (MQ=255)

                                                                                                                                                   ttACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGTAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTCAACAGTa                                                                                                                                                                                                                                                        <  6:1327425/149‑1 (MQ=255)

                                                                                                                                                   ttACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTaaa                                                                                                                                                                                                                                                             <  6:956215/144‑1 (MQ=255)

                                                                                                                                                   ttACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCAACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGt                                                                                                                          >  2:599703/1‑275 (MQ=255)

                                                                                                                                                       cTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCAc                                                                                                                                                                                                                                                                              <  5:721107/123‑1 (MQ=255)

                                                                                                                                                       cTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAAt                                                                                                                                                                                                                                                    >  5:996410/1‑149 (MQ=255)

                                                                                                                                                             tGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGga                                                                                                                                                                                                                                              <  5:768231/149‑1 (MQ=255)

                                                                                                                                                              ggCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGgat                                                                                                                                                                                                                                             <  7:561190/149‑1 (MQ=255)

                                                                                                                                                                 gcTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGgatgat                                                                                                                                                                                                                                          >  6:646048/1‑149 (MQ=255)

                                                                                                                                                                   ttGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATaa                                                                                                                                                                                                                                        >  5:1488254/1‑149 (MQ=255)

                                                                                                                                                                    tGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCa                                                                                                                                                                                                                                                                           >  4:479194/1‑113 (MQ=255)

                                                                                                                                                                     tggCTGGTACTTTGTTACGGAATGAAAAACAGTTGGCCGTGGGCGGGCGG‑GGGTAAGGTGAAGACCGCGCAACCGGAGCAAGAAAACAAAAATAAAGCGTTAGGTGCACTTACGCTTTTGCTTAAACAGGAAAATGGCCGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCa                                                                                                        <  4:374086/274‑1 (MQ=255)

                                                                                                                                                                     gggCTGTGATTAGGTTAAGCAATGAAAAACAG‑TGGCCGGTGGCTGGCGGCGCGGAGTGCAAAGAAGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTTAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCTGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCa                                                                                                        <  4:239977/275‑1 (MQ=255)

                                                                                                                                                                       gCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTcc                                                                                                                                                                                                                                    <  7:1050463/149‑1 (MQ=255)

                                                                                                                                                                       gCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTcc                                                                                                                                                                                                                                    >  5:908108/1‑149 (MQ=255)

                                                                                                                                                                        cTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCt                                                                                                                                                                                                                                                                       <  2:479198/113‑1 (MQ=255)

                                                                                                                                                                         tGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCgg                                                                                                                                                                                                                                  >  8:743881/1‑149 (MQ=255)

                                                                                                                                                                                 ggTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCt                                                                                                                                                                                                                                                                 <  6:1301329/110‑1 (MQ=255)

                                                                                                                                                                                 ggTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCt                                                                                                                                                                                                                                                                 >  5:1301329/1‑110 (MQ=255)

                                                                                                                                                                                     aCGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTg                                                                                                                                                                                                                      <  6:726209/149‑1 (MQ=255)

                                                                                                                                                                                      cGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGt                                                                                                                                                                                                                     <  5:391322/149‑1 (MQ=255)

                                                                                                                                                                                       gCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGGAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTc                                                                                                                                                                                                                    >  6:251521/1‑149 (MQ=255)

                                                                                                                                                                                       gCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTc                                                                                                                                                                                                                    >  5:813762/1‑149 (MQ=255)

                                                                                                                                                                                       gCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTc                                                                                                                                                                                                                    >  6:65509/1‑149 (MQ=255)

                                                                                                                                                                                        cAATGACAAACAG‑TGGCCAGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAGAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACTGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCa                                                                                                                                                                                                                   <  8:813473/149‑1 (MQ=255)

                                                                                                                                                                                         aaTGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCg                                                                                                                                                                                                                                   <  5:1167538/132‑1 (MQ=255)

                                                                                                                                                                                         aaTGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAg                                                                                                                                                                                                                  >  5:879362/1‑149 (MQ=255)

                                                                                                                                                                                          aTGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAgg                                                                                                                                                                                                                 >  6:403020/1‑149 (MQ=255)

                                                                                                                                                                                             aaaaaCAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGAtt                                                                                                                                                                                                              >  8:1155440/1‑149 (MQ=255)

                                                                                                                                                                                                  cAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGt                                                                                                                                                                                                         <  6:323579/149‑1 (MQ=255)

                                                                                                                                                                                                  cAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGt                                                                                                                                                                                                         <  8:60222/149‑1 (MQ=255)

                                                                                                                                                                                                  cAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGt                                                                                                                                                                                                         <  6:1404320/149‑1 (MQ=255)

                                                                                                                                                                                                        gCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCa                                                                                                                                                                                                                                                                           >  4:546763/1‑78 (MQ=255)

                                                                                                                                                                                                        gCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGGTACCATCATTAATCAGGAAAATGACCCCTGTTTCGATAGCCGACATCAg                                                                      >  2:344249/1‑275 (MQ=255)

                                                                                                                                                                                                           tgtgGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGgatga                                                                                                                                                                                                                                           <  1:58900/107‑1 (MQ=255)

                                                                                                                                                                                                           tgtgGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGgatga                                                                                                                                                                                                                                           >  3:58900/1‑107 (MQ=255)

                                                                                                                                                                                                            gtgGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGAc                                                                                                                                                                                                <  6:1183302/149‑1 (MQ=255)

                                                                                                                                                                                                                       cgcgTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGAt                                                                                                                                                                                                               >  2:387750/1‑123 (MQ=255)

                                                                                                                                                                                                                       cgcgTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTa                                                                                                                                                                                     <  6:813762/149‑1 (MQ=255)

                                                                                                                                                                                                                        gcgTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTAt                                                                                                                                                                                    >  7:469745/1‑149 (MQ=255)

                                                                                                                                                                                                                             gTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCAc                                                                                                                                                                               <  5:1618963/149‑1 (MQ=255)

                                                                                                                                                                                                                              tGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCAcc                                                                                                                                                                              >  5:335273/1‑149 (MQ=255)

                                                                                                                                                                                                                               gCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCa                                                                                                                                                                             <  5:1345824/149‑1 (MQ=255)

                                                                                                                                                                                                                               gCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCa                                                                                                                                                                             >  6:524176/1‑149 (MQ=255)

                                                                                                                                                                                                                                 aaaGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGa                                                                                                                                                                           >  6:70455/1‑149 (MQ=255)

                                                                                                                                                                                                                                  aaGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGaa                                                                                                                                                                          <  6:892274/149‑1 (MQ=255)

                                                                                                                                                                                                                                  aaGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGGGCGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGAATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGAAAATGACCCCTGTTTCGATAGCCGACATCAGGCACAGCATAACAGGTTCGTTCGAg                                             <  2:694569/274‑1 (MQ=255)

                                                                                                                                                                                                                                   aGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTc                                                                                                                                                                                                                                     >  6:1494590/1‑89 (MQ=255)

                                                                                                                                                                                                                                    gACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATg                                                                                                                                                                        >  8:1623380/1‑149 (MQ=255)

                                                                                                                                                                                                                                    gACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATg                                                                                                                                                                        >  8:1091260/1‑149 (MQ=255)

                                                                                                                                                                                                                                    gACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATg                                                                                                                                                                        >  7:489904/1‑149 (MQ=255)

                                                                                                                                                                                                                                        gcgcACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACACATGGCGGCATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTCTTTTGCAAATCCCACACACCACGGTCGATACCATCATTAATCAGGAAACTGACCCCTGTTTCGATAGCCGAACTCAGCAACCGCCAAACAGGTTCGTCACAGGTTTAAc                                      >  4:257698/1‑275 (MQ=255)

                                                                                                                                                                                                                                          gcACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTAACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCATGGTCGATACCATCATTAATCAGGAAAATGACCCCTGTTTCGCTAGCCGACATCAGGCACAGCATAAGAGGTTCTTTCGAGgtgcaccc                                     >  2:568601/1‑269 (MQ=255)

                                                                                                                                                                                                                                             ccGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGc                                                                                                                                                               >  5:470116/1‑149 (MQ=255)

                                                                                                                                                                                                                                              cGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGAAAATGACCCCTGTTTCGATAGCCGACATCAGGCACAGCATAACAGGTTCGTTCGAGGTGTAACCCACTTtcc                             >  1:102048/1‑278 (MQ=255)

                                                                                                                                                                                                                                                     aGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGc                                                                                                                                                       >  8:23919/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     aGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGc                                                                                                                                                       <  5:499005/149‑1 (MQ=255)

                                                                                                                                                                                                                                                             aaaaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTAt                                                                                                                                                                                    >  7:1617672/1‑112 (MQ=255)

                                                                                                                                                                                                                                                             aaaaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTAtttt                                                                                                                                               <  7:870990/149‑1 (MQ=255)

                                                                                                                                                                                                                                                             aaaaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTAtttt                                                                                                                                               >  8:33255/1‑149 (MQ=255)

                                                                                                                                                                                                                                                             aaaaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTAtttt                                                                                                                                               >  5:1281063/1‑149 (MQ=255)

                                                                                                                                                                                                                                                              aaaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTg                                                                                                                                              <  6:501810/149‑1 (MQ=255)

                                                                                                                                                                                                                                                               aaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGc                                                                                                                                                                                         >  7:1687385/1‑105 (MQ=255)

                                                                                                                                                                                                                                                               aaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGc                                                                                                                                             >  5:662866/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                    aaGCGTGAGTGCCCCTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTg                                                                                                                                              <  7:1175832/143‑1 (MQ=255)

                                                                                                                                                                                                                                                                     aGCGTGAGTGGCACTCCAGCTTTCGCTTAAACAGTAAAATGCCGGATGTTAATTCCTGTTTTTTTATTTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGAAAATGACCCCTGTTTCGATAGCCGACATCAGGCACAGCATAACAGGttcgttcg                                               <  3:24477/237‑1 (MQ=255)

                                                                                                                                                                                                                                                                     aGCGTGAGTGCCCCTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAAc                                                                                                                                                                                                    <  8:1258344/88‑1 (MQ=39)

                                                                                                                                                                                                                                                                     aGCGTGAGTGCCCCTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATcc                                                                                                                                       >  7:1399653/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     aGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGt                                                                                                                                                                                                         <  5:1402021/83‑1 (MQ=255)

                                                                                                                                                                                                                                                                      gCGTGAGTGCCCCTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATccc                                                                                                                                      >  6:481626/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                         tGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAAcc                                                                                                                                                                                                   >  5:1654696/1‑85 (MQ=255)

                                                                                                                                                                                                                                                                         tGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAAcc                                                                                                                                                                                                   <  6:1654696/85‑1 (MQ=255)

                                                                                                                                                                                                                                                                             gtgCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGGTTCCGGTGGAACCGACCTATGGCGGTATTTCACCAGGATGTCATTCTGCCGTTCTGCTTTATTTTGCGAATCCCACAGAcc                                                                                                                               >  6:612105/3‑149 (MQ=255)

                                                                                                                                                                                                                                                                             gggCACTCAGGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCcaca                                                                                                                                   >  5:1539310/1‑145 (MQ=255)

                                                                                                                                                                                                                                                                             gggCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGAcc                                                                                                                               >  8:1334063/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                              ggCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCa                                                                                                                              <  6:113200/149‑1 (MQ=255)

CATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAG‑TGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGAAAATGACCCCTGTTTCGATAGCCGACATCAGGCACAGCATAACAGGTTCGTTCGAGGTGTAACCCACTTCCCCTTCAAGCAAGCGCTGGTAGTCAATAA  >  NC_000913/1100521‑1101062

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