New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3579387	106 (1.030)	3 (0.030)	3/286	12.2	3.1%	coding (259/696 nt)	yhhW	quercetinase activity in vitro
	?	NC_000913	= 3579470	86 (0.900)					coding (176/696 nt)	yhhW	quercetinase activity in vitro
Rejected: Position hash score below cutoff.

GCATCGAAGCGACGCTGTTCATAACGCGGCGTAATACCGTTTTCTTCAGGCATGATCCAGATCTGATACAGATGCAGACGCTCGGTGCTGCTTGGGTTGTACTCTGAGTGACGAATACCCGTACCAGCACTCATAATCTGGAACTCACCCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3579236‑3579387 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑actcacccgcATAAAGATATGGAAATTTTGACCTACGTGCTGGAAGGTACTGTTGAGCATCAGGACAGCATGGGCAATAAAGAGCAGGTTCCGGCGGGTGAGTTCCAGATTATGAGTGCTGGTACGGGTATTCGTCACTCAGAGTACAACCCAAGCAGCACCGAGCGTCTGCATCTGTATCAGATCTGGATCATGCCTGAAGAAAAC  <  NC_000913/3579470‑3579274                                                                                                                                                                                                                                                                                                                                                               GCATCGAAGCGACGCTGTTCATAACGCGGCGTAATACCGTTTTCTTCAGGCATGATCCAGATCTGATACAGATGCAGACGCTCGGTGCTGCTTGGGTTGTACTCTGAGTGACGAATACCCGTACCAGCACTCATAATCTGGAACTCACCC                                                                                                                                                                                                         >  3:78533/1‑150                 GTTCATAACGCGGCGTAATACCGTTTTCTTCAGGCATGATCCAGATCTGATACAGATGTAGACGCTCGGTGCTTCTTGGGTTGTACTCTGAGTGACGAATACCCGTACCAGCACTCATCATCTGGAACTCAC                                                                                                                                                                                                           >  2:77369/1‑132                  TTCATAATGCGGCGTAATACCGTTTTCTTCAGGCATGATCCAGATCTGATACAGATGCAGACGCTCGGTGCTGCTTGGGTTGTACTCTGAGTGACGAATACCCGTACCAGCACTCATAATCTGGAACTCACC                                                                                                                                                                                                          <  2:438153/132‑1                  TTCATAACGCGGCGTAATACCGTTTTCTTCAGGCATGATCCAGATCTGATACAGATGCAGACGCTCGGTGCTGCTTGGGTTGTACTCTGAGTGACGAATACCCGTACCAGCACTCATAATCTGGAACTCACC                                                                                                                                                                                                          <  2:436459/132‑1                                                                       GATGCAGACGCTCGGTGCTGCTTGGGTTGTACTCTGAGTGACGAATACCCGTACCAGCACTCATAATCTGGAACTCACCCGCATAAAGATATGGAAATTTTGACCTACGTGCTGGAAGGTACTGTTGAGCATCAGGACAGCATGGGCAATAAAGAGCAGGTTCCGGCGGGTGAGTTCCAGATTATGAGTGCTGGTACGGG                                                                                 >  4:773337/1‑200                                                                                              GGGTTGTACTCTGAGTGACGAATACCCGTACCAGCACTCATAATCTGGAACTCACCCGCATAAAGATATGGAAATTTTGACCTACGTGCTGGAAGGTACTGTTGAGCATCAGGACAGCATGGGCAATAAAGAGCAGGTTCCGGCGGGTGAGTTCCAGATTATGAGTGCTGGTACGGGTATTCGTAACTCAGAGTACAACC                                                          >  4:480892/1‑200                                                                                                          GAGTGACGAATACCCGTACCAGCACTCATAATCTGGAACTCACCCGCATAAAGATATGGAAATTTTGACCTACGTGCTGGAAGGTACTGTTGAGCATCAGGACAGCATGGGCAATAAAGAGCAGGTTCCGGCGGGTGAGTTTCAGATTATGAGTGCTGGTACGGGTATTCGTCACTCAGAGTACAACCCCAGCAGCACCG                                              >  4:122021/1‑200                                                                                                                                                CTCACCCGCATAAAGATATGGAAATTTTGACCTACGTGCTGGAAGGTACTGTTGAGCATCAGGACAGCATGGGCAATAAAGAGCAGGTTCCGGCGGGTGAGTTCCAGATTATGAGTGCTGGTACGGGTATTCGTCACTCAG                                                                   >  1:547350/1‑141                                                                                                                                                 TCACCCGCATAAAGATATGGAAATTTTGACCTACGTGCTGGAAGGTACTGTTGAGCATCAGGACAGCATGGGCAATAAAGAGCAGGTTCCGGCGGGTGAGTTCCAGATTATGAGTGCTGGTACGGGTATTCGTCACTCAGAGTACAACCCAAGCAGCACCGAGCGTCTGCATCTGTATCAGATCTGGATCATGCCTGAAG       <  4:653590/200‑1                                                                                                                                                 TCACCCGCATAAAGATATGGAAATTTTGACCTACGTGCTGGAAGGTACTGTTGAGCATCAGGACAGCATGGGCAATAAAGAGCAGGTTCCGGCGGGTGAGTTCCAGATTATGAGTGCTGGTACGGGTATTCGTCACTCAGAGTACAACCC                                                         <  3:575947/150‑1                                                                                                                                                 TCAACCGCATAAAGATATGGAAATTTTGACCTACGTGCTGGAAGGTACTGTTGAGCATCAGGACAGCATGGGCAATAAAGAGCAGGTTCCGGCGGGTGAGTTCCAGATTATGAGTGCTGGTACGGGTATTCGTCACTCAGAGTACAACCCAAGCAGCACCGAGCGTCTGCATCTGTATCAGATCTGGATCATGCCTGAAG       <  4:170200/200‑1                                                                                                                                                  CACCCGCATAAAGATATGGAAATTTTGACCTACGTGCTGGAAGGTACTGTTGAGCATCAGGACAGCATGGGCAATAAAGAGCAGGTTCCGGCGGGTGAGTTCCAGATTATGAGTGCTGGTACGGGTATTCGT                                                                          >  2:330195/1‑132                                                                                                                                                    CCCGCATAAAGATATGGAAATTTTGACCTACGTGCTGGAAGGTACTGTTGAGCATCAGGACAGCATGGGCAATAAAGAGCAGGTTCCGGCGGGTGAGTTCCAGATTATGAGTGCTGGTACGGGTATTCGTCACTCAGAGTA                                                               >  1:559447/1‑141                                                                                                                                                      CGCATAAAGATATGGAAATTTTGACCTACGTGCTGGAAGGTACTGTTGAGCATCAGGACAGCATGGGCAATAAAGAGCAGGTTCCGGCGGGTGAGTTCCAGATTATGAGTGCTGGTACGGGTATTCGTCACTCAGAGTACAACCCAAGCAGCACCGAGCGTCTGCATCTGTATCAGATCTGGATCATGCCTGAAGAAAAC  <  4:683661/200‑1                                                                                                                                                      CGCATAAAGATATGGAAATTTTGACCTACGTGCTGGAAGGTACTGTTGAGCATCAGGACAGCATGGGCAATAAAGAGCAGGTTCCGGCGGGTGAGTTCCAGATTATGAGTGCTGGTACGGGTATTCGTCACTCAGAGTACA                                                             >  1:33352/1‑141                                                                                                                                                       GCATAAAGATATGGAAATTTTGACCTACGTGCTGGAAGGTATTGTTGAGCATCAGGACAGCATGGGCAATAAAGAGCAGGTTCCGGCGGGTGAGTTCCAGATTATGAGTGCTGGTAAGGGTATTCGTCACTC                                                                     >  2:556460/1‑132                                                                                                                                                       GCATAAAGATATGGAAATTTTGACCTACGTGCTGGAAGGTACTGTTGAGCATCAGGACAGCATGGGCAATAAAGAGCAGGTTCCGGCGGGTGAGTTCCAGATTATGAGTGCTGGTACGGGTATTCGTCACTCAGAGTACAA                                                            >  1:415291/1‑141                                                                                                                                                       GCATAAAGATATGGAAATTTTGACCTACGTGCTGGAAGGTACTGTTGAGCATCAGGACAGCATGGGCAATAAAGAGCAGGTTCCGGCGGGTGAGTTCCAGATTATGAGTGCTGGTACGGGTATTCGTCACTC                                                                     >  2:629669/1‑132                                                                                                                                                                                                                                                                                                                                                               GCATCGAAGCGACGCTGTTCATAACGCGGCGTAATACCGTTTTCTTCAGGCATGATCCAGATCTGATACAGATGCAGACGCTCGGTGCTGCTTGGGTTGTACTCTGAGTGACGAATACCCGTACCAGCACTCATAATCTGGAACTCACCCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3579236‑3579387 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑actcacccgcATAAAGATATGGAAATTTTGACCTACGTGCTGGAAGGTACTGTTGAGCATCAGGACAGCATGGGCAATAAAGAGCAGGTTCCGGCGGGTGAGTTCCAGATTATGAGTGCTGGTACGGGTATTCGTCACTCAGAGTACAACCCAAGCAGCACCGAGCGTCTGCATCTGTATCAGATCTGGATCATGCCTGAAGAAAAC  <  NC_000913/3579470‑3579274

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.