BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	475,971	C→G	8.9%	intergenic (+20/+11)	ybaY → / ← ybaZ	outer membrane lipoprotein/excision repair protein, alkyltransferase‑like protein ATL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	475,971	0	C	G	8.9%	868.5 / 61.4	256	intergenic (+20/+11)	ybaY/ybaZ	outer membrane lipoprotein/excision repair protein, alkyltransferase‑like protein ATL
Reads supporting (aligned to +/- strand): ref base C (119/114); new base G (12/11); total (131/125)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.47e-01

Only 100 of 260 total aligned reads displayed.

GATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGGAGACGCTTTAACACACCGCCCACCTGGCGCGCGGCGCGGGGCGATCCCGCCAGTTTCGCCACATCACCGTAAGTGGTGACATAGCCTTCGGGA  >  NC_000913/475703‑476240

              cgcgcgTATTCTGTTGAGTGCGGCGATCACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCCCCAACTCAAGTGAATCCGTCTTCGGCAGTTCCCGCTCCCACGCAATATTAAGCGACGTTAACCCTCTCCCCCCGGGAAGGGTtag                                                                                                                                                                                                                                                                 >  4:192906/1‑269 (MQ=255)

                  cgTATTATGTTGAGTGCGGCGATCACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTagtagt                                                                                                                                                                                                                                                             <  2:192906/269‑1 (MQ=255)

                         tGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGGAGGGGTTAGTAGTTCCAGCGATAACg                                                                                                                                                                                                                                                >  2:306487/1‑275 (MQ=255)

                                                                tatcccccgagcccGTCCAGCCTGTGACNANCCCGGGCGGCACTAAAGCCGCCCTGACATGGGGGCCGGAACAGCANGCCGCCGTGCCGGGTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGNATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCGCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACAcc                                                                                                                                                                                                         <  3:37561/263‑1 (MQ=255)

                                                                       cGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGc                                                                                                                                                                                                   >  2:268069/1‑274 (MQ=255)

                                                                                                  gggCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTCATGCCTGTCGCTGACGCTGTaaa                                                                                                                                                                       >  2:74632/1‑275 (MQ=255)

                                                                                                                             ggTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGga                                                                                                                                                                                                                                                                          <  5:333756/149‑1 (MQ=255)

                                                                                                                             ggTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGga                                                                                                                                                                                                                                                                          <  8:852650/149‑1 (MQ=255)

                                                                                                                                 ccGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAgg                                                                                                                                                                                                                                                                      <  8:555403/149‑1 (MQ=255)

                                                                                                                                 ccGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAgg                                                                                                                                                                                                                                                                      <  8:122935/149‑1 (MQ=255)

                                                                                                                                 ccGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAgg                                                                                                                                                                                                                                                                      <  7:1105178/149‑1 (MQ=255)

                                                                                                                                 ccGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAgg                                                                                                                                                                                                                                                                      <  6:591394/149‑1 (MQ=255)

                                                                                                                                   ggTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAGGGt                                                                                                                                                                                                                                                                    <  8:1334374/149‑1 (MQ=255)

                                                                                                                                   ggTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGt                                                                                                                                                                                                                                                                    <  7:285955/149‑1 (MQ=255)

                                                                                                                                   ggTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGa                                                                                                                                      >  4:70413/1‑275 (MQ=255)

                                                                                                                                    gTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAGGGtt                                                                                                                                                                                                                                                                   <  6:941575/149‑1 (MQ=255)

                                                                                                                                    gTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGtt                                                                                                                                                                                                                                                                   <  5:708064/149‑1 (MQ=255)

                                                                                                                                     tACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTta                                                                                                                                                                                                                                                                  >  7:233686/1‑149 (MQ=255)

                                                                                                                                        aGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagta                                                                                                                                                                                                                                                               <  5:8542/149‑1 (MQ=255)

                                                                                                                                        aGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGGTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagca                                                                                                                                                                                                                                                               >  8:1252316/1‑147 (MQ=255)

                                                                                                                                        aGCAAACCGCCGTGCCGGCTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTCCCACTTCACCAACTCAGGTGACTCCGTCTTCGGCGGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGACGGCtagta                                                                                                                                                                                                                                                               >  8:329393/1‑149 (MQ=255)

                                                                                                                                         gCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagtag                                                                                                                                                                                                                                                              >  8:899861/1‑149 (MQ=255)

                                                                                                                                          cAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCCAGTCGATTTGCCCGCTTCCCTATCCCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTACATCCGGTCCGGTTAGCGAACTTGTGCCGTGGCGACTAacca                                                                                                                               >  4:88974/1‑275 (MQ=255)

                                                                                                                                             aCCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAaccacc                                                                                                                             <  4:109998/274‑1 (MQ=255)

                                                                                                                                             aCCGCCGTGCCGGGTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCGTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAaccacc                                                                                                                             <  2:128015/274‑1 (MQ=255)

                                                                                                                                               cgccgTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTCAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTCACCACCCgg                                                                                                                          >  2:339748/1‑275 (MQ=255)

                                                                                                                                                  cgTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAAACCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTg                                                                                                                        <  1:37008/274‑1 (MQ=255)

                                                                                                                                                   gTGCCGGGTCAGGCCAGCGGGGGCGCAACGACTACCGTACCTTCGACTTAACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGc                                                                                                                       <  4:268069/274‑1 (MQ=255)

                                                                                                                                                    tGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTGGACCGCTTCCCGATACCATCACACCTTCTGCAAGTAATGCCTGTCGCTGACGCTGTAAATCAGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGTTGcc                                                                                                                      >  2:156924/1‑274 (MQ=255)

                                                                                                                                                     gCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATaa                                                                                                                                                                                                                                                  <  6:80282/149‑1 (MQ=255)

                                                                                                                                                      ccGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAAc                                                                                                                                                                                                                                                 <  6:450925/149‑1 (MQ=255)

                                                                                                                                                              ggCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGTTATCGCTGGaac                                                                                                                                                                                                                                         <  6:1004200/149‑2 (MQ=255)

                                                                                                                                                              ggCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAg                                                                                                                                                                                                                                         >  6:255434/1‑149 (MQ=255)

                                                                                                                                                                   gCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATCAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGAt                                                                                                                                                                                                                                    <  5:837628/149‑1 (MQ=255)

                                                                                                                                                                           gcAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGagag                                                                                                                                                                                                                                                                       >  5:55162/1‑106 (MQ=255)

                                                                                                                                                                           gcAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGc                                                                                                                                                                                                                            >  7:727013/1‑149 (MQ=255)

                                                                                                                                                                           gcAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGc                                                                                                                                                                                                                            >  8:1199194/1‑149 (MQ=255)

                                                                                                                                                                            cAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCt                                                                                                                                                                                                                           <  7:440272/149‑1 (MQ=255)

                                                                                                                                                                             aaCGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCtt                                                                                                                                                                                                                          >  7:853806/1‑149 (MQ=255)

                                                                                                                                                                              aCGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTc                                                                                                                                                                                                                         >  7:1416623/1‑149 (MQ=255)

                                                                                                                                                                                gACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTccc                                                                                                                                                                                                                       <  5:145357/149‑1 (MQ=255)

                                                                                                                                                                                   tACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGAt                                                                                                                                                                                                                    >  7:1163555/1‑149 (MQ=255)

                                                                                                                                                                                     ccGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagt                                                                                                                                                                                                                                                                <  5:42931/103‑1 (MQ=255)

                                                                                                                                                                                        tACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCAt                                                                                                                                                                                                               >  7:1036388/1‑149 (MQ=255)

                                                                                                                                                                                        tACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCAt                                                                                                                                                                                                               <  6:358628/149‑1 (MQ=255)

                                                                                                                                                                                            ttCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATcaca                                                                                                                                                                                                           >  7:587804/1‑149 (MQ=255)

                                                                                                                                                                                            ttCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATcaca                                                                                                                                                                                                           <  7:752265/149‑1 (MQ=255)

                                                                                                                                                                                             tCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATAcc                                                                                                                                                                                                                 >  4:346178/1‑142 (MQ=255)

                                                                                                                                                                                              cGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGAt                                                                                                                                                                                                                    <  4:47666/138‑1 (MQ=255)

                                                                                                                                                                                                aCTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCtt                                                                                                                                                                                                       >  6:1350659/1‑149 (MQ=255)

                                                                                                                                                                                                aCTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCtt                                                                                                                                                                                                       >  6:1079484/1‑149 (MQ=255)

                                                                                                                                                                                                  ttCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCt                                                                                                                                                                                                     >  8:1122844/1‑149 (MQ=255)

                                                                                                                                                                                                  ttCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCt                                                                                                                                                                                                     >  6:11634/1‑149 (MQ=255)

                                                                                                                                                                                                   tCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGcc                                                                                                                                                                                                                                                                              <  3:49690/75‑1 (MQ=255)

                                                                                                                                                                                                    cACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGc                                                                                                                                                                                                   <  5:36586/149‑1 (MQ=255)

                                                                                                                                                                                                       cAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAgg                                                                                                                                                                                                                                                                      >  5:1446821/1‑79 (MQ=255)

                                                                                                                                                                                                        aaCTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGt                                                                                                                                                                                               >  5:1107987/1‑149 (MQ=255)

                                                                                                                                                                                                        aaCTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGt                                                                                                                                                                                               >  7:835664/1‑149 (MQ=255)

                                                                                                                                                                                                          cTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTaa                                                                                                                                                                                             <  5:224460/149‑1 (MQ=255)

                                                                                                                                                                                                             aGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGc                                                                                                                                                                                          <  5:255434/149‑1 (MQ=255)

                                                                                                                                                                                                              ggTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGcc                                                                                                                                                                                         >  6:1466832/1‑149 (MQ=255)

                                                                                                                                                                                                               gTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCt                                                                                                                                                                                        >  8:1410352/1‑149 (MQ=255)

                                                                                                                                                                                                                   aTCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCt                                                                                                                                                                                        >  5:389281/1‑145 (MQ=255)

                                                                                                                                                                                                                   aTCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCt                                                                                                                                                                                        <  6:389281/145‑1 (MQ=255)

                                                                                                                                                                                                                       gTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGa                                                                                                                                                                                >  6:501491/1‑149 (MQ=255)

                                                                                                                                                                                                                       gTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGa                                                                                                                                                                                <  8:1175085/149‑1 (MQ=255)

                                                                                                                                                                                                                        tCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGAc                                                                                                                                                                               >  6:1299959/1‑149 (MQ=255)

                                                                                                                                                                                                                             ggCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGt                                                                                                                                                                          <  5:1087471/149‑1 (MQ=255)

                                                                                                                                                                                                                              gCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTa                                                                                                                                                                         >  5:41982/1‑149 (MQ=255)

                                                                                                                                                                                                                                 gTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAt                                                                                                                                                                      >  8:707097/1‑149 (MQ=255)

                                                                                                                                                                                                                                  ttCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtc                                                                                                                                                                     <  6:66662/149‑1 (MQ=255)

                                                                                                                                                                                                                                        gtgcaaagCAATAAGAAGCGACGTTAACCATCTCCGCCCGGAGAGGGTTAGAAGTTCCAGCGATAACGCTGCAAGTAGATGTGCACGCTTCCCGATACCATCACACCTGCTGCCAGGAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAAACACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGGAGACGCTTTAACACACCGCCCACCTGGCGCGCGGCGCGGGGCGATCCCGCCAGTTTCGCcac                                 <  1:73061/268‑1 (MQ=255)

                                                                                                                                                                                                                                             aCGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTa                                                                                                                                                          <  6:197298/149‑1 (MQ=255)

                                                                                                                                                                                                                                              cGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGGTACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAg                                                                                                                                                         >  8:783893/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  atatTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGaa                                                                                                                                                     <  8:1400573/149‑1 (MQ=255)

                                                                                                                                                                                                                                                   tatTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGaaa                                                                                                                                                    <  8:39549/149‑1 (MQ=255)

                                                                                                                                                                                                                                                    atTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAAt                                                                                                                                                   >  8:243368/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     ttAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAAtt                                                                                                                                                  >  6:1383803/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     acgacacaCGTCAACCCTCCCCGCCCGGAGAGGGTTATTAGTTCCAGCGATCACGCTGCACGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCGGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCTGTCCGTTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGGAGACGCTTTAACACACCGCCCACCTGGCGCGCGGCGCGGGGCGATCCCGCCAGTTTCGCCACATCACCGTAAgtg                    <  1:6658/268‑1 (MQ=255)

                                                                                                                                                                                                                                                        aGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCCGGGTAAGGTGCTGCCTTCAGGGAGCCGCTTTAACACACCTCCCACCTGGCGCGCGGCGCGGGGCGACCCCGCCAGTTTCCCCACTTCCCCGTAAggggg                  >  4:182870/1‑273 (MQ=255)

                                                                                                                                                                                                                                                        aGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGGAGACGCTTTAACAAACCGCCCACCTGGCGCGCGGCGCGGGGCGATCCCGCCAGTTTCGCCACATCACCGTAAgtggtc                 >  2:174541/1‑274 (MQ=255)

                                                                                                                                                                                                                                                           actactAACCCTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCg                                                                                                                                            >  7:1266659/6‑149 (MQ=255)

                                                                                                                                                                                                                                                            aCGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGt                                                                                                                                           >  5:739153/1‑149 (MQ=255)

                                                                                                                                                                                                                                                             tactAACCCTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTACTGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTg                                                                                                                                          >  6:1123037/4‑149 (MQ=255)

                                                                                                                                                                                                                                                             cGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCAGGTGCCTGGGTAAGGTGCTGCCTTCAGGGAGACGCTTTAAAACACCGCCCACCTGGCGCGCGGCGCGGGGCGATCCCGCCAGTTTCGCCACATCCCCGTAAGTGGTGACAt             >  1:48738/1‑274 (MQ=255)

                                                                                                                                                                                                                                                               ttAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGc                                                                                                                                        <  7:955277/149‑1 (MQ=255)

                                                                                                                                                                                                                                                               ttAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGt                                                                                                                         >  4:279812/1‑164 (MQ=255)

                                                                                                                                                                                                                                                                tAACCCTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCg                                                                                                                                       >  8:1021486/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                tAACCCTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCg                                                                                                                                       >  7:1388910/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                 aaCCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGa                                                                                                                                      <  5:906527/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                   cccTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGAtt                                                                                                                                    <  6:739153/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                   cccTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGc                                                                                                                       <  2:363721/162‑1 (MQ=255)

                                                                                                                                                                                                                                                                   cccTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCa                                                                                                                     <  2:279812/164‑1 (MQ=255)

                                                                                                                                                                                                                                                                    ccTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTa                                                                                                                                   >  7:1332950/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     ctctCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccggt                                                                                                                                                                 >  8:1046933/1‑118 (MQ=255)

                                                                                                                                                                                                                                                                     ctctCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTaa                                                                                                                                  <  5:1079484/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                     ctctCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTaa                                                                                                                                  >  5:1215490/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     ctctCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTaa                                                                                                                                  <  5:1350659/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                     ctctCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTaa                                                                                                                                  <  7:1016429/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                     ctctCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTaa                                                                                                                                  <  6:821989/149‑1 (MQ=255)

GATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGGAGACGCTTTAACACACCGCCCACCTGGCGCGCGGCGCGGGGCGATCCCGCCAGTTTCGCCACATCACCGTAAGTGGTGACATAGCCTTCGGGA  >  NC_000913/475703‑476240

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