BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	68,678	T→C	A457A (GCA→GCG)	araB ←	ribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	68,678	0	T	C	97.6%	387.0 / NA	128	A457A (GCA→GCG)	araB	ribulokinase
Reads supporting (aligned to +/- strand): ref base T (1/1); major base C (58/65); minor base A (2/1); total (61/67)
Fisher's exact test for biased strand distribution p-value = 6.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

GATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGCTCGCTGCACGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGAGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGTGCCATCACGTTATTAACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCC  >  NC_000913/68431‑68921

gATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGcc                                                                                                                                                                                                                                                     >  1:240572/1‑251 (MQ=255)

   aGCGGAGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATc                                                                                                                                                                                                                                                  <  2:706426/251‑1 (MQ=255)

      ggcgATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGCTGCCGCCCAGCGCCCTCACg                                                                                                                                                                                                                                               >  1:174677/1‑251 (MQ=255)

      ggcgATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACg                                                                                                                                                                                                                                               <  1:305464/251‑1 (MQ=255)

         gATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGtta                                                                                                                                                                                                                                            <  1:962990/251‑1 (MQ=255)

            aaaGCTGTTCAAAGCGTTGGGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCACACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCCGCACGTCGCAGCAGGCCTGAAAAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCCACACGttttt                                                                                                                                                                                                                                          >  2:652197/1‑247 (MQ=255)

            aaaGCTGTTCAAAGCGTTGAGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGAAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGCGGCAAAAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGAGCGATGCCGCCCAGCGCCATCACGttattg                                                                                                                                                                                                                                         >  2:626271/1‑250 (MQ=255)

             aaGCTGTTCAAAGCGTTGTGCCTTTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGttattga                                                                                                                                                                                                                                        <  2:768888/251‑3 (MQ=255)

                  gtgcAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCGCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAAGCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTGTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACgg                                                                                                                                                                                                                                     <  1:859672/246‑1 (MQ=255)

                  gTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCg                                                                                                                                                                                                                                   >  1:723363/1‑251 (MQ=255)

                     cAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGAt                                                                                                                                                                                                                                 <  2:558497/250‑1 (MQ=255)

                     cAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATc                                                                                                                                                                                                                                <  1:724736/251‑1 (MQ=255)

                        aGCGTTGCGCCTGTTCGCTGCGCGGTTGCAGGGTTTGCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGTCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGCCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATcccc                                                                                                                                                                                                                             <  1:530425/251‑1 (MQ=255)

                         gcgttggtccTGTTCGCTGCGCGGTTGCAGGGGTTTCGCTTCCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCTACGGCAGCAAAAATCGCCTCACCGAGCTCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGGCGCAGCAGTCCTGCATAATGACCTTGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCt                                                                                                                                                                                                                            <  1:887335/243‑1 (MQ=255)

                           gTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGCTGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCCCGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGCTCCCCTgc                                                                                                                                                                                                                          >  2:355020/1‑250 (MQ=255)

                              gtgCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGCTGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCg                                                                                                                                                                                                                       >  2:242556/1‑251 (MQ=255)

                              gtgCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATc                                                                                                                                                                                                                        <  1:291223/250‑1 (MQ=255)

                                   tgttCGCTGCGAGGTTGCAGGGTTTTCTCTAACGCACTGGCCATTTTTTGCTGGGCTGGTGGGATGTATGCGTGCACTTTCGCGGCGACGGCAGCAACAATCGCCGCACCGAGCGCACAGCAATGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGCCGCAGCAGGCCTGCATAAAGACCTGGGTTTTCCGCGCGATGCGGCCCAGCGCCATCACGTTATTGACCGCGATCCCCTGCTCggcac                                                                                                                                                                                                                   >  2:122442/4‑247 (MQ=255)

                                     cTCGCTGCACGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGAGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGTGCCATCACGTTATTAACGGCGATCCCCTGATCGGTAA‑AGc                                                                                                                                                                                                               >  1:286342/1‑251 (MQ=255)

                                       cGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCa                                                                                                                                                                                                              <  1:770721/250‑1 (MQ=255)

                                          tacgcggatGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAACCCTCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCAGCCCAGCTCCCTCCCGCTATTGCCGGCGATCCCCTGATCAGTAA‑Aacacacc                                                                                                                                                                                                          >  2:727382/9‑248 (MQ=255)

                                            cgcgttgccAGGGTTTTCTCTTCCGCACTGGCCATTTTTTGCTGGGCTGAGTGGATGTCTGCGTGCACTTTCGCGGCGACTGCAGCAAAAATGGCCGCTCCGAGCGCACAGCACTGGTCAGAGGCATCAACTTGCAGCGGGCGATTCAGCACGTCGCTGCAGGCCTGCGTAATGGCCGGGTTTTTTCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCAt                                                                                                                                                                                                        <  2:813700/243‑1 (MQ=255)

                                             gcGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATg                                                                                                                                                                                                       <  1:429245/250‑1 (MQ=255)

                                             ccGGTTGTAGGGTTTTTTCTACCTCACTGGCCATTTTTTGCTGGGATGATGGGATGTCGGCGTGCCCTTACGCGTCGACGGCAGCAAAAATCGCCGCACCGAGCGCTCAGCACTGGTAAGAGGCACCAATTTGCTGCGGGCGATGCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCTCCGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATg                                                                                                                                                                                                       <  2:514186/250‑1 (MQ=255)

                                               ggTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGGCGCAGCCGGCCTGCATAATGACCTGGTCTTTCCGCGCGATGCCCACCAGCGCCATCACGTTATTGACAGCGATCCCCTGATCGGTCA‑AGACCTCCAttat                                                                                                                                                                                                     >  2:251746/1‑248 (MQ=255)

                                                    cAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTgcgc                                                                                                                                                                                                >  1:991721/1‑251 (MQ=255)

                                                        gTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTgc                                                                                                                                                                                            <  1:369897/251‑1 (MQ=255)

                                                          tttCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTgcgc                                                                                                                                                                                          >  1:478184/1‑251 (MQ=255)

                                                                 aCCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCaaa                                                                                                                                                                                      <  2:240572/248‑1 (MQ=255)

                                                                    gaaCTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGCCGGCAGCAAAACTCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCCACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCCTCAAGTTATTGACGGCGCTCCCCTGATCGGTCA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCggg                                                                                                                                                                                >  2:124415/3‑250 (MQ=255)

                                                                    gCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGCCGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAACGGCAACAATTTGCAGCGGGCGATTCCGCACGTCGCAGCCGGCCTGCCTAATGCCCTGGTTTTTCCGCGCGCTGCCGACCAGCGCCATCACGTTATTGACGGCGCTCCCCTGTTAGATAA‑AGCCATCCATGATTGCGCGTGCGCCAAAGGCggg                                                                                                                                                                                >  2:596853/1‑250 (MQ=255)

                                                                        tGGCCATTTTTTGCTGGGCTGATGGGATGGCTGCGTGCACTTTCGAGGCGACGGCGGCACAAATCGACGCACCGAGCGCACAGCCCTGGTCAGAGGCAAAAATTTGCAGCGGGCGATTCAGCAGGTCGCAGCAGGCCTGAATCATGACCTGGGTTGTACGCGCGCTGCCGCCCAGAGCAGTCACGTTCTTGACGTCGATACCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCACAGGCGGtcgg                                                                                                                                                                             >  2:354189/1‑247 (MQ=255)

                                                                         ggCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCCCTTTCGCGGCGCCGGCAGCAACACTCGCCGCACCGCGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGACGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCGGCCCAGAGCCATCACGTTATTGACGGCGATCCACTGATCGGTAACACCA‑TCCATGATTGCGCGCGCGCCaacgtgggtggccg                                                                                                                                                                           >  2:211720/1‑239 (MQ=255)

                                                                         ggCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGgc                                                                                                                                                                             >  1:356010/1‑249 (MQ=255)

                                                                         ggCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGGAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGgc                                                                                                                                                                             >  2:797329/1‑249 (MQ=255)

                                                                           ccATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGgcagca                                                                                                                                                                         <  1:472739/251‑1 (MQ=255)

                                                                                 tttGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCaa                                                                                                                                                                    <  2:723363/250‑1 (MQ=255)

                                                                                  ttGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAAc                                                                                                                                                                  >  1:444672/1‑251 (MQ=255)

                                                                                  ttGCTGGGCTGAGGGGATGTCTGCGTGCAATTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGGGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCaaa                                                                                                                                                                   <  2:224799/250‑1 (MQ=255)

                                                                                    gctgggctgaTGGGATGTCTTCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAAcc                                                                                                                                                                 <  1:992335/250‑1 (MQ=255)

                                                                                     ctgggctgaTGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGAGGCAGCAATCAAAccgc                                                                                                                                                               >  1:868099/5‑251 (MQ=255)

                                                                                            gaTGGGAGGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAACCGCCGCACCGAGCGCACAGCACGGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAAca‑g                                                                                                                                                       <  2:868099/251‑1 (MQ=255)

                                                                                              tGGGATGTCTGCGTGCACTTTCGCAGCGCCGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCCGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGCGGCAGCCATCAAACCGCCGAAaa‑gc                                                                                                                                                      >  2:558238/1‑250 (MQ=255)

                                                                                                    gTCTGCGTGCACTTTCGCGGCGACGGCAGAAAAAATCGCCGCCCCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAg                                                                                                                                               <  1:472000/251‑1 (MQ=255)

                                                                                                      cTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACAGCCGAACA‑GCAGCGGAGCg                                                                                                                                             >  2:340776/1‑251 (MQ=255)

                                                                                                          gTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCgg                                                                                                                                         >  1:431346/1‑251 (MQ=255)

                                                                                                            gcccttgcgcgGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTa                                                                                                                                       <  1:391813/244‑1 (MQ=255)

                                                                                                              acgttCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCTCGTTATTGACGGc                                                                                                                                                                                                                                    <  2:687875/155‑1 (MQ=255)

                                                                                                                 ttCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCg                                                                                                                                    <  1:893263/249‑1 (MQ=255)

                                                                                                                 tgcggggcgACGGCAGCAAAAAGCCCCGCCCCGAGCGCCCAGCCCGGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCCCGTCGCAGCAGGCCTGCATAATGACCTGGTGTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCg                                                                                                                                    <  2:440756/244‑1 (MQ=255)

                                                                                                                 gtCGCGGCGCAGGCAGCAAAAAGCGCCGCACCGAGCGCACAGCACTGTGCAGAGGCAACAAGTTGCAGCGGGCGATTCAGCACGTTGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTAGTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCg                                                                                                                                    <  2:511731/248‑1 (MQ=255)

                                                                                                                   cgcgGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAgg                                                                                                                                 <  1:165524/250‑1 (MQ=255)

                                                                                                                    gcggcgACGGCAGCAAAAATCGCCGCCCCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATGGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGtt                                                                                                                               <  1:63141/251‑1 (MQ=255)

                                                                                                                    gcggcgACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGtt                                                                                                                               <  1:603229/251‑1 (MQ=255)

                                                                                                                     ggggggcGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGGCAGAGGCAACAATTTGCAGCGGGCGATTCAGCCCGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGGGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTa                                                                                                                              <  2:480244/245‑1 (MQ=255)

                                                                                                                     cggcgACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGCTTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTCA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCCATCAAACCGCCGAACA‑GCAGCGGGGCGTCGGGAGCgggggt                                                                                                                               >  2:912124/1‑245 (MQ=255)

                                                                                                                           cGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGAGAGGTTa                                                                                                                              >  1:942333/1‑245 (MQ=255)

                                                                                                                              cagcaAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGt                                                                                                                      >  1:117405/1‑250 (MQ=255)

                                                                                                                                 aaaaaaTCGCCGCACCGAGCGCACAGCAGTGGACAGAGGCAACAATTTGCAGGGGGCGATTCAGCACGCCGCAGCAGGCCTGCATAAGGCCCTGGTGTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGAGCCCCTGATCGGTAA‑AGCACGCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATc                                                                                                                  <  2:219912/250‑1 (MQ=255)

                                                                                                                                  aaaaaTGGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCCGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCa                                                                                                                 <  2:163609/251‑1 (MQ=255)

                                                                                                                                   aaaaTCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCAc                                                                                                                >  1:950527/1‑251 (MQ=255)

                                                                                                                                     aaTCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTGGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCAccc                                                                                                              <  1:114375/251‑1 (MQ=255)

                                                                                                                                                 gAGCGCACAGCACGGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGtt                                                                                                 <  1:857941/251‑1 (MQ=255)

                                                                                                                                                  ggcgacaaGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGAGTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCAATGATTGCGCGTGCGCCAAAGGCGGTGGAAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTg                                                                                                <  2:733598/244‑1 (MQ=255)

                                                                                                                                                       ccaGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGTTTCAGCACGGCGCAGCAGGCCCGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAg                                                                                           <  2:361901/250‑1 (MQ=255)

                                                                                                                                                       acaGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATACCCTGCTCGGTAA‑AGAACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCCGCGGAGCGTCGGGAGCGAGGTTAAGATCGGAAATCAACC‑CTTTCCGGCGTTGGTTa                                                                                            >  2:715376/1‑250 (MQ=255)

                                                                                                                                                            aCTGGTCAGAGGCAACAATTTGCAGCGGGCGCTTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGttt                                                                                      >  1:376478/1‑250 (MQ=255)

                                                                                                                                                               ggTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGtttggc                                                                                   >  1:638661/1‑247 (MQ=255)

                                                                                                                                                                gTCAGAGGCAACAATTGGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCTTTGGTTAGCGtttggcg                                                                                  <  1:548495/251‑6 (MQ=255)

                                                                                                                                                                 tCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGtttggcg                                                                                  >  1:350875/1‑245 (MQ=255)

                                                                                                                                                                         aacaaTTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑CCAGCGGAGCGTCGGTAGCGAGGTTAAGAACGGGAATCACCC‑CTTTCAGGCGTTTGTTAGCGtttggcgaaggactac                                                                         >  1:685024/1‑237 (MQ=255)

                                                                                                                                                                          acaaTTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAAGGGCCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACTTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGGGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGAGCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGtttggcgagcgacgacc                                                                        <  2:143746/251‑16 (MQ=255)

                                                                                                                                                                          acaaTTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGACTGCATAATGACCTGGTTTTTCCGCGCGATGCCCCCCAGCGCCATCACGTTATTGACGGCGATACCCTGATCGGTAA‑AGCCATCCATGATTGCGCGTGCGCCACAGGCGGTGGCCGAAATCAAACC‑CCACACAAGCATCGAAGCGTCGGTAGATGGGGTAAGATCGGTAAT‑ACCCACTTTCAGGCGGTGGTTCGCGtttggcgagcttagccc                                                                        >  2:484852/1‑236 (MQ=255)

                                                                                                                                                                           cattttGCCGCGGGCTATCCAGCACGTCGCGTCGGGCCAGCATATTTACCTGTTTTTTCCGCGCGATTCCGCCCATCGCCTGCCCTTTATTGCCGGCTATCCCCTGATCGTTAA‑AGCACCCCTTGATTGCGCGTGCGCCAAGGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGATCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTGTCAGGCGTTGGTTAGCGtttggcgagcgacgaccg                                                                       <  2:762485/248‑17 (MQ=255)

                                                                                                                                                                             aTTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGtttggcgagcgacgaccgtt                                                                     <  1:339354/251‑19 (MQ=255)

                                                                                                                                                                                  cAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAAcc                                                                >  1:299422/1‑251 (MQ=255)

                                                                                                                                                                                    gCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAAcc                                                                >  1:311519/1‑249 (MQ=255)

                                                                                                                                                                                    gCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGAGGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAg                                                              <  1:468091/251‑1 (MQ=255)

                                                                                                                                                                                     cGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTTGGCGAGAGACGACCGTTAAACCAGt                                                             >  1:62061/1‑251 (MQ=255)

                                                                                                                                                                                          gATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAATCGAGc                                                        >  1:552691/1‑251 (MQ=255)

                                                                                                                                                                                                gCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCCCGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCCCTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCaccacc                                                  <  1:150458/251‑1 (MQ=255)

                                                                                                                                                                                                gCACGCCGCAGCGGGCCTACATAATTACTGGGGTGTTCCGCGGCATGCCGCCCCGCGCCATCACTTTATTCCCGGCGTTCCCCTGATCGTTAA‑AGCACCCCATATTTGCGCTGGCGCCAAAGTGGGTGGGAGCAATAAAACCGCCGAACA‑GCAGTGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGGTAAACCAGTCGTGcacctcc                                                  <  2:706012/251‑4 (MQ=255)

                                                                                                                                                                                                 cACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGGCGAGCACCACCg                                                 <  1:195945/251‑1 (MQ=255)

                                                                                                                                                                                                   cGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGGAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGACGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGGTAGCGTTTGGAGAGGGACGACCGTTAAACACGTCGATCACCACCGGc                                               >  1:791024/1‑251 (MQ=255)

                                                                                                                                                                                                    gTCGCAGCAGGCCTGCCTAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTCA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTAGGTAGGGAGGGTAAGCTCGGTAATCACCC‑CTTTAAGGCGTTGGTTAGCGCTTGGCGAGCGACGGCCCATAAACCCGTCGAGCACCACCGGcg                                              >  2:946119/1‑250 (MQ=255)

                                                                                                                                                                                                       gcagcaGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGt                                           <  1:326062/251‑1 (MQ=255)

                                                                                                                                                                                                         agcagGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACC‑GAAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGGAATCCACC‑CTTTCCGGCGTTGGTTAGCGTTTGGCGGGCGGCGACCGTTAACCCAGACGAGCcccacccgccggggg                                         >  2:850158/1‑242 (MQ=255)

                                                                                                                                                                                                            gggCCTGCATAATGACCTGGTTGTCCCGTGCGTAGCCGCCCAAAGCCATCACGTTTTTTACGGCGATCCCCTGAACGGTAA‑AGCCCTCCAGGATTGCGCGTCCGCCAATGGCTGTTGCTGCAATTAATCCGCCGAACA‑GCAGCGGCGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTTGGCTAGCGACGACCGTTAAACCAGTCGAGCAACTCCGGCAGGTGaacc                                      <  2:809632/250‑4 (MQ=255)

                                                                                                                                                                                                                cTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGGGATCCag                                    >  1:619327/1‑249 (MQ=255)

                                                                                                                                                                                                                           ccTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCCCTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAACAATCAACACGACGGACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCCGTAATCACCC‑CTTTCAGGCGATGGGTCGCGTTTGGCGAGCGAAGACCGTTAAACCAGTCGCGCACCACCCGCCGGTGGTCCAGAGACGGttttttt                       >  2:322751/1‑250 (MQ=255)

                                                                                                                                                                                                                              ggTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCTCTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGc                     <  1:784809/250‑1 (MQ=255)

                                                                                                                                                                                                                              ggTTTTTCCCCGCGAGGCCGCCCAGCGCCAACCCGTTATGGGCGGCGATCCCCGGATCGGTAA‑AGCACCCCAGGATTGCGCGGGCGCCAAAGGCGGGGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATGTTTGGcc                    <  2:117405/251‑1 (MQ=255)

                                                                                                                                                                                                                              ggTTTGTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCGGACCGGTAA‑AGCACTCCATGATGGCGGGGGCGCCAAAGGCGGGGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGcc                    <  2:299422/251‑1 (MQ=255)

                                                                                                                                                                                                                               gTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACC‑GCAGCGGAGCGGCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGGTAAACCAGTCGAGCACCACCGGCAGGGGGTTCAGAGACGGATTTTTTGccc                   >  2:184386/1‑251 (MQ=255)

                                                                                                                                                                                                                                 ttttCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGGCGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCAAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCa                  <  1:391614/250‑1 (MQ=255)

                                                                                                                                                                                                                                  tttCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGGGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCAt                 <  2:733407/250‑1 (MQ=255)

                                                                                                                                                                                                                                       gcgcgATGCCGCCAAGCGCCATAGCGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACCCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGAGTTTTGGCCCATGCTTCg           <  2:963955/251‑1 (MQ=255)

                                                                                                                                                                                                                                        cgcgATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGAACCACCGGCAAGTGATCCAGAAACGGCTTTTTGTCCCATGCTTCgg          >  1:180146/1‑251 (MQ=255)

                                                                                                                                                                                                                                           gATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCCAGCTTCGGTcc       >  1:120781/1‑250 (MQ=255)

                                                                                                                                                                                                                                              gccgccCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAgc     <  1:850158/250‑1 (MQ=255)

GATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGCTCGCTGCACGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGAGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGTGCCATCACGTTATTAACGGCGATCCCCTGATCGGTAA‑AGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACA‑GCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCC‑CTTTCAGGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCC  >  NC_000913/68431‑68921

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