BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,032,771	A→G	11.5%	E94G (GAA→GGA)	yigZ →	UPF0029 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,032,771	0	A	G	11.5%	353.9 / 3.0	115	E94G (GAA→GGA)	yigZ	UPF0029 family protein
Reads supporting (aligned to +/- strand): ref base A (22/76); new base G (15/2); total (37/78)
Fisher's exact test for biased strand distribution p-value = 3.34e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 115 total aligned reads displayed.

ATCAAAAAGAGCCGCTTCATAACGATGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTAACCGGCATTGAAGCGTTGCTGGGGCAGTGTGACGGCAAAATTATCAACAGTGATTATCAGGCATTCGTT  >  NC_000913/4032539‑4033003

aTCAAAAAGAGCCGCTTCATAACGATGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCCGCTAATGGGGGGCGGGGTCGGGGGAAAtt                                                                                                                                                                                                                                            >  2:819510/1‑238 (MQ=255)

aTCAAAAAGAGCCGCTTCATAACGATGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACAGCAGGGAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGGC‑GGGGGAAttt                                                                                                                                                                                                                                           >  2:72228/1‑237 (MQ=255)

 tCAAAAAGAGCCGCTTCATAACGATGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTAc                                                                                                                                                                                                                                          <  2:1162221/238‑1 (MQ=255)

 tCAAAAAGAGCCGCTTCATAACGATGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTAc                                                                                                                                                                                                                                          <  2:1257472/238‑1 (MQ=255)

 tCAAAAAGAGCCGCTTCATAACGATGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGGGCGGGCAGAACACCCCGATGCCCGGCTCCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTAc                                                                                                                                                                                                                                          <  2:800557/238‑1 (MQ=255)

  cAAAAAGAGCCGCTTCATAACGATGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCg                                                                                                                                                                                                                                        <  2:947285/239‑1 (MQ=255)

          gCCGCTTCATAACGATGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTgg                                                                                                                                                                                                                                  <  2:1045378/237‑1 (MQ=255)

          gCCGCTTCATAACGATGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAAACGATGGCCGCCCCGCCTATGGGGAGGGGCGGC‑GGGGGAATTACCGcgtggggt                                                                                                                                                                                                                                >  2:211038/1‑232 (MQ=255)

            cGCTTCATAACGATGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAg                                                                                                                                                                                                                               <  1:1122272/238‑1 (MQ=255)

            cGCTTCAGAACGATGTTGGCGCATGCCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAg                                                                                                                                                                                                                               <  2:1323844/238‑1 (MQ=255)

             gCTTCATAACGATGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTg                                                                                                                                                                                                                             <  2:202905/239‑1 (MQ=255)

               ttCATAACGATGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTgc                                                                                                                                                                                                                            <  2:1172302/238‑1 (MQ=255)

                     aCGATGTTGGCGCAGACCGGTGGCGTTGAGGCGGGGAAAGCGTTGGTGGAATCGGTGCGGGCAGAGCACCCCGATGCCCGGCGCCAGTGCGTGGCGTGGGTCGCGGGTGCGCCGGGTGATTCTCAACAGCTGGGTTGCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGctaata                                                                                                                                                                                                                     <  2:66014/239‑1 (MQ=255)

                         tggtGGCGAATACCGAAGGCGGTGAGGCGGCGAAAGCGTTTGTGGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGGGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTTTCTGACGAAGGGGAGGCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTAcgg                                                                                                                                                                                                                  <  2:13838/235‑1 (MQ=255)

                         tGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCATAACACCCCGATGCCCGGCTCCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTAcgg                                                                                                                                                                                                                  <  1:594993/238‑1 (MQ=255)

                         tGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTGTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTAcgg                                                                                                                                                                                                                  <  2:1362642/238‑1 (MQ=255)

                         tGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTAcgg                                                                                                                                                                                                                  <  2:592333/238‑1 (MQ=255)

                         tGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAAGCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTAcgg                                                                                                                                                                                                                  <  1:1386915/238‑1 (MQ=255)

                          gTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTAcggcg                                                                                                                                                                                                                <  2:790111/239‑1 (MQ=255)

                             ggCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCCCTGGGGTAGTGCGCTACTAcggcgg                                                                                                                                                                                                               >  1:38875/1‑237 (MQ=255)

                                gcATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑Cata                                                                                                                                                                                                          <  1:387729/238‑1 (MQ=255)

                                  aTACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTg                                                                                                                                                                                                      <  1:1038855/239‑1 (MQ=255)

                                   tACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTAcgg                                                                                                                                                                                                                  <  2:728036/228‑1 (MQ=255)

                                   tACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTAcgg                                                                                                                                                                                                                  >  1:728036/1‑228 (MQ=255)

                                    aCCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCt                                                                                                                                                                                                    <  1:42309/239‑1 (MQ=255)

                                      cGAGGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTg                                                                                                                                                                                                  <  1:942907/239‑1 (MQ=255)

                                        aTGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGAAAAACGCTAGTCGCCCCGAGCACTGGCGCCGGCGGC‑GGGGGGAATATCGCTGTGGTGGTGCGatactcctcggtgccattgcttgg                                                                                                                                                                                                   >  2:278298/1‑219 (MQ=255)

                                             ggtGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCgg                                                                                                                                                                                           <  1:413528/237‑1 (MQ=255)

                                             gTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCg                                                                                                                                                                                            <  2:605287/238‑1 (MQ=255)

                                             gTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGAACCATTGCGTGGCGTGGGTCGCTGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCg                                                                                                                                                                                            <  2:856787/238‑1 (MQ=255)

                                                gAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTgg                                                                                                                                                                                        <  2:1283011/239‑1 (MQ=255)

                                                  ggcggcGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCCCCATTGCGTGGCGGGGGTCGCGGGTGCGCCGGATGCTTCTCAACAGCTGGGTTTCTCTGACGACGGGGCGCCGGCGGGCACGGCAGGGAAACCGTCGCCCGCCCCGCATAAGGGCAGCGGCCTC‑GGGGGAAATGACGGTTTGGGTGGGCGCCACTACGGCGGCCCAA‑TGGCTTGGCGCCgggggg                                                                                                                                                                                       >  2:530423/1‑239 (MQ=255)

                                                   gcggcgAAAGCGTGTGTGGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCGGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAAGGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGt                                                                                                                                                                                      <  1:908178/238‑1 (MQ=255)

                                                             cGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAAACGATGCTCGCCCAGCTAATGGGCAGCGGCGGC‑GGGGGAAATTCCCCTGTGGTAGGGGGCTACTCCGGCGG‑CATATTTGTTTGGAACCGGGGGGttgtggaaagg                                                                                                                                                                            >  2:1174522/1‑230 (MQ=255)

                                                                 tgttgAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGGC‑GGGGAAATTACCGCTGTGGGAGTGCGCTACTACGGCGG‑CCTA‑TTGCTTGGCACCGGTGGGGTAGTGAAAGCGTATg                                                                                                                                                                       >  1:548246/1‑239 (MQ=255)

                                                                                  cAGACCACCCCGCTGCCCGGTACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCa                                                                                                                                                      <  2:53783/238‑1 (MQ=255)

                                                                                          cccGATGCACGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCCGCGCCCGCCCAAGTAAAGGGCACCGGGGGC‑GGGGGAAATACCGCTTGGGGAGTGCGCCACATCGGCGG‑GACATTTGCTGGGCCCCGG‑GGGGTAGGGAAAGCGTGCTGG‑AGAGGCGTGGAACAGGCGCTgg                                                                                                                                              >  2:859418/1‑237 (MQ=255)

                                                                                           tagATGCCCGACGCCCTGGGGTGGAGTGGGTGGGGGTTGGGCCGGATGATTCTCAGCAGCTGGGTTTACCTGACGACCGTGGGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCTGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTGCGGCGG‑CATA‑TTGCTTGGCTCCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTgcgc                                                                                                                                            <  2:790221/237‑1 (MQ=255)

                                                                                           ccGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCATGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTgcgc                                                                                                                                            <  2:932960/239‑1 (MQ=255)

                                                                                           ccGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCCGCTAATGGGGAGCGGCGCC‑GGGGGAATTACCGCTGTGGTAGTGGGCGACTACGGCGG‑CCTA‑TTGCTTGGCACCGGTGGGGTAGGGTAAGCGTGTggc‑ggcgggcgggaccaagcggcgcc                                                                                                                                             >  2:534689/1‑220 (MQ=255)

                                                                                           ccGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTAACGCTGTGGGAGTGCGCTACTACGGCGG‑CATA‑TTTCTTGGCACCGGGGGGGTAGTGAAAGCGGATGGC‑GGCGGCGGGGATGAGgcgcgggg                                                                                                                                             >  2:1399893/1‑234 (MQ=255)

                                                                                               tGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑Cata                                                                                                                                                                                                          >  1:1076612/1‑175 (MQ=255)

                                                                                               tGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCAGCTAATGGGGAGCGGCGTC‑GGGGGAATTACCGCTGTGGTAGGGCGCTACTAACGCGG‑CATA‑TTGCTTTGCCCCGGGGGGTGAGTGGAAGAGGATGGC‑GGCGGCGTGGAACAGGCGCTGCGCCCGc                                                                                                                                        >  2:1362461/1‑239 (MQ=255)

                                                                                                  ccGGCTCCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTaa                                                                                                                                     <  1:1399893/239‑1 (MQ=255)

                                                                                                  ccGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTaa                                                                                                                                     <  1:211038/239‑1 (MQ=255)

                                                                                                  ccGGCACCATTGCGTGGCGTGGGCCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGAAGGGGAGCCGGCGGTAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTaa                                                                                                                                     <  2:120286/239‑1 (MQ=255)

                                                                                                   cGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAAc                                                                                                                                    <  1:198081/239‑1 (MQ=255)

                                                                                                    ggCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGGTGCTCGCCCAGCCAATGGGCAGCGGCGGC‑GGGGGAATTACCCCTGTGGGAGGGCGCTACTACCGCCG‑CCTA‑TTGCTTGGCACCGGGGGGGGAGTGGAAACGGAAGGC‑GGCGGCGGGGATCCGGCGCCGCCGCAGcccacg                                                                                                                                   >  2:208681/1‑234 (MQ=255)

                                                                                                                tggcgtggGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAg                                                                                                                       >  1:684367/1‑239 (MQ=255)

                                                                                                                tggcgtggGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGGGGTGCGCTACTACGGCGG‑CCAA‑TTGGTTGGCAACGGGGGGGTAGGGAAAGGGGAAGGCGGGCGGCAGGAA‑CCGGGGCGCGGCCAGCCAACGAGCCAACGccca                                                                                                                       >  2:964352/1‑236 (MQ=255)

                                                                                                                 ggcgtggGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAg                                                                                                                       <  1:289182/238‑1 (MQ=255)

                                                                                                                 ggcgtggGTCGCGGGGGCGCGGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAg                                                                                                                       <  1:691953/238‑1 (MQ=255)

                                                                                                                   cgtggGTCGCGGGTGCTCCGTATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGAc                                                                                                                     <  1:606753/238‑1 (MQ=255)

                                                                                                                   cgtggGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCCGCTAATGGGCAGCGGCGGC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CCTA‑TTGCTTGGCCCCGGTGGGTTAGTGAAAGCGGTTGGG‑GGCGGGGGGAATCAGGCGCTGCGCCAGCTCAAGACCCCAccccacac                                                                                                                     >  2:918303/1‑231 (MQ=255)

                                                                                                                   cgtggGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCCATGCTCGCCCAGCTAATGGGCAGCGGCGGC‑GGGGGAATTACACCTGTGGTAGTGCGCTACTACGGCGG‑GAGA‑TTGCTTGGCACCGGGGGGGTAGAGAAAACGTATGGC‑GGGGGCGTGAAACACGCGCTGGGccccccacaagccacacgccagac                                                                                                                     >  2:1229815/1‑216 (MQ=255)

                                                                                                                     tggGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGGTGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CCTA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGc                                                                                                                   >  1:112475/1‑238 (MQ=255)

                                                                                                                          cgcgGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑a                                                                                                             <  1:1253375/238‑1 (MQ=255)

                                                                                                                          cgcgGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑a                                                                                                             <  2:224026/238‑1 (MQ=255)

                                                                                                                          cgcgGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCGCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑a                                                                                                             <  1:1072869/238‑1 (MQ=255)

                                                                                                                          cgcgGGTGCGCCGGATGATTCTCAACAGCTGGGGTTCTCTGACGGCGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑a                                                                                                             <  2:900965/238‑1 (MQ=255)

                                                                                                                          cgcgGGTGCGCCGGATGATTCGCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑a                                                                                                             <  2:262201/238‑1 (MQ=255)

                                                                                                                               gTGCGCCGGATGCTTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGa                                                                                                        <  1:1008123/238‑1 (MQ=255)

                                                                                                                                         tgatTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGc                                                                                            <  1:979099/239‑1 (MQ=255)

                                                                                                                                            ttCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAAt                                                                                         <  2:717496/239‑1 (MQ=255)

                                                                                                                                             tctcAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACAGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAAtg                                                                                        <  1:72228/239‑1 (MQ=255)

                                                                                                                                               tcAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCCGCCAATGGGCAGCGGCGGC‑GGGGGAATTACCGCCGTGGGAGTGCGCTACTACGGCGG‑CCTA‑TTGCTTGGCACCCGGGGGGTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGGGGCGCGCCCGCTAACGACCCAACGCCAGGCGCCCTT‑AAACGGATAT‑ACATTGCAAtgtg                                                                                      >  2:66113/1‑239 (MQ=255)

                                                                                                                                                 aaCAGCGTGGTTTATCGTACGAGAGGGAGCCGGCGTGAACGGTGGGTAAACCGATGCTCGCCCGGCTAATGGTCAGCGGCGTC‑GGGGGAATTACCTGTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGTC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCGAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGa                                                                                     <  2:197809/238‑1 (MQ=255)

                                                                                                                                                     gCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCAGCTAATGGGCAGCGGGGTC‑GGGGGAATTACCGCTGTGGGAGTGCGCTACTACGGCGG‑GATA‑TTGCTTGGCACCGGGGGGGTTGTGAAAGCGGATGGC‑GGGGGCGGGGATCAGGGGCTGGGCCCGCTAACAACCCCCCGCAAGACCCCCTT‑AACCGAATAT‑AATTTGCCATTTGagaatc                                                                                >  2:953294/1‑234 (MQ=255)

                                                                                                                                                      cTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCCATGTGCATatca                                                                               >  1:1209736/1‑239 (MQ=255)

                                                                                                                                                       tGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATatc                                                                                <  1:66113/237‑1 (MQ=255)

                                                                                                                                                        gggTGTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATatcatc                                                                             <  1:278298/239‑1 (MQ=255)

                                                                                                                                                           tttCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑g                                                                          <  1:83246/238‑1 (MQ=255)

                                                                                                                                                               tctGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTAAc                                                                     <  1:1304311/239‑1 (MQ=255)

                                                                                                                                                                ctGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGGAGTGCGCTACTACGGCGG‑CCTA‑TTGCTTGGCACCGGGGGGGTAGGGAAAGCGGATTGC‑GGCGGGGGGAAAAAAGGGGCGGGCCCGGCAACGACCCAAACCCAGACGCCATT‑AACCCGATAA‑ACTTTTCAATGTGAaaatacacaggtaaa                                                                      >  2:1304311/1‑227 (MQ=255)

                                                                                                                                                                         gggAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGGAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CCTA‑TTGGTTGGCACCGGGGGGGTAGTGAAAGCGTATGGG‑GGGGGGGGGGATTAGGCGCCGCCGCCGCCAACGACCCCACCCCAAACGCCAAT‑TACCGGAAATAACTTTGCAA‑GGGAAAATCATCACGTTAACCggatttaa                                                            >  2:7393/1‑233 (MQ=255)

                                                                                                                                                                           gAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTAACCGGCATTGAAg                                                          <  2:1209736/238‑1 (MQ=255)

                                                                                                                                                                            aGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTAACCGGCATTGAAGc                                                         <  1:1355493/238‑1 (MQ=255)

                                                                                                                                                                            aGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTAACCGGCATTGAAGc                                                         <  2:684367/238‑1 (MQ=255)

                                                                                                                                                                                 gcggGAACGGCAGGGAAAACGATGCTCGCCCAGCTAATGGGCAGCGGCGGC‑GGGGGAATTACCGCCGTGGTAGTGCGCTACTACGGCGG‑GCTA‑TTGCTTTGGACCGGGGGGGTAGTGGAAGCGGATGGCGGGGGGGGGGATTCAGG‑GCCGCGCCCGCCAAACACCCCCAGCCAAACGCCCTTAAACCGAA‑AA‑ACCTTGTCATGTGAATATAATCA‑GTTAAACCGCATTTTAGCGtt                                                      >  2:1107766/1‑236 (MQ=255)

                                                                                                                                                                                               gaaaCCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGGGCGCTACTACGGCGG‑CATA‑TTGCTTGGGACCGGAGGGTTAGTGAAAGCGGATGGC‑GGCGGGGGGAATCAGGGGCTGGGCCAACCAACGACCCAACGCGCGACGCCATT‑AACCGAAAAT‑AACTTGCAATGTGAATATCATAA‑ATTAACCGGCACTGAAGCGTTGCTCGGGCAGTGTGACg                                     >  2:449355/2‑239 (MQ=255)

                                                                                                                                                                                                  aCCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTAACCGGCATTGAAGCGTTGCTGGGGCAGTGTGACGGCa                                  <  1:1174522/239‑1 (MQ=255)

                                                                                                                                                                                                  aCCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTAACCGGCATTGAAGCGTTGCTGGGGCAGTGTGACGGCa                                  <  1:918303/239‑1 (MQ=255)

                                                                                                                                                                                                  aCCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTAACCGGCATTGAAGCGTTGCTGGGGCAGTGTGACGGCa                                  <  1:1151918/239‑1 (MQ=255)

                                                                                                                                                                                                    cGATGCTCGCCCAGCTAATGGGCGGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTAACCGGCATTGAAGCGTTGCTGGGGCAGTGTGACGGCaaa                                <  2:1238133/239‑1 (MQ=255)

                                                                                                                                                                                                            gCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTAACCGGCATTGAAGCGTTGCTGGGGCAGTGTGACGGCAAAATTATCaa                        <  1:7393/239‑1 (MQ=255)

                                                                                                                                                                                                            gCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTAACCGGCATTGAAGCGTTGCTGGGGCAGTGTGACGGCAAAATTATCaa                        <  1:1362461/239‑1 (MQ=255)

                                                                                                                                                                                                            gCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTAACCGGCATTGAAGCGTTGCTGGGGCAGTGTGACGGCAAAATTATCaa                        <  1:953294/239‑1 (MQ=255)

                                                                                                                                                                                                            gCCCAGCGAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTAACCGGCATTGAAGCGTTGCTGGGGCAGTGTGACGGCAAAATTATCaa                        <  1:1210577/239‑1 (MQ=255)

                                                                                                                                                                                                             cccAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGGGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAAAGCAAGACCCCCTT‑AACCGAATAT‑ACTTTGCCATGTGAATATCATCA‑GTTAACCGGCATTGACACGTTGCTGGGGCAGTTTGACGGCAAAATTAtaca                        >  2:974059/1‑235 (MQ=255)

                                                                                                                                                                                                              ccAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTAACCGGCATTGAAGCGTTGCTGGGGCAGTGTGACGGCAAAATTATCAACa                      <  1:819510/239‑1 (MQ=255)

                                                                                                                                                                                                              ccAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGGTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCCTT‑AACCGAATAT‑ACTTTGCAATGAGAATATCATCA‑GTTAACCGGCATTGAAGCGTTGCTGGGGCAGTGGGACGGCAAAATTATCAACa                      >  1:1260145/1‑239 (MQ=255)

                                                                                                                                                                                                               cgggtAATGGGCAGCGGGGTC‑GGGGAAATTTCCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTAACCGGCATTGAAGCGTTGCTGGGGCAGTGTGACGGCAAAATTATCAACAg                     <  2:1260145/235‑1 (MQ=255)

                                                                                                                                                                                                                 gCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTAACCGGCATTGAAGCGTTGCTGGGGCAGTGTGACGGCAAAATTATCAACAGTg                   >  1:106204/1‑239 (MQ=255)

                                                                                                                                                                                                                          cAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTAACCGGCATTGAAGCGTTGCTGGGGCAGTGTGACGGCAAAATTATCAACAGTGATTATCAgg          >  1:1345925/1‑239 (MQ=255)

                                                                                                                                                                                                                            gcggcgTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGCCCCAACGCACGACGCCATT‑AACCGCATAT‑ACTTTGCACTGTGAATATCATCA‑GTTCACCGGCATTGACGCGTTCCTGGGGCAGGTTGACGGCAAACTTACCAACCGTGATTATCAgg          >  2:102823/1‑237 (MQ=255)

                                                                                                                                                                                                                            gcggcgTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTAACCGGCATTGAAGCGTTGCTGGGGCAGTGTGACGGCAAAATTATCAACAGTGATTATCAGGCa        >  1:553100/1‑239 (MQ=255)

                                                                                                                                                                                                                            gcgcgctC‑GGGTGAATTACCGGTGTGGTAGTTCGCTGCTCCGGCGG‑CATA‑ATGCTTGGCACCGGTGGGTTACTGAAAGCGTACGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTAACCGGCATTGAAGCGTTGCTGGGGCAGTGTGACGGCAAAATTATCAACAGTGATTATCAGGc         <  2:1224569/232‑1 (MQ=255)

                                                                                                                                                                                                                              ggcggc‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTAACCGGCATTGAAGCGTTGCTGGGGCAGTGTGACGGCAAAATTATCAACAGTGATTATCAGGCAt       <  2:1375195/238‑1 (MQ=255)

                                                                                                                                                                                                                                  tC‑GGGGAAATTACCGCTGTGGTAGAGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTAACCGGCATTGAAGCGGTGCTGGGGCAGTGTGACGGCAAAATTATCAACAGTGATTATCAGGCATTCGtt  >  2:1126453/1‑239 (MQ=255)

                                                                                                                                                                                                                                    ctgggAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTACCCGGCATTGAAGCGTTGCTGGGGCAGTGTGACGGCAAAATTATCAACAGTGATTATCAGGCATTCGtt  <  1:638195/236‑1 (MQ=255)

ATCAAAAAGAGCCGCTTCATAACGATGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTC‑GGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGG‑CATA‑TTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGC‑GGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATT‑AACCGAATAT‑ACTTTGCAATGTGAATATCATCA‑GTTAACCGGCATTGAAGCGTTGCTGGGGCAGTGTGACGGCAAAATTATCAACAGTGATTATCAGGCATTCGTT  >  NC_000913/4032539‑4033003

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