BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,822,679	A→C	5.1%	intergenic (‑40/+29)	recX ← / ← recA	regulatory protein for RecA/DNA recombination and repair protein; ssDNA‑dependent ATPase; synaptase; ssDNA and dsDNA binding protein; ATP‑dependent homologous DNA strand exchanger; recombinase A; LexA autocleavage cofactor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,822,679	0	A	C	5.1%	755.5 / 25.5	215	intergenic (‑40/+29)	recX/recA	regulatory protein for RecA/DNA recombination and repair protein; ssDNA‑dependent ATPase; synaptase; ssDNA and dsDNA binding protein; ATP‑dependent homologous DNA strand exchanger; recombinase A; LexA autocleavage cofactor
Reads supporting (aligned to +/- strand): ref base A (97/107); new base C (7/4); total (104/111)
Fisher's exact test for biased strand distribution p-value = 3.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.47e-01

Only 100 of 216 total aligned reads displayed.

AACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCACCTTTGTAGCTGTACCACGCGCCTGCTTTCTCGATCAGCTTCTCTTTTACGCCCAGGTCAACCAGTT  >  NC_000913/2822415‑2822949

aacaaATCGGCCGTCATCCAGATAGCCATGTACATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCtt                                                                                                                                                                                                                                                                      <  2:136925/275‑1 (MQ=255)

                                          ggCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTTCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAACAACGCAAAAGGGCCGCAGATGCGCCCCTTGTGTATCaaa                                                                                                                                                                                                                                                            >  2:146663/1‑243 (MQ=255)

                                              aTAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGc                                                                                                                                                                                                                                                                             <  4:69492/222‑1 (MQ=255)

                                                          tAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGtt                                                                                                                                                                                                                                 <  2:120487/254‑1 (MQ=255)

                                                                 ctGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATTCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCCCCATCATCTCCaa                                                                                                                                                                                                     >  4:317353/1‑274 (MQ=255)

                                                                                   ctttggggCCATTTTTGACCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGt                                                                                                                                                                                   <  2:324199/270‑1 (MQ=255)

                                                                                             aTTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGGGATCGCGCACCGCCAGAAAGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCt                                                                                                                                                                         <  4:249160/275‑1 (MQ=255)

                                                                                                    cccATATTCGGGGCCGCGAGTTTACGTCGCAGTTCTTGCTACCTGTGAGCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCagcagc                                                                                                                                                                  <  4:243818/275‑1 (MQ=255)

                                                                                                    cccATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCagcagc                                                                                                                                                                  <  2:317353/275‑1 (MQ=255)

                                                                                                    cccATAATCGGTGCCGCGAGTTTACGTCGCAGGTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGAGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCagcagc                                                                                                                                                                  <  2:318942/275‑1 (MQ=255)

                                                                                                        tAATCGGTGCCGCGAGTTTACGTCGAAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCAGGGCGACGGGCTGTTGATTCTGCCATGGCATATCCTTACAACTTAAAAACGAACAAGGGCAGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTAGCATTGACAGAGAAATCCGGCGGTGAGTTCGGGTTGCTCAGAAGCaa                                                                                                                                                                >  4:46006/1‑273 (MQ=255)

                                                                                                              gTGCCGCGAGATTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGAATAAAAATCTTCGTTAGTTTCTGCTACGCCTTAGCTTTCATCTACAGAGACATCCGGCGTTGCGTTCGGGTTGCTCATCAGCACCTCACg                                                                                                                                                          >  2:96159/1‑273 (MQ=255)

                                                                                                                      aGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATg                                                                                                                                                                                                                                                                              <  5:429534/149‑1 (MQ=255)

                                                                                                                             gTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCt                                                                                                                                                                                                                                                                       >  5:535700/1‑149 (MQ=255)

                                                                                                                             gTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCt                                                                                                                                                                                                                                                                       >  8:1221982/1‑149 (MQ=255)

                                                                                                                               cGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTtg                                                                                                                                                                                                                                                                     >  5:982727/1‑149 (MQ=255)

                                                                                                                               cGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGACGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTtg                                                                                                                                                                                                                                                                     >  7:1414930/1‑149 (MQ=255)

                                                                                                                                gCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTtgt                                                                                                                                                                                                                                                                    <  8:509936/149‑1 (MQ=255)

                                                                                                                                gCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTtgt                                                                                                                                                                                                                                                                    <  8:486162/149‑1 (MQ=255)

                                                                                                                                 cAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGGATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTtgtg                                                                                                                                                                                                                                                                   <  7:1104015/149‑1 (MQ=255)

                                                                                                                                 cAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACa                                                                                                                                                                                                                                                          >  4:31297/1‑158 (MQ=255)

                                                                                                                                     tCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATc                                                                                                                                                                                                                                                               >  5:553444/1‑149 (MQ=255)

                                                                                                                                       ttGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTa                                                                                                                                                                                                                                                                 >  6:538211/1‑145 (MQ=255)

                                                                                                                                        tGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCaaa                                                                                                                                                                                                                                                            <  8:317143/149‑1 (MQ=255)

                                                                                                                                        tGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCaaa                                                                                                                                                                                                                                                            <  5:1470949/149‑1 (MQ=255)

                                                                                                                                           tccctaggaccgcgAACCGCTATAATGCATACCGAACGAACCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTcc                                                                                                                           <  1:37206/265‑1 (MQ=255)

                                                                                                                                           tCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACaa                                                                                                                                                                                                                                                         >  8:1132571/1‑149 (MQ=255)

                                                                                                                                           tCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACaa                                                                                                                                                                                                                                                         >  6:905941/1‑149 (MQ=255)

                                                                                                                                           tCACTGTGAGCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCGGCAGATGCGACCCTTGTGTATCAAACaa                                                                                                                                                                                                                                                         >  8:1254981/1‑149 (MQ=255)

                                                                                                                                             aCTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGTCGCATCTGCGGCCCttt                                                                                                                                                                                                                                                                     <  8:676937/135‑2 (MQ=255)

                                                                                                                                             aCTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGa                                                                                                                                                                                                                                                       >  6:528228/1‑149 (MQ=255)

                                                                                                                                             aCTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGa                                                                                                                                                                                                                                                       >  7:940918/1‑149 (MQ=255)

                                                                                                                                                       cgcACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTaaaaa                                                                                                                                                                                                                                             <  5:247133/149‑1 (MQ=255)

                                                                                                                                                         cACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGAc                                                                                                                                                                                                                                                                          >  2:304798/1‑118 (MQ=255)

                                                                                                                                                            cGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTc                                                                                                                                                                                                                                        <  8:878231/149‑1 (MQ=255)

                                                                                                                                                            cGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTc                                                                                                                                                                                                                                        <  5:1121650/149‑1 (MQ=255)

                                                                                                                                                             gCCGGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCg                                                                                                                                                                                                                                       >  7:1044529/1‑149 (MQ=255)

                                                                                                                                                             gCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCg                                                                                                                                                                                                                                       >  6:288696/1‑149 (MQ=255)

                                                                                                                                                              ccAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGt                                                                                                                                                                                                                                      <  6:1090167/149‑1 (MQ=255)

                                                                                                                                                                aGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTa                                                                                                                                                                                                                                    <  5:1318092/149‑1 (MQ=255)

                                                                                                                                                                        tACCGCACGACCCAACAGGCGGGGATATTCCGGGCGACGGGATGTTGATTCTGTCATTGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCAt                                                                                               <  2:226090/274‑1 (MQ=255)

                                                                                                                                                                           cGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTAc                                                                                                                                                                                                                         >  5:286616/1‑149 (MQ=255)

                                                                                                                                                                           cGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTAc                                                                                                                                                                                                                         >  5:732891/1‑149 (MQ=255)

                                                                                                                                                                              aCGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGcc                                                                                                                                                                                                                      >  5:179843/1‑149 (MQ=255)

                                                                                                                                                                                gATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCtt                                                                                                                                                                                                                    >  7:116615/1‑149 (MQ=255)

                                                                                                                                                                                gATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCtt                                                                                                                                                                                                                    >  7:1151539/1‑149 (MQ=255)

                                                                                                                                                                                 agccaacggccttgCATATTCCGGTCGAGNNTATGTTTGTTCTTGCATTGCATATCCGTTCNACTTAAAGAAGCATACGTGCAGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGGTTNGGCTACGCCTTCGCTAACAGCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAGCTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTAcc                                                                                     <  3:47097/262‑1 (MQ=255)

                                                                                                                                                                                 aTCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTtgt                                                                                                                                                                                                                                                                    <  5:627173/100‑1 (MQ=255)

                                                                                                                                                                                 aTCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTc                                                                                                                                                                                                                   <  6:928115/149‑1 (MQ=255)

                                                                                                                                                                                     aaCAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTa                                                                                                                                                                                                               >  5:1446399/1‑149 (MQ=255)

                                                                                                                                                                                            gAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTTAcctaact                                                                          >  4:270384/1‑270 (MQ=255)

                                                                                                                                                                                              gCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACa                                                                                                                                                                                                      <  6:982727/149‑1 (MQ=255)

                                                                                                                                                                                                atatGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACaga                                                                                                                                                                                                    <  6:647018/149‑1 (MQ=255)

                                                                                                                                                                                                atatGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACaga                                                                                                                                                                                                    >  7:1352559/1‑149 (MQ=255)

                                                                                                                                                                                                       gggCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATcc                                                                                                                                                                                             <  6:736881/149‑1 (MQ=255)

                                                                                                                                                                                                        ggCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCACCTTTg                                                               >  2:164705/1‑274 (MQ=255)

                                                                                                                                                                                                           gACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCg                                                                                                                                                                                         >  7:1326844/1‑149 (MQ=255)

                                                                                                                                                                                                             cGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGtt                                                                                                                                                                                       <  5:1014606/149‑1 (MQ=255)

                                                                                                                                                                                                             cGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGtt                                                                                                                                                                                       <  7:1446232/149‑1 (MQ=255)

                                                                                                                                                                                                              gggATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCt                                                                                                                                                                                                                >  8:133605/1‑123 (MQ=255)

                                                                                                                                                                                                              gggATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCATTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTg                                                                                                                                                                                      <  7:460040/149‑1 (MQ=255)

                                                                                                                                                                                                                  tGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACt                                                                                                                                                              <  4:14288/169‑1 (MQ=255)

                                                                                                                                                                                                                    ttGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGcc                                                                                                                                                                                                                      <  7:10169/111‑1 (MQ=255)

                                                                                                                                                                                                                    ttGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCg                                                                                                                                                                                >  5:777157/1‑149 (MQ=255)

                                                                                                                                                                                                                    ttGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCg                                                                                                                                                                                <  5:652092/149‑1 (MQ=255)

                                                                                                                                                                                                                     tGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCgg                                                                                                                                                                               <  8:1018526/149‑1 (MQ=255)

                                                                                                                                                                                                                     tGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTCCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCgg                                                                                                                                                                               <  7:778931/149‑1 (MQ=255)

                                                                                                                                                                                                                      gATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCggg                                                                                                                                                                              <  7:653713/149‑1 (MQ=255)

                                                                                                                                                                                                                        ttCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGtt                                                                                                                                                                            >  6:1149635/1‑149 (MQ=255)

                                                                                                                                                                                                                        ttCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGtt                                                                                                                                                                            <  6:1044653/149‑1 (MQ=255)

                                                                                                                                                                                                                          cTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGc                                                                                                                                                                          <  6:553444/149‑1 (MQ=255)

                                                                                                                                                                                                                          cTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGc                                                                                                                                                                          >  5:294346/1‑149 (MQ=255)

                                                                                                                                                                                                                           tGTAATGGCATATCCCTACAACCTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCCAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCACCTTTGTAGCTGTACCACGCGCCTg                                            <  4:321871/274‑1 (MQ=255)

                                                                                                                                                                                                                             tCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTca                                                                                                                                                                       >  6:911369/1‑149 (MQ=255)

                                                                                                                                                                                                                             tCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTca                                                                                                                                                                       >  8:1355922/1‑149 (MQ=255)

                                                                                                                                                                                                                             tCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTca                                                                                                                                                                       >  8:1176474/1‑149 (MQ=255)

                                                                                                                                                                                                                               aTGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTcagc                                                                                                                                                                     >  8:749499/1‑149 (MQ=255)

                                                                                                                                                                                                                               aTGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTcagc                                                                                                                                                                     <  5:326173/149‑1 (MQ=255)

                                                                                                                                                                                                                                tGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCgg                                                                                                                                                                               >  5:910109/1‑138 (MQ=255)

                                                                                                                                                                                                                                  gCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAg                                                                                                                                                                                                                                                        <  4:382207/63‑1 (MQ=255)

                                                                                                                                                                                                                                  gCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTc                                                                                                                                                                                                                   >  8:562830/1‑100 (MQ=255)

                                                                                                                                                                                                                                  gCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTc                                                                                                                                                                                                                   <  7:562830/100‑1 (MQ=255)

                                                                                                                                                                                                                                     tatCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAAc                                                                                                                                                               <  8:1326844/149‑1 (MQ=255)

                                                                                                                                                                                                                                      atCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTc                                                                                                                                                                                 >  5:803501/1‑130 (MQ=255)

                                                                                                                                                                                                                                       tCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTc                                                                                                                                                             <  7:1000837/149‑1 (MQ=255)

                                                                                                                                                                                                                                           tACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGt                                                                                                                                                         >  6:1360864/1‑149 (MQ=255)

                                                                                                                                                                                                                                            aCAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTg                                                                                  >  2:253315/1‑219 (MQ=255)

                                                                                                                                                                                                                                                 ttAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTttc                                                                                                                                                   >  8:162587/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   aaaaaaGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTcagca                                                                                                                                                                    <  6:294807/130‑1 (MQ=255)

                                                                                                                                                                                                                                                   aaaaaaGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTttctt                                                                                                                                                 <  6:1464418/149‑1 (MQ=255)

                                                                                                                                                                                                                                                     aaaaGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTtcttct                                                                                                                                               >  7:1192898/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     aaaaGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGAAGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCACCTTTGTAGCTGTACCACGCGCCTGCTTTCTCGATCAGCTTCTCTTTTACGc                 <  1:48675/275‑1 (MQ=255)

                                                                                                                                                                                                                                                           aaaaGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTc                                                                                                                                                             >  5:689415/1‑129 (MQ=255)

                                                                                                                                                                                                                                                           aaaaGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTc                                                                                                                                                             <  6:689415/129‑1 (MQ=255)

                                                                                                                                                                                                                                                           aaaaGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGAtct                                                                                                                                         >  6:1221123/1‑149 (MQ=255)

                                                                                                                                                                                                                                                           aaaaGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGAtct                                                                                                                                         >  7:252557/1‑149 (MQ=255)

                                                                                                                                                                                                                                                           aaaaGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGAtct                                                                                                                                         >  5:191636/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                ggCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTc                                                                                                                                    <  6:186701/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                   cGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTcgcg                                                                                                                                 <  6:1118130/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                    gCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCACCTTTGTAGCTGTACCACGCGCCTGCTTTCTCGATCAGCTTCTCTTTTACGCCCAGGTCAACCAGtt  <  2:53577/275‑1 (MQ=255)

AACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCACCTTTGTAGCTGTACCACGCGCCTGCTTTCTCGATCAGCTTCTCTTTTACGCCCAGGTCAACCAGTT  >  NC_000913/2822415‑2822949

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