BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,351,590	C→T	V25I (GTC→ATC)	dcuS ←	sensor histidine kinase DcuS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,351,590	0	C	T	100.0%	487.8 / NA	138	V25I (GTC→ATC)	dcuS	sensor histidine kinase DcuS
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (55/83); total (55/83)

Only 100 of 138 total aligned reads displayed.

CGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTG‑CGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGG  >  NC_000913/4351343‑4351835

 gAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTG‑CGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCAtt                                                                                                                                                                                                                                                   <  1:1191014/251‑1 (MQ=255)

 cgacatgacgttcttccgatctgACAGTAAACAGCAGTGCGTTGCGTTTG‑CGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCAtt                                                                                                                                                                                                                                                   <  2:959785/229‑1 (MQ=255)

  aaGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTG‑CGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCAtt                                                                                                                                                                                                                                                   >  1:963563/1‑250 (MQ=255)

  aaGACGTTGCATATCGGTAACGACAATAAACAGCCGCTCGTTGC‑TTTGTCGTACGGCTTCCGCGCTGGCCTCGATGCCTCTCTGCGGCGGTTTTTCCTGCAAGCCCTGACGGATTTCCGGCGAGTTGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCTCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGAGGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTa                                                                                                                                                                                                                                                  <  2:650245/251‑1 (MQ=255)

       tttGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTG‑CGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTa                                                                                                                                                                                                                                                  >  2:1191993/1‑246 (MQ=255)

         tGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTG‑CGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGCGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACAAATAGCACCGCGAACAGTACCGCACTGATCATTAAGATCAc                                                                                                                                                                                                                                           >  2:633187/1‑251 (MQ=255)

          gCATATCGGTAACGACAATAAACAGCAGATCGTTGCGGTTG‑CGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAGGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACt                                                                                                                                                                                                                                          <  2:263057/251‑1 (MQ=255)

           cATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTG‑CGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACtg                                                                                                                                                                                                                                         <  1:27844/251‑1 (MQ=255)

              atCGGTAACGACAATAAACAGCAGATCGTTGCGTTTG‑CGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACtgtg                                                                                                                                                                                                                                       <  1:839722/250‑1 (MQ=255)

                   tAACGACAATAAACAGCAGATCGTTGCGTTTG‑CGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGATAACAGTACCGCACTGATCATTAAGATCACTGTGGTACt                                                                                                                                                                                                                                  >  2:26714/1‑250 (MQ=255)

                      cGACAATAAACAGCAGATCGTTGCGTTTG‑CGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAAt                                                                                                                                                                                                                              <  2:963563/251‑1 (MQ=255)

                       gACAATAAACAGCAGATCGTTGCGTTTG‑CGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCAtt                                                                                                                                                                                                                                                   >  1:959785/1‑229 (MQ=255)

                       gACAATAAACAGCAGATCGTTGCGTTTG‑CGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACtgt                                                                                                                                                                                                                                        <  2:909463/240‑1 (MQ=255)

                       gACAATAAACAGCAGATCGTTGCGTTTG‑CGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACtgt                                                                                                                                                                                                                                        >  1:909463/1‑240 (MQ=255)

                            aaaaCAGGAGATCGTTGCGTTTG‑CGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATc                                                                                                                                                                                                                        <  2:201545/250‑1 (MQ=255)

                             aaaCAGCAGATCGTTGCGTTTG‑CGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCg                                                                                                                                                                                                                       <  1:196801/251‑1 (MQ=255)

                                 agcagATCGTTGCGTTTG‑CGTACGGCTTCCGCGATGGCCTGGATGCCCCTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCCTCTCGTGTCATATCACTGATTTGCGAGCAGTAAATCAGATGCACCACCAATAGCACCGCGAACAGTACCGCACTGATCCTTAAGATAACTGTGGTACTCAATTTCATCGGACg                                                                                                                                                                                                                   >  2:880422/1‑251 (MQ=255)

                                    agATCGTTGCGTTTG‑CGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGttt                                                                                                                                                                                                                >  2:417896/1‑251 (MQ=255)

                                    agATCGTTGCGTTTG‑CGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGttt                                                                                                                                                                                                                <  2:446946/251‑1 (MQ=255)

                                     gATCGTTGCGTTTG‑CGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTg                                                                                                                                                                                                               <  1:654212/251‑1 (MQ=255)

                                        cGTTGCGTTTG‑CGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATGTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCg                                                                                                                                                                                                             <  2:261867/250‑1 (MQ=255)

                                           tgcgtttg‑cgtACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCACGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAAc                                                                                                                                                                                                         >  2:542719/1‑251 (MQ=255)

                                            gcgtttg‑cgtACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTc                                                                                                                                                                                                                           <  2:808567/232‑1 (MQ=255)

                                            gcgtttg‑cgtACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTc                                                                                                                                                                                                                           >  1:808567/1‑232 (MQ=255)

                                            gcgtttg‑cgtACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAAc                                                                                                                                                                                                         >  1:948589/1‑250 (MQ=255)

                                               tttg‑cgtACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGc                                                                                                                                                                                                     >  1:1022862/1‑251 (MQ=255)

                                                 tccgtgctggCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACAt                                                                                                                                                                                                       <  2:1022489/240‑1 (MQ=255)

                                                      tACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAgg                                                                                                                                                                                               >  2:444982/1‑251 (MQ=255)

                                                       ccGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAggg                                                                                                                                                                                              <  2:1222471/250‑1 (MQ=255)

                                                       aCGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAgg                                                                                                                                                                                               <  1:132412/250‑1 (MQ=255)

                                                        cGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGc                                                                                                                                                                                             <  2:575093/251‑1 (MQ=255)

                                                         ggCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACAt                                                                                                                                                                                                       >  1:1022489/1‑240 (MQ=255)

                                                         ggCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGc                                                                                                                                                                                             >  2:1161688/1‑250 (MQ=255)

                                                                 cgATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGt                                                                                                                                                                                    >  2:390517/1‑251 (MQ=255)

                                                                  gATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGtc                                                                                                                                                                                   <  1:26714/251‑1 (MQ=255)

                                                                    tGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTg                                                                                                                                                                                                               >  1:379607/1‑221 (MQ=255)

                                                                     ggCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGtctc                                                                                                                                                                                 >  1:252039/1‑250 (MQ=255)

                                                                         tGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCaatga                                                                                                                                                                                        <  1:201982/239‑1 (MQ=255)

                                                                          ggATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAg                                                                                                                                                                            <  1:391579/250‑1 (MQ=255)

                                                                             tGCCGCTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCTCCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCtt                                                                                                                                                                         <  1:1108290/250‑1 (MQ=255)

                                                                             tGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCaatgaat                                                                                                                                                                                      <  1:1058367/237‑1 (MQ=255)

                                                                             tGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCaatgaat                                                                                                                                                                                      >  2:1058367/1‑237 (MQ=255)

                                                                               ccACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCt                                                                                                                                                                      <  1:1161688/251‑1 (MQ=255)

                                                                                    ctcCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTtgt                                                                                                                                                                   <  1:429369/249‑1 (MQ=255)

                                                                                     tcCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTtgtgtg                                                                                                                                                                >  1:154364/1‑251 (MQ=255)

                                                                                      ccTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTtgtgtg                                                                                                                                                                <  2:1022862/250‑1 (MQ=255)

                                                                                      ccTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTtgtgt                                                                                                                                                                 <  1:648985/249‑1 (MQ=255)

                                                                                         gCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACa                                                                                                                                                             <  1:1130662/250‑1 (MQ=255)

                                                                                            gTTTGTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAAtt                                                                                                                                                         <  1:633187/251‑1 (MQ=255)

                                                                                                     gCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGc                                                                                                                                                 <  1:542719/250‑1 (MQ=255)

                                                                                                      cAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCAt                                                                                                                                               >  1:709685/1‑251 (MQ=255)

                                                                                                        aGCCCTGTCGGATTTCCGGCGAGTCGGCGAGGGTACTCGCCACTGGCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCTAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTCCCGCACTGATCATTAAGATCACTGTGGTACTCAATCTCATCGGACGCTTGCGTAACATGCGGTAGGGCAATGAATGTCtcatca                                                                                                                                                                             <  1:1203781/219‑1 (MQ=255)

                                                                                                        aGCCCTGACGGATTTCCGGCGACTCGACGGGGGTACGAACCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGGCATATCACTGATTTGCGAGAAGTAAAACAGATGCAACACCAATAGCACCGAGAACAGTACGGCCCTGATCATTAAGATCACTGCGGTACTCAATTTCATCCGACGTTTGCGGCACATGAGGGAGGGCAATGAATGtaccacca                                                                                                                                                                             >  2:1203781/1‑212 (MQ=255)

                                                                                                              gACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTAtctc                                                                                                                                         <  1:390517/249‑1 (MQ=255)

                                                                                                                    tttCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCtt                                                                                                                                                     >  1:140189/1‑231 (MQ=255)

                                                                                                                    tttCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCtt                                                                                                                                                     <  2:140189/231‑1 (MQ=255)

                                                                                                                                cGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCa                                                                                                                      <  2:396784/250‑1 (MQ=255)

                                                                                                                                 ggCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCaa                                                                                                                     >  1:610198/1‑250 (MQ=255)

                                                                                                                                      gggTCCGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAg                                                                                                               <  2:1071980/251‑1 (MQ=255)

                                                                                                                                      gggTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACa                                                                                                                                                             >  2:78254/1‑205 (MQ=255)

                                                                                                                                      gggTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACa                                                                                                                                                             <  1:78254/205‑1 (MQ=255)

                                                                                                                                            gcgcCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAg                                                                                                                                  <  2:298791/226‑1 (MQ=255)

                                                                                                                                                cACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAAATTAAAGGGAGgaagaaaaat                                                                                                     >  2:1105386/1‑247 (MQ=255)

                                                                                                                                                     ccAGTGCCGTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATGTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCt                                                                                                <  1:807346/251‑1 (MQ=255)

                                                                                                                                                      cAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGCGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCa                                                                                                                                                >  1:70646/1‑202 (MQ=255)

                                                                                                                                                      cAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCTTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCa                                                                                                                                                <  2:70646/202‑1 (MQ=255)

                                                                                                                                                                 tgGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGaa                                                                                     <  2:709685/250‑1 (MQ=255)

                                                                                                                                                                       gtccATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCAACGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTGCCGTGTGTGACAAATTTCTTATGCAGTATCTCTGATGAGGCGGTTAATTCAAAGGGAGGAAGAATGATTGGCTATATAGGGGAAGAGActc                                                                              <  2:412552/249‑1 (MQ=255)

                                                                                                                                                                       gCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGCGGAAGAATGATTGGCTATATAGGGGAAGGGActc                                                                              >  1:630659/1‑251 (MQ=255)

                                                                                                                                                                                  gTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTTGCTATATAGGGGAAGAGACTCTGGCAACGGa                                                                    >  1:1190136/1‑250 (MQ=255)

                                                                                                                                                                                    cATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTAGCTCTGATGAGGCGGGTAATGCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAAc                                                                 <  1:1214965/251‑1 (MQ=255)

                                                                                                                                                                                          aCTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACtgtg                                                                                                                                                                                                                                       <  1:60833/79‑1 (MQ=255)

                                                                                                                                                                                             gATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTAt                                                                                                                                            >  1:62308/1‑167 (MQ=255)

                                                                                                                                                                                             gATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTAt                                                                                                                                            <  2:62308/167‑1 (MQ=255)

                                                                                                                                                                                               tttGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGt                                                                                                                            <  1:1188515/181‑1 (MQ=255)

                                                                                                                                                                                                     agaAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGaa                                                <  1:711314/251‑1 (MQ=255)

                                                                                                                                                                                                         gTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTTAAAGGGGGTAAGAATGGTTGGGTCTATATGGGGAGAGAATCTGGCAACGGAAACTGCCAGTGCGGTATGAAGAtt                                            >  1:410888/1‑251 (MQ=255)

                                                                                                                                                                                                           aaaTCAGATGCACCACCAATAGCACCGAGACCAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAAGGAGGAAGAATGATTGGTTATATTGGGGGGGAGACCGTGGCAACCGAAACttcactagcttgttttagaattc                                          >  1:330155/1‑229 (MQ=255)

                                                                                                                                                                                                            aaTCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCg                                         <  1:298067/251‑1 (MQ=255)

                                                                                                                                                                                                               cAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCggg                                       <  2:610198/250‑1 (MQ=255)

                                                                                                                                                                                                                    gCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTGCTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTTTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCACCGGAAACTGCCAGTGCTTTATGAAGATTCCGGGGCTATg                                 <  2:1218657/251‑1 (MQ=255)

                                                                                                                                                                                                                        caccaATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCg                                                                                                                               <  2:763327/153‑1 (MQ=255)

                                                                                                                                                                                                                        caccaATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCg                                                                                                                               >  1:763327/1‑153 (MQ=255)

                                                                                                                                                                                                                        caccaATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTta                             <  1:1159560/251‑1 (MQ=255)

                                                                                                                                                                                                                         accaATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTTGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCAGGGGCTATGCTtat                            <  1:537959/251‑1 (MQ=255)

                                                                                                                                                                                                                            attAGCACCGAGACCAGTACCGCACTGCTCAGTAAGATGACTGTGTTTCTCAGTTTCATCGTGCGTTTGCGTAACATGCGGTAGGGCGATGAATGTCTCATCAGCCTCCTTGTGTGACAAAGTTATTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATGGGCTATATAGGCGACGAGACTCTGGCATCTGAAGCTGCCAGGGCTGTCTGAAGATTCCGGGTCTaggtgtatagc                         <  2:813045/249‑11 (MQ=255)

                                                                                                                                                                                                                            aaTAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGc                         <  1:442251/251‑1 (MQ=255)

                                                                                                                                                                                                                             aTAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAAAATGATTGGCTCTATAAGGGAAGAGACTCTGGCAAAGGAAACTGCCAGTGGTGGATGAAGATTCCGGGGCTATGCTTATAGCg                        >  1:567219/1‑251 (MQ=255)

                                                                                                                                                                                                                               atcACCGAGAACAGTACCGCACTTATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGa                       <  2:453814/248‑1 (MQ=255)

                                                                                                                                                                                                                                 cACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCtt                              <  2:882716/241‑1 (MQ=255)

                                                                                                                                                                                                                                 cACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTCAAAGGGAGTAAGAATGATTGGCTTTATAGGGGGAGAGACTCTGGCAACGGAAAATGCCCGGGCTGTATGAAGATTTCGGGGCTATGCtt                              >  1:882716/1‑241 (MQ=255)

                                                                                                                                                                                                                                     gagaACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTc                                                                                                                       <  1:893216/148‑1 (MQ=255)

                                                                                                                                                                                                                                     gagaACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCAt                <  1:1111551/251‑1 (MQ=255)

                                                                                                                                                                                                                                        aaCAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGc                         <  2:387398/239‑1 (MQ=255)

                                                                                                                                                                                                                                          cAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACtga           <  1:555765/251‑1 (MQ=255)

                                                                                                                                                                                                                                          cAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACtga           <  2:345039/251‑1 (MQ=255)

                                                                                                                                                                                                                                            gTACCGCACTGATCATTAAGATAACTGTGGTACTCCATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTCTCTCTGATGAGGCGGACAATTAACAGGGAGTAAAAATGAATGGCTATATATGGGGAGAGAATCTGGGAACCGGAACTGCCAGTGCTTTATAAAGAATCCCGGGCTATGCTCATAGAGATAATCATACtgatg         >  2:157845/1‑251 (MQ=255)

                                                                                                                                                                                                                                                  cACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGGGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTAGGAAGATTGCGGGGCTATGCTTATAGCGATAATCATACTGATGagagag   <  2:562166/251‑1 (MQ=255)

                                                                                                                                                                                                                                                   aCTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCCAAGGGAGGAAGAATGATTGGCTATATATGGGGAGAGACTCTTGCAACGGAAAATGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTTTAGCGATAATCATACTGATGAGAGAgg  >  1:747583/1‑251 (MQ=255)

                                                                                                                                                                                                                                                   aCTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTCAAAGGGAGTAAGAATGATTGGGTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGGGCTGTATGAAGATTCCGGGGCTATGCTTATAGAGGTAATCATACTGATGAGAGAgg  >  1:870969/1‑251 (MQ=255)

CGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTG‑CGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGG  >  NC_000913/4351343‑4351835

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