BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,877,256	A→G	7.4%	L70P (CTC→CCC)	yidB ←	DUF937 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,877,256	0	A	G	7.4%	464.7 / 2.3	164	L70P (CTC→CCC)	yidB	DUF937 family protein
Reads supporting (aligned to +/- strand): ref base A (78/67); new base G (3/15); total (81/83)
Fisher's exact test for biased strand distribution p-value = 4.74e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 166 total aligned reads displayed.

CCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑CAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTG  >  NC_000913/3877018‑3877493

cccggtggAAAGGGGGGATGCAACGCGCGTTGCGTGCGCCTTTGTTGGTTAGCGGAAGAGTTCCCCTTTCAGCAGTTCCATGCCTTCGGAAAGCAGATCGTTATTTGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGGTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGa                                                                                                                                                                                                                                                  <  2:349455/233‑1 (MQ=255)

   gAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAgcg                                                                                                                                                                                                                                               <  1:151500/239‑1 (MQ=255)

           aGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCCTTTGT‑GCCGAGCGCCGATTcc                                                                                                                                                                                                                                       >  1:545919/1‑239 (MQ=255)

                 aTGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCCCGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGGGCCGATTCCAGTTg                                                                                                                                                                                                                                  >  1:830397/1‑238 (MQ=255)

                  tGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTTT‑GCCGAGCGCCGATTCCAGTTGc                                                                                                                                                                                                                                 >  1:200304/1‑238 (MQ=255)

                  gtcaacggcggTTGGGTGCTCCCTTTTTGCTTAGCGGAAGGTTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGc                                                                                                                                                                                                                                 <  1:41080/228‑1 (MQ=255)

                       cgcgcgTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTTGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCa                                                                                                                                                                                                                           <  1:465601/239‑1 (MQ=255)

                       cgcgcgTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGTCCGcc                                                                                                                                                                                                                            >  1:847636/1‑238 (MQ=255)

                         cgctTTGGGTTCTCCCTTTTTGGTTTGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTAGTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACt                                                                                                                                                                                                                         <  2:830397/239‑1 (MQ=255)

                         cgcgTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACt                                                                                                                                                                                                                         >  1:28106/1‑239 (MQ=255)

                                   gCTCCCTTTTTTCTTAGCGGAAGATTTTCCCTTTCAGCATTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCgattga                                                                                                                                                                                                               <  2:471105/239‑1 (MQ=255)

                                    cttccTTTTTGCTTAGCCGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCgattga                                                                                                                                                                                                               <  1:359686/235‑1 (MQ=255)

                                        cTTTTTGCTTAGCGGAAGTGTTTCCCTTTCAGCAGTTCCATGCCGTCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGAATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGcc                                                                                                                                                                                                          <  1:944096/239‑1 (MQ=255)

                                        cTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGcc                                                                                                                                                                                                          <  1:380743/239‑1 (MQ=255)

                                        cTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGcc                                                                                                                                                                                                          <  1:365439/239‑1 (MQ=255)

                                            ttGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGt                                                                                                                                                                                                      >  1:188813/1‑239 (MQ=255)

                                               ttgtccgtaggattttCCCGTTAAGCTTTTCCCAGCCTGCGGGACGGAGATTGTTATTGGCTTGGGGTGACACTTCACCTTGCCGTGAGAGCGCAACAATAATCTTCGGTAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGa                                                                                                                                                                                                   <  2:987340/226‑1 (MQ=255)

                                               cTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGCCACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGATCGCCACTAACCGATTGATTGCCCTGTTGa                                                                                                                                                                                                   >  1:40060/1‑239 (MQ=255)

                                                 tAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTAACCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTATTGAtt                                                                                                                                                                                                 >  1:757442/1‑239 (MQ=255)

                                                  tggggAAGGTTTTCCCTTTCAGGCTTTTCATGCCGGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCCCCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTa                                                                                                                                                                                                <  2:446350/236‑1 (MQ=255)

                                                     gggagggTTCCCCTTTCCGGCGTTCCATGCCTGCGGAAAGCCGGTCGTTATTGGCTTGTGGTGACCCTTCACCTTGGGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGGAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACt                                                                                                                                                                                              <  2:766240/232‑1 (MQ=255)

                                                     ggATGAGTTTCCCTTTCAGCAGTTCCCTGGCTGCGGAAAGCATACCGTTATTGGCCTGGGGTGACCCTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACt                                                                                                                                                                                              <  2:774810/238‑1 (MQ=255)

                                                             ttCCTTTTCAGCCGTTCCATGCCTGGGGAAAGCAGATCGTTATTGGCTTGGGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGt                                                                                                                                                                                    <  2:28049/239‑1 (MQ=255)

                                                             ttCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGTTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGt                                                                                                                                                                                    <  1:196880/239‑1 (MQ=255)

                                                             ttCCCTTTCAGCAGTTCCAATCCTGCGGAAAGCAGTTCGTTATTGGTTTGGGGTGACCCTTCACCTTGCGGTTAGAGCGCATCAATTATCTTCGGCAATTTTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGt                                                                                                                                                                                    <  2:368159/239‑1 (MQ=255)

                                                                        cagTTCCATGCCTGCGGAAAGCCGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGCCACCGCATTTGT‑GCCGAGCGCCCATTCCCGTTGCTCGCCCCTAACCGCTTGCTTGCCCTGTTGATTACTCCGCC‑AGGTTGTGGGAATg                                                                                                                                                                          >  1:739986/1‑238 (MQ=255)

                                                                           ttCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGGAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGTTTGCCCAGTTGATTTCTCAGCC‑AGGtttgaggaatagcc                                                                                                                                                                       >  1:914874/1‑226 (MQ=255)

                                                                           ttCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGGGAATGGcc                                                                                                                                                                       >  1:340822/1‑238 (MQ=255)

                                                                           ttCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGcc                                                                                                                                                                       >  1:700642/1‑238 (MQ=255)

                                                                            tcccttcctggtgaaagcgattCGTTTTTGGCTTGTGGGGACACTTCCCCTTGCGGTGAGAGCGGATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGcccc                                                                                                                                                                     <  2:481551/218‑1 (MQ=255)

                                                                                     gCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAaccgc                                                                                                                                                             >  1:973182/1‑238 (MQ=255)

                                                                                      cGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCCCTAACCGATTGATTGCCCTTTTGATTACTCAACC‑AGGTTGAGAGAAAGGCCCCCAAcccgcc                                                                                                                                                            >  2:729420/1‑238 (MQ=255)

                                                                                           aGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCCCCCCttt                                                                                                                                                       >  1:1017416/1‑238 (MQ=255)

                                                                                                aTCGTTATTGGCTTGTGGTGACACTTCACCTTTCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCACCACTTTTCAGtt                                                                                                                                                 >  1:775403/1‑239 (MQ=255)

                                                                                                          ttttgggttgcctcTTCCCTTTGGGGGGGGAGCGCTTATATAACTTTGGGAAATTGTTCT‑GCCAGTAAACGGGAAGCTGTTCTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTTATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCcagcag                                                                                                                                       <  2:140252/226‑1 (MQ=255)

                                                                                                          gCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCCAAGCCGCCACTTTGCCGTTTTTCacgcag                                                                                                                                       >  1:442825/1‑239 (MQ=255)

                                                                                                           gttttggtagttcTTCACCTGTCGGTGAGAGCGCTCAAAAAATCTGCGGCAAGTGTTCT‑GCCAGTAAACTGGAAGGTGTTCTGGTATCCACGCCAAGTTTTTGCCCGGGATCGGACACCGCATTTGG‑GCCGGCCGCCGATTCCAGTTTCTCGCCACTAACCGATTGATTGCCCTGTTGATTCCTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCTACTTTGCAGTTTTTCcagcag                                                                                                                                       <  2:573524/226‑1 (MQ=255)

                                                                                                              gtggtgACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCCCTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGGGGGAATGGCCCCTAAACCCCCCCTTTGCAGTTTTTCCAGCAGCAcc                                                                                                                                   >  2:885005/1‑239 (MQ=255)

                                                                                                              gtggtgACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCCCTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCCAACCCCCCACTTTGCCGTTTTTCcaacccacc                                                                                                                                    >  2:554657/1‑231 (MQ=255)

                                                                                                              gtggtgACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAgcagca                                                                                                                                     >  1:237044/1‑237 (MQ=255)

                                                                                                              gtggtgACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAgcagca                                                                                                                                     >  1:1025749/1‑237 (MQ=255)

                                                                                                               tggtgccctttCCCCTTGGGGTGGGAGCGCATTAATAATCTGCGGAAATTTTCCT‑GCCAGTAAACGGGAAGCTGTCCTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACCCCGCTTTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCt                                                                                                                                  <  2:1025749/230‑1 (MQ=255)

                                                                                                               tggtgACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGCGAGAATGGCCCCCAAGCCCCCCCTTTTCAGTTTTTCCCgaagca                                                                                                                                     >  2:4523/1‑236 (MQ=255)

                                                                                                               ggttgccctttccctttGGGGGGAGGGCGTCTCAATTATTTTCGGAATTTTTTCT‑GCCATTAAACTGGAAGCTGTACTGGTACCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCt                                                                                                                                  <  2:188813/224‑1 (MQ=255)

                                                                                                                    acacTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAAt                                                                                                                              >  1:136750/1‑238 (MQ=255)

                                                                                                                    acacTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGACTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCCCTTTGCAGTTTTTCCAGCAGCACCTGAAt                                                                                                                              >  1:19662/1‑238 (MQ=255)

                                                                                                                      acTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTTCAGTTTTTCCAGCAGCACCTGAAtcc                                                                                                                            >  2:1063449/1‑236 (MQ=255)

                                                                                                                       cTTCCCTTTGGGGGGAGAGCGCTTCAATTATTCTGGGAATTTGTTCT‑GCCGGTAACCGGGAAGCTGTTCTGGTACCCGCGCAAAGTGTTGGCCCGAGATAGGACCCCGCATTTGT‑GCCGAGCGCGGATTCCAGTTGCGCGCCACTAACCGATGGATTGCCCTTTGGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATg‑ccg                                                                                                                         <  2:545919/239‑1 (MQ=255)

                                                                                                                            ccTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCATCACCTGAATG‑CCGCCCt                                                                                                                     >  1:907952/1‑238 (MQ=255)

                                                                                                                                gCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGCTTGATTGCCCTCTTGATTACTCAGCC‑CGGTTGAGCGAATGGCCCCTAAGCCGCCCCTTTTCAATTTTTCCAGCACGACCtacagaccgcccctccc                                                                                                                  >  2:213702/1‑222 (MQ=255)

                                                                                                                                      agagCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCCCTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CGCCCCTGCTCCCCAAc                                                                                                           >  1:889730/1‑238 (MQ=255)

                                                                                                                                       gaggcatgcaaatattttCGGCATTTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGTATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAAcc                                                                                                          <  1:792341/222‑1 (MQ=255)

                                                                                                                                       gagCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGt                                                                                                                                                                                    <  1:866329/165‑1 (MQ=255)

                                                                                                                                       gagCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAAcc                                                                                                          <  1:623242/238‑1 (MQ=255)

                                                                                                                                       gacactataattattCTTCGGCTATT‑TTCTGGCCATTCACCTGGACGCTGTCTTTGTATCCCCGCCAAGTTTTTGCCCGAGATCGGACCCCGCATTTGG‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAAcc                                                                                                          <  2:136750/224‑1 (MQ=255)

                                                                                                                                       aacggcttCAATAACCTTCGGCAATTGTTTG‑GCCATTAAACGGGAAGCTGTACTGGTTTCCCAGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAAcc                                                                                                          <  1:956224/231‑1 (MQ=255)

                                                                                                                                           gcataaaatacttttggCAATTTTTCG‑GCAGGTAAACGGGAGCCGGTACTGGCCCCCACGCACGTTTTTTCCCCGAGATCGGACACCGCATTTGG‑GCCGAGCGCCGATTCCAGTTGCTCGCCTCTAACCGCTGGATGGCCCTGTTGATATCTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCACAACCCAAc                                                                                                      <  2:866150/223‑1 (MQ=255)

                                                                                                                                           gcATCAATAATCTCCGGGAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCGCGCCAAGTTTTTGCCCGAGATCGGACCCCGCATTTGT‑GCCGGGCGCCAATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAAc                                                                                                      <  1:213702/238‑1 (MQ=255)

                                                                                                                                                attattgtgcggaattttttCG‑GCCAGTAAACGGGAAGGTGTACTGGTACCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTaaaa                                                                                                <  2:35077/221‑1 (MQ=255)

                                                                                                                                                aatattttccGGAAATTGTTCG‑CCCAGTAAACTGGAAGCTGTACGGGTATCCCGGCCAGGTTTTTGCCCGAATCCGGACCCCGCATTTGT‑GCCGGGCGCCGATCCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCGAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTaaaa                                                                                                <  2:700642/230‑1 (MQ=255)

                                                                                                                                                  taatCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCt                                                                                                                                  >  2:297064/1‑204 (MQ=255)

                                                                                                                                                  taatCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCt                                                                                                                                  <  1:297064/204‑1 (MQ=255)

                                                                                                                                                         cGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGCGAGAATGGCCCCTAAGCCGCCACTTTGAAGGTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCAAACTTAAAATAGccgcat                                                                                       >  1:582396/1‑235 (MQ=255)

                                                                                                                                                         cGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGGTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCCACTTAAAATAGCCTgct                                                                                       >  1:1003013/1‑237 (MQ=255)

                                                                                                                                                         cGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTAGTCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACATGAATG‑CCCCCCTCGCCCCCAACCCAACTTAAAAAACCCTGaa                                                                                       >  2:205407/1‑238 (MQ=255)

                                                                                                                                                          ggCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAAGT‑CCCCCCTGCCCCCCAACCCAACTTAAAAAAGCCtaa                                                                                        >  2:940074/1‑235 (MQ=255)

                                                                                                                                                               ttGTTCG‑GCCGTTAACCTGGAAGTTGTACTGGTATCCACGCCAAGGTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGTCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCTTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTcc                                                                                 <  1:429753/239‑1 (MQ=255)

                                                                                                                                                                                      gCTGTACTGGTATCCACGCCAAG‑TTTTGCCCGGGGTCCGGCACCGCATTCGG‑GCCCGGCGCCGATTCCAGGTTCTCCGCACTAACCGATTGATTTCCCTTTTGATTAATCAGCCAAGGTGGAG‑GAAAGGGCCCTTAGCCGCCACTTTGCCGGTTTTCCCACAGCCCCTAGATGCCCCCCCTCCCCCCCAACCCACTTTAAAAAACCCTAAAATTTCCCCCCCCCCCCTTTCGAAAAgg                                                            >  2:1013027/1‑238 (MQ=255)

                                                                                                                                                                                                  tCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGcc                                                                                                                                                                                                          <  1:827018/86‑1 (MQ=255)

                                                                                                                                                                                                  tCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGcc                                                                                                                                                                                                          >  2:827018/1‑86 (MQ=255)

                                                                                                                                                                                                   ccACGCCAAGTTTTGGCCAGAGATCGGACCCCGCATTTGT‑GCCGAGCGCCGATTCACGTTGCTCGCCACAAACCGATTGATGGCCCTGTTGATTCCTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTtca                                              <  1:1052992/239‑1 (MQ=255)

                                                                                                                                                                                                    cACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGCATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACCtcat                                             >  2:906864/1‑239 (MQ=255)

                                                                                                                                                                                                    cACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGACTG‑CCGCCCGGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCCCATCGCCTTTCAGACAGGCACCGACAACTtcat                                             >  1:1042523/1‑239 (MQ=255)

                                                                                                                                                                                                    cACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAACATAGCCTGATATTTCCCCGCATCGCCTTTCAGCAAGGCACCGACAACTtca                                              >  1:572577/1‑238 (MQ=255)

                                                                                                                                                                                                    cACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCCCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTtcat                                             >  2:71704/1‑239 (MQ=255)

                                                                                                                                                                                                       gCCAAGTTTTTGCCCGAGATCGGACCCCGGATTTGT‑GCCGAGCTCCGATTCCAGTGGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑Caa                                         <  1:993227/239‑1 (MQ=255)

                                                                                                                                                                                                        ccatttttttGCCCGAGTCGGGACCCCGCTTTTGT‑GCCGGGCCTCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACCCAGCC‑AGGTTGAAAGATTGGCCCCTAAGCCGCCACTTTGCAG‑TTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑Caaaa                                       <  2:1044096/234‑1 (MQ=255)

                                                                                                                                                                                                        ccacttttttGCCCGGGACCGGCACCCGTATTTGT‑GCCGAGGGCAGATTCCAGTGTGTCGCCCCTAACTGTTTGATGGCCCGTTTGATGACCCAGCC‑AGGTTGAGAGAATTGCCCCTAAGCCGCCACTTGGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑Caaa                                        <  2:324873/234‑1 (MQ=255)

                                                                                                                                                                                                        ccAGGTTTTTCCCCGAGCGCGGACCCCGGATTTGT‑GCCGTGCGCCGATCCCAGTTGCGCGCCACTAACCGATGGATTCCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑Caaa                                        <  2:889730/239‑1 (MQ=255)

                                                                                                                                                                                                        ccAATTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑Caaa                                        <  1:345333/239‑1 (MQ=255)

                                                                                                                                                                                                         cAAGTTTTTGCCCGAGATCGGACCCCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑Caaa                                        <  1:1058350/238‑1 (MQ=255)

                                                                                                                                                                                                           aatttttGCCCGGGATCGGACCCCGGATTTGG‑GCCGGGCGCCGATTCCAGTTGCCCCCCCATAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑CAAAAAg                                     <  1:1013027/237‑1 (MQ=255)

                                                                                                                                                                                                              ttgtGCCCGAGATCGGACCCCGCATTTGG‑GCCGAGCGCCGATTCCAGTGGCTCGCCATTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑CAAAAAGc                                    <  1:255881/234‑1 (MQ=255)

                                                                                                                                                                                                                ttGCCCGGGATCGGACCCCGCATTTGG‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑CAAAAAGCCCCa                                <  1:940074/239‑1 (MQ=255)

                                                                                                                                                                                                                ttGCCCGGGATCCGACCCCGCATTTGG‑GCCGAGCGCCGATTCCTGTTGCTCGCCACTAACGGATTGATTGCCCTTTTAATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑CAAAAAGCCCCa                                <  1:476949/239‑1 (MQ=255)

                                                                                                                                                                                                                ttGCCCGGCACCGGACCCCGTTTTTGT‑GCCGGGCGCCGTTTCAGGTTGCGCCCCCCTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAAAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑CAAAAAGCCCCa                                <  2:973182/239‑1 (MQ=255)

                                                                                                                                                                                                                 tGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCCCCCTGCTCCTCAACCCAACTTAAAATAGCCGAGAATTTCCCCCCACCCCCTTTCAAAAAAGAACACAAAACTTAATAAAAAAAGCCCCa                                >  2:417869/1‑239 (MQ=255)

                                                                                                                                                                                                                     cggggccggcccccctttttGT‑GCCGGGGCCCGATTCCAGTTGCCCCCCAATAACCAATTGGTTGCCGTTTTGATTCCTCAGCA‑AGGTTGAGAGAATGCCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑CAAAAAGCCCCATGATa                           <  2:442825/222‑1 (MQ=255)

                                                                                                                                                                                                                     cGAGATCGGACACCGAATTTGT‑GCCGAGCGCCGATTCCAGTTGCCCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑CAAAAAGCCCCATGATa                           <  1:703558/239‑1 (MQ=255)

                                                                                                                                                                                                                       ggctcggccccgctttttGG‑GCCGGGCGCCTATTCCATTTGCCCGCCACTACCCGATTGATTCCCCTTTTGATTACTCAGCC‑AGGTTGAGAGAATGCCCCATAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑CAAAAAGCCCCATGATAAt                         <  2:228884/224‑1 (MQ=255)

                                                                                                                                                                                                                       agaTCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCCCATCCCCTTTCAGAAAGGCACCGACAACTTCAT‑CAAAAACCCCCATgcaaat                         >  2:871398/1‑234 (MQ=255)

                                                                                                                                                                                                                        aaTCGGACACCGTATTGGT‑GCCGAGCCCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑CAAAAAGCCCCATGATAAt                         <  1:765367/237‑1 (MQ=255)

                                                                                                                                                                                                                         ttCGGAC‑CCGCCTTTGTGGCCGGGCGCCTTTTCCAGGTGCTCGCCACTAACCGGTTTGTTGCCTTTTTGATTCACCAGCC‑AGGTTGGGAAAATGGCCCCGAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑CAAAAAGCCCCATGATAATc                        <  2:658211/237‑1 (MQ=255)

                                                                                                                                                                                                                                cccGGATTTGG‑GCCGAGCGCCGTTTCCAGTTGCCCGCCACTAACCGATTGATTGCCCTGTTGATTACCCAGCC‑AGGTTGAAAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑CAAAAAGCCCCATGATAATCACCTGTaa                <  1:424733/238‑1 (MQ=255)

                                                                                                                                                                                                                                 cccccttttggcccGGGCGCCGTTTCCAGTGGGTCGCCATTACCCGATTCATTGCCCGTGTGATTCCTCAGCC‑AGGTTGAGAGATTGGCCCCTAAGCCGCCCCTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑CAAAAAGCCCCATGATAATCACCTGTaa                <  2:747578/227‑1 (MQ=255)

                                                                                                                                                                                                                                  cGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCa                                                                                                                                                    >  1:450052/1‑107 (MQ=255)

                                                                                                                                                                                                                                  cGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑CAAAAAGCCCCATGTTAATCACCTGTaaaa              >  1:48745/1‑238 (MQ=255)

                                                                                                                                                                                                                                  cGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑CAAAAAGCCCCAGGAAAAACACCCGaaaaa              >  2:684333/1‑238 (MQ=255)

                                                                                                                                                                                                                                  cGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCCCATCGCCTTTCAGAAAGGCACCGACAACTtcat‑ca                                          >  1:452217/1‑211 (MQ=255)

                                                                                                                                                                                                                                      tttGT‑GCCGAGGCGCGATGCCAGGTGAGCCCAAATAACCGATAGATTGCCCTGTTGATTACTCAGCC‑GGGTTGACAAAATGGCCCTTAAACCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑CAAAAAGCCCCATGATAATCACCTGTAAAGCGtt          <  1:152107/239‑1 (MQ=255)

CCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGTTATTGGCTTGTGGTGACACTTCACCTTGCGGTGAGAGCGCATCAATAATCTTCGGCAATTGTTCT‑GCCAGTAAACTGGAAGCTGTACTGGTATCCACGCCAAGTTTTTGCCCGAGATCGGACACCGCATTTGT‑GCCGAGCGCCGATTCCAGTTGCTCGCCACTAACCGATTGATTGCCCTGTTGATTACTCAGCC‑AGGTTGAGAGAATGGCCCCTAAGCCGCCACTTTGCAGTTTTTCCAGCAGCACCTGAATG‑CCGCCCTGCTCCTCAACCCAACTTAAAATAGCCTGATATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCAT‑CAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTG  >  NC_000913/3877018‑3877493

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