BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	891,829	A→G	15.6%	E75G (GAG→GGG)	nfsA →	nitroreductase NfsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	891,829	0	A	G	15.6%	214.9 / 6.1	119	E75G (GAG→GGG)	nfsA	nitroreductase NfsA
Reads supporting (aligned to +/- strand): ref base A (42/52); new base G (17/4); total (60/59)
Fisher's exact test for biased strand distribution p-value = 3.31e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

AAGAGAAAAAGATAATGACGCCAACCATTGAACTTATTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATCGGCGGCCTGCGCAATAATATTGAAGCGGTGACGAAACTGCTTAAATTACCGCAGCATGTTCTGCCGCTGT  >  NC_007779/891592‑892062

aaGAGAAAAAGATAATGACGCCAACCATTGAACTTATTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGa                                                                                                                                                                                                                                              <  1:518969/238‑1 (MQ=255)

   agaAAAAGATAATGACGCCAACCATTGAACTTATTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGt                                                                                                                                                                                                                                            >  1:1215002/1‑237 (MQ=255)

   agaAAAAGATAATGACGCCAACCATTGAACTTATTCGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCCAAAACAAGTAGAGCAAGAGGAGGAGttt                                                                                                                                                                                                                                          >  2:518969/1‑238 (MQ=255)

    gaAAAAGATAATGACGCCAACCATTGAACTTATTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTc                                                                                                                                                                                                                                          <  1:210665/238‑1 (MQ=255)

      aaaaGATAATGACGCCAACCATTGAACTTATTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑Tgg                                                                                                                                                                                                                                      <  2:717910/239‑1 (MQ=255)

             aaTGACGCCAACCATTGAACTTATTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGACCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCCAAAAAACGAAGGGCAAGGCGCGGGGGTC‑TTGGTGttgt                                                                                                                                                                                                                               >  2:334930/1‑237 (MQ=255)

              aTGACGCCAACCATTGAACTTATTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑Tgg                                                                                                                                                                                                                                      <  2:354373/231‑1 (MQ=255)

               tGACGCCAACCATTGAACTTATTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTACCGCAAGCAGCGGAGTGC‑TGGGTGGTcggg                                                                                                                                                                                                                             >  2:1231969/1‑236 (MQ=255)

                  cGCCAACCATTGAACTTATTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGc                                                                                                                                                                                                                           <  1:442706/238‑1 (MQ=255)

                          aTTTAACTTATTTGTGGCCCTCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACCCGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTaa                                                                                                                                                                                                                  <  1:1177577/239‑1 (MQ=255)

                           ttGAACTTATTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTaa                                                                                                                                                                                                                  >  1:225040/1‑238 (MQ=255)

                               ctttttttGGGGCCATGGCTCCAGTCGCCTTTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGGGGCGATTATTAACGGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGAAAAAGCTTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTACCCCATGCGGGGGAGTTC‑TGGGTCTTCTGTGCCGACTTTAACCGcc                                                                                                                                                                                                             <  2:1803995/234‑1 (MQ=255)

                               aCTTATTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAGCCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGcc                                                                                                                                                                                                             <  1:1492757/239‑1 (MQ=255)

                               aCTTATTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGcc                                                                                                                                                                                                             <  1:1402895/239‑1 (MQ=255)

                               aCTTATTGGTGGCCACCGCTCCATTCGCCATTTCACTGATGATCCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTCCCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGCCCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGcc                                                                                                                                                                                                             <  1:432297/239‑1 (MQ=255)

                                 ttATTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACTGACAAAGCGTTACGTGAAGACCTGGTGACGCTGACCGGCGGGCCAAAACACGTAGCGCAAGCGGCGGAGTTT‑CGGCGGCTGTGGGCCGACTTTAACCACCAt                                                                                                                                                                                                           >  2:1249052/1‑239 (MQ=255)

                                   aTTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGGGGTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCCttt                                                                                                                                                                                                         >  1:1168705/1‑239 (MQ=255)

                                   aTTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCCttt                                                                                                                                                                                                         >  1:536449/1‑239 (MQ=255)

                                   aTTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGGTCTGTGCCGACTTTAACCGCCAttt                                                                                                                                                                                                         >  1:1538456/1‑239 (MQ=255)

                                   aTTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGGTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCAttt                                                                                                                                                                                                         >  1:200235/1‑239 (MQ=255)

                                   aTTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACAAGTAGCGCAAGCGGAGGGGGTA‑TGGGGGGTTTGTGCCGGCTTTTAACGCCCttt                                                                                                                                                                                                         >  2:1177577/1‑239 (MQ=255)

                                         ggCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCATTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACTTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGAt                                                                                                                                                                                                   <  1:479390/239‑1 (MQ=255)

                                          gCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGGGGCGGGGGTCTTGGGGTTTCTG‑GCCGACCTTAACCGCCCATTAAAGaa                                                                                                                                                                                                   >  2:1941036/1‑237 (MQ=255)

                                          gCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAAAACCTAGCGGAACAGGCGGAGGTC‑TGGGGGGTCTGTGCCGACCTTAATCGccccttaacgta                                                                                                                                                                                                   >  2:645076/1‑228 (MQ=255)

                                             aTCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGGGGTC‑TGGGGGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGt                                                                                                                                                                                               >  1:1379947/1‑239 (MQ=255)

                                              tCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTCCCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACCCGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTc                                                                                                                                                                                              <  1:1891696/239‑1 (MQ=255)

                                                      ttCGCCATTTCACTGATGAACCCATTTCCGAAGCGTAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGGAAAAACACGTAGCGCAAGCGGCGGGGGTC‑TGGGTGTTCTGGGCCGACTTTAACCGCCATTTACAGGACTGTCCGGATGc                                                                                                                                                                                       >  2:573074/1‑238 (MQ=255)

                                                       aaaacttttCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGGGTTCCGTGCCGACTTTTACCCCCATTTACAGATCCGTCCGGAGGcg                                                                                                                                                                                      >  2:319093/7‑237 (MQ=255)

                                                                  cTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGGTT‑TGGGTGTTCTGTGCCGACTTTATCCGCCCTTTATCGATCTGTCCGGATGCTCAGCTCGGcccg                                                                                                                                                                          >  2:1626715/1‑237 (MQ=255)

                                                                   tgatgaACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCACCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGCGGTCTGTGCCGACTTTAACCGCCATTTTCAGATCTGTTCGGGTGCTCAGCTCGGCCTgg                                                                                                                                                                         >  1:1832335/1‑239 (MQ=255)

                                                                     atgaACCCATTTCCGAAGCGCAGCGTGAGGCGATTTTTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGGGGTT‑TGGGGGGTCTGGGCCGGCCTTAACCGCCCTTTTCACGTCTTTCCGGGTTCCCAGCCCGGCCTGGC‑g                                                                                                                                                                      >  2:529230/1‑239 (MQ=255)

                                                                       gaACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGGTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGa                                                                                                                                                                    >  1:158171/1‑239 (MQ=255)

                                                                          cccATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCCCTCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGCACTGGTGACGCTGACCGGCGGGCAAAAACACCCAGCGCAAGCGGCGGAGCTC‑TGGGGGGTCTGTTCCGGCATTATCCGGCATTTACACACATGGTCCGCTGCTGAGCCGCGCCTGGC‑GGgaca                                                                                                                                                                 >  2:173493/1‑239 (MQ=255)

                                                                                     aGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATGCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGCCCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGGGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCg                                                                                                                                                      <  1:2153840/239‑1 (MQ=255)

                                                                                       cgaaGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTGGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGGGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCg                                                                                                                                                      <  1:1626715/234‑1 (MQ=255)

                                                                                                 ggCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGAt                                                                                                                                           <  1:573074/238‑1 (MQ=255)

                                                                                                 ggCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCACTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCCAGCGGCGGAGTTC‑TGGGGGTTCTGTGCCGACTTTAACCGCCCTTTACAGAGCTGTCCGGATGCTCAGCTCTGCCCGGC‑GGAACAACTGTTGCCCCGGGTCGTGGATa                                                                                                                                          >  2:2086808/1‑239 (MQ=255)

                                                                                                      ttattaACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGc                                                                                                                                                                             <  1:41764/200‑1 (MQ=255)

                                                                                                      ttattaACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGaac                                                                                                                                                                  <  1:255/210‑1 (MQ=255)

                                                                                                      ttattaACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAACGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAAAAACGAACGCCAGGGGCGGATGTC‑TTGGAGTTATTGGGCGACTTTTACAACCCATTACCGATATTTCCGGAATCTCCACTCCGCCCGGCGGGAGCAAGTTTTG‑TCGGTGTCggttgaccagga                                                                                                                                     >  2:1407180/1‑228 (MQ=255)

                                                                                                               gcgcCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGACCAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGgc                                                                                                                             >  1:1620690/1‑238 (MQ=255)

                                                                                                                cgcCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGGGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCCTTTACAGAACTGTCCGGATGCTCAGCCCCGCCTGGC‑GGAACAACTGTTGCTCTGTGTCGTTTGTACGGCActgattggc                                                                                                                            >  2:1770757/1‑234 (MQ=255)

                                                                                                                cgcCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGGGTTC‑TGGGTGTTCTGGGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGgcgc                                                                                                                           >  1:523794/1‑239 (MQ=255)

                                                                                                                cgcCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTT‑TGGGTGTTCTGTGCCGACTTTAAACGCCATTTACCGATCTGTCCGGAGGCTCAGCTCGGCCTGGC‑GGGACAACTGTTGcc                                                                                                                                                        >  2:255/1‑209 (MQ=255)

                                                                                                                cgcCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATgggg                                                                                                                            >  1:1265850/1‑236 (MQ=255)

                                                                                                                cgcCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGCTC‑TGGGCGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCCCGGTGTCGTTGATACGGCAATGATGgcgc                                                                                                                           >  1:132399/1‑239 (MQ=255)

                                                                                                                cgcCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGAGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAAACGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGGGTCGTTGGTACGGCCATGATGgcgc                                                                                                                           >  1:1878741/1‑239 (MQ=255)

                                                                                                                 gcCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGGGTTCTGGGCCGACTTTAACCGCCATTTACCGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCa                                                                                                                          >  1:296650/1‑239 (MQ=255)

                                                                                                                 gcCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGAGGCGGAGGTC‑TGGGGGTTCTGTGCCGACTTTAACCCCCACTTACACATACGTCCGGATGCTCAGCTCGGCCTGGC‑GGCACAACTGTTTCTCGGTGTCGTTGATACGGGAatgatgtggcg                                                                                                                          >  2:217439/1‑234 (MQ=255)

                                                                                                                 ccccGGGCGACGCCCAGTCCCAGTTTTTTGCGGTTCGGTAGCATTATTCGCAGGCCCACAAAAGCGTTCCGTGAAGACCGGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGCCTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCa                                                                                                                          <  2:2058170/238‑1 (MQ=255)

                                                                                                                   ccGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAACCACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAGCCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAg                                                                                                                         <  1:1942141/238‑1 (MQ=255)

                                                                                                                      tggaccGTCCAGTTCCGGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCAACAAAGCGTTACGTGAAGAACGGGGGACGCGGACCGGCGGGAAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGATTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAAt                                                                                                                      <  2:1168705/233‑1 (MQ=255)

                                                                                                                             tccagttccagtTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGGGGTC‑TGGGTGGTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTTCGGATGCTCCGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGGGGCGGTGTTACCGCCATGTTGGCGCAggagacaataaa                                                                                                              >  2:101343/1‑228 (MQ=255)

                                                                                                                             tccagttccagtTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGgcg                                                                                                                            <  1:484583/225‑1 (MQ=255)

                                                                                                                                     cagttTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGGGGTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCCTGC‑GGGACAAATGTTGCTCGGGGTCGTTGCTACGGCAATGATGGCGCAGCACGCATTAATcgccgctg                                                                                                      >  2:807353/1‑237 (MQ=255)

                                                                                                                                       tttttttGAAGTGCAGTAGCTTTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGAAAAAACCCGTAGCGCAAGCGGCGGATTTT‑TGGGTGTTCTGTGCCGCCTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGa                                                                                                     <  2:280249/237‑1 (MQ=255)

                                                                                                                                       gttTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGGGTTC‑TGGGGGTTCTGTGCCGACTTTAACCGCCCTTTACAGAGCTGTCCGGGTGCTCAGCCCGGCCCGGC‑GGAACAACTTGTGCCAGGTGTCGTTGGTACGGACATGCAGGGGCAGAAAGCATTAATCGCAGCGGac                                                                                                    >  2:746215/1‑238 (MQ=255)

                                                                                                                                                  tgcagtAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGGGTTCTGTGCCGACTTTAACCGCCCTTTACAGATCTGTCCGGGTGCCCCGCCCGGGCTGGC‑GGAACAACTGTTGCTCCGGGTCGTTGATACGGCAATTATGGCGCAGAAAACATTAATCGCAGCGGG‑ACCCCTGGGGtt                                                                                        >  2:2181443/1‑239 (MQ=255)

                                                                                                                                                     agtGGCATTATTCGCATTACCGACAAAGCGTGACGTGAAGAGCTGGTGACGGTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTTTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTgg                                                                                      <  2:536449/238‑1 (MQ=255)

                                                                                                                                                        ggaattattCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGGCGCTGCCCGGCGGGCAAAAACACGTGGCGCAGGCGGGGGAGTTC‑TGGGTGTTCTGTGCCGTCTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCg                                                                                   <  2:1791818/235‑1 (MQ=255)

                                                                                                                                                        ggCATTATTCGCATTACCGACAAAGCGTTACGTGAAGACCTGGTGCCGCTGACCGGCGGGAAAAAACACGTAGGGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAAAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTggg                                                                                     <  2:523794/235‑1 (MQ=255)

                                                                                                                                                        atcATTATTCGCATTACCGACAAAGCGTTACGTGAAGAGCTGGTGCCGCTGACCGGCGGTCAAAACCACGTCGCGCAGGCGGCGGTGTTC‑TGGGTGTTCTGTGCCGCCTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTggg                                                                                     <  2:1605002/234‑1 (MQ=255)

                                                                                                                                                        aGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGGGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTTTTCTGTGCCGACTTTAACCGCCATTTTCAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGc                                                                                    <  1:1722929/237‑1 (MQ=255)

                                                                                                                                                        aGCATTATGCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGGCGCTGGCCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTggg                                                                                     <  1:746215/236‑1 (MQ=255)

                                                                                                                                                        aGCATAATTCGCATTACCGAAAAAGCGTTACGTGAAGAACTGGTGACGATGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑GGGGTTTTCTGGGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTggg                                                                                     <  2:158171/236‑1 (MQ=255)

                                                                                                                                                         gCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCggg                                                                                 >  1:1550519/1‑239 (MQ=255)

                                                                                                                                                          cATTATTCGCATTCCCAACAAAGCGTTACGTGAAGAACTGGTGACGCTGCCGGGGGGGCAAAAACAGTTAGCGCAAGCGGCGGAGTTC‑GGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCggg                                                                                 <  2:296650/238‑1 (MQ=255)

                                                                                                                                                           attattCGCATTCCCGACAAAGGGTTACGTGAAGACCTGGTGCCGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCggg                                                                                 <  1:709492/237‑1 (MQ=255)

                                                                                                                                                            ttattCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACCAGTAACGCAAGCGGGGGGGGTC‑TGGGGGGTCTGTGCCGACCTTAAAACCCATTTACAGGTCTGTCTGGAGGCCCAGCTCGGCCCGGC‑GGGCGAACCGTTGCTCGGGGTCGTTGGAACGGGCATGATTGCGGAGAATGCCTTTATCGACGGCGA‑ATCGCTGGGGTTggggggggga                                                                              >  2:1628946/1‑237 (MQ=255)

                                                                                                                                                            ttattCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGGGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCCGCGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCggggga                                                                              >  2:709492/1‑237 (MQ=255)

                                                                                                                                                            ttattCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAAAACGTAGCGCAAGCGGCGGGGGTC‑TGGGGGTTCTGTTGCGACTTTAACCCCCATTTACAGATCTGGCCCGGTGCTCAGCTCGGCCTGGC‑GGAAACACTCTTGCCCGGGGTTGTTGGTACGGGCATTATGGGCCCGCAAGCATTAAACGGAGCCGACATGGCTGGATTTGGGAgggga                                                                               >  2:1522396/1‑238 (MQ=255)

                                                                                                                                                              attCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGtata                                                                            <  1:319093/239‑1 (MQ=255)

                                                                                                                                                                 cGCATTACCGACAAAGCGTTACGTGAAGAACGGGTGACGCTGACCGGCGGGCAAAACCACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTggg                                                                                     <  1:1795394/227‑1 (MQ=255)

                                                                                                                                                                   cATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGcc                                                                                                                                                                                                             >  2:1987554/1‑107 (MQ=255)

                                                                                                                                                                    aTTACCGACAAAGCGTTACGTGAAGACCTGGTGACGCTGCCCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATcg                                                                       <  1:334930/238‑1 (MQ=255)

                                                                                                                                                                    aTTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATcg                                                                       >  1:1441139/1‑238 (MQ=255)

                                                                                                                                                                    aTTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGGTC‑TGGGGGGTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGGGTCGTTGATACGGCAATGATGGCGCAGAAAGCATTAATCGCCGCGGA‑AGTGGCGGGGGTGGGCGGGGGATTGATcgg                                                                      >  2:272512/1‑239 (MQ=255)

                                                                                                                                                                           aaaaggCGTTCCGGGAAAAATGGGGGAGGCTGGCGGGGGGGAAAAACCCCGTGGCGCAGGCGGCGGGGTTC‑TGGGTGTTTTGTGCCGCCTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATCGGCGGCCTg                                                               <  2:1215002/234‑1 (MQ=255)

                                                                                                                                                                           aCAAAGCGTTCTTAGTAGAGCTGGTGAGGCTGACCGGCGGGCAAAACCACGTACCGAAAGCGGCGGAGTTC‑TGGGTATTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATCGGCGGCCTg                                                               <  1:1359488/239‑1 (MQ=255)

                                                                                                                                                                           aCAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGGTC‑TGGGTGGTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGGACCACTGTTGCTCGGTGTCGGTGATATGGCGATGATGGCGCCGAAGGCATTAATCGCAGCGGA‑ATCGCTGGGGGTGGGCGGGGgatatacttgttcttat                                                               >  2:1795394/1‑228 (MQ=255)

                                                                                                                                                                             aaggCGTTACGTGAAGAACTGGTGACGTGGGCCGGCGGGCAAAAACTCGTGGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGTCTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATCGGCGGCCTgcg                                                             <  2:291806/236‑1 (MQ=255)

                                                                                                                                                                             aaaGCGTTGCGTGAAGACCTGGTGCCGCTGGCCGGCGGGCAAAAACCCGTAGCGCAAGCGGCGGAGTTC‑TGGTTGTTCTGTGCCGACTTTACCCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATCGGCGGCCTgcg                                                             <  1:2136187/239‑1 (MQ=255)

                                                                                                                                                                                 cgtaccGTGAAGAACTGGTGACGCTGACCGGCGGGAAAAAACACTTAGCGCAAGCGGCGGAGTTG‑TGGGTGTTCTGTGCGGATTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATCGGCGGCCTGCGCaa                                                          <  1:645076/233‑1 (MQ=255)

                                                                                                                                                                                 cgtaacggaaagaaCGGGGGAGGCTGACGGGGGGGCAAAACCACGTGGGCCAGGCGGGGGAGTTC‑GGGGTGTTCTGGGCCGCCTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATCGGCGGCCTGCGCaa                                                          <  2:1265850/229‑1 (MQ=255)

                                                                                                                                                                                  gTTACTTAAGGACCTGGTCACGCGGACGGGGGGGCAAAAACACGTAGCCCAGGCGGCGGGGTTC‑TGGGTGTTCTGTGCCGACTTGAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATCGGCGGCCTGCGCaat                                                         <  2:132399/238‑1 (MQ=255)

                                                                                                                                                                                               cgggtggcggtggcCGGCGGGCAAAAACACGTGGCGCAGGGGGCGGAGTTG‑TGGGGTTTTTGTCCCGCCTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATCGGCGGCCTGCGCAATAATATTGAAGCgg                                            <  2:1379947/236‑1 (MQ=255)

                                                                                                                                                                                                  gTGGCGCTGCCCGGCGGGCAAACACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGATTTTAGCCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGa                                                                                                                                            <  1:1602197/140‑1 (MQ=255)

                                                                                                                                                                                                       gCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATCGGCGGCCTGCGCAATAATATTGAAGCGGTGACGaa                                     >  1:98926/1‑237 (MQ=255)

                                                                                                                                                                                                          gACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTTACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGGGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGAGTGGGCGGGGTATATATCGGCGGGCCGCGCAATAAAATTTGAACGGTGACCGAACTg                                 >  2:1940475/1‑238 (MQ=255)

                                                                                                                                                                                                          gACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATCGGCGGCCTGCGCAATAATATTGAAGCGGTGACGAAACt                                  >  1:63241/1‑237 (MQ=255)

                                                                                                                                                                                                                   gCAAAAACAAGAAGCGCAAGCGGCGGAGTTC‑GGGGGGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATCGGCGGCCTGCGCAATAATATTGAAGCGGTGACGAAACTGCTTAAATTa                        <  1:1130069/238‑1 (MQ=255)

                                                                                                                                                                                                                       aaaCACGTAGCGCAAGCGGCGGAGTTC‑TGGTTTTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATCGGCGGCCTGCGCAATAATATTGAAGCGGTGACGAAACTGCTTAAATTACCgca                   <  1:101343/239‑1 (MQ=255)

                                                                                                                                                                                                                              tAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGc                                                                                    >  1:947736/1‑167 (MQ=255)

                                                                                                                                                                                                                                 cgcAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATCGGCGGCCTGCGCAATAATATTGAAGCGGAGACGAACCCTCTTAAATTACCGCACCAAGGTCt          >  2:169760/1‑238 (MQ=255)

                                                                                                                                                                                                                                    aggggcgggATTTC‑TGGGGTTTCTGTGCGGACTTAAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATCGGCGGCCTGCGCAATAATATTGAAGCGGTGACGAAACTGCTTAAATTACCGCAGCATGTTCTGcc       <  2:1550519/231‑1 (MQ=255)

                                                                                                                                                                                                                                    aaGCGGCGGTGGTC‑TGGGTTTTCTGTCCTGCCTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATCGGCGGCCTGCGCAATAATATTGAAGCGGTGACGAAACTGCTTAAATTACCGCAGCATGTTCTGcc       <  2:1257410/238‑1 (MQ=255)

                                                                                                                                                                                                                                    aaGCGGCGGCGTTC‑TGGGTGTTGTGTGCCGACTTTAGCCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATCGGCGGCCTGCGCAATAATATTGAAGCGGTGACGAAACTGCTTAAATTACCGCAGCATGTTCTGcc       <  1:2086808/238‑1 (MQ=255)

                                                                                                                                                                                                                                     aGCGGCGGGGTTC‑TGGGTGTTCTGGGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATCGGCGGCCTGCGCAATAATATTGAAGCGGTGACGAAACTGCTTAAATTACCGCAGCATGTTCTGCCGc     <  1:1770757/239‑1 (MQ=255)

                                                                                                                                                                                                                                         gcggAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATCGGCGGCCTGCGCAATAATATTGAAGCGGTGACGAAACTGCTTAAATTACCGCAGCATGTTCTGCCGCTGt  <  1:807353/238‑1 (MQ=255)

                                                                                                                                                                                                                                         gcggAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATCGGCGGCCTGCGCAATAATATTGAAGCGGTGACGAAACTGCTTAAATTACCGCAGCATGTTCTGCCGCTGt  <  1:529230/238‑1 (MQ=255)

AAGAGAAAAAGATAATGACGCCAACCATTGAACTTATTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTC‑TGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGC‑GGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGA‑ATCGCTGGGATTGGGCGGGGTATATATCGGCGGCCTGCGCAATAATATTGAAGCGGTGACGAAACTGCTTAAATTACCGCAGCATGTTCTGCCGCTGT  >  NC_007779/891592‑892062

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