BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,603,724	A→C	5.4%	T236P (ACC→CCC)	lsrC →	autoinducer 2 import system permease protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,603,724	0	A	C	5.4%	566.1 / 18.2	170	T236P (ACC→CCC)	lsrC	autoinducer 2 import system permease protein
Reads supporting (aligned to +/- strand): ref base A (64/95); new base C (3/7); total (68/102)
Fisher's exact test for biased strand distribution p-value = 7.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.62e-01

Only 100 of 173 total aligned reads displayed.

GCCGAACTGAAACAGCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGCAAAAATATGCCCGCTTTATGACGCCAC  >  NC_000913/1603454‑1603997

               ctctCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTCTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTGGCTTCGCAGATTGGTTTTATCCCCAACCAGACAGGTACCGGGCTGgaga                                                                                                                                                                                                                                                                 >  1:15533/1‑274 (MQ=255)

                    cGCCCCGCTTCTGCTTTGCGTTTCAGCAATTGGTTGGTTGACTATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAg                                                                                                                                                                                                                                                           <  2:220711/275‑1 (MQ=255)

                                      cGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCCGATTGGTTTTATCCCCAACCAGACCGGGACCGGGCTGGAGATGAAGACAATTGAAGCCTGCGTga                                                                                                                                                                                                                                         >  2:340459/1‑274 (MQ=255)

                                                        gttgACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGACTTGGGTTTGGTTCGCGATTTGGTTTTATCCACAACCAGCCCGGTACCGGGCGGGGGGTTAAAGCCAANNNACCCTTGGTGCTgggggggctgttgtttc                                                                                                                                                                                                                       >  3:195572/1‑264 (MQ=255)

                                                          tgACGATCATTTGGGTGGCATTTCTGGCCTGGCTGCTGGCAAAGACTGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTc                                                                                                                                                                                                                     <  2:182197/275‑1 (MQ=255)

                                                             cGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGGATTTTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAAGTGCAGCCTTCGTGCTGGGCGGCATTAGTTTGATCgg                                                                                                                                                                                                                   >  2:439834/1‑274 (MQ=255)

                                                                       ggtggCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGt                                                                                                                                                                                                         <  2:202063/274‑1 (MQ=255)

                                                                                       tGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGt                                                                                                                                                                                        >  2:170238/1‑275 (MQ=255)

                                                                                         gctgctGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCCGATTGGTTTTATCCCCAACCAGGCCGGTACCGGGCTGGGGGTGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTTCCCGGTTGTTCCGGTGCGATCATTGGGGCGGtta                                                                                                                                                                                      >  4:92536/1‑273 (MQ=255)

                                                                                         gctgctGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGAcc                                                                                                                                                                                                                                                                                 >  2:747236/1‑184 (MQ=255)

                                                                                                        ggCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTc                                                                                                                                                                        <  4:170238/274‑1 (MQ=255)

                                                                                                              tGGACGCAGGTTTTAGGCAACGGGCGATAATTTACAGGGGGCGCGTAAACTGGGCGTTCGTACTGAAGCCAGTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACg                                                                                                                                                                  <  4:190937/274‑1 (MQ=255)

                                                                                                              gggcgtcggTTTTTATGCAACGGGCGTTAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATGGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACg                                                                                                                                                                  <  4:508220/266‑1 (MQ=255)

                                                                                                                  cgaaGTTGTTATGCCACGGGCGCTAATTTACAGGGCGCTCGTAAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGGAATTGCAGGCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGAt                                                                                                                                                             <  2:532183/272‑1 (MQ=255)

                                                                                                                       tttttATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTACCTGACGCAGATCGa                                                                                                                                                          <  1:199473/273‑1 (MQ=255)

                                                                                                                           tATGCCCCGGGCGATAAATTACATGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATGCGCATTTGTGCATTTTCTTTGAACGGCTGCATGGCGGCACTGGCGGGAATTTTGTTTGCTTCGCAGATTGGTTTTCTCCCCATCCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGt                                                                                                                                                    <  1:51126/275‑1 (MQ=255)

                                                                                                                           tATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGt                                                                                                                                                    <  2:92536/275‑1 (MQ=255)

                                                                                                                              gcaaaggggGATAATTTACAGGGCGCTCGGTAACTGGGCGTTCGGACTGGAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTAc                                                                                                                                                  <  4:357144/269‑1 (MQ=255)

                                                                                                                               cgaCGGGCGATAAGTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGGACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTAc                                                                                                                                                  <  1:207065/271‑1 (MQ=255)

                                                                                                                                   gggCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCg                                                                                                                                                                                                                                                                          <  7:241072/149‑1 (MQ=255)

                                                                                                                                   gggCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCg                                                                                                                                                                                                                                                                          <  6:733401/149‑1 (MQ=255)

                                                                                                                                       gATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCTCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCt                                                                                                                                                                                                                                                                      >  7:989629/1‑149 (MQ=255)

                                                                                                                                       gATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCt                                                                                                                                                                                                                                                                      >  5:529861/1‑149 (MQ=255)

                                                                                                                                         tAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTgg                                                                                                                                                                                                                                                                    <  5:116731/149‑1 (MQ=255)

                                                                                                                                          aaTTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCg                                                                                                                                                                                                                                                                                >  4:91535/1‑136 (MQ=255)

                                                                                                                                          aaTTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGga                                                                                                                                                                                                                                                                   <  8:30196/149‑1 (MQ=255)

                                                                                                                                          aaTTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGga                                                                                                                                                                                                                                                                   <  7:69807/149‑1 (MQ=255)

                                                                                                                                          aaTTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGga                                                                                                                                                                                                                                                                   <  5:27147/149‑1 (MQ=255)

                                                                                                                                          aaTTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCGTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGga                                                                                                                                                                                                                                                                   <  6:849591/149‑1 (MQ=255)

                                                                                                                                          aaTTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCATCGCAGATTGGTTTTAGCCCCAACCAGCCCGGTACCGGTCTggt                                                                                                                                                                                                                                                                   <  6:324116/149‑2 (MQ=255)

                                                                                                                                            tttACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGCCCGGTACCGGTCTggt                                                                                                                                                                                                                                                                   <  8:745240/147‑2 (MQ=255)

                                                                                                                                               aCAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGCCCGGTACCGGTCTggttgggg                                                                                                                                                                                                                                                              <  8:657602/149‑7 (MQ=255)

                                                                                                                                                 aGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGaaa                                                                                                                                                                                                                                                            >  7:736473/1‑149 (MQ=255)

                                                                                                                                                  gggCGCTCGGCAACTGGGCGTTCGTATTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGACGGCACTGGCGGGAATTGTGTTTGCTTCGCCGATTGGGTTTATCCCCAACCTGACCGGTACCGGGCTGGAGGTGAAAg                                                                                                                                                                                                                                                           >  6:40108/1‑149 (MQ=255)

                                                                                                                                                     cggtCGTCAACTGGTCGATCGTACTGAAGACATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCa                                                                                                                           <  2:374256/271‑1 (MQ=255)

                                                                                                                                                     cgcTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCt                                                                                                                                                                                                                                                                      <  6:51482/135‑1 (MQ=255)

                                                                                                                                                     cgcTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCAtt                                                                                                                                                                                                >  2:151994/1‑205 (MQ=255)

                                                                                                                                                       cTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCg                                                                                                                                                                                                                                                                          <  2:155144/129‑1 (MQ=255)

                                                                                                                                                        tCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTg                                                                                                                                                                                                                                                     >  7:870967/1‑149 (MQ=255)

                                                                                                                                                         cGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTg                                                                                                                                                                                            <  4:151994/205‑1 (MQ=255)

                                                                                                                                                          gTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCa                                                                                                                                                                                                                                                   <  6:572567/149‑1 (MQ=255)

                                                                                                                                                            cAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGc                                                                                                                                                                                                                                                 <  7:559991/149‑1 (MQ=255)

                                                                                                                                                                    cGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGc                                                                                                                                                                                                                                         <  5:428722/149‑1 (MQ=255)

                                                                                                                                                                    cGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACGCGTCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCCTGGGGGAATGATTTTa                                                                                                           >  1:160640/1‑275 (MQ=255)

                                                                                                                                                                      ttCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTg                                                                                                                                                                                                                                       >  6:850046/1‑149 (MQ=255)

                                                                                                                                                                         gTACTGGAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTTTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATcgc                                                                                                       <  1:195572/274‑1 (MQ=255)

                                                                                                                                                                             tGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCt                                                                                                                                                                                                                                                                      >  7:46070/1‑111 (MQ=255)

                                                                                                                                                                             tGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCt                                                                                                                                                                                                                                                                      <  8:46070/111‑1 (MQ=255)

                                                                                                                                                                             tGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCa                                                                                                                                                                                                                                >  5:754391/1‑149 (MQ=255)

                                                                                                                                                                             tGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCa                                                                                                                                                                                                                                <  7:590503/149‑1 (MQ=255)

                                                                                                                                                                              gAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCAt                                                                                                                                                                                                                               >  8:822935/1‑149 (MQ=255)

                                                                                                                                                                                  ccATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGt                                                                                                                                                                                                                           <  5:807174/149‑1 (MQ=255)

                                                                                                                                                                                   cATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGtt                                                                                                                                                                                                                          <  6:529861/149‑1 (MQ=255)

                                                                                                                                                                                   cATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGtt                                                                                                                                                                                                                          <  5:385523/149‑1 (MQ=255)

                                                                                                                                                                                    aTTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGttt                                                                                                                                                                                                                         >  7:709965/1‑149 (MQ=255)

                                                                                                                                                                                       cGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCt                                                                                                                                                                                                                      <  7:125888/149‑1 (MQ=255)

                                                                                                                                                                                        gCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTc                                                                                                                                                                                                                     <  8:555094/149‑1 (MQ=255)

                                                                                                                                                                                           ttGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCggt                                                                                                                                                                                                                  <  5:762486/149‑1 (MQ=255)

                                                                                                                                                                                            tgtgGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTTGGCGGCATTAGTTTGCTCGGTGGTTTCGGTGCGATCan                                                                                                                                                                                                 >  4:742261/1‑164 (MQ=255)

                                                                                                                                                                                                   ttttCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTAGCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCgg                                                                                                                                                                                                          <  6:456168/149‑1 (MQ=255)

                                                                                                                                                                                                    tttCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGt                                                                                                                                                                                                         <  7:403949/149‑1 (MQ=255)

                                                                                                                                                                                                         ttGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGAt                                                                                                                                                                                                                                                                <  7:789543/89‑1 (MQ=255)

                                                                                                                                                                                                             aCGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCAtt                                                                                                                                                                                                <  6:745329/149‑1 (MQ=255)

                                                                                                                                                                                                              cGGCTGCATGGCGACACTGGCGGGAAGTGTGGGTGCTTCGCAGATTGGTTGTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCAGTAGTTTGCTCGGTGGTTCCGGTGCGATCATTg                                                                                                                                                                                               <  6:341456/149‑1 (MQ=255)

                                                                                                                                                                                                                  tGCAGGGCGGGATGGGCGGGAATTGGGTTTGGTTCGCAGATTGGTTTTATCCCCAACCAGAGCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCt                                                             <  1:5792/275‑1 (MQ=255)

                                                                                                                                                                                                                   gCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCg                                                                                                                                                                                          <  7:1054766/149‑1 (MQ=255)

                                                                                                                                                                                                                           gcACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCgg                                                                                                                                                                                  >  8:2393/1‑149 (MQ=255)

                                                                                                                                                                                                                               tGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGgc                                                                                                                                                                                 <  6:925326/146‑1 (MQ=255)

                                                                                                                                                                                                                               tGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGgc                                                                                                                                                                                 >  5:925326/1‑146 (MQ=255)

                                                                                                                                                                                                                                ggCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCAt                                                                                                                                                                             <  8:585139/149‑1 (MQ=255)

                                                                                                                                                                                                                                ggCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCAt                                                                                                                                                                             <  5:630749/149‑1 (MQ=255)

                                                                                                                                                                                                                                   gggAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTa                                                                                                                                                                                       >  5:36541/1‑136 (MQ=255)

                                                                                                                                                                                                                                   gggAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTa                                                                                                                                                                                       <  6:36541/136‑1 (MQ=255)

                                                                                                                                                                                                                                    ggAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGtt                                                                                                                                                                         <  8:769940/149‑1 (MQ=255)

                                                                                                                                                                                                                                     gAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCggtg                                                                                                                                                                                                                 >  5:58621/1‑108 (MQ=255)

                                                                                                                                                                                                                                     gAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCggtg                                                                                                                                                                                                                 <  6:58621/108‑1 (MQ=255)

                                                                                                                                                                                                                                      aaTTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTcc                                                                                                                                                                       <  8:717189/149‑1 (MQ=255)

                                                                                                                                                                                                                                       aTTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCt                                                                                                                                                                      >  7:122201/1‑149 (MQ=255)

                                                                                                                                                                                                                                               tGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGAAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGa                                                                                                                                                              <  6:1042965/149‑1 (MQ=255)

                                                                                                                                                                                                                                                gCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCACATGGTTCCTGACGCAGAt                                                                                                                                                             <  6:207833/149‑1 (MQ=255)

                                                                                                                                                                                                                                                   tCGCAGATTGGTTTTATCCCCAACCAGACCGGTCCCGGGCTGGAGATGAACGCAATTGCAGCCTGCGTGCTTGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACggcggc                             <  2:433660/274‑1 (MQ=255)

                                                                                                                                                                                                                                                     gCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATa                                                                                                                                                        >  6:589043/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     gCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATa                                                                                                                                                        >  7:353060/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     gCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATa                                                                                                                                                        <  8:353060/149‑1 (MQ=255)

                                                                                                                                                                                                                                                         aTTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGt                                                                                                                                                    >  6:56154/1‑149 (MQ=255)

                                                                                                                                                                                                                                                         aTTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATAGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGt                                                                                                                                                    >  7:924890/1‑149 (MQ=255)

                                                                                                                                                                                                                                                              ttttATCCCCAACCAGACCGTTACCGGGCTGGGGATGAAAGGAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTgg                                                                                                                                               >  8:207355/1‑149 (MQ=255)

                                                                                                                                                                                                                                                              ttttATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGAt                                                                                                                                                             >  6:769060/1‑135 (MQ=255)

                                                                                                                                                                                                                                                              ttttATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTgg                                                                                                                                               >  6:670166/1‑149 (MQ=255)

                                                                                                                                                                                                                                                              ttttATCCCCAACCAGACCGGGACCGGGCCGGAGAGGAAAGCAATTGCAGCCTGCTTGCTGGGCGGCAATAGTTGGCTCTGTGGTACCGTTGCGATCATTGGTGCAGTACTCCGCGCATAGTACCTGACGCAGATCGATAGCGTACTgg                                                                                                                                               >  8:670815/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                ttATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTg                                                                                                                                             <  8:457041/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                  aTCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCt                                                                                                                                           <  5:566146/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                  aTCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCt                                                                                                                                           >  6:567143/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                      ccAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTg                                                                                                                                       >  7:497598/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                      ccAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTg                                                                                                                                       >  7:328264/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                      ccAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTg                                                                                                                                       >  8:731525/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                       cAACCAGACCAGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTgc                                                                                                                                      <  8:924890/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                             gACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTTCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGCAAAAATATGCCCGCTTTATGACGCCa   >  2:56739/1‑274 (MQ=255)

                                                                                                                                                                                                                                                                             caCCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTc                                                                                                                                <  7:315729/148‑1 (MQ=255)

                                                                                                                                                                                                                                                                              aCCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGCAAAAATATGCCCGCTTTATGACGCCAc  >  2:379936/1‑274 (MQ=255)

GCCGAACTGAAACAGCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGCAAAAATATGCCCGCTTTATGACGCCAC  >  NC_000913/1603454‑1603997

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