BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	475,972	C→G	5.9%	intergenic (+21/+10)	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,972	0	C	G	5.9%	736.2 / 28.4	217	intergenic (+21/+10)	ybaY/ybaZ	outer membrane lipoprotein/excision repair protein, alkyltransferase‑like protein ATL
Reads supporting (aligned to +/- strand): ref base C (94/109); new base G (5/8); total (100/117)
Fisher's exact test for biased strand distribution p-value = 7.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.21e-02

Only 100 of 221 total aligned reads displayed.

AGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGGAGACGCTTTAACACACCGCCCACCTGGCGCGCGGCGCGGGGCGATCCCGCCAGTTTCGCCACATCACCGTAAGTGGTGACATAGCCTTCGGGA  >  NC_000913/475702‑476240

agtgttttaGCCGAACTCGCGTTATCTTTTGAGGTCGGCTATTACCATGAATGACAAACTGGTCTTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGga                                                                                                                                                                                                                                                                          <  1:185494/267‑1 (MQ=255)

 gtggtgaaGCCTAACGCGTGTATTATGGTAGGTGCGGCGAATACATGGAATGACAAACTGGTATTTTACGCCGATACCGGTCAGCCGGTGAACAACCAGGGCGGAACTAAAGCCGACCTGACTTGGGTGCGGGTACAGCAAAGGGCAGTGCCGGTGCAGGCCAGTGGTGCCGCAACGACGACCGTACCTTCGACTTCACCAACTCAGGTGTAGCCGTCTTCGGCAGTTCCAGCTGGAACGCAATATGAA‑GCGACGTTAACCCTCTCCGCCCGgag                                                                                                                                                                                                                                                                         <  4:472425/268‑1 (MQ=255)

    gTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACACACTGGTATTTACCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGAACTGACATTGGTGCAGGTACAGCAAACCGCCGTGACGGTTCAGGCCAGCGGGTGCGCAACGACTACCGTACCTTCGACGACACCAACTCAGGTGAAGCGGACGTCGTCAGTTCCAGCCCCTACACCATATTAACGCGAC‑TACACCCTCCCCGCCCGGAgggg                                                                                                                                                                                                                                                                      >  4:183381/1‑275 (MQ=255)

              aCGCGCGTATTCTGTGGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGGACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGtt                                                                                                                                                                                                                                                            <  1:210717/275‑1 (MQ=255)

               cgcgcgTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGtt                                                                                                                                                                                                                                                            >  2:199873/1‑274 (MQ=255)

                                                      cAAACTGGTATGTATCACCCATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTGTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGa                                                                                                                                                                                                                     <  2:288742/274‑1 (MQ=255)

                                                              tATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACGCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCCATACCTTccc                                                                                                                                                                                                            >  4:136394/1‑273 (MQ=255)

                                                                   aTCACCGATCCGGTCCAGGCGGGGCTCAACCATGGAGCAAGTAAAGCGGACCTGATATTGGGGCCGGTACAGAAAACCCCCGTCCCGGTTCAGGCCAGCGGGGGCGCAAGGACTACCGTACCTTCGACTTCCCCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCGCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCt                                                                                                                                                                                                        <  4:355861/274‑1 (MQ=255)

                                                                     cACCGATACCGGTCAGCCGGGGATCAGCCGGGGCGGAACTAAAGCCTCCCTGAGATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCt                                                                                                                                                                                                     <  4:199873/275‑1 (MQ=255)

                                                                          aTACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAGTATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGc                                                                                                                                                                                                   >  1:95095/1‑272 (MQ=255)

                                                                                          gATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATACCGCTGCAAGCCGATTTGCCCGCTTCCCGATACCANNNCCCCTTCTCCCAGTAATGCCGTTCGCtc                                                                                                                                                                                 >  3:195084/1‑273 (MQ=255)

                                                                                                         aaCTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTg                                                                                                                                                                                                                                             >  4:650910/1‑199 (MQ=255)

                                                                                                         aaCTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccggt                                                                                                                                                                 >  2:284759/1‑275 (MQ=255)

                                                                                                                                 gCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGagag                                                                                                                                                                                                                                                                       >  8:975295/1‑149 (MQ=255)

                                                                                                                                 gCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGagag                                                                                                                                                                                                                                                                       <  5:850707/149‑1 (MQ=255)

                                                                                                                                  ccGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAgg                                                                                                                                                                                                                                                                      <  6:328359/149‑1 (MQ=255)

                                                                                                                                  ccGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAgg                                                                                                                                                                                                                                                                      <  5:871306/149‑1 (MQ=255)

                                                                                                                                  ccGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAgg                                                                                                                                                                                                                                                                      <  5:996556/149‑1 (MQ=255)

                                                                                                                                   cGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGGGCGGAGAggg                                                                                                                                                                                                                                                                     <  6:1093990/149‑1 (MQ=255)

                                                                                                                                   cGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAggg                                                                                                                                                                                                                                                                     >  7:744384/1‑149 (MQ=255)

                                                                                                                                         aGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCCCTGACGCTGTAAATCCGGTCCGGTTAGCGCAATTTGGCCCTGGCGATTAac                                                                                                                                 >  4:303982/1‑275 (MQ=255)

                                                                                                                                           cAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTagtagt                                                                                                                                                                                                                                                             <  5:654464/149‑1 (MQ=255)

                                                                                                                                                cgccgTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCgg                                                                                                                          >  2:26803/1‑275 (MQ=255)

                                                                                                                                                     tGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCa                                                                                                                     <  2:232448/275‑1 (MQ=255)

                                                                                                                                                     ggCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCa                                                                                                                     <  2:645309/274‑1 (MQ=255)

                                                                                                                                                     ggCCGGTTCAGGCCAGCGGTGGCCCAACGACTACCTTACCTTCGACTTCACCAACTCAGGTGAATCCGACTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGcc                                                                                                                      <  4:438541/273‑1 (MQ=255)

                                                                                                                                                          gttacggCCAGCGGTGGCGCAACGACTACCTTACCTTCGACTTCCCCAACTCAGCTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTa                                                                                                                <  2:557116/270‑1 (MQ=255)

                                                                                                                                                          gTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGTGAggg                                                                                                                                                                                                                                                                     >  2:485027/1‑126 (MQ=255)

                                                                                                                                                            tCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGc                                                                                                                                                                                                                                            <  6:875150/149‑1 (MQ=255)

                                                                                                                                                              aGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGTGAGGGTtag                                                                                                                                                                                                                                                                 <  4:485027/126‑1 (MQ=255)

                                                                                                                                                                gCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGtgc                                                                                                          <  4:26803/275‑1 (MQ=255)

                                                                                                                                                                    ggggTGGGGCAACGACTACCGTACCTTCCACTTCACCAACTCAGGGGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGtgctgc                                                                                                       <  4:284759/272‑1 (MQ=255)

                                                                                                                                                                      ggTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGAttt                                                                                                                                                                                                                                  <  7:816761/149‑1 (MQ=255)

                                                                                                                                                                        tGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCTCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGc                                                                                                                                                                                                                                <  8:488295/149‑1 (MQ=255)

                                                                                                                                                                        tGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGc                                                                                                                                                                                                                                <  7:975295/149‑1 (MQ=255)

                                                                                                                                                                            gcAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGc                                                                                                                                                                                                                            >  7:908996/1‑149 (MQ=255)

                                                                                                                                                                            gcAACGACTACCGTAACTTCGACTTCACCAACTAAGGTGAATCGGTATTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCACGAAGAGGGTTAGTAGTTCCAGCGATACCGCTGCAAGGCGATTTGCCCGc                                                                                                                                                                                                                            >  6:758589/1‑149 (MQ=255)

                                                                                                                                                                               aCGACTACCTGACCTGCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCAAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGCGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGga                                                                                              <  4:66516/272‑1 (MQ=255)

                                                                                                                                                                               aCGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGgaga                                                                                            <  1:117159/274‑1 (MQ=255)

                                                                                                                                                                                    tACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGAt                                                                                                                                                                                                                    <  6:742760/149‑1 (MQ=255)

                                                                                                                                                                                      ccGTACCTTCGACTTCACCAACGCAGGTGCATCCGTCTTCGGCAGTTCCTGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTGCCAGCGATAACTCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGGAGACGCTTTa                                                                                     <  2:183381/274‑1 (MQ=255)

                                                                                                                                                                                           ccTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATca                                                                                                                                                                                                             <  7:410177/149‑1 (MQ=255)

                                                                                                                                                                                            cTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGAt                                                                                                                                                                                                                    >  5:302994/1‑141 (MQ=255)

                                                                                                                                                                                                gACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGGAGACGCTTTAACACACCGccc                                                                          >  2:370452/1‑275 (MQ=255)

                                                                                                                                                                                                  cgtCACCAACTCAGGTGAGTCCGTCTTCGGCAGGTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGc                                                                                                                                                                                                                            <  3:37758/125‑1 (MQ=255)

                                                                                                                                                                                                   ttCACCAACTCCGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCt                                                                                                                                                                                                     >  7:613904/1‑149 (MQ=255)

                                                                                                                                                                                                     cACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATAcc                                                                                                                                                                                                                 >  2:187665/1‑135 (MQ=255)

                                                                                                                                                                                                      aCCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGcc                                                                                                                                                                                                  >  5:77816/1‑149 (MQ=255)

                                                                                                                                                                                                           cTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTaa                                                                                                                                                                                             >  6:474986/1‑149 (MQ=255)

                                                                                                                                                                                                              aGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGGGCGgaga                                                                                                                                                                                                                                                                        <  8:128638/71‑1 (MQ=39)

                                                                                                                                                                                                              aGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGGGCGgaga                                                                                                                                                                                                                                                                        >  7:128638/1‑71 (MQ=39)

                                                                                                                                                                                                               ggTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGcc                                                                                                                                                                                         <  8:613904/149‑1 (MQ=255)

                                                                                                                                                                                                                    aTCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCg                                                                                                                                                                                    >  8:253472/1‑149 (MQ=255)

                                                                                                                                                                                                                     tCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCTGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGGAGACGCTTTAACACACCGCCCACCTGgcgcg                                                                <  2:577262/264‑1 (MQ=255)

                                                                                                                                                                                                                     tCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGc                                                                                                                                                                                   >  5:791458/1‑149 (MQ=255)

                                                                                                                                                                                                                     tCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGc                                                                                                                                                                                   <  7:728923/149‑1 (MQ=255)

                                                                                                                                                                                                                     tCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGc                                                                                                                                                                                   >  5:1080373/1‑149 (MQ=255)

                                                                                                                                                                                                                      ccGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGAtt                                                                                                                                                                                                                                   >  1:40638/1‑100 (MQ=255)

                                                                                                                                                                                                                      ccGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCt                                                                                                                                                                                  <  6:593313/149‑1 (MQ=255)

                                                                                                                                                                                                                      acGTCTTCGGCAGTTCCAACTCCTAGACAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCt                                                                                                                                                                                  >  6:392549/2‑149 (MQ=255)

                                                                                                                                                                                                                       cGTCTTCGGCAATTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGc                                                                                                                                                                                                                                               >  7:613151/1‑87 (MQ=255)

                                                                                                                                                                                                                        gTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGa                                                                                                                                                                                <  6:1070305/149‑1 (MQ=255)

                                                                                                                                                                                                                         tCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGAc                                                                                                                                                                               >  8:756316/1‑149 (MQ=255)

                                                                                                                                                                                                                               gCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGAttt                                                                                                                                                                                                                                  <  6:862481/92‑1 (MQ=255)

                                                                                                                                                                                                                               gCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCAGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTa                                                                                                                                                                         >  6:803065/1‑149 (MQ=255)

                                                                                                                                                                                                                                 aGTTCCAGCTCCTACGCAATATTTA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATATCGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTaaa                                                                                                                                                                       <  5:721246/149‑1 (MQ=255)

                                                                                                                                                                                                                                   ttCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCAt                                                                                                                                                                                                               <  5:386432/107‑1 (MQ=255)

                                                                                                                                                                                                                                     ccAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccg                                                                                                                                                                   <  5:876679/149‑1 (MQ=255)

                                                                                                                                                                                                                                      cAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCATTAATGCCTGTCGCTGACGCTGTAAAtccgg                                                                                                                                                                  >  8:642271/1‑149 (MQ=255)

                                                                                                                                                                                                                                       aGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccggt                                                                                                                                                                 <  6:380109/149‑1 (MQ=255)

                                                                                                                                                                                                                                       aGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGc                                                                                                                       <  2:651550/191‑1 (MQ=255)

                                                                                                                                                                                                                                          tCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATAcc                                                                                                                                                                                                                 >  5:642629/1‑98 (MQ=255)

                                                                                                                                                                                                                                              aCGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTa                                                                                                                                                          <  7:1022434/149‑1 (MQ=255)

                                                                                                                                                                                                                                                gCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGc                                                                                                                                                        <  7:554141/149‑1 (MQ=255)

                                                                                                                                                                                                                                                 cAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGc                                                                                                                                                        <  5:1061651/148‑1 (MQ=255)

                                                                                                                                                                                                                                                  aaTATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGa                                                                                                                                                      >  5:267289/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    tatTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGtt                                                                                                                                                           >  4:242566/1‑142 (MQ=255)

                                                                                                                                                                                                                                                    tatTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGaaa                                                                                                                                                    <  8:554836/149‑1 (MQ=255)

                                                                                                                                                                                                                                                      ttAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCACGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAAtt                                                                                                                                                  >  8:159345/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        aa‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTGATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATtgt                                                                                                                                                <  7:469312/149‑1 (MQ=255)

                                                                                                                                                                                                                                                        aa‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGgaga                                                                                            <  2:703333/201‑1 (MQ=255)

                                                                                                                                                                                                                                                         a‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATtgtg                                                                                                                                               >  5:1028581/1‑149 (MQ=255)

                                                                                                                                                                                                                                                           gCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGc                                                                                                                                              <  7:303644/149‑1 (MQ=255)

                                                                                                                                                                                                                                                             gACGTTAACCCTCGCCGCCCGGAGAGGGTTAGTAGTTCCAGCCATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGGAGACGCTTTAATACACCGCCCACCTGGCGCGCGGCGCGGGGCGATCCCGCCAGTTTCGCCACATCACCGTAAGTGGTGACa              <  2:136394/275‑1 (MQ=255)

                                                                                                                                                                                                                                                              aCGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGt                                                                                                                                           <  8:541418/149‑1 (MQ=255)

                                                                                                                                                                                                                                                               cGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTg                                                                                                                                                                                                                                 >  3:100692/1‑62 (MQ=255)

                                                                                                                                                                                                                                                               cGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCt                                                                                                                                                                                                        >  4:34233/1‑87 (MQ=255)

                                                                                                                                                                                                                                                                 ttAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGc                                                                                                                                        <  6:18176/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                  tacacccctcCGCCCGGGGGTGGTTACTNGNTCCCCCGATAGCGCCGCACGTCGACTTGCCCGCTTCCCGATGCCGNAACACCTTCTGCCAGGAATGCCTGTCGCGGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGNCGATTAACCCACCGGTGACAGGTTAAGGTGCTGCCTTCAGGGAGACGCTTTAACACACCGCCCGCCTGGCGCGCGGCGCGGGGCGATCCCGCCAGTTTCGCCACATCACCGTAAGTGGTGACATAGcc         <  3:55707/268‑1 (MQ=255)

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

                                                                                                                                                                                                                                                                   aaCCCTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTACCGAAATTGTGCCGTGTCGa                                                                                                                                      >  8:78890/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                   aaCCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCg                                                                                                                                                                                                                             <  1:100692/62‑1 (MQ=255)

                                                                                                                                                                                                                                                                   aaCCCCCTCGGCCCGGAGAGGGTTAGGTGTGCCAGCGATAACGTTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAAGGCCTGTAGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCCCCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGga                                                                                              <  4:730317/189‑1 (MQ=255)

                                                                                                                                                                                                                                                                    aCCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACATTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGTAATTGTGCCGTGGCGAt                                                                                                                                     >  8:305027/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     cccTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGAtt                                                                                                                                    >  5:436053/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     cccTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGAtt                                                                                                                                    >  8:78133/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                        tctcCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccg                                                                                                                                                                   <  7:638985/115‑1 (MQ=255)

                                                                                                                                                                                                                                                                         ctcCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAgg                                                                                                >  2:250895/1‑181 (MQ=255)

                                                                                                                                                                                                                                                                         ctcCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGGAGACGCTTTAACACACCGCCCACCTGGCGCGCGGCCCGGGGCGATCCCGCCCGTTTCGCCACATCACCTGAAGTGGGGACATAGCCCTAgggg  >  4:583957/1‑274 (MQ=255)

                                                                                                                                                                                                                                                                             gCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAaccacc                                                                                                                             <  2:327961/148‑1 (MQ=255)

AGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAA‑GCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGTGCCAGGGTAAGGTGCTGCCTTCAGGGAGACGCTTTAACACACCGCCCACCTGGCGCGCGGCGCGGGGCGATCCCGCCAGTTTCGCCACATCACCGTAAGTGGTGACATAGCCTTCGGGA  >  NC_000913/475702‑476240

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