BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,468,768	C→T	W362* (TGG→TGA)	chiA ←	endochitinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,468,768	0	C	T	98.7%	272.8 / NA	75	W362* (TGG→TGA)	chiA	endochitinase
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (35/38); minor base G (0/1); total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCACCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGG‑TCACCGGTG‑TTTCCGGAGTTGGATCAAC‑CGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCT  >  NC_000913/3468529‑3469005

gCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGGACTCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCtt                                                                                                                                                                                                                                          <  1:661361/251‑1 (MQ=255)

 cGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTg                                                                                                                                                                                                                                         >  2:223573/1‑251 (MQ=255)

  gTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATACTCTCCGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGt                                                                                                                                                                                                                                        >  1:425250/1‑251 (MQ=255)

     ttGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTaa                                                                                                                                                                                                                                      <  1:284887/250‑1 (MQ=255)

           gACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTAc                                                                                                                                                                                                                               <  1:235869/251‑1 (MQ=255)

             cAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCa                                                                                                                                                                                                                             <  1:272218/251‑1 (MQ=255)

               gAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCAtt                                                                                                                                                                                                                           >  2:400571/1‑251 (MQ=255)

                      tataATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAAt                                                                                                                                                                                                                     <  2:189095/250‑1 (MQ=255)

                              gCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGt                                                                                                                                                                                                             <  1:608306/250‑1 (MQ=255)

                                   gTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAt‑tt                                                                                                                                                                                                      <  1:218915/251‑1 (MQ=255)

                                   gTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAt‑tt                                                                                                                                                                                                      >  1:22913/1‑251 (MQ=255)

                                   gTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAt‑tt                                                                                                                                                                                                      <  1:324200/251‑1 (MQ=255)

                                   gTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGAAAACCATCCGCAACAATAACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAt‑tt                                                                                                                                                                                                      <  1:324173/251‑1 (MQ=255)

                                          tctacgatgccatcacagcTGTACGCTTTGCTGGCCTGGAAACACTCCGCACCAATGACGCCGCCGTTATCCCTTTCACAGGTAGTCGGGCTACCGTACTGACTCAATTCCGCAGCGGTCGCTGTACGCTCGAGACGCCATGGGTTATTGGAGATTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTAACATTGAAAATACCATGGTAATT‑TTCGCTAcc                                                                                                                                                                                               >  2:91268/18‑251 (MQ=255)

                                             aCGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCt                                                                                                                                                                                              <  1:343404/249‑1 (MQ=255)

                                              cGATGCTGTCCGCGCTGTCCGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCtt                                                                                                                                                                                             <  2:8288/249‑1 (MQ=255)

                                              cGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCtt                                                                                                                                                                                             >  1:8288/1‑249 (MQ=255)

                                               gATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCt                                                                                                                                                                                          >  2:512475/1‑251 (MQ=255)

                                                 tGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTgg                                                                                                                                                                                        >  1:611094/1‑251 (MQ=255)

                                                    tGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCg                                                                                                                                                                                     >  1:624101/1‑251 (MQ=255)

                                                     gTCCGCGCTGTCCGCTTTGCGGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCGGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGcc                                                                                                                                                                                      <  2:611094/249‑1 (MQ=255)

                                                     gTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAAGGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGcc                                                                                                                                                                                      <  1:96660/249‑1 (MQ=255)

                                                       ccGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCt                                                                                                                                                                                   <  1:486832/250‑1 (MQ=255)

                                                          cgcTGTACGCTTTGCGGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCAc                                                                                                                                                                               <  1:400571/251‑1 (MQ=255)

                                                           gcTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACg                                                                                                                                                                              >  1:197333/1‑251 (MQ=255)

                                                              gTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCtt                                                                                                                                                                           <  2:537908/251‑1 (MQ=255)

                                                              gTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACAGTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCt                                                                                                                                                                            <  1:486254/250‑1 (MQ=255)

                                                                aCGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGc                                                                                                                                                                         >  1:123676/1‑251 (MQ=255)

                                                                aCGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGc                                                                                                                                                                         <  1:682134/251‑1 (MQ=255)

                                                                    tttGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACAtt                                                                                                                                                                                                                                 <  1:72861/192‑1 (MQ=255)

                                                                    tttGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCAAAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACAtt                                                                                                                                                                                                                                 >  2:72861/1‑192 (MQ=255)

                                                                          ggCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAAt                                                                                                                                                               <  1:425789/251‑1 (MQ=255)

                                                                           gCCTGGAACCCATCCGCAACAAGGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATc                                                                                                                                                              <  2:573002/251‑1 (MQ=255)

                                                                                gAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTa                                                                                                                                                          >  1:135494/1‑250 (MQ=255)

                                                                                 aaaCCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAAc                                                                                                                                                        <  2:729310/251‑1 (MQ=255)

                                                                                   cccagccGCAACAATGACGCCGCCGTTAGCAAGTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATGACCATGGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCTACGCTTGCAATGAGTAATCGTTAACGc                                                                                                                                                      <  1:474637/246‑1 (MQ=255)

                                                                                           caacaaTGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCt                                                                                                                                              <  1:87868/251‑1 (MQ=255)

                                                                                            aacaaTGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTg                                                                                                                                             >  2:133542/1‑251 (MQ=255)

                                                                                              caaTGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑gt                                                                                                                                          >  2:716483/1‑251 (MQ=255)

                                                                                               aaTGACGCCGCAGATATCATTTTAACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCCGCGGGAGCTGTACGCTCGAGACGCCATGGGTTATTGGGGTCTTCAGCGCTGGTATCCCGTGGGCAATCTTTAGACCCACCCCATCAGGCGTTCTTGTAATCATTACCATTGAAAAAAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCCCGACGCTTGCAATGAGTAGTCGTTACAGCTATCCGCTG‑gtg                                                                                                                                         >  2:671939/1‑251 (MQ=255)

                                                                                                aTGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑gtgt                                                                                                                                        >  1:400484/1‑251 (MQ=255)

                                                                                                         ccgTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCACGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGCACATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTc                                                                                                                              >  1:621481/1‑251 (MQ=255)

                                                                                                         ccgTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAAGTCCGCAGCGGTAGCTGTACGCTCGAGACGCCAGGGGTTATTGGAGTTTTCAGCGCTGGTGCCCCGTGGGCAATCTTTATACCCAACCCATCAGGCGTTCTTGTAAACATTTCCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTc                                                                                                                              <  2:202243/251‑1 (MQ=255)

                                                                                                           gTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATg                                                                                                                                                                     <  1:297587/212‑1 (MQ=255)

                                                                                                           gTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGAGGTAGCTGTACGCTCGAGCCGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATg                                                                                                                                                                     >  2:297587/1‑212 (MQ=255)

                                                                                                             tATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCt                                                                                                                           <  1:133542/250‑1 (MQ=255)

                                                                                                             tATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTg                                                                                                                          >  1:35687/1‑251 (MQ=255)

                                                                                                                cAATTTCACAGGTAGTCGGGTTACCGTCCTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGGGTAATCGTTAACGCTATCCGCTG‑GCGTTG‑AGGGCTCTCTGttt                                                                                                                       >  2:572096/1‑250 (MQ=255)

                                                                                                                 aaTTTCACAGGTCGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGCCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTAt                                                                                                                      >  1:585696/1‑251 (MQ=255)

                                                                                                                 aaTTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTAt                                                                                                                      >  1:579655/1‑251 (MQ=255)

                                                                                                                    ttCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTATcc                                                                                                                    <  1:223573/250‑1 (MQ=255)

                                                                                                                      cacaGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGGTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTATCCgg                                                                                                                  <  1:275695/250‑1 (MQ=255)

                                                                                                                               tCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGG‑Tc                                                                                                       >  2:185388/1‑251 (MQ=255)

                                                                                                                                cGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGG‑TCa                                                                                                      <  1:673869/251‑1 (MQ=255)

                                                                                                                                 gggTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑GTCGCTACCTTCCGGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGG‑TCAc                                                                                                     <  2:624101/251‑1 (MQ=255)

                                                                                                                                    ttACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGG‑TCACCgg                                                                                                  >  2:162553/1‑251 (MQ=255)

                                                                                                                                     tACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGG‑TCACCGgt                                                                                                 >  1:503793/1‑251 (MQ=255)

                                                                                                                                              gACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGGTAGCCGCTG‑GTGTTG‑ATGGCGCTCTGGTATCCGGAGTTGGGG‑TTACCGGTG‑TTTCCgg                                                                                        <  2:305925/250‑1 (MQ=255)

                                                                                                                                               aCTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGG‑TCACCGGTG‑TTTCCgg                                                                                        >  1:573039/1‑249 (MQ=255)

                                                                                                                                               aCTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGG‑TCACCGGTG‑TTTCCGGAg                                                                                      <  2:22913/251‑1 (MQ=255)

                                                                                                                                                 ccgagtCCGCAGCGGTAGCGATACGCTCGAGACGCCATGGGTTAGTGGAGTGTTCAGCGCTGGTTCCCCGTGGGCAATCGGTAGACCCAACCCATCAGTCGGTCTTGTAAACATTACCATTGAAAATAACAGGGTAAAG‑TTCGCTACCTTCCTGGCCGCTCCACTCTTGCAATGAGTAGTCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTAACCGGAGTTGGGG‑TCACCGGTG‑TTTCCGGAGtt                                                                                    <  2:66905/246‑1 (MQ=255)

                                                                                                                                                      tCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGG‑TCACCGGTG‑TTTCCGGAGTTGGAt                                                                                <  1:708235/250‑1 (MQ=255)

                                                                                                                                                           aGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGG‑TCACCGGTG‑TTTCCGGAGTTGGATCAAC‑Cg                                                                         >  1:39880/1‑251 (MQ=255)

                                                                                                                                                           aGCGGTAGCTGTACGAACGAGGCGCCATGGGATATTAGCGTTTTCCGCGCTGGAACCCCGTGGGCAATCTTTAGACCCAACACATCAGGCGTTCTTGTAAACATTACCATTGACAATAACCTGGTAAAT‑TTAGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAATAATCGTTAACGCTCTCCGCTG‑GTGTTG‑AGGGCTCGCTGTTAACGGGAGTTGGGG‑CCAACGGTG‑TTTCCGGAGTTGGCTAAAC‑Cg                                                                         >  2:63647/1‑251 (MQ=255)

                                                                                                                                                                         gCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGG‑TCACCGGTG‑TTTCCGGAGTTGGATCAAC‑CGGAGTCACAGCACCg                                                           <  2:503793/251‑1 (MQ=255)

                                                                                                                                                                                      tGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGG‑TCACCGGTG‑TTTCCGGAGTTGGATCAAC‑CGGAGTCACAGCACCGCCATTATTATCCg                                              >  1:37889/1‑251 (MQ=255)

                                                                                                                                                                                          ttATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGG‑TCACCGGTG‑TTTCCGGAGTTGGATCAAC‑CGGAGTCACAGCACCGCCATTATTATCCGGCtt                                          >  2:616464/1‑251 (MQ=255)

                                                                                                                                                                                            aTTGGAGTTTTCAGCGCTGGTACCACGTGGGCAATCTTTAAACACAACCCATCAGGAGGTCTTGTACACATCACCATTGAAAATAACATGGTAACT‑TTCGCTACCTTCCTGGCCGATCCACGCTTGCACTGCGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGG‑TCACCGGTG‑TTTCCGGAGTTGGCTCAAC‑CGGAGTCACAGCACCGCCATTATTATCCGGCTTa                                         >  2:72596/1‑250 (MQ=255)

                                                                                                                                                                                                 ggTTTTCAGCGTGGGTACCCGGGGGGAAATCTTAAGACCCAACCCAGCAGGCGTTCTGGAAAACATTACCATTGAAAATAATATGG‑AAATGTTCGCTGCCTTCCGGCCCGCCCCCCGCTTGAAAGGAGTAAGTGTAAAGGCTATCGGC‑G‑GTTTTGGAGGGGTTGCGGTTAGCGGGAGGGGGGG‑TC‑CCGGTGGTTTCGGGAGT‑GGAGAAACGCGTAGTCACACCACCGCCATTATTATCCGGCTTAACGGaa                                   <  1:698683/250‑1 (MQ=255)

                                                                                                                                                                                                              ggtccccggTGGGCATTTTTTAAACCCACCCCATCAGGCGTTTTTGTAAACATTCCCATTAAAAATAACATGGTAAAT‑TTCGCTCCCTTCCTGGCCGCCCCACGTTTGCAATGAGTAATCGTTAACGCTACCCGCTG‑GTGTTG‑ATGGCTCGCTGTTATCCGGAGTGGGGG‑TCCCCGGTG‑TTTCCGGAGTTGGATCACC‑CGGAGTCACACCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGt                      <  1:206797/243‑1 (MQ=255)

                                                                                                                                                                                                               gtcccccggggacaattttaacacccaccccaACAGGCTTTTTTATAAACATTACCATTAAAAAAAACATGTTAAAT‑TCCCCTCCCTCCCTGGCCCCCCCCCTCTTGCAATGAGTAATGGTTAACGCGTTCCGCGG‑GTGTTG‑TGGGCTCGCTGTTATCGGGGGTTGGGG‑TCACCGTTG‑TTTCGGGGGTTGGACCAAC‑CGGAGCCCCAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGtt                     <  2:579655/229‑1 (MQ=255)

                                                                                                                                                                                                                 cccccGGGGGCAATCTTAAGACCCACCCCTTGAGGCTTTTTTGAAAACTTTCCCATAAAAAAAATCAGGGAAATT‑TTCCCTTCCTTCTTGCCCGCTCCACGCTTGCAATGATTAAGCTTCAACGCT‑TCCCGTGTGTTTTG‑AGGGCTCGTTTTTA‑CCGGGGTTGGGGGTCCCCGTTG‑TTTCCGGAGTTGGATCACC‑CGGGGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTAc                   <  1:320474/250‑1 (MQ=255)

                                                                                                                                                                                                                            aaTCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGG‑TCACCGGTG‑TTTCCGGAGTTGGATCAAC‑CGGCGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTg        >  2:193047/1‑251 (MQ=255)

                                                                                                                                                                                                                            aaTCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGG‑TCACCGGTG‑TTTCCGGAGTTGGATCAAC‑CGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTg        >  1:534076/1‑251 (MQ=255)

                                                                                                                                                                                                                                  tAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCCTTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGG‑TCACCGGTG‑TTTCCGGAGTTGGATCAAC‑CGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCt  >  1:424838/1‑251 (MQ=255)

GCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCACCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAT‑TTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTG‑GTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGG‑TCACCGGTG‑TTTCCGGAGTTGGATCAAC‑CGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCT  >  NC_000913/3468529‑3469005

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