BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,708,595	G→A	5.4%	intergenic (‑890/+21)	dppA ← / ← proK	dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor/tRNA‑Pro

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,708,595	0	G	A	5.4%	530.6 / 12.5	163	intergenic (‑890/+21)	dppA/proK	dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor/tRNA‑Pro
Reads supporting (aligned to +/- strand): ref base G (82/72); new base A (4/5); total (86/77)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-02

Only 100 of 169 total aligned reads displayed.

ATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGGGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATA  >  NC_000913/3708327‑3708864

aTGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCCCATGAGCTGCAGATGCAAAAAAGCCTGATCGTTTAGc                                                                                                                                                                                                                                                                         >  2:176945/1‑275 (MQ=255)

                         gCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCTGGCCGTTGCTGCCGCCCTTTTCTTTCGCCTCGCTCCAGTCCAACCTTCATTGCCCCCGCGGTTCCCATCCCCATGAGCCGCACATGCCAAAAAGCCTGCACGTTGAGGCGGCATTTCCAATTTGGTAGGttaa                                                                                                                                                                                                                                                >  4:124600/1‑272 (MQ=18)

                           cTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCCCGCTTCGCCCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAAGTTGGTCGGTGATa                                                                                                                                                                                                                                               >  2:497732/1‑274 (MQ=255)

                                   agtcgaACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGt                                                                                                                                                                                                                                                                              >  4:510341/1‑235 (MQ=255)

                                       gaACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAg                                                                                                                                                                                                                                                                          <  2:510341/235‑1 (MQ=255)

                                                                  tCCCCATGAGCTGCACATGAGAAAGAGCCTGCTCGTTGAGCAGGCTTTTAGAATGTGGTCGGGGATAGAGTATTACTCGCCACTTCGTGGCTCGCCCTGGGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTg                                                                                                                                                                                                       <  1:62038/275‑1 (MQ=37)

                                                                            cTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAgg                                                                                                                                                                                             <  2:450414/275‑1 (MQ=255)

                                                                                                        aGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGAGTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAGAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCa                                                                                                                                                                 <  2:162808/275‑1 (MQ=255)

                                                                                                                  tCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCCAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTTGTCTGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCTGGGCGCCTCTTACTgc                                                                                                                                                        >  2:422788/1‑274 (MQ=255)

                                                                                                                          ggTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGtt                                                                                                                                                                                                                                                                             <  5:910490/149‑1 (MQ=255)

                                                                                                                             cGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAg                                                                                                                                                                                                                                                                          <  5:270631/149‑1 (MQ=255)

                                                                                                                                tGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAgg                                                                                                                                                                                                                                                                      >  2:207551/1‑150 (MQ=255)

                                                                                                                                  aTAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGc                                                                                                                                                                                                                                                                     <  7:753466/149‑1 (MQ=255)

                                                                                                                                    agagGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCtt                                                                                                                                                                                                                                                                   <  7:527736/149‑1 (MQ=255)

                                                                                                                                    agagGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCAACCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCtc                                                                                                                                                                                                                                                                   <  6:935348/149‑2 (MQ=255)

                                                                                                                                      agGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCTCTCGAGTCGAACCTCCCTGCCCCCGGAGGATCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCtttt                                                                                                                                                                                                                                                                 >  5:52172/1‑149 (MQ=255)

                                                                                                                                      agGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTACTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGACGCAAAAAAGCCTGCTCAACGAGCAGGCtttt                                                                                                                                                                                                                                                                 <  8:662666/149‑1 (MQ=25)

                                                                                                                                          ttACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCg                                                                                                                                                                        >  4:192864/1‑234 (MQ=255)

                                                                                                                                              tCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGGGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTg                                                                                                                                                                                                       <  2:12238/199‑1 (MQ=255)

                                                                                                                                              tCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTg                                                                                                                                                                                                                                                         <  7:253171/149‑1 (MQ=255)

                                                                                                                                              tCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTg                                                                                                                                                                                                                                                         >  7:887604/1‑149 (MQ=255)

                                                                                                                                              tCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGgc                                                                                                                                                                    <  2:192864/234‑1 (MQ=255)

                                                                                                                                              tCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCGCGTTGAGCAGGCTTTTCGAATTTg                                                                                                                                                                                                                                                         >  6:828676/1‑149 (MQ=255)

                                                                                                                                              tCGCCACTGCGTGGCTCGCCCGCCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTg                                                                                                                                                                                                                                                         <  8:505334/149‑1 (MQ=255)

                                                                                                                                                    cTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGt                                                                                                                                                                                                                                                   <  6:27646/149‑1 (MQ=255)

                                                                                                                                                        gTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATa                                                                                                                                                                                                                                               >  8:104768/1‑149 (MQ=255)

                                                                                                                                                        gTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATa                                                                                                                                                                                                                                               >  7:585239/1‑149 (MQ=255)

                                                                                                                                                          ggCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATaga                                                                                                                                                                                                                                             <  5:647943/149‑1 (MQ=255)

                                                                                                                                                             tCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCAACGAGCAGGCtttt                                                                                                                                                                                                                                                                 <  8:155182/126‑1 (MQ=255)

                                                                                                                                                               gCCCTGCGGGCCGTTGCTGGCGCAGCGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGAtt                                                                                                                                                                                                                                        <  6:719796/149‑1 (MQ=255)

                                                                                                                                                               gCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGCATTTGGTCGGTGAGAGAGGCTTCGAACCTCCGACCCCTTCGTCCCCCACGAAATCCCCTCCCAGGCTGCCNCAATCCCCCAATGCGGGGCCACACTTCCTGCGCAGATACACCCTCCTCACTCCCTTAACAACAAAACACCTTTTTAt                                                                                                          >  3:160620/1‑275 (MQ=21)

                                                                                                                                                               gCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGAt                                                                                                          >  2:336823/1‑275 (MQ=255)

                                                                                                                                                                cccTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGCGCTGCTGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTc                                                                                                                                                                                                                                       >  6:341187/1‑149 (MQ=255)

                                                                                                                                                                 ccTGCGGGCCGTTGCTGGGTCAACGTTCTCTCGCTTCGCTCGAGTCTAACCTCCCTGCCCCCGGGGGTTCTAATCCTCATGGGCTACAGATGCAAAAAAGCGTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCg                                                                                                        <  4:497732/275‑1 (MQ=255)

                                                                                                                                                                  cTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGa                                                                                                                                                                                                                                     >  8:1000464/1‑149 (MQ=255)

                                                                                                                                                                   tGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAg                                                                                                                                                                                                                                                                          <  2:133828/111‑1 (MQ=255)

                                                                                                                                                                    gCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAAc                                                                                                                                                                                                                                   <  7:655377/149‑1 (MQ=255)

                                                                                                                                                                        gccctgcgtGGCGCCTCGTTCGCTCGCTTCGCTTGCGTCCACCCTCCCTCCCCCCCGAGGTTCTCATCCTCATGTGCTGCAGATGCAAAAACGCCTGCTCGTTGCGCAGGCTTTTCGAATTTGGTCGGTGATAGATGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGaa                                                                                                 <  1:102012/267‑1 (MQ=37)

                                                                                                                                                                         ccGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCg                                                                                                                                                                                                                              <  6:620895/149‑1 (MQ=255)

                                                                                                                                                                               cTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATGGAGGATTCGAACCTCCGACCCCt                                                                                                                                                                                                                        <  7:1000464/149‑1 (MQ=255)

                                                                                                                                                                                 ggCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTc                                                                                                                                                                                                                      <  5:127736/149‑1 (MQ=255)

                                                                                                                                                                                   cgtaccGTTCTGTCGCTTTGGTTGGGTCGAACCTCCCTGGACCCGGTGGTTCTCATCGTCATGAGCTGGAGATGGAAAAAAGCCTGGTCGTTGGTCAGGCTTTTCGAAGTTGGTCGGTGGTAGAGGATTCGAACCTCCGACCCCTTCTGCCCGAAGGAAGTGCGCTACCAGGCTGCGCCAAGCACCGAATGCGGGGCGACTCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAgg                                                                                       <  4:629403/269‑1 (MQ=255)

                                                                                                                                                                                    gcAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTc                                                                                                                                                                                                                   >  8:32223/1‑149 (MQ=255)

                                                                                                                                                                                      aaCGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGaa                                                                                                                                                                                                                                    <  4:144358/130‑1 (MQ=255)

                                                                                                                                                                                      aaCGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTccc                                                                                                                                                                                                                 >  8:316462/1‑149 (MQ=255)

                                                                                                                                                                                         gTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCcgaa                                                                                                                                                                                                              <  6:1027586/149‑1 (MQ=255)

                                                                                                                                                                                             tctcGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCctaacgaa                                                                                                                                                                                                          <  7:104768/149‑1 (MQ=255)

                                                                                                                                                                                                  cttcgctCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTgcg                                                                                                                                                                                                     >  5:452854/1‑149 (MQ=255)

                                                                                                                                                                                                  cttcgctCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCCCCGACCCCTTCGTCCcgaacga                                                                                                                                                                                                           >  1:81467/1‑143 (MQ=255)

                                                                                                                                                                                                      gctCNGGTCGAACCTCCCTGCCCCCGGTGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTg                                                                                                                                                                                                       <  3:81467/143‑1 (MQ=255)

                                                                                                                                                                                                       gaggatgagAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCACTCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTGGCCGCTGCct                                                                  >  4:262564/9‑275 (MQ=255)

                                                                                                                                                                                                       ctCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCCCCGAATGCAGGGCACATCTTACTTCGCAGCTACGCCCTCCTCAATCCCTTAATAGCAAAAATCCTTTTGATCGGCGAAAACGTTAGCAGGCCACCTCGTTAGCCCCCGccc                                                                  >  4:11573/1‑274 (MQ=25)

                                                                                                                                                                                                       ctCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCgggg                                                                                                                                                                     >  1:67894/1‑176 (MQ=255)

                                                                                                                                                                                                           aGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGc                                                                                                                                                                                            <  7:940205/149‑1 (MQ=255)

                                                                                                                                                                                                               gAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcg                                                                                                                                                                                        >  8:834944/1‑149 (MQ=255)

                                                                                                                                                                                                                aaCCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcgc                                                                                                                                                                                       >  6:26138/1‑149 (MQ=255)

                                                                                                                                                                                                                aaCCTCCCTGCCCCCAGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgc                                                                                                                                                                                         <  5:983587/147‑1 (MQ=255)

                                                                                                                                                                                                                  ccTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCa                                                                                                                                                                                     >  8:456889/1‑149 (MQ=255)

                                                                                                                                                                                                                   cTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCaa                                                                                                                                                                                    >  6:859074/1‑149 (MQ=255)

                                                                                                                                                                                                                    tCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAAt                                                                                                                                                                                   >  8:763161/1‑149 (MQ=255)

                                                                                                                                                                                                                          cccccGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGcc                                                                                                                                             <  2:55281/181‑1 (MQ=255)

                                                                                                                                                                                                                            cccGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAAt                                                                                                                                                                           <  7:495302/149‑1 (MQ=255)

                                                                                                                                                                                                                             ccGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATg                                                                                                                                                                          >  8:110411/1‑149 (MQ=255)

                                                                                                                                                                                                                              cGGAGGTTCTCCTCCTCATGGGCGGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCGTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGGGATGTGacac                                           <  4:422788/275‑1 (MQ=255)

                                                                                                                                                                                                                                gAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGt                                                                                                                                        >  2:204961/1‑180 (MQ=255)

                                                                                                                                                                                                                                 aGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCggg                                                                                                                                                                      <  5:292538/149‑1 (MQ=255)

                                                                                                                                                                                                                                    ttCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAAt                                                                                                                                    <  4:204961/180‑1 (MQ=255)

                                                                                                                                                                                                                                      ctcATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCa                                                                                                                                                                 >  6:914414/1‑149 (MQ=255)

                                                                                                                                                                                                                                       tcATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCAt                                                                                                                                                                >  8:549085/1‑149 (MQ=255)

                                                                                                                                                                                                                                         aTCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGaa                                                                                                                                                                            <  7:452749/135‑1 (MQ=255)

                                                                                                                                                                                                                                           ccTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTa                                                                                                                                                            >  7:757469/1‑149 (MQ=255)

                                                                                                                                                                                                                                           ccTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTa                                                                                                                                                            <  7:274846/149‑1 (MQ=255)

                                                                                                                                                                                                                                            cTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTAc                                                                                                                                                           <  7:316462/149‑1 (MQ=255)

                                                                                                                                                                                                                                            cTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTAc                                                                                                                                                           >  5:397099/1‑149 (MQ=255)

                                                                                                                                                                                                                                             tCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACt                                                                                                                                                          <  8:897258/149‑1 (MQ=255)

                                                                                                                                                                                                                                             tCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACt                                                                                                                                                          >  5:1069779/1‑149 (MQ=255)

                                                                                                                                                                                                                                             tCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACt                                                                                                                                                          <  5:859074/149‑1 (MQ=255)

                                                                                                                                                                                                                                              cATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTg                                                                                                                                                         >  8:477760/1‑149 (MQ=255)

                                                                                                                                                                                                                                               aTGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCt                                                                                                                                           >  2:741861/1‑162 (MQ=255)

                                                                                                                                                                                                                                                   gCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTc                                                                                                                                       <  4:741861/162‑1 (MQ=255)

                                                                                                                                                                                                                                                    cTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGGGATGTGACACCAGTTGTTGTTTTCGTTAATTc                     >  4:472364/1‑275 (MQ=255)

                                                                                                                                                                                                                                                       cAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATAc                                                                                                                                                <  7:800523/149‑1 (MQ=255)

                                                                                                                                                                                                                                                       cAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGGGATGTGACACCAGTTGTTGTTTTCGTTAATTcca                   <  2:472364/274‑1 (MQ=255)

                                                                                                                                                                                                                                                          aTGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGcc                                                                                                                                             >  6:431371/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          aTGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGcc                                                                                                                                             >  5:866674/1‑149 (MQ=255)

                                                                                                                                                                                                                                                              aaaaaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCg                                                                                                                                         >  7:122213/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               aaaaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTTCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGt                                                                                                                                        >  6:761351/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               aaaaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGGGATGTGACACAAGTTGTTGTTTTCGTGAATTCCACCAGCCGGt          >  2:352464/1‑275 (MQ=255)

                                                                                                                                                                                                                                                                aaaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTc                                                                                                                                       <  6:452854/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                 aaaGCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCa                                                                                                                                      >  8:513349/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                 aaaGCCTGATTGTTGAGGAGGCTTTGCGAATTTGGTCGGTGTTAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAAGGCCTTTTGATCGGCGGGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGGGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTg         <  4:434862/274‑1 (MQ=255)

                                                                                                                                                                                                                                                                  aaGCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCaa                                                                                                                                     >  8:39573/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     ccTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATcc                                                                                                                                  >  7:546741/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                       tGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCATTCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCt                                                                                                                                <  7:1044401/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                        gCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGGGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGtata  <  4:336823/274‑1 (MQ=255)

                                                                                                                                                                                                                                                                          tCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTaa                                                                                                                             <  8:1034343/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                            gTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGGTTCGAACCTCCGACGCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATa                                                                                                                           >  7:687523/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                            gTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATa                                                                                                                           >  8:650714/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                            gTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATa                                                                                                                           >  6:144898/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                            gTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCCTCAATCCCTTAATa                                                                                                                           <  5:828676/149‑1 (MQ=255)

ATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGCAAAATGCCTTTTGATCGGCGAGAAAGTCAGCAGGCCGCTTAGTTAGCCGCTGCCTCTTTTGCCTGCGGGATGTGACACCAGTTGTTGTTTTCGTTAATTCCACCATCCGGTGAAGTATA  >  NC_000913/3708327‑3708864

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