New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 3456325	29 (0.580)	3 (0.060)	3/316	NT	8.6%	coding (138/480 nt)	eco	Ecotin
	?	AE015451	= 3456502	37 (0.790)					intergenic (‑40/+16)	eco/hbdH	Ecotin/3‑hydroxybutyrate dehydrogenase

CACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGGTAGGGCGCATAAGGCCTCCTTGCTGATAAACGAATGCGCGACCGTCGCCAGGCCCAGGCCGCCCGTCACTGGGCCAACCAGCCACCATCGACGTTC  >  AE015451/3456112‑3456548                                                                                                                                                                                                                      |                                                                                                                                                                                                                                caccGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCggg                                                                                                                                                                                                           <  1:543741/236‑1 (MQ=255)   ccGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAg                                                                                                                                                                                                                            >  2:809382/1‑217 (MQ=255)                           tGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTgg                                                                                                                                                                                                                        <  2:934223/197‑1 (MQ=255)                                  cGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGGGAAGCCTTTTTCCGCCTCGGGGTAGGGcgcg                                                                                                                                                                           >  1:159223/1‑235 (MQ=255)                                  cGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCAGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCttt                                                                                                                                                                                                 <  1:151551/213‑1 (MQ=255)                                                        gccgggccgCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGGGGATCACCCGCCGGGGGAAGCCTTTTTCCGCCTCGGGGGAGGGCGCGAt                                                                                                                                                                         >  1:85207/1‑215 (MQ=255)                                                                            cAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTgcggcgg                                                                                                                                      >  1:350306/1‑230 (MQ=255)                                                                               gCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAgg                                                                                                                              <  1:779246/235‑1 (MQ=255)                                                                                 ggTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCt                                                                                                                            >  1:799065/1‑235 (MQ=255)                                                                                        aTTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAgg                                                                                                                              <  2:457553/226‑1 (MQ=255)                                                                                                cccAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTgcggcggcg                                                                                                                                    >  2:588303/1‑212 (MQ=255)                                                                                                       ttCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTgcggcggct                                                                                                                                    >  2:340773/1‑204 (MQ=255)                                                                                                                  gtgCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGGTTATAGGCCGACTCGTCGGCCTGGGCCGGCCGGTGGATAACCTGCCGGGTGAAGCCGTTTTCCGCCTCGGGGTCGGGCGCGATGTc                                                                                                                                                                      >  2:17176/1‑160 (MQ=255)                                                                                                                     cTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCtggcgacggtcgcgcattcgtttatcagcaaggaggccttatgcgccctaccccaatgaccgcgattctggcc                                                                                                                                                        <  2:596631‑M1/171‑74 (MQ=255)                                                                                                                     cTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCtggcgacggtcgcgcattcgtttatcagcaaggaggccttatgcgccctaccccaatgaccgcgattctggcc                                                                                                                                                        >  1:596634‑M1/1‑98 (MQ=255)                                                                                                                      ttCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAAt                                                                                                                >  1:791406/1‑210 (MQ=255)                                                                                                                      ttCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGGTAGGGCGCATa                                                                                       >  1:391693/1‑235 (MQ=255)                                                                                                                      ttCCAGGTTGCCGCCCAGGCGCTGACGATTACAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGGTAGGGCGCATaa                                                                                      >  1:530365/1‑236 (MQ=255)                                                                                                                          aGGTTGCCTGCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTAGAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAAcc                                                                                                                                                                                                                  <  2:59771/108‑1 (MQ=255)                                                                                                                                   cccAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGAt                                                                                                                                                                                                                      >  1:226141/1‑95 (MQ=255)                                                                                                                                   cccAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGAt                                                                                                                                                                                                                      <  2:226139/95‑1 (MQ=255)                                                                                                                                         cgcTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGGt                                                                                                 <  1:320000/206‑1 (MQ=255)                                                                                                                                              aCGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGc                                                                                                                                                                                           >  2:435987/1‑111 (MQ=255)                                                                                                                                              aCGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGc                                                                                                                                                                                           <  1:435990/111‑1 (MQ=255)                                                                                                                                                           aCCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAgg                                                                                                                              >  2:336895/1‑159 (MQ=255)                                                                                                                                                           aCCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGGTAGGGCGCATaa                                                                                      <  2:530362/199‑1 (MQ=255)                                                                                                                                                             cTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATcgcg                                                                                                            >  2:641094/1‑175 (MQ=255)                                                                                                                                                                cAGCGTCTTGCCAGCGAGGATTTCCTGTTTATAGGCCGACTCGTCGGCCTGGGCCtggcgacggtcgcgcattcgtttatcagcaaggaggccttatgcgccctaccccaatgaccgcgat                                                                                                                                                               <  2:61081‑M1/121‑67 (MQ=255)                                                                                                                                                                   cGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGCCTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGGTAGGGCGCATAAGGCCTCCTTGCTGATAAACGAATGCGCGa                                                         <  2:845421/220‑1 (MQ=255)                                                                                                                                                                          ccAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTgg                                                                                                    <  2:857989/170‑1 (MQ=255)                                                                                                                                                                                             tataGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGGTAGGGCGCa                                                                                         <  2:364844/162‑1 (MQ=255)                                                                                                                                                                                                        tcgtcgGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGGTAGGGCGCATAAGGCCTCCTTGCTGATAAACGAATGCGCGACCGTCGCCAGGCCCAGGCCGCCCGTCACTGGGCCAACCAGCCACCATCGa       >  2:692891/1‑233 (MQ=255)                                                                                                                                                                                                          gtcgGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGGTAGGGCGCATAAGGCCTCCTTGCTGATAAACGAATGCGCGACCGTCGCCAGGCCCAGGCCGCCCGTCACTGGGCCAACCAGCCACCATCGACGTTc  >  1:440685/1‑236 (MQ=255)                                                                                                                                                                                                             ggcctgggcctggcgacggtcgcgcattcgtttatcagcaaggaggccttatgcgccctaccccaatgaccgcgattctggccctgagcctggccgccgccgcacctgcgcgggccgccagcctgaaggacatcgcgccctaccccgaggcggaaaaaggctt                                                                        >  1:128898‑M1/1‑10 (MQ=255)                                                                                                                                                                                                                  gggcctggcgacggtcgcgcattcgattatcagcaaggaggccttatgcgccctaccccaatgaccgcgattctggccctgagcctggccgccgcc                                                                                                                                      >  2:735462‑M1/1‑5 (MQ=255)                                                                                                                                                                                                                      |                                                                                                                                                                                                                                CACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGGTAGGGCGCATAAGGCCTCCTTGCTGATAAACGAATGCGCGACCGTCGCCAGGCCCAGGCCGCCCGTCACTGGGCCAACCAGCCACCATCGACGTTC  >  AE015451/3456112‑3456548

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