BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,307,294	A→T	7.5%	L169Q (CTG→CAG)	alsK ←	D‑allose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,307,294	0	A	T	7.5%	635.0 / 30.4	187	L169Q (CTG→CAG)	alsK	D‑allose kinase
Reads supporting (aligned to +/- strand): ref base A (88/84); new base T (6/8); total (94/93)
Fisher's exact test for biased strand distribution p-value = 5.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-03

Only 100 of 194 total aligned reads displayed.

CGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTT  >  NC_000913/4307022‑4307571

cGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAgg                                                                                                                                                                                                                                                                                      <  4:378787/275‑1 (MQ=255)

                gAAGGCGGGCCTCTCCCTCCCGCCACCGCCCGGGATCACCGCTCCGGGATCGAACAGATTAATGCTGGTGGCAATGGGCCGGGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGt                                                                                                                                                                                                                                                                       <  4:724114/274‑1 (MQ=255)

                  aGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCAAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGGTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTc                                                                                                                                                                                                                                                                     <  2:125378/274‑1 (MQ=255)

                            atCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCCAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGGTCGTACAAGCGTCGTAGCGCCATTCCAGAGCAATGGGTTGCCAAGCACCCAGAATGTCCCTAAGCGCAGCGTTGAGTGATATTCCCNNNGGGCAGCGGGCCCAGATGGACTCCCACA‑CCgc                                                                                                                                                                                                                                                      >  3:213341/1‑277 (MQ=255)

                                                      ccGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGCAAAGATCGCGCAATGGGTAATTTCGGTGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGTTTGCCACACGCGCAGTGTTTGGCCATATATCCCAGGGGGATATGACCCATTTCGCCTG‑CACACCCGTGTTAACCCCTCCACGGCGCAcc                                                                                                                                                                                                                                >  2:195213/1‑274 (MQ=255)

                                                      ccGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTTCACCCGTCCACGGCGCACCGTTc                                                                                                                                                                                                                            >  1:87391/1‑278 (MQ=255)

                                                                    aaCAGAGTAAGTCTGGGGGAAAGGGCCCGTGCCGCGGTTTCAAGCAGACTGTGGAGGAAAGGTGCGTGTGCCGCATGGACGAAAAGATCGCGAAATGGGTAATGTCGGGGCTGTTGGTCGTACCAGCGTCTTAGCGCCATGCCAGAGCAATTGGTGTCCAGGCACCCAGGATTGCCACACGCGCAGTGGTGGGTCAGATCTCCCAGGGGGATATGCCCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGCTCATCCACACTGCg                                                                                                                                                                                                                <  1:76143/276‑1 (MQ=255)

                                                                            agtccTGGTGGCAATGTCCCGTGCCTGGTTTTCAAGCCGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGGCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATccc                                                                                                                                                                                                     <  1:213341/275‑1 (MQ=255)

                                                                             aTGCTGGGGGAGATGGCCCGTGCCGCGTTTTCAAGCAGACTGTGGACGAAAGGGGCGTTTTCCGCATGGAGGAAAAGACCGGGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGGGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCa                                                                                                                                                                                                         <  4:448669/274‑1 (MQ=255)

                                                                                   gtggCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCCGTTTTTGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGTACCCGTCCACGGCGCACCGTTCATCCACACTCTGAACCCCATCCCCg                                                                                                                                                                                                   >  2:374958/1‑274 (MQ=255)

                                                                                   gtggCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGCTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGCAATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCCAACCGTTCCTCCACACTGCGAACCCCATCCCCGt                                                                                                                                                                                                  >  2:475919/1‑275 (MQ=255)

                                                                                                 tGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCAGTACCGAGACAGgccg                                                                                                                                                                                    >  2:335881/1‑275 (MQ=255)

                                                                                                                   aCTCTGGACGACAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCATCGTCTTAGCGCCATTCCCGAGCAATTGGTTACCAGGCACCCAGGACTACCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCACCTGCTCCA‑CCGTGGGCACCCTTCCACGACGCCGCGGTCAGCCACACTCCTAACCCCATCCACGTACCGAGCTAGGCCCCCAGAACCAtttgatccg                                                                                                                                                                  >  4:679623/1‑270 (MQ=255)

                                                                                                                       tGGACGAAAGGGGCGTTTTCCGCATGGCCGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTCTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCCTCCACACTGCGAACCCCATCCCCGTACCCAGACAGCCCCCCAGAACCCGTTTTTGCGTaa                                                                                                                                                               >  4:41051/1‑274 (MQ=255)

                                                                                                                        ggACGAAAGTGGGTGTTTCCGCATGGACGAAAAGGTCGCCCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATGGGTTTCCAGGCACCCAGGATTGCCACACGCGAAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAgg                                                                                                                                                             <  4:585183/275‑1 (MQ=255)

                                                                                                                           cGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCa                                                                                                                                                                                                                                                 <  2:431692/188‑1 (MQ=255)

                                                                                                                             aaaGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAg                                                                                                                                                                                                                                                                                       <  6:791543/149‑1 (MQ=255)

                                                                                                                               aGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAgg                                                                                                                                                                                                                                                                                      <  8:237859/148‑1 (MQ=255)

                                                                                                                               aGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCacac                                                                                                                                                                                                                    >  4:470052/1‑213 (MQ=255)

                                                                                                                                ggggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAggg                                                                                                                                                                                                                                                                                     <  8:89276/148‑1 (MQ=255)

                                                                                                                                 gggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAggggg                                                                                                                                                                                                                                                                                   >  6:61268/1‑149 (MQ=255)

                                                                                                                                  ggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGa                                                                                                                                                                                                                                                                                  <  5:175289/149‑1 (MQ=255)

                                                                                                                                   gCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGa                                                                                                                                                                                                                                                                                  <  7:77667/148‑1 (MQ=255)

                                                                                                                                      ttttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATg                                                                                                                                                                                                                                                                              >  5:354493/1‑149 (MQ=255)

                                                                                                                                        ttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGa                                                                                                                                                                                                                                                                             <  8:798862/148‑1 (MQ=255)

                                                                                                                                         tCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAcc                                                                                                                                                                                                                                                                           <  6:513888/149‑1 (MQ=255)

                                                                                                                                         tCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAcc                                                                                                                                                                                                                                                                           <  5:84708/149‑1 (MQ=255)

                                                                                                                                                   cGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGcc                                                                                                                                                                                                                                                                  <  8:317220/148‑1 (MQ=255)

                                                                                                                                                    gAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGTTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTg                                                                                                                                                                                                                                                                >  5:1016565/1‑149 (MQ=255)

                                                                                                                                                            cgcgCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGt                                                                                                                                                                 >  4:526850/1‑235 (MQ=255)

                                                                                                                                                             gcgcAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCgt                                                                                                                                                                                                                                                      <  5:271495/149‑1 (MQ=255)

                                                                                                                                                             gcgcAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCg                                                                                                                                                                                                                                                       <  7:149112/148‑1 (MQ=255)

                                                                                                                                                                cAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGtgtg                                                                                                                                                                                                                                                   >  6:281536/1‑149 (MQ=255)

                                                                                                                                                                    gggTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCacac                                                                                                                                                                                                                    >  2:453204/1‑176 (MQ=255)

                                                                                                                                                                     ggTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCAccc                                                                                                                                                                                                                                              >  6:185420/1‑149 (MQ=255)

                                                                                                                                                                      gTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGa                                                                                                                                                                                                                                                                             >  6:424199/1‑118 (MQ=255)

                                                                                                                                                                       tAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGatat                                                                                                                                                                                                                                                                               >  7:259359/1‑115 (MQ=255)

                                                                                                                                                                       tAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGatat                                                                                                                                                                                                                                                                               <  8:259359/115‑1 (MQ=255)

                                                                                                                                                                               gggCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCCCCCCTCCACGGCGCACCGTTCATCCCCACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGACCCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTTCCGGTTCCCGTCTCTGGAAAACTcc                                                                                                       >  4:319571/1‑273 (MQ=255)

                                                                                                                                                                               gggCTGTTGTTCGGACCAGCGTCTTAGCGCCATTCAAGAGCAATTGGTTTCCAGGCACCCAGTATGGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCaa                                                                                                      <  2:202707/275‑1 (MQ=255)

                                                                                                                                                                               gggCTGTTCTTCGTACCAGCGTCTTAGCGCCATTCCGGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCGGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCaa                                                                                                      <  2:509237/275‑1 (MQ=255)

                                                                                                                                                                                ggCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGgcg                                                                                                                                                                                                                                    >  8:449163/1‑148 (MQ=255)

                                                                                                                                                                                   tgttgtTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCa                                                                                                                                                                                                                                                 <  2:588614/132‑1 (MQ=255)

                                                                                                                                                                                    gttgttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAc                                                                                                                                                                                                                                                                            <  5:584735/105‑1 (MQ=39)

                                                                                                                                                                                          cGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCAt                                                                                                                                                                                                                          <  7:733934/148‑1 (MQ=255)

                                                                                                                                                                                                aGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCacac                                                                                                                                                                                                                    >  8:1022768/1‑148 (MQ=255)

                                                                                                                                                                                                aGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACt                                                                                                                                                                                                                   >  5:498277/1‑149 (MQ=255)

                                                                                                                                                                                                aGCGTCTTAGCGCCATGCCAGAGTAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGGAGTGTTGGGGCATATCTCCCAGGGGGATATGAGCCAGTTCGCCTGCCACA‑CAGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACt                                                                                                                                                                                                                   >  6:693975/1‑149 (MQ=255)

                                                                                                                                                                                                     cTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGa                                                                                                                                                                                                               <  7:761213/148‑1 (MQ=255)

                                                                                                                                                                                                      ttAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAc                                                                                                                                                                                                             >  6:134836/1‑149 (MQ=255)

                                                                                                                                                                                                            ccATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATc                                                                                                                                                                                                       >  6:1061733/1‑149 (MQ=255)

                                                                                                                                                                                                            ccATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGGGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTa                                                                         <  2:679623/275‑1 (MQ=255)

                                                                                                                                                                                                             cATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATc                                                                                                                                                                                                       >  7:808539/1‑148 (MQ=255)

                                                                                                                                                                                                               ttCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCAGA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCTTATCCCCGTACCGAGATAGGCCGCCCGAACCAGTTTTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAAGTCAACCGGACAATTCAGTGTATTTTAGAGCTTATCg                                                                      >  2:500140/1‑275 (MQ=255)

                                                                                                                                                                                                                 ccAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCgaga                                                                <  1:132078/279‑1 (MQ=255)

                                                                                                                                                                                                                 ccAGAACAATCGGTTTCCGGTAACCCAGGATTGCACCACGCGTAGTGTTGGGTCATATCCCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCAACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGc                                                                    <  2:84665/275‑1 (MQ=255)

                                                                                                                                                                                                                    gagCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTAc                                                                                                                                                                                                >  7:807948/1‑148 (MQ=255)

                                                                                                                                                                                                                     agCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCacac                                                                                                                                                                                                                    >  2:539097/1‑127 (MQ=255)

                                                                                                                                                                                                                        aaTTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCgaga                                                                                                                                                                                           >  5:1002951/1‑149 (MQ=255)

                                                                                                                                                                                                                         aTTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCCCCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTGTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGCAAACTCAACCGGACCATTCAGCGTCTTTTCGAGCTTATCGGCGAGATCata                                                           >  1:187399/1‑276 (MQ=255)

                                                                                                                                                                                                                          ttGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGgcgc                                                                                                                                                                                                                                   >  4:783104/1‑107 (MQ=255)

                                                                                                                                                                                                                          ttGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGgcg                                                                                                                                                                                                                                    <  1:35317/106‑1 (MQ=255)

                                                                                                                                                                                                                          ttGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGAt                                                                                                                                                                                          >  8:508613/1‑148 (MQ=255)

                                                                                                                                                                                                                           tGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATa                                                                                                                                                                                         <  8:718318/148‑1 (MQ=255)

                                                                                                                                                                                                                            ggTTTCCAGGCACCCAGGCTTGCCACACGCGCAGTGTTGGGTCTTATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAgg                                                                                                                                                                                       <  6:843649/149‑1 (MQ=255)

                                                                                                                                                                                                                                  cAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCa                                                                                                                                                                                 >  5:462316/1‑149 (MQ=255)

                                                                                                                                                                                                                                        cccAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAAcc                                                                                                                                                                            <  7:523729/148‑1 (MQ=255)

                                                                                                                                                                                                                                        cccAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCa                                                                                                                                                                           >  6:675969/1‑149 (MQ=255)

                                                                                                                                                                                                                                        cccAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGt                                                                                                                                                                 >  4:569769/1‑159 (MQ=255)

                                                                                                                                                                                                                                          cAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgt                                                                                                                                                                         >  6:235546/1‑149 (MQ=255)

                                                                                                                                                                                                                                            ggATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAgg                                                                                                                                                             <  2:569769/159‑1 (MQ=255)

                                                                                                                                                                                                                                              aTTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttgtt                                                                                                                                                                     <  5:390874/149‑1 (MQ=255)

                                                                                                                                                                                                                                              aTTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttgtt                                                                                                                                                                     <  5:670006/149‑1 (MQ=255)

                                                                                                                                                                                                                                              aTTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttgt                                                                                                                                                                      >  7:433519/1‑148 (MQ=255)

                                                                                                                                                                                                                                                    acacGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTa                                                                                                                                                                >  7:948481/1‑148 (MQ=255)

                                                                                                                                                                                                                                                     cacGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTaa                                                                                                                                                               >  8:580298/1‑148 (MQ=255)

                                                                                                                                                                                                                                                      acGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAgg                                                                                                                                                             <  6:1002951/149‑1 (MQ=255)

                                                                                                                                                                                                                                                      acGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAg                                                                                                                                                              <  7:478705/148‑1 (MQ=255)

                                                                                                                                                                                                                                                        gcgtcgtgtTCTGTCATATCTCCCAGGGGGATATGACCCAGTTCTCCTGCCTCA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACCGCGAAACCCATCCCCGTACCTAGATAGGCCGCCAGAACCAGTTGTTTCGGAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTTGCTGGAAAACTCAACCGGACAATTCAGCTTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGt                          <  1:124087/273‑1 (MQ=255)

                                                                                                                                                                                                                                                           cAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGtt                                                                                                                                                        <  6:36168/149‑1 (MQ=255)

                                                                                                                                                                                                                                                              tgtTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGtttt                                                                                                                                                      <  7:322287/148‑1 (MQ=255)

                                                                                                                                                                                                                                                                 tGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTcta                                                                                                                                                   <  8:802129/148‑1 (MQ=255)

                                                                                                                                                                                                                                                                  gggTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctact                                                                                                                                                 >  5:593067/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                  gggTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctac                                                                                                                                                  <  8:343505/148‑1 (MQ=255)

                                                                                                                                                                                                                                                                   ggTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACTTCGCGGGAAAACTCAACCGACCAATTCAGCGTATTTTCGAGCTTAACGGCCAGATCATATAAATCCCCCCCTGTTAACGGCAGGTTAGGCCTatgaaaaa                  >  4:61840/1‑268 (MQ=255)

                                                                                                                                                                                                                                                                   ggTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGCGATCATATAAATCAGCCGCTGTTAACGGCAGGTTAGGCGTTGAaataa                  >  2:230747/1‑275 (MQ=255)

                                                                                                                                                                                                                                                                      cATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgt                                                                                                                                                                         <  7:668107/121‑1 (MQ=255)

                                                                                                                                                                                                                                                                      cATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgt                                                                                                                                                                         >  8:668107/1‑121 (MQ=255)

                                                                                                                                                                                                                                                                        tatCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTc                                                                                                                                            >  7:245776/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                        tatCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTg             >  2:97526/1‑275 (MQ=255)

                                                                                                                                                                                                                                                                         atCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTccc                                                                                                                                          >  6:715949/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                         atCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTccc                                                                                                                                          >  6:966183/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                          tctcCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCa                                                                                                                                         <  5:244236/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                          tctcCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGc            >  4:128305/1‑274 (MQ=255)

                                                                                                                                                                                                                                                                              ccAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAg                                                                                                                                                                                        <  6:498731/98‑1 (MQ=255)

                                                                                                                                                                                                                                                                              ccAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGGTCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGCAATAATGGTGCggtgt       >  4:19866/1‑271 (MQ=255)

                                                                                                                                                                                                                                                                               cAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGgaga                                                                                                                                    >  5:221703/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                                aGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGgaga                                                                                                                                    >  8:103058/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                                aGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGagag                                                                                                                                   >  5:350017/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                                aGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGagag                                                                                                                                   >  5:916003/1‑149 (MQ=255)

CGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACA‑CCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTT  >  NC_000913/4307022‑4307571

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