BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,395,454	T→C	5.0%	C463C (TGT→TGC)	nnr →	bifunctional NAD(P)H‑hydrate repair enzyme; C‑terminal domain ADP‑dependent (S)‑NAD(P)H‑hydrate dehydratase and N‑terminal domain NAD(P)H‑hydrate epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,395,454	0	T	C	5.0%	571.5 / 10.3	172	C463C (TGT→TGC)	nnr	bifunctional NAD(P)H‑hydrate repair enzyme; C‑terminal domain ADP‑dependent (S)‑NAD(P)H‑hydrate dehydratase and N‑terminal domain NAD(P)H‑hydrate epimerase
Reads supporting (aligned to +/- strand): ref base T (77/84); new base C (6/3); total (83/89)
Fisher's exact test for biased strand distribution p-value = 3.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.22e-01

Only 100 of 179 total aligned reads displayed.

GCCGCACGGTTGTTAGGCTGTTCCGTCGCTGAAATTGAAAGTGACCGCTTACATTGCGCCAAACGTCTGGTACAACGTTATGGCGGCGTAGCGGTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGG  >  NC_000913/4395200‑4395723

                                                                                     gcgTAGCGGTGCTGAAAGGTTCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCCTTGGTATTATTGGCGCATTTCTTGGCCAAAAACTGTCGCAGTATGATGCAGCCCGTGCAGCCTGTGTCGCGCACGGTGCCGCAGCTGACGTACTGGCCGCGCGCTTTCGATATCCCTGGGTGCGTGCCCACGCTCTCTTGTCCACCCTACCGGGTATtttt                                                                                                                                                                      >  4:167020/1‑275 (MQ=255)

                                                                                             gTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCCGCCTGTGCAGGCTGTGTCGCGCACGGTGGGGCAGCTGACGTAGTGGCCGCGCGTGTTGGAACGCGCGGGATGCTGGCAAACGATCTCTTTGCCCCCCGACAGAGACACGTTAACCCGGa                                                                                                                                                              >  4:77684/1‑275 (MQ=255)

                                                                                                  gAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGa                                                                                                                                                         <  2:77684/275‑1 (MQ=255)

                                                                                                          ccGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACgg                                                                                                                                                                                                                                                           >  4:155233/1‑169 (MQ=255)

                                                                                                           cGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCtgtg                                                                                                                                                                                                                                                                              >  8:681200/1‑149 (MQ=255)

                                                                                                            ggAATCGTGTTCGGCGCCCATCCGGACGCTTTAGGCAGTATTGATGCCGGAAATGCA‑GCATGGCGAGCGGCGGCATTGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGAAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATaaaga                                                                                                                                               <  2:323928/274‑3 (MQ=255)

                                                                                                            ggAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGc                                                                                                                                                                                                                                                                             >  8:163319/1‑149 (MQ=255)

                                                                                                              aaCCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCg                                                                                                                                                                                                                                                       <  2:155233/169‑1 (MQ=255)

                                                                                                                ccGTGGTCGCCGCCCACCCTGACGGTTTAGGCATTATTGTTGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCACCCTGTGCAGGc                                                                                                                                                                                                                                                                         <  7:1280086/149‑1 (MQ=255)

                                                                                                                  gTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGGGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGCACAGGCtg                                                                                                                                                                                                                                                                       <  8:161974/149‑1 (MQ=255)

                                                                                                                  gTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCtg                                                                                                                                                                                                                                                                       <  6:208006/149‑1 (MQ=255)

                                                                                                                  gTGGTCGCCCCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCtg                                                                                                                                                                                                                                                                       >  6:976951/1‑149 (MQ=255)

                                                                                                                   tGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTCCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCtgt                                                                                                                                                                                                                                                                      <  5:985154/149‑1 (MQ=255)

                                                                                                                    ggTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCtgtg                                                                                                                                                                                                                                                                     >  5:1176473/1‑149 (MQ=255)

                                                                                                                      tCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTc                                                                                                                                                                                                                                                                   <  8:1263921/149‑1 (MQ=255)

                                                                                                                      tCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTc                                                                                                                                                                                                                                                                   <  6:378546/149‑1 (MQ=255)

                                                                                                                        gccgccCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTcgc                                                                                                                                                                                                                                                                 <  7:316774/149‑1 (MQ=255)

                                                                                                                        gccgccCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCCGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTcgc                                                                                                                                                                                                                                                                 <  7:1015041/149‑1 (MQ=255)

                                                                                                                           gccCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCa                                                                                                                                                                                                                                                              >  8:1217776/1‑149 (MQ=255)

                                                                                                                             ccATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGCGCAGGCTGTGTCGCGCACg                                                                                                                                                                                                                                                            >  6:309894/1‑149 (MQ=255)

                                                                                                                                 ccTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCa                                                                                                                                                                                                                                                              >  6:1062770/1‑143 (MQ=255)

                                                                                                                                 ccTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCa                                                                                                                                                                                                                                                              <  5:1062770/143‑1 (MQ=255)

                                                                                                                                   tGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCgg                                                                                                                                                                                                                                                      <  5:578801/149‑1 (MQ=255)

                                                                                                                                    gACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCATAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGc                                                                                                                                                                                                                                                     >  8:892705/1‑149 (MQ=255)

                                                                                                                                      cGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAg                                                                                                                                                                                                                                                   <  7:99216/149‑1 (MQ=255)

                                                                                                                                        cTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCt                                                                                                                                                                                                                                                 <  6:189380/149‑1 (MQ=255)

                                                                                                                                         tttAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGAGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTg                                                                                                                                                                                                                                                >  5:1000295/1‑149 (MQ=255)

                                                                                                                                          ttAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGa                                                                                                                                                                                                                                               >  6:72029/1‑149 (MQ=255)

                                                                                                                                          ttAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGCTGAATCGAGTAATTCCGCTCCCtg                                                                                                                 >  2:160051/1‑275 (MQ=255)

                                                                                                                                           tAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCtga                                                                                                                >  2:240985/1‑275 (MQ=255)

                                                                                                                                              gCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTa                                                                                                                                                                                                                                           <  7:649204/149‑1 (MQ=255)

                                                                                                                                              gCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTa                                                                                                                                                                                                                                           >  8:1249230/1‑149 (MQ=255)

                                                                                                                                              gCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTa                                                                                                                                                                                                                                           >  7:157497/1‑149 (MQ=255)

                                                                                                                                              gCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTa                                                                                                                                                                                                                                           >  7:413957/1‑149 (MQ=255)

                                                                                                                                               cATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTAc                                                                                                                                                                                                                                          <  5:194864/149‑1 (MQ=255)

                                                                                                                                                  tattGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCAAGATGAATCGAGTAATTCCGCTCCCTGATGAGCa                                                                                                          >  2:301555/1‑274 (MQ=255)

                                                                                                                                                  tattGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGCGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTgg                                                                                                                                                                                                                                       <  5:309894/149‑1 (MQ=255)

                                                                                                                                                    ttGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTggcg                                                                                                                                                                                                                                     >  7:1297182/1‑149 (MQ=255)

                                                                                                                                                     tgcggCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATTTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGc                                                                                                       <  4:301555/270‑1 (MQ=255)

                                                                                                                                                            ggAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGttt                                                                                                                                                                                                                             >  7:1076247/1‑149 (MQ=255)

                                                                                                                                                            ggAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGttt                                                                                                                                                                                                                             >  8:1035670/1‑149 (MQ=255)

                                                                                                                                                             gAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGtttt                                                                                                                                                                                                                            >  8:1417827/1‑149 (MQ=255)

                                                                                                                                                             gAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGtt                                                                                                                                                                                                                              >  6:1215112/1‑147 (MQ=255)

                                                                                                                                                              aaaTGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTg                                                                                                                                                                                                                           <  5:407533/149‑1 (MQ=255)

                                                                                                                                                              aaaTGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCCGCGCGTTTTGGAACGCGCGGGATCCTGGCACCCGATCTCTTTTCCACGCTACAGCGTATTCTTAACCCGGAAGTGGCTGATAAAAACCCTGATGAATCGAGGTATTCCCCCCCCCTATGAGCatgacacaataaca                                                                                             >  4:301477/1‑262 (MQ=255)

                                                                                                                                                                 tGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGaa                                                                                                                                                                                                                        >  5:62287/1‑149 (MQ=255)

                                                                                                                                                                   cAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAAcg                                                                                                                                                                                                                      >  7:1121371/1‑149 (MQ=255)

                                                                                                                                                                    aGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAAcgc                                                                                                                                                                                                                     >  8:835624/1‑149 (MQ=255)

                                                                                                                                                                    aGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAAcgc                                                                                                                                                                                                                     >  8:15170/1‑149 (MQ=255)

                                                                                                                                                                      gCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAAcgcgc                                                                                                                                                                                                                   >  6:1295500/1‑149 (MQ=255)

                                                                                                                                                                      gCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAAcgcgc                                                                                                                                                                                                                   >  5:671114/1‑149 (MQ=255)

                                                                                                                                                                             gAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGc                                                                                                                                                                                                            <  7:637835/149‑1 (MQ=255)

                                                                                                                                                                             gAGCGGCGGCATGGGCGATGTGCCCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGc                                                                                                                                                                                                            <  6:542082/149‑1 (MQ=255)

                                                                                                                                                                              aGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACAt                                                                                                  <  4:201624/254‑1 (MQ=255)

                                                                                                                                                                                cggcggCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATa                                                                                                                                                   >  2:105899/1‑203 (MQ=255)

                                                                                                                                                                                    ggcATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATaaaaa                                                                                                                                               <  4:105899/203‑1 (MQ=255)

                                                                                                                                                                                        tGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGAt                                                                                                                                                                                                 <  6:317756/149‑1 (MQ=255)

                                                                                                                                                                                        tGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGAt                                                                                                                                                                                                 <  5:762868/149‑1 (MQ=255)

                                                                                                                                                                                            cGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATctct                                                                                                                                                                                             <  5:122480/149‑1 (MQ=255)

                                                                                                                                                                                            cGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATctct                                                                                                                                                                                             <  5:72029/149‑1 (MQ=255)

                                                                                                                                                                                                gtgCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTc                                                                                                                                                                                         <  6:167023/149‑1 (MQ=255)

                                                                                                                                                                                                 tgCTCTCTGGTATTCTTGGCGCATCGCTTGGGCAAAAACTGTCGCCGTATGACGAAGCCTGTGCAGGCTGTGTCGCGCATGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGAGGGATGCTGGCAACGGATCTCTTTTcc                                                                                                                                                                                        >  6:821388/1‑149 (MQ=255)

                                                                                                                                                                                                 tgCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTcc                                                                                                                                                                                        <  7:1417827/149‑1 (MQ=255)

                                                                                                                                                                                                     ctctGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTGGCGGTATGATGCAGCCTGTGCAGGCTGAGTCACGCGCGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACTCGCGGGATGCTGGCAACCGATCTCTTGTCcgcgg                                                                                                                                                                                    <  6:1172487/149‑2 (MQ=255)

                                                                                                                                                                                                        tGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTAc                                                                                                                                                                                 >  8:324752/1‑149 (MQ=255)

                                                                                                                                                                                                         ggTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACa                                                                                                                                                                                >  5:272119/1‑149 (MQ=255)

                                                                                                                                                                                                           tattatTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGc                                                                                                                                                                              >  7:1162450/1‑149 (MQ=255)

                                                                                                                                                                                                           tattatTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGc                                                                                                                                                                              <  7:324752/149‑1 (MQ=255)

                                                                                                                                                                                                             ttattGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGt                                                                                                                                                                            >  5:1310007/1‑149 (MQ=255)

                                                                                                                                                                                                              tattGGGGTATTGCTTGGGGAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAatctat                                             <  4:240985/275‑1 (MQ=255)

                                                                                                                                                                                                                     gcATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTaa                                                                                                                                                                    >  7:382838/1‑149 (MQ=255)

                                                                                                                                                                                                                      cATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGTTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCggg                                                                                <  1:11364/232‑1 (MQ=255)

                                                                                                                                                                                                                        ttGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCgg                                                                                                                                                                                                                 <  8:192517/101‑1 (MQ=255)

                                                                                                                                                                                                                        ttGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAAccc                                                                                                                                                                 <  7:7167/149‑1 (MQ=255)

                                                                                                                                                                                                                        ttGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAatctatcgggacggcc                                   >  2:279581/1‑266 (MQ=255)

                                                                                                                                                                                                                              gggCAAAATCTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGt                                                                                                                                                           <  5:1420142/149‑1 (MQ=255)

                                                                                                                                                                                                                                gCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGa                                                                                                                                                         >  8:463927/1‑149 (MQ=255)

                                                                                                                                                                                                                                gCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGa                                                                                                                                                         <  7:835624/149‑1 (MQ=255)

                                                                                                                                                                                                                                  aaaaaCTGTCGCCGTATGATGCAGCCTGTGCTGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGgcg                          <  2:167020/274‑1 (MQ=255)

                                                                                                                                                                                                                                  aaaaaCTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACt                                                                                                                                                       >  6:1355220/1‑149 (MQ=255)

                                                                                                                                                                                                                                   aaaaCTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTg                                                                                                                                                      >  5:470116/1‑149 (MQ=255)

                                                                                                                                                                                                                                    aaaCTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGa                                                                                                                                                     >  6:605078/1‑149 (MQ=255)

                                                                                                                                                                                                                                      aCTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATa                                                                                                                                                   >  8:953425/1‑149 (MQ=255)

                                                                                                                                                                                                                                          tCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATaaaaa                                                                                                                                               <  5:517409/149‑1 (MQ=255)

                                                                                                                                                                                                                                          tCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATaaaaa                                                                                                                                               <  6:554000/149‑1 (MQ=255)

                                                                                                                                                                                                                                          tCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTaaaa                 >  2:2150/1‑275 (MQ=255)

                                                                                                                                                                                                                                            gCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGTTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAAcc                                                                                                                                             <  7:172655/149‑1 (MQ=255)

                                                                                                                                                                                                                                             ccGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCa                                                                                                                                            <  5:742073/149‑1 (MQ=255)

                                                                                                                                                                                                                                                tATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCatga                                                                                                                                         <  8:413957/149‑1 (MQ=255)

                                                                                                                                                                                                                                                 gcgacacAGCCTGCACAGGCTGCGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCat                                                                                                                                           >  6:295680/7‑146 (MQ=255)

                                                                                                                                                                                                                                                 atgatgCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCatgat                                                                                                                                        <  7:993196/149‑1 (MQ=255)

                                                                                                                                                                                                                                                   gacacAGCCTGCACAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCAtgatga                                                                                                                                      >  6:710368/5‑149 (MQ=255)

                                                                                                                                                                                                                                                     cacAGCCTGCACAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAAt                                                                                                                                    >  7:1015938/3‑149 (MQ=255)

                                                                                                                                                                                                                                                      gCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATc                                                                                                                                   <  5:904635/149‑1 (MQ=255)

                                                                                                                                                                                                                                                        aGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGa                                                                                                                                 <  8:1297182/149‑1 (MQ=255)

                                                                                                                                                                                                                                                         gCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGACCGCGCGGGATGCTGGCAACCGATCTCTTTTCCGCGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAATCCATGATGAACCGAg                                                                                                                                <  8:1076247/149‑1 (MQ=255)

                                                                                                                                                                                                                                                         gCCTGGGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGAGGGAGGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCggg  <  2:301477/275‑1 (MQ=255)

                                                                                                                                                                                                                                                          ccTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGgcg                                                           <  4:273852/217‑1 (MQ=255)

                                                                                                                                                                                                                                                             gtgCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAAt                                                                                                                            >  7:208848/1‑149 (MQ=255)

                                                                                                                                                                                                                                                              tgCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAatct                                               <  1:45718/225‑1 (MQ=255)

GCCGCACGGTTGTTAGGCTGTTCCGTCGCTGAAATTGAAAGTGACCGCTTACATTGCGCCAAACGTCTGGTACAACGTTATGGCGGCGTAGCGGTGCTGAAAGGTGCCGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGG  >  NC_000913/4395200‑4395723

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