BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,770,203	T→A	9.0%	intergenic (+17/‑169)	ydiK → / → rprA	UPF0118 family inner membrane protein/sRNA antisense activator for RpoS translation, Hfq‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,770,203	0	T	A	9.0%	512.3 / 30.7	157	intergenic (+17/‑169)	ydiK/rprA	UPF0118 family inner membrane protein/sRNA antisense activator for RpoS translation, Hfq‑dependent
Reads supporting (aligned to +/- strand): ref base T (75/67); new base A (6/8); total (82/75)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

Only 100 of 160 total aligned reads displayed.

ACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCC  >  NC_000913/1769941‑1770459

aCGCTGGATAACGTCATCTGCCCAATGTTAATTCGCCTGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGa                                                                                                                                                                                                                                                       <  2:324870/274‑1 (MQ=255)

              cATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGAGGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGc                                                                                                                                                                                                                                         <  2:36225/274‑1 (MQ=255)

                tCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAacca                                                                                                                                                                                                                                             >  4:135658/1‑268 (MQ=255)

                 ccGCCCAATGTTAGTTCGCTTGGGTGCGTATTTACCGTTTGTCCTGATTCTTTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAAt                                                                                                                                                                                                                                     <  4:177597/275‑1 (MQ=255)

                    cccAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGc                                                                                                                                                                                                                                         <  2:135658/268‑1 (MQ=255)

                                      agggccagATTTACCGTTGATCCTGATTTTTCTTGGGGTTATTGGTGGTTTGATTGCTTTCGGGATGAGCGGTCTGTTTATTGGTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTcc                                                                                                                                                                                                                 <  4:217049/267‑1 (MQ=255)

                                                         aTCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGctc                                                                                                                                                                                                                     >  2:217049/1‑251 (MQ=255)

                                                                               ttGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATTTTTGCTTGCCTTTCCCCATTTTTTCCTGCTTAACCCCATTTACCTggctataccct                                                                                                                                                                        >  4:291086/1‑265 (MQ=255)

                                                                                   tggtTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCa                                                                                                                                                                   >  4:34577/1‑275 (MQ=255)

                                                                                   tggtTTGATTGCTTTCGGGAGGATCGGTCTTTTTATTGTTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTTGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAg                                                                                                                                                                                     <  2:291086/257‑1 (MQ=255)

                                                                                             gCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTAc                                                                                                                                                         <  2:34577/275‑1 (MQ=255)

                                                                                             gCTTTAGGGATGATCCGTGTGTGTATTGGTCCGGTGCTGTTAGACGTTTCCGGGCGTCTGTTTGCAGCGTGGGTGGAAGAAGTCCCGACGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCAACAGGCGGCTCTGCCGCCGGATTGTTAACCACCGCTAATTGATATTTGGTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTAc                                                                                                                                                         <  2:196063/275‑1 (MQ=255)

                                                                                              ttttCGGGGTGGTCGGTTTGTGTATTGGTCCGGTGCTGTTAGCGGTTTCTGGGGGTCGGGTGGCGTGGTGGGTGGAAGAAGGCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACa                                                                                                                                                        <  4:213995/274‑1 (MQ=255)

                                                                                               tttCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCACCCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGATAAACCGAATAAGTAATTTCTCATCAGGCGGCTATGCCGCCTGATTGTTAACCACCGCTAATTAATATTTTCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCGTTCGCAAATTCtttcctg                                                                                                                                                       >  4:336465/1‑271 (MQ=255)

                                                                                                                tggttATTGGTCCGGTTCGGTTAGCCGTTGCCGTGCGTCTGTTTGCCGAGTGGGTAGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATgaga                                                                                                                                      <  2:232171/272‑1 (MQ=255)

                                                                                                                  tttATTGGTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGct                                                                                                                                                                                                                                                                  <  5:762817/149‑1 (MQ=255)

                                                                                                                           ccGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCt                                                                                                                                                                                                                                                         <  6:325479/149‑1 (MQ=255)

                                                                                                                            cGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTg                                                                                                                                                                                                                                                        <  6:927374/149‑1 (MQ=255)

                                                                                                                             ggTGCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGa                                                                                                                         <  2:48295/275‑1 (MQ=255)

                                                                                                                              gTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAt                                                                                                                                                                                                                                                      <  5:589019/149‑1 (MQ=255)

                                                                                                                               ttCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGAtt                                                                                                                                                                                                                                                     <  7:1380036/149‑1 (MQ=255)

                                                                                                                                tCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAtga                                                                                                                                                                                                                                                    <  5:279744/149‑3 (MQ=255)

                                                                                                                                  tGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGtt                                                                                                                                                                                                                                                  <  7:231223/149‑1 (MQ=255)

                                                                                                                                    ttAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCaaa                                                                                                                   >  2:274149/1‑274 (MQ=255)

                                                                                                                                         cGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAATAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTc                                                                                                              <  2:237014/274‑1 (MQ=255)

                                                                                                                                          gTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCg                                                                                                                                                                                                                                          >  8:1426130/1‑149 (MQ=255)

                                                                                                                                           tttCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGc                                                                                                                                                                                                                                         >  7:454856/1‑149 (MQ=255)

                                                                                                                                             tccttggtgttGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGACGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTa                                                                                                          <  4:28352/264‑1 (MQ=255)

                                                                                                                                               cTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAAt                                                                                                                                                                                                                                     >  6:723520/1‑149 (MQ=255)

                                                                                                                                                tGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAAtt                                                                                                                                                                                                                                    >  8:1328797/1‑149 (MQ=255)

                                                                                                                                                 ggCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATGTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTg                                                                                                                                                                                                                                   >  5:435482/1‑149 (MQ=255)

                                                                                                                                                  gCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGa                                                                                                                                                                                                                                  <  7:1426130/149‑1 (MQ=255)

                                                                                                                                                   cGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGat                                                                                                                                                                                                                                 <  8:850920/149‑1 (MQ=255)

                                                                                                                                                   cGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGat                                                                                                                                                                                                                                 >  8:685699/1‑149 (MQ=255)

                                                                                                                                                    gTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGata                                                                                                                                                                                                                                >  5:620499/1‑149 (MQ=255)

                                                                                                                                                        gTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTaa                                                                                                                                                                                                                                      >  7:941423/1‑139 (MQ=255)

                                                                                                                                                         tttGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgc                                                                                                                                                                                                                           >  5:487414/1‑149 (MQ=255)

                                                                                                                                                         tttGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgc                                                                                                                                                                                                                           <  5:723520/149‑1 (MQ=255)

                                                                                                                                                          ttGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgct                                                                                                                                                                                                                          <  5:85614/149‑1 (MQ=255)

                                                                                                                                                           tGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGTACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgctt                                                                                                                                                                                                                         <  5:1213643/149‑1 (MQ=255)

                                                                                                                                                           tGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgctt                                                                                                                                                                                                                         <  8:8175/149‑1 (MQ=255)

                                                                                                                                                            gCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgcttg                                                                                                                                                                                                                        <  8:454856/149‑1 (MQ=255)

                                                                                                                                                            gCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTc                                                                                          >  4:13906/1‑275 (MQ=255)

                                                                                                                                                              cgcgTGGGTGGAAGAAGTCGCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccgc                                                                                                                                                                                                                                                           >  7:856042/1‑112 (MQ=255)

                                                                                                                                                              cgcgTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccgc                                                                                                                                                                                                                                                           <  8:856042/112‑1 (MQ=255)

                                                                                                                                                                  tgggtggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTc                                                                                                                                                                                                                  >  8:777031/1‑149 (MQ=255)

                                                                                                                                                                   gggtggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTcc                                                                                                                                                                                                                 <  7:1221666/149‑1 (MQ=255)

                                                                                                                                                                        gaagaaGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCata                                                                                                                                                                                                            >  5:927398/1‑149 (MQ=255)

                                                                                                                                                                        gaagaaGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCata                                                                                                                                                                                                            <  5:305778/149‑1 (MQ=255)

                                                                                                                                                                              gTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTcc                                                                                                                                                                                                      >  5:744780/1‑149 (MQ=255)

                                                                                                                                                                               tCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCt                                                                                                                                                                                                     <  7:368298/149‑1 (MQ=255)

                                                                                                                                                                               tCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCt                                                                                                                                                                                                     <  6:487414/149‑1 (MQ=255)

                                                                                                                                                                                 ccgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGc                                                                                                                                                                                                   >  5:448303/1‑149 (MQ=255)

                                                                                                                                                                                   gccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCa                                                                                         <  2:307804/253‑1 (MQ=255)

                                                                                                                                                                                         gACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGcc                                                                                                                                                                                           >  8:665984/1‑149 (MQ=255)

                                                                                                                                                                                            ggACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCAtt                                                                                                                                                                                        >  6:1223049/1‑149 (MQ=255)

                                                                                                                                                                                            ggACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCAtt                                                                                                                                                                                        >  6:1049664/1‑149 (MQ=255)

                                                                                                                                                                                              aCCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCa                                                                                                                                                                                      >  6:345958/1‑149 (MQ=255)

                                                                                                                                                                                              aCCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCa                                                                                                                                                                                      <  7:756305/149‑1 (MQ=255)

                                                                                                                                                                                                cAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGc                                                                                                                                                                                    <  7:705229/149‑1 (MQ=255)

                                                                                                                                                                                                cAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGc                                                                                                                                                                                    <  8:626968/149‑1 (MQ=255)

                                                                                                                                                                                                     ggAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgta                                                                                                                                                                               <  7:638735/149‑1 (MQ=255)

                                                                                                                                                                                                            aTTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACGACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTg                                                                                                                                                                        >  8:1062670/1‑149 (MQ=255)

                                                                                                                                                                                                             ttCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGt                                                                                                                                                                       <  5:1144251/149‑1 (MQ=255)

                                                                                                                                                                                                                tcGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGc                                                                                                                                                                    >  6:552160/1‑149 (MQ=255)

                                                                                                                                                                                                                tcGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGc                                                                                                                                                                    >  8:1455219/1‑149 (MQ=255)

                                                                                                                                                                                                                tcGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGc                                                                                                                                                                    >  8:753062/1‑149 (MQ=255)

                                                                                                                                                                                                                  gaagaaCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCaa                                                                                                                                                                  >  8:1329806/1‑149 (MQ=255)

                                                                                                                                                                                                                    agaaCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAAt                                                                                                                                                                >  7:1213282/1‑149 (MQ=255)

                                                                                                                                                                                                                    agaaCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAAt                                                                                                                                                                >  7:779175/1‑149 (MQ=255)

                                                                                                                                                                                                                        cTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCtt                                                                                                                                                            <  7:1329806/149‑1 (MQ=255)

                                                                                                                                                                                                                          tGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTa                                                                                                                                                          <  6:661839/149‑1 (MQ=255)

                                                                                                                                                                                                                           ggCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTAc                                                                                                                                                         >  5:704130/1‑149 (MQ=255)

                                                                                                                                                                                                                                  tCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCTGCTGGTAGTACCTGTCGCAAATTCTTTACAGttttt                                                                                                                                                  >  7:973602/1‑149 (MQ=255)

                                                                                                                                                                                                                                    gagaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTaa                                                                                                                                                >  7:404292/1‑149 (MQ=255)

                                                                                                                                                                                                                                     agaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTaaa                                                                                                                                               <  7:247923/149‑1 (MQ=255)

                                                                                                                                                                                                                                     agaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTaaa                                                                                                                                               <  7:665984/149‑1 (MQ=255)

                                                                                                                                                                                                                                          ccGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAAt                                                                                                                                          <  6:538124/149‑1 (MQ=255)

                                                                                                                                                                                                                                            gAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAAc                                                                          >  4:382771/1‑211 (MQ=255)

                                                                                                                                                                                                                                               tAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGAc                                                                                                                                     <  6:927398/149‑1 (MQ=255)

                                                                                                                                                                                                                                                aaGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACg                                                                                                                                    >  6:184516/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    aaTTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTc                                                                                                                                                                      >  3:24055/1‑111 (MQ=255)

                                                                                                                                                                                                                                                     aTTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCt                                                                                                                               >  7:736237/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     aTTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTNTACAGTTNTTANACTAATNAGACGAATCTGATCGACGCAAAAAGTCCGTa                                                                                                          >  3:19940/1‑170 (MQ=255)

                                                                                                                                                                                                                                                      tttCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATGCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTcc  <  2:336465/273‑1 (MQ=255)

                                                                                                                                                                                                                                                        tctcATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCaa                                                                                                                                                                  <  1:24055/111‑1 (MQ=255)

                                                                                                                                                                                                                                                         ctcATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATc                                                                                                                           <  6:141229/149‑1 (MQ=255)

                                                                                                                                                                                                                                                         ctcATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGcc                                                                                                      <  1:19940/170‑1 (MQ=255)

                                                                                                                                                                                                                                                           aatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCtt                                                                                                                                                            >  5:298027/2‑114 (MQ=255)

                                                                                                                                                                                                                                                            atcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGAc                                                                                                                        >  5:1421410/1‑149 (MQ=255)

                                                                                                                                                                                                                                                            atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCaaa                                                                                                                                                                 <  6:297120/108‑1 (MQ=255)

                                                                                                                                                                                                                                                            atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATgaga                                                                                                                                      >  8:1262751/1‑135 (MQ=255)

                                                                                                                                                                                                                                                            atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGAc                                                                                                                        >  6:680001/1‑149 (MQ=255)

                                                                                                                                                                                                                                                             tcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGACATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATga                                                                                                                                        >  8:957505/1‑132 (MQ=255)

                                                                                                                                                                                                                                                               aGGCGGCTCTGCCGCCTGATTGTTAACTACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCa                                                                                                                     >  7:975312/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                ggcggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCaa                                                                                                                    >  8:366166/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                  cggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCaaaa                                                                                                                  >  5:27034/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                   ggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCaaaaa                                                                                                                 <  8:1162069/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                      tctGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTc                                                                                                              >  6:127452/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                      tctGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTc                                                                                                              <  5:683298/149‑1 (MQ=255)

ACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCC  >  NC_000913/1769941‑1770459

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