New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 3797840	34 (1.470)	4 (0.180)	4/238	NT	11.4%	coding (197/1449 nt)	vdh	vanillin dehydrogenase
	?	NC_002947	= 3797929	30 (1.360)					coding (108/1449 nt)	vdh	vanillin dehydrogenase

GTCATGGCTGCAGCCTCGCGCAGCATGTTGGCGGCCAACTTCACGTTGAAGCCATACCAGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCAT  >  NC_002947/3797693‑3797987                                                                                                                                                    |                                                                                                                                                    gTCATGGCTGCAGCCTCGCGCAGCATGTTGGCGGCCAACTTCACGTTGAAGCCGTACCAGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTgcgc                                                                                                                                                   >  2:256786/1‑150 (MQ=255)        cTGCAGCCTCGCGCAGCATGTTGGCGGCCAACTTCACGTTGAAGCCATACCAGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGc                                                                                                                                            >  1:423292/1‑150 (MQ=255)             gCCTCGCGCAGCATGTTGGCGGCCAACTTCACGTTGAAGCCATACCAGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGt                                                                                                                                       <  1:402112/150‑1 (MQ=255)                          tGTTGGCGGCCAACTTCACGTTGAAGCCATACCAGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCATGCTGCTGCGCCGCTCGCCCGGTGCCAGGGCGGccc                                                                                                                          <  1:273845/150‑1 (MQ=255)                           gTTGGCGGCCAACTTCACGTTGAAGCCAGACCAGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCa                                                                                                                         >  1:545381/1‑150 (MQ=255)                                        ttCACGTTGAAGCCATACCAGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCa                                                                                                                         >  2:300020/1‑137 (MQ=255)                                        ttCACGTTGAAGCCATACCAGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCa                                                                                                                         <  1:300020/137‑1 (MQ=255)                                             gTTGAAGCCATACCAGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCt                                                                                                       >  2:219080/1‑150 (MQ=255)                                               tGAAGCCATACCAGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTgg                                                                                                     <  1:29249/150‑1 (MQ=255)                                                 aaGCCATACCAGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGgca                                                                                                   <  2:376653/150‑1 (MQ=255)                                                   gCCATACCAGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCcgccgccacactggccgatgccgatgccgcggtggctgctgccagcgcggcg                                                                                                 >  1:516775‑M1/1‑98 (MQ=255)                                                         ccAGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCAc                                                                                           >  1:373671/1‑150 (MQ=255)                                                          cAGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCAcc                                                                                          >  1:330881/1‑150 (MQ=255)                                                             ttGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACcgcg                                                                                       >  2:235972/1‑150 (MQ=255)                                                                 ccATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGcggcat                                                                                   <  2:8487/150‑1 (MQ=255)                                                                   aTGGCCCCGGTTTCACGGGCGGCGGCGATGAACTTGGCGGTCATCGCCTGCAACAACTAGGAGCCTGCTAAGAAGCGGCTGCACCGCTCGCCCGGCGCCAGGGAGGCCCAGGCCGGAAACCCCGCGCTGGCAGCAGCCACCGcggcctcc                                                                                 >  1:537779/1‑149 (MQ=255)                                                                   aTGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGACCCTCGCCTGCAGCATATCAGCGCCTGCCGGTAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGcggcatcg                                                                                 >  1:141353/1‑150 (MQ=255)                                                                        cccGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCggcatcggcatc                                                                            <  2:545381/150‑1 (MQ=255)                                                                        cccGGTTTCACCGGCGGCGGCGATGAACTCGGCGACCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCggcatcggcatc                                                                            >  2:502050/1‑150 (MQ=255)                                                                                           gcgATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGgca                                                         >  2:423292‑M2/1‑147 (MQ=255)                                                                                                 aaCTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGc                                                   <  2:564290/150‑1 (MQ=255)                                                                                                 aaCTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGc                                                   <  1:256786/150‑1 (MQ=255)                                                                                                    tCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGac                                                >  1:373371/1‑150 (MQ=255)                                                                                                            ccTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCcgccgccacactggccgatgccgatgccgcggtggctgctgccagcgcggcgtttccggcctgggccgccctggcaccgggcga                                                                 >  2:165977‑M1/1‑41 (MQ=255)                                                                                                            ccTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCcgccgccacactggccgatgccgatgccgcggtggctgctgccagcgcggcgtttccggcctgggccgccctggcaccgggcga                                                                 <  1:165977‑M1/125‑85 (MQ=255)                                                                                                               cGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGcc                                     <  1:219080/150‑1 (MQ=255)                                                                                                                      aaCAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGgcagccgcttgctggcaggcgctgatctgt                              <  2:516775‑M2/150‑31 (MQ=255)                                                                                                                       aCAGATCAGCGCCTGCCAGCAAGCGGCTGCcgccgccacactggccgatgccgatgccgcggtggctgctgccagcgcggcgtttccggcctgggccgccctggcaccgggcgagcggcgcagccgcttgctggcaggcgctgatctgtt                             <  2:516775‑M1/150‑121 (MQ=255)                                                                                                                       aCAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGt                                                                                                                                       <  2:592262/44‑1 (MQ=255)                                                                                                                       aCAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGt                                                                                                                                       >  1:592262/1‑44 (MQ=255)                                                                                                                         aGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCggg                           <  1:502050/150‑1 (MQ=255)                                                                                                                               gcgcCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGcggcatcg                                                                                 <  1:115293/90‑1 (MQ=255)                                                                                                                               gcgcCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGcggcatcg                                                                                 >  2:115293/1‑90 (MQ=255)                                                                                                                                       cAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGcc                                                                                                                           >  2:504967/1‑40 (MQ=255)                                                                                                                                       cAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGcc                                                                                                                           <  1:504967/40‑1 (MQ=255)                                                                                                                                         gCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAgg           >  1:19267/1‑150 (MQ=255)                                                                                                                                            ggCGGCTGCGCCGCTCGCGCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTGACAGCGCTCGAAGGTTc        <  2:537779/149‑1 (MQ=255)                                                                                                                                             gCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCg       >  2:540914/1‑150 (MQ=255)                                                                                                                                             gCGGCTGCGCCGCTCGCCCGGTGCCAAGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCg       >  1:401320/1‑150 (MQ=255)                                                                                                                                              cGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGt      <  2:373371/150‑1 (MQ=255)                                                                                                                                              cGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGt      <  1:235972/150‑1 (MQ=255)                                                                                                                                                gCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTcc    >  2:551031/1‑150 (MQ=255)                                                                                                                                                  tGCcgccgccacactggccgatgccgatgccgcggtggctgctgccagcgcggcgtttccggcctgggccgccctggcaccgggcgagcggcgcagccgcttgctggcaggcgctgatctgttgcaggcgagggccgccgagttcatcgc  >  2:423292‑M1/1‑3 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GTCATGGCTGCAGCCTCGCGCAGCATGTTGGCGGCCAACTTCACGTTGAAGCCATACCAGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCAT  >  NC_002947/3797693‑3797987

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ 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.