BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,822,684	A→G	5.1%	intergenic (‑45/+24)	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,684	0	A	G	5.1%	722.0 / 18.1	217	intergenic (‑45/+24)	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 (95/110); new base G (7/4); total (103/114)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.46e-01

Only 100 of 218 total aligned reads displayed.

AACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCACCTTTGTAGCTGTACCACGCGCCTGCTTTCTCGATCAGCTTCTCTTTTACGCCCAGGTCAACCAGTTCGCCGTAG  >  NC_000913/2822415‑2822957

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

                               ccATGGTACCAGGCAATAACGCGCTCGTAGTCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCAACGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCg                                                                                                                                                                                                                                               <  2:318530/274‑1 (MQ=255)

                                        cAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCCGATGCGACCCTTTTGTATCAAACAAGACTATTAAAAACCTTCGTTAGTTTc                                                                                                                                                                                                                                       >  4:318942/1‑274 (MQ=255)

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

                                                      cTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGtt                                                                                                                                                                                                                                             >  4:120487/1‑254 (MQ=255)

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

                                                                       tAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTCCAGAGAAAt                                                                                                                                                                                                       >  2:7402/1‑275 (MQ=255)

                                                                         gCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATcc                                                                                                                                                                                                     >  2:321871/1‑275 (MQ=255)

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

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

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

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

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

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

                                                                                                                               cGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTtg                                                                                                                                                                                                                                                                             <  8:1254232/149‑1 (MQ=255)

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

                                                                                                                                gCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTtgt                                                                                                                                                                                                                                                                            <  5:802817/149‑1 (MQ=255)

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

                                                                                                                                gCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTtgt                                                                                                                                                                                                                                                                            <  6:1175286/149‑1 (MQ=255)

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

                                                                                                                                 cAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATCTCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAAAAATCCGGCGTTGAGTTCGGGTTCCTCAGCTGCACCTTACCTACTttcttccggaccccttt                                                                                                                                             >  4:313528/1‑264 (MQ=255)

                                                                                                                                     tCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATc                                                                                                                                                                                                                                                                       <  6:919636/149‑1 (MQ=255)

                                                                                                                                     tCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGAc                                                                                                                                                                                                                                                              <  2:31297/158‑1 (MQ=255)

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

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

                                                                                                                                           tCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACaa                                                                                                                                                                                                                                                                 >  7:904842/1‑149 (MQ=255)

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

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

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

                                                                                                                                             aCTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGa                                                                                                                                                                                                                                                               <  6:219988/149‑1 (MQ=255)

                                                                                                                                             aCTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGa                                                                                                                                                                                                                                                               >  5:150179/1‑149 (MQ=255)

                                                                                                                                               tgtgATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACg                                                                                                                                                                                                                                                             >  5:928115/1‑149 (MQ=255)

                                                                                                                                                 tgATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGAt                                                                                                                                                                                                                                                           <  6:1110104/149‑1 (MQ=255)

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

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

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

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

                                                                                                                                                             gCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCtt                                                                                                                                                                                                                                                                              <  4:304798/118‑1 (MQ=255)

                                                                                                                                                                aGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTa                                                                                                                                                                                                                                            <  7:71748/149‑1 (MQ=255)

                                                                                                                                                                        tACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGc                                                                                                                                                                                                                                    <  6:507471/149‑1 (MQ=255)

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

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

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

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

                                                                                                                                                                                gATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCtt                                                                                                                                                                                                                            <  6:1478577/149‑1 (MQ=255)

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

                                                                                                                                                                                 aTCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTtgt                                                                                                                                                                                                                                                                            >  6:627173/1‑100 (MQ=255)

                                                                                                                                                                                     aaCAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTa                                                                                                                                                                                                                       >  6:1021925/1‑149 (MQ=255)

                                                                                                                                                                                     aaCAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTa                                                                                                                                                                                                                       >  7:1232377/1‑149 (MQ=255)

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

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

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

                                                                                                                                                                                                 tatGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGTCGCATCTGCGGCCCtt                                                                                                                                                                                                                                                                              >  8:471119/1‑82 (MQ=38)

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

                                                                                                                                                                                                        ggCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCg                                                                                                                                                                                                    >  6:648605/1‑149 (MQ=255)

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

                                                                                                                                                                                                         gCGACGGGATGTTTATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCACCtt                                                                         <  2:245546/271‑1 (MQ=255)

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

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

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

                                                                                                                                                                                                              gggATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCt                                                                                                                                                                                                                        <  7:133605/123‑1 (MQ=255)

                                                                                                                                                                                                              gggATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTg                                                                                                                                                                                              <  6:535700/149‑1 (MQ=255)

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

                                                                                                                                                                                                               ggATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGa                                                                                                                                                                                             <  6:1339104/149‑1 (MQ=255)

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

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

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

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

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

                                                                                                                                                                                                                         tCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTg                                                                                                                                                                                   >  6:54647/1‑149 (MQ=255)

                                                                                                                                                                                                                            gTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTc                                                                                                                                                                                >  5:350512/1‑149 (MQ=255)

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

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

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

                                                                                                                                                                                                                                tGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCgg                                                                                                                                                                                       <  6:910109/138‑1 (MQ=255)

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

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

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

                                                                                                                                                                                                                                  gCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCagcagc                                                                                                                                                                          <  5:54647/149‑1 (MQ=255)

                                                                                                                                                                                                                                   cATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCACCTTTGTAGCTGTACCACGCGCCTGCTTTCTCGa                                           <  2:270384/275‑1 (MQ=255)

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

                                                                                                                                                                                                                                       tCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGcc                                                                                                                                                                                                                              >  3:39428/1‑92 (MQ=255)

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

                                                                                                                                                                                                                                          ttACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTCTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACg                                                                                                                                                                  >  6:1054216/1‑149 (MQ=255)

                                                                                                                                                                                                                                            aCAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTCTCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATTTACAGATAAATCCGGCGTTGAGTTCGGGTTGCTCCTCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTAAGCCAGGCAGTAGCATTAGCTTTACCCTGACCGATCTTCCCACTTTTGTAGCTGTACCACGCCCCTGCTTTCTCGATCAtcttct                                  >  4:90140/1‑275 (MQ=255)

                                                                                                                                                                                                                                                cTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACttt                                                                                                                                                            >  5:1464418/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 ttAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAACCAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAAc                                                                                                                                                                       <  5:290022/137‑1 (MQ=255)

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

                                                                                                                                                                                                                                                 ttAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTttc                                                                                                                                                           >  5:1372684/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   aaaaaaGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTcagca                                                                                                                                                                            >  5:294807/1‑130 (MQ=255)

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

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

                                                                                                                                                                                                                                                      aaaGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTtctc                                                                                                                                                      <  6:939008/149‑1 (MQ=255)

                                                                                                                                                                                                                                                           acaaGGGTCGCATCTGCGGCCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGAtct                                                                                                                                                 >  5:189173/3‑149 (MQ=255)

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

                                                                                                                                                                                                                                                           aaaaGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGAtct                                                                                                                                                 <  6:833098/149‑1 (MQ=255)

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

                                                                                                                                                                                                                                                                      aGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGtt                                                                                                                                      <  5:237912/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                         tGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTcc                                                                                                                                   <  5:905941/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                           cGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCACCTTTGTAGCTGTACCACGCGCCTGCTTTCTCGATCAGCTTCTCTTTTACGCCCAGGTCAACCAGTTCGCCGTa   <  2:46006/275‑1 (MQ=255)

AACAAATCGGCTGTCATCGAGATAGCCATGTTCATGGCACCAGGCAATAACGCGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCACCTTTGTAGCTGTACCACGCGCCTGCTTTCTCGATCAGCTTCTCTTTTACGCCCAGGTCAACCAGTTCGCCGTAG  >  NC_000913/2822415‑2822957

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: `‑`