BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	Exported	522,467	2 bp→CGG	coding (2546‑2547/2592 nt) coding (377‑378/462 nt)	BCHEFKBB_00548 → BCHEFKBB_00549 ←	restriction endonuclease EcoPI subunit R[NGO_0546] hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	522,466	1	.	C	97.8%	357.3 / 2.0	90	coding (2545/2592 nt) coding (379/462 nt)	BCHEFKBB_00548 BCHEFKBB_00549	restriction endonuclease EcoPI subunit R[NGO_0546] hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (2/0); new base C (34/54); total (36/54)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.09e-01
*	Exported	522,467	0	T	G	100.0%	365.0 / NA	90	L849R (CTA→CGG) V126G (GTA→GGC)	BCHEFKBB_00548 BCHEFKBB_00549	restriction endonuclease EcoPI subunit R[NGO_0546] hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); new base G (36/54); total (36/54)
*	Exported	522,468	0	A	G	100.0%	362.6 / NA	90	L849R (CTA→CGG) V126G (GTA→GGC)	BCHEFKBB_00548 BCHEFKBB_00549	restriction endonuclease EcoPI subunit R[NGO_0546] hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (36/54); total (36/54)

Only 100 of 130 total aligned reads displayed.

GGATTACCGGCTTGTCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCAGTATATGCC‑G‑C‑TACTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTATTCGCCCGACTT  >  Exported/522342‑522618

|||

ggATTACCGGCTTGTCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTa                                                                                                                                >  1:102250/1‑150 (MQ=255)

ggATTACCGGCTTGTCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTa                                                                                                                                >  1:878762/1‑150 (MQ=255)

ggATTACCGGCTTGTCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTa                                                                                                                                >  1:550405/1‑150 (MQ=255)

ggATTACCGGCTTGTCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTa                                                                                                                                >  1:193673/1‑150 (MQ=255)

ggATTACCGGCTTGTCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTa                                                                                                                                >  2:802417/1‑150 (MQ=255)

ggATTACCGGCTTGTCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTa                                                                                                                                >  1:138814/1‑150 (MQ=255)

 gATTACCGGCTTGTCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACg                                                                                                                              >  2:159158/1‑151 (MQ=255)

        ggCTTGTCGGCAGCAATATCCATCCGTCCAAATTTAGCAATAAAGACGGCAAACCGCGCGCGGTGAAACAAGCAGATTTGGGCTGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGAGTTTCTACGAttggg                                                                                                                        >  2:1496896/1‑147 (MQ=255)

          cTTGTCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGTAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTGCGATTCGGata                                                                                                                     >  2:1282287/1‑151 (MQ=255)

          cTTGTCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGata                                                                                                                     <  2:326282/151‑1 (MQ=255)

              tCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAg                                                                                                                  <  2:391433/150‑1 (MQ=255)

              tCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGa                                                                                                                 <  2:44905/151‑1 (MQ=255)

               cGGCAGCACACTCCATCCGACCAAATTTACCAATGAAGACGGCAAACCGCGCGCGGTGACAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGaa                                                                                                                <  2:1455165/151‑1 (MQ=255)

               cGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGaa                                                                                                                <  1:302243/151‑1 (MQ=255)

               cGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGaa                                                                                                                <  1:219974/151‑1 (MQ=255)

               cGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGaa                                                                                                                <  1:1369791/151‑1 (MQ=255)

               cGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGaa                                                                                                                <  1:1270371/151‑1 (MQ=255)

                    gcaAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGa                                                                                                           >  1:907426/1‑151 (MQ=255)

                    gaaaaaTCCATCCGACCAAATTTACCAATAAATACGTCAAACCGCGAGCTGTGAAAAAACCAGATTTGGTCAGATTAGAAGATACGGAGCACCGG‑‑CCTGCC‑GAC‑GGCTATCTCTTCGGCGAGAATTTCTACGAGTCGGATATAGAACaata                                                                                                           >  2:1002951/3‑148 (MQ=255)

                      aatctccatccGACCATATTTACCAATAAAGACGGCCAAGCGCGCGCGGTGATAAAAGCCGATTTGGGCAGATTTGAAGAGACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGATCATGat                                                                                                          <  2:930930/146‑2 (MQ=255)

                       aaaTCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGaaaa                                                                                                        >  2:1066017/1‑151 (MQ=255)

                        aaTCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAAt                                                                                                       <  2:1351690/151‑1 (MQ=255)

                        aaTCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAAt                                                                                                       <  2:570320/151‑1 (MQ=255)

                         atccatccGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATg                                                                                                      <  1:821792/151‑1 (MQ=255)

                                  aCCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCcaac                                                                                              <  1:585784/150‑1 (MQ=255)

                                      aaTTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCa                                                                                         >  1:123899/1‑151 (MQ=255)

                                         ttACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAAt                                                                                      <  1:1421460/151‑1 (MQ=255)

                                         ttACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAAt                                                                                      <  1:1421452/151‑1 (MQ=255)

                                           aCCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTg                                                                                    <  1:859535/151‑1 (MQ=255)

                                           aCCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTg                                                                                    <  1:283184/151‑1 (MQ=255)

                                           aCCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTg                                                                                    <  1:283171/151‑1 (MQ=255)

                                               aTAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAgg                                                                                <  2:171157/151‑1 (MQ=255)

                                               aTAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAgg                                                                                <  2:1187064/151‑1 (MQ=255)

                                                tAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGc                                                                               <  1:76692/151‑1 (MQ=255)

                                                 aaaGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCg                                                                              >  1:1072976/1‑151 (MQ=255)

                                                 aaaGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCg                                                                              >  1:1073048/1‑151 (MQ=255)

                                                 aaaGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCg                                                                              >  1:1073049/1‑151 (MQ=255)

                                                      cGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAAt                                                                          >  2:478838/1‑150 (MQ=255)

                                                      cGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATc                                                                         <  1:1066352/151‑1 (MQ=255)

                                                      cGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATc                                                                         <  1:629783/151‑1 (MQ=255)

                                                      cGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATc                                                                         <  1:629786/151‑1 (MQ=255)

                                                       ggCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCg                                                                        >  2:637455/1‑151 (MQ=255)

                                                        gCAAACCGCGCGCTGTTAAAAAAGCAGATGTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCg                                                                        <  1:1316716/150‑1 (MQ=255)

                                                        gCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCg                                                                        <  1:238897/150‑1 (MQ=255)

                                                        gCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCg                                                                        <  2:482096/150‑1 (MQ=255)

                                                        gCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCg                                                                        <  1:791379/150‑1 (MQ=255)

                                                        gCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCg                                                                        <  1:1067939/150‑1 (MQ=255)

                                                          aaaCCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTAt                                                                     <  2:1431340/151‑1 (MQ=255)

                                                            aCCGCGCGCGGTGAACACAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTAttt                                                                   <  2:1532525/151‑1 (MQ=255)

                                                            aCCGCGCGCGGTGAAAACAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTAttt                                                                   <  2:1532548/151‑1 (MQ=255)

                                                                cgcgcgGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTAcc                                                                <  1:106882/150‑1 (MQ=255)

                                                                       tGAAAACAGCAGATTTGGGCAGATTTGCAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATGTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATAc                                                         <  2:814449/150‑1 (MQ=255)

                                                                       tGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATAc                                                         <  1:763073/150‑1 (MQ=255)

                                                                       tGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATAc                                                         <  1:786303/150‑1 (MQ=255)

                                                                       tGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATAc                                                         <  1:981115/150‑1 (MQ=255)

                                                                        gAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCg                                                       >  2:1469343/1‑151 (MQ=255)

                                                                        gAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCg                                                       >  1:1617400/1‑151 (MQ=255)

                                                                                  gATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCAC‑C‑GGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCa                                             >  1:148108/1‑151 (MQ=255)

                                                                                  gATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCa                                             >  1:267894/1‑151 (MQ=255)

                                                                                  gATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCa                                             >  1:1420862/1‑151 (MQ=255)

                                                                                  gATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGAATTTACCAAAATACCGAGAAACTCCa                                             >  1:581495/1‑151 (MQ=255)

                                                                                  gATTTGGGCAGATTTGAAGATACGGAGCACC‑G‑CCCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAAAACCGAGAAACTCCa                                             >  1:581499/1‑151 (MQ=255)

                                                                                   aTTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTc                                               <  1:1369818/148‑1 (MQ=255)

                                                                                   aTTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTc                                               <  2:1552080/148‑1 (MQ=255)

                                                                                   aTTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTc                                               <  1:242409/148‑1 (MQ=255)

                                                                                   aTTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTc                                               <  1:242511/148‑1 (MQ=255)

                                                                                   aTTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTc                                               <  1:834792/148‑1 (MQ=255)

                                                                                   aTTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTc                                               <  1:834786/148‑1 (MQ=255)

                                                                                       gggCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGTCGTAATCTTATTTACAAAAATACCGAGAAACTCCAACaaa                                        >  1:1630058/1‑151 (MQ=255)

                                                                                       gggCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCC‑G‑CCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCaaa                                        >  1:1630057/1‑151 (MQ=255)

                                                                                                                   ccggcctgccgccggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACg                >  2:1625727/16‑151 (MQ=255)

                                                                                                                     ggcctgccgccggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGt               >  2:664676/14‑150 (MQ=255)

                                                                                                                     ggcctgccgccggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTa              >  2:70939/14‑151 (MQ=255)

                                                                                                                     ggcctgccgccggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTa              >  2:964580/14‑151 (MQ=255)

                                                                                                                     ggcctgccgccggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTa              >  2:964596/14‑151 (MQ=255)

                                                                                                                     ggcctgccgccggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTa              >  2:1211545/14‑151 (MQ=255)

                                                                                                                       cctgccgccggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTTATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTAtt            <  1:16103/140‑1 (MQ=255)

                                                                                                                       cctgccgccggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTAtt            <  2:669603/140‑1 (MQ=255)

                                                                                                                       cctgccgccggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTAtt            <  1:16102/140‑1 (MQ=255)

                                                                                                                        ctgccgccggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTATTc           <  2:1339824/141‑1 (MQ=255)

                                                                                                                          gccgccggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTATTCg          >  1:364487/9‑150 (MQ=255)

                                                                                                                           ccgccggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTATTCGcc        <  2:1634255/144‑1 (MQ=255)

                                                                                                                           ccgccggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTATTCGcc        <  1:578845/144‑1 (MQ=255)

                                                                                                                           ccgccggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTATTCGcc        <  1:578841/144‑1 (MQ=255)

                                                                                                                           ccgccggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTATTCGcc        <  1:55919/144‑1 (MQ=255)

                                                                                                                           ccgccggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTATTCGcc        <  2:550239/144‑1 (MQ=255)

                                                                                                                           ccgccggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTATTCGc         >  1:1346154/8‑150 (MQ=255)

                                                                                                                           ccgccggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTATTCGc         >  1:154858/8‑150 (MQ=255)

                                                                                                                           ccgccggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATAGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTATTCGcc        <  2:1376191/144‑1 (MQ=255)

                                                                                                                           ccgccggctaTCTCTTCGGAGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTATTCGcc        <  1:967745/144‑1 (MQ=255)

                                                                                                                                ggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTAtt            <  2:156156/140‑1 (MQ=255)

                                                                                                                                ggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTAtt            >  1:156209/3‑142 (MQ=255)

                                                                                                                                ggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTAtt            <  2:156140/140‑1 (MQ=255)

                                                                                                                                ggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTAtt            >  1:156193/3‑142 (MQ=255)

                                                                                                                                ggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTATTCGCCCGa     >  2:1643046/3‑149 (MQ=255)

                                                                                                                                ggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTATTCGCCCGa     >  1:336693/3‑149 (MQ=255)

                                                                                                                                ggctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTATTCGCCCGa     >  1:621869/3‑149 (MQ=255)

                                                                                                                                 gctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTATTCGCCCGACtt  >  1:1056986/2‑151 (MQ=255)

                                                                                                                                 gctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTATTCGCCCGACtt  >  1:1056952/2‑151 (MQ=255)

                                                                                                                                 gctaTCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTATTCGCCCGACtt  >  1:1056795/2‑151 (MQ=255)

|||

GGATTACCGGCTTGTCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCAGTATATGCC‑G‑C‑TACTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAATCCCTGTTGCCGGCGGCGGCACGTATTCGCCCGACTT  >  Exported/522342‑522618

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