BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	694,721	C→T	G79S (GGC→AGC)	miaB ←	isopentenyl‑adenosine A37 tRNA methylthiolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	694,721	0	C	T	99.2%	495.4 / NA	131	G79S (GGC→AGC)	miaB	isopentenyl‑adenosine A37 tRNA methylthiolase
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (63/67); minor base G (0/1); total (63/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 131 total aligned reads displayed.

CAATATTTATTGCAGCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTTTGGTCATCGACTTGCTCTT  >  NC_000913/694474‑694967

cAATATTTATTGCAGCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGAt                                                                                                                                                                                                                                                     >  1:79336/1‑251 (MQ=255)

    atttattGCAGCCTTCCCTGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATc                                                                                                                                                                                                                                                 >  2:160266/1‑251 (MQ=255)

    atttattGCAGCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATc                                                                                                                                                                                                                                                 >  1:246467/1‑251 (MQ=255)

     tttagtGCAGCCTTCCATGATGGGGACAAACGCGGTCGGCTCGTCGGCTCGCGGTTCCGGTAGACGGGCAAACTGCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGTCATAGGGGGGGCGCTGGCGAATGTGCTCGCCTTCTGGCGATGCCACGCAGCCACCGACGCTGATaagc                                                                                                                                                                                                                                                 <  2:356935/250‑3 (MQ=255)

       tattGCAGCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAg                                                                                                                                                                                                                                               <  2:192875/250‑1 (MQ=255)

        attGCAGCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGt                                                                                                                                                                                                                                             >  1:683248/1‑251 (MQ=255)

          tGCAGCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTTCGATGCCACGCAGCCACCGACGCTGATAATCAGGTc                                                                                                                                                                                                                                            <  1:360180/250‑1 (MQ=255)

            cAGCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACTGACGCTGATAAGCAGGTCTgg                                                                                                                                                                                                                                         <  1:294509/251‑1 (MQ=255)

            cAGCCTGCCATGAGGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTgg                                                                                                                                                                                                                                         <  2:683248/251‑1 (MQ=255)

              gCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGAt                                                                                                                                                                                                                                       <  1:901933/251‑1 (MQ=255)

               ccTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGAtt                                                                                                                                                                                                                                      >  1:662748/1‑251 (MQ=255)

               ccTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGAtt                                                                                                                                                                                                                                      <  2:246467/251‑1 (MQ=255)

                cTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGAttc                                                                                                                                                                                                                                     >  1:332342/1‑251 (MQ=255)

                cTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGAttc                                                                                                                                                                                                                                     >  1:12887/1‑251 (MQ=255)

                  tCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGAttctt                                                                                                                                                                                                                                   >  1:893325/1‑251 (MQ=255)

                       gatgGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCtt                                                                                                                                                                                                                              <  2:301977/251‑1 (MQ=255)

                                      ggTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGCTCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGc                                                                                                                                                                                                               >  1:754274/1‑251 (MQ=255)

                                      ggTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGc                                                                                                                                                                                                               >  2:449695/1‑251 (MQ=255)

                                           gCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGACGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGCATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTACCCGTTTCCAGCGAccc                                                                                                                                                                                                          >  2:160978/1‑251 (MQ=255)

                                           gCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAAGAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACTCAGCCACCGACGCTTATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGAccc                                                                                                                                                                                                          <  2:79336/251‑1 (MQ=255)

                                            atcaggggCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTGCCAGCGACCCa                                                                                                                                                                                                         <  1:722760/245‑1 (MQ=255)

                                                 cGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGACACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTCACCGTTTCCAGCGACCCAACTgg                                                                                                                                                                                                    >  2:713858/1‑250 (MQ=255)

                                                  ggCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGAt                                                                                                                                                                                                   <  1:449695/251‑1 (MQ=255)

                                                       gcgGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAg                                                                                                                                                                                              <  1:694617/251‑1 (MQ=255)

                                                         ggTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGAc                                                                                                                                                                                            <  1:921869/251‑1 (MQ=255)

                                                         ggTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGAc                                                                                                                                                                                            >  1:575267/1‑251 (MQ=255)

                                                          gTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACt                                                                                                                                                                                           <  2:736518/251‑1 (MQ=255)

                                                               ggTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTc                                                                                                                                                                                      <  1:13584/251‑1 (MQ=255)

                                                                   gACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCCGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTg                                                                                                                                                                                   >  2:452505/1‑250 (MQ=255)

                                                                   gACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTg                                                                                                                                                                                   >  2:361542/1‑250 (MQ=255)

                                                                      ggTCAAACTGCCCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGcc                                                                                                                                                                               <  2:323327/251‑1 (MQ=255)

                                                                        gtcgcacttccGATCTCCGGGAAGCGGATATCTACAACCGGGCGGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGAGAATCAGGTCTGGATTCTTCTCTGTTAACAGTTTCCAGCTACCCAACTGATGGAAGACTTGTTCCTGAGCCtt                                                                                                                                                                             <  2:502215/241‑1 (MQ=255)

                                                                         cAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAACATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTGCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCtt                                                                                                                                                                             <  1:284050/250‑1 (MQ=255)

                                                                             cTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTctct                                                                                                                                                                                                                               >  1:261369/1‑196 (MQ=255)

                                                                                  cGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATATCCGGCAGACGGTGCCGCGTTTGCGCCCCACCAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGCTGCCACGCAGCCACCGCCGCTGATAATCAGGCCTGGATTCTTCTCTTTTAACAGTTTCCAGCGAACCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGAtt                                                                                                                                                                   >  1:941442/1‑251 (MQ=255)

                                                                                       tCCGGCAAGCTGATATCTACAACCGGGCTGCGGCCGCAACGCACGGAGTTGATAATCTCCGGCAGACGGTGCCGCGTTTGCCGCCAAAAAATACTAACGACATACCGGGCGCGCTGGCGAATGAGCTCGCCTTCTTGCGCTGCAACGCAGCCACCGACGCTGATAATTAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACCTTTTCCTGAGCATTCTCGCGGATTgtgc                                                                                                                                                               >  2:577674/1‑247 (MQ=255)

                                                                                        ccGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCACCTGATGGAAGACTTTTTCCTGAGCCtt                                                                                                                                                                             >  1:502215/1‑235 (MQ=255)

                                                                                          gggAAGCTGCTATCTACCACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATATCCGGCAGACGGTGCAGCGTTTGCGGCCAAAACATAATATCGACATAGTGGGCGCGCTGGCGAATGGGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTc                                                                                                                                                                                                                   >  2:623937/1‑195 (MQ=255)

                                                                                          gggAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCGGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTc                                                                                                                                                                                                                   <  1:623937/195‑1 (MQ=255)

                                                                                           ggAAGCTGCTATCTACAACCGGGCTGCGGTCGCAACGCACGGCGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAAAAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGCGCCTTCTCGCggc                                                                                                                                                                     >  2:866650/1‑239 (MQ=255)

                                                                                           ggAAGCGGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGa                                                                                                                                                                     <  1:866650/240‑1 (MQ=255)

                                                                                              aGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACGGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTc                                                                                                                                                       <  1:160978/251‑1 (MQ=255)

                                                                                              aGCTGATATCTACAACCGGGCGGCAGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCTCCGACACTGATAATCAGGTCTGGATTCTTCTCGTTTAACAGTTTCCAGCGACCCAACTGATGGAGGACTTTTTCCTTAGCCTTCTCGCGGATTGAGCAGGTGTTc                                                                                                                                                       <  2:423255/251‑1 (MQ=255)

                                                                                               gCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGACCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTc                                                                                                                                                       >  2:393694/1‑250 (MQ=255)

                                                                                                       cTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCagcagcagc                                                                                                                                              >  2:94605/1‑251 (MQ=255)

                                                                                                         aCAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGcac                                                                                                                                            <  1:361542/251‑1 (MQ=255)

                                                                                                            aCCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACAGAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTGTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATc                                                                                                                                         <  2:893325/251‑1 (MQ=255)

                                                                                                             ccGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATcc                                                                                                                                        <  2:662748/251‑1 (MQ=255)

                                                                                                                    gCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGcttctt                                                                                                                                 <  1:452505/251‑1 (MQ=255)

                                                                                                                     cgggcGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAgg                                                                                                                                                            <  2:81181/219‑1 (MQ=255)

                                                                                                                     cGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAgg                                                                                                                                                            >  1:81181/1‑223 (MQ=255)

                                                                                                                          gCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGCCATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGCTGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGcc                                                                                                                            >  2:363618/1‑250 (MQ=255)

                                                                                                                          gCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTg                                                                                                                                                                                                                                          <  2:8559/140‑1 (MQ=255)

                                                                                                                          gCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGAGGCCCAACAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTTAGCAGCAGCACATCCGCTCCTTCCGccc                                                                                                                           >  2:868546/1‑250 (MQ=255)

                                                                                                                                 ccgGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGt                                                                                                                    <  2:465788/250‑1 (MQ=255)

                                                                                                                                      gTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGtt                                                                                                               <  1:393694/251‑1 (MQ=255)

                                                                                                                                         gATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGAt                                                                                                            <  1:38538/251‑1 (MQ=255)

                                                                                                                                             atcTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGcc                                                                                                        >  2:729450/1‑251 (MQ=255)

                                                                                                                                                tcCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCGCGCCGTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCGGGAGTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACGTTTTCCTGAGCCTGCGCGCTGATTGAGCAGGTGTTc                                                                                                                                                       <  2:590686/201‑1 (MQ=255)

                                                                                                                                                tcCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTCCTCGCCTTCTTGCGCTGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTc                                                                                                                                                       >  1:590686/1‑201 (MQ=255)

                                                                                                                                                    gCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCgtgggtg                                                                                                 <  1:94605/251‑1 (MQ=255)

                                                                                                                                                     cAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCgtgggtgg                                                                                                >  2:699141/1‑251 (MQ=255)

                                                                                                                                                            tGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGAATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGc                                                                                                                             >  2:769171/1‑215 (MQ=255)

                                                                                                                                                            tGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGAATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGc                                                                                                                             <  1:769171/215‑1 (MQ=255)

                                                                                                                                                                gCGTTTGCGGCACAAAAATAATATCGACAGAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGGCTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGGTCAGCAGCAGCACATCCGCTTCTTCGGCCACGTCGGTCAGTTGATGGCCGTGGGTGGCATCcagca                                                                                       <  2:87013/249‑1 (MQ=255)

                                                                                                                                                                  gTTTGAGGCCCAAAAATTATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGACGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCCGCAGCACCTCCGCTTCTTCCGCAACGTCGGTCAGTTGATAGCCGTGGGTGGCGTCcagc                                                                                        >  2:840954/1‑246 (MQ=255)

                                                                                                                                                                    ttGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCcagc                                                                                        >  2:121837/1‑244 (MQ=255)

                                                                                                                                                                     tGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGCCTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGc                                                                                >  2:309966/1‑251 (MQ=255)

                                                                                                                                                                       cGGCCCCAAAATAATATCGACATAGTGGGCGCGATGGCGAATGTGCTCGCCTCCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGCCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGcc                                                                               >  2:444717/1‑250 (MQ=255)

                                                                                                                                                                        ggCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCAt                                                                             <  1:121837/251‑1 (MQ=255)

                                                                                                                                                                           ccAAAAATAATAGCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACGGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGGTCAGCAGCAGCACATCCGCTTCTGCCGCCACGTCGGTCAGTGGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCtt                                                                          <  2:939810/251‑1 (MQ=255)

                                                                                                                                                                            cAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTc                                                                         >  2:339517/1‑251 (MQ=255)

                                                                                                                                                                               aaaTAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCAt                                                                             <  1:782440/244‑1 (MQ=255)

                                                                                                                                                                               aaaTAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCAt                                                                             >  2:782440/1‑244 (MQ=255)

                                                                                                                                                                               aaaTAATATCGACATAGCGGGCGCGCTGGCGAATGTGCTCGCATTCTTGCGATGCCACGCAGCCACCGACGCTGATACTCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCgtgggtgg                                                                                                >  1:842371/1‑225 (MQ=255)

                                                                                                                                                                               aaaTAATATCGACAGAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCGTCTCGCGGATTGAGCAGGGGTTCAGCATCAGCACATCCGCTTCTTCCGCCACGTCGGTGAGTTGATAGCCgtgggtgg                                                                                                <  2:842371/225‑1 (MQ=255)

                                                                                                                                                                                 ataataTCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATg                                                                     >  2:733226/1‑250 (MQ=255)

                                                                                                                                                                                  aaataTCGCCAGAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTGGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGAGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGaa                                                                   <  2:222736/250‑1 (MQ=255)

                                                                                                                                                                                    atatCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATc                                                                 <  1:729450/251‑1 (MQ=255)

                                                                                                                                                                                    atatCGACATAGTGGGCGCGCTGGCGAAGGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCGTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCTGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATc                                                                 <  1:840954/251‑1 (MQ=255)

                                                                                                                                                                                      atCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGt                                                               >  2:81057/1‑251 (MQ=255)

                                                                                                                                                                                              aGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCttcttc                                                                                                                                >  2:771962/1‑178 (MQ=255)

                                                                                                                                                                                              aGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCttcttc                                                                                                                                <  1:771962/178‑1 (MQ=255)

                                                                                                                                                                                                  ggCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCt                                                   <  2:168470/251‑1 (MQ=255)

                                                                                                                                                                                                             gAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCAc                                                                                                                          >  1:628395/1‑169 (MQ=255)

                                                                                                                                                                                                             gAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGGGTTCAGCAGCAGCACATCCGCTTCTTCCGCCAc                                                                                                                          <  2:628395/169‑1 (MQ=255)

                                                                                                                                                                                                               aTGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGCTTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTACTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGGGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCACCAGGt                                      >  2:355558/1‑251 (MQ=255)

                                                                                                                                                                                                               aTGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGt                                                                                                                <  1:824415/177‑1 (MQ=255)

                                                                                                                                                                                                                       gCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAatat                              <  1:362570/251‑1 (MQ=255)

                                                                                                                                                                                                                               gcgcagCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGGCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTGAATATGGAGttt                       <  2:575267/245‑1 (MQ=255)

                                                                                                                                                                                                                                 gATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGtttttt                    <  1:713858/251‑1 (MQ=255)

                                                                                                                                                                                                                                    gCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTTTGGt                 <  1:733226/251‑1 (MQ=255)

                                                                                                                                                                                                                                           aGCCACCGACGCTGATCATCAGGTCTGGAATCTTCTCTTTTTACAGTTTCCAGAGACCCAACTAATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGCGTTTTTTGGTCATCGAc          >  1:923181/1‑251 (MQ=255)

                                                                                                                                                                                                                                           aGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCGTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATc                                                                 <  2:366644/196‑1 (MQ=255)

                                                                                                                                                                                                                                            gCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGc                                                                                                                                      >  1:728164/1‑126 (MQ=255)

                                                                                                                                                                                                                                             ccACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTTTg                   <  1:863702/240‑1 (MQ=255)

                                                                                                                                                                                                                                              cACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCCGCGACACAAATGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTTTGGTCAttgccttg       >  2:429498/1‑244 (MQ=255)

                                                                                                                                                                                                                                              cACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCGGACAGCCCCAGGTTTTAATATGGAGTTTTTTGGTCATCGACTTg       <  1:868546/251‑1 (MQ=255)

                                                                                                                                                                                                                                                  gACGCTGATAATCAGTGCGGGATTCTTCTTTTTTAATCGTTTCCAGCGACCCAACTGATGGAAGCCTTTTTCCTGAGCCTCCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACACCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCAGCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTGTGGTCATCGACTTGctct   <  2:772284/251‑1 (MQ=255)

                                                                                                                                                                                                                                                   aCGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTTTGGTCATCGACTTGCTCtt  >  1:243682/1‑251 (MQ=255)

CAATATTTATTGCAGCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTGCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTTTGGTCATCGACTTGCTCTT  >  NC_000913/694474‑694967

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