New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	5234910 =	16 (0.400)	3 (0.080)	3/318	NT	12.0%	coding (244/2325 nt)	fecA	outer membrane ferric citrate porin
	?	AE015451	5235001 =	29 (0.770)					coding (335/2325 nt)	fecA	outer membrane ferric citrate porin

CTGCTGCTCACCACCACCCTCGGCCTGGCCATGCCCGCCGCCCAGGCCGAAACCCGCAGCTACCACATCGCCGCCGGTTCTCTTGAAGACGCCCTCAACCAGTTCGGCCGCGAAAGCGGTGCACTGATCTCGTTCGGCTCGCAACTGACCCAGGGCATGAGCACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTAC  >  AE015451/5234700‑5235124                                                                                                                                                                                                                   |                                                                                                                                                                                                                       agttcagcgccaggtcatggctgccggtaccgttgttgtcgggggcgttgaccccggggatgcggttgagcacttcgcgggcggtggtagcaccgttgcgttcgaattcctcacggcgcactacgtcgcgggcgccggggtgctcgaatacgttgtcctgctgcgcctcggccagccagtcgccgaccacggtggaggcgcctagctcgaCTGCTGCGCGGCAGCGGCCTGCaa                                                                                                                                                                                                 >  1:475361‑M1/211‑234 (MQ=255)       ctCACCACCACCCTCGGCCTGGCCATGCCCGCCGCCCAGGCCGAAACCCGCAGCTACCACATCGCCGCCGGTTCTCTTGAAGACGCCCTCAACCAGTTCGGCCGCGAAAGCGGTGCACTGATCTCGTTCGGCTCGCAACTGACCCAGGGCATGAGCACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCCCGACGCCctgcgg                                                                                                                                                                                                                   >  1:180697/1‑208 (MQ=255)               caccCTCGGCCTGGCCATGCCCGCCGCCCAGGCCGAAACCCGCAGCTACCACATCGCCGCCGGTTCTCTTGAAGACGCCCTCAACCAGTTCGGCCGCGAAAGCGGTGCACTGATCTCGTTCGGCTCGCAACTGACCCAGGGCATGAGCACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCTCGACGCCctgct                                                                                                                                                                                                                    <  2:630272/201‑1 (MQ=255)                               aTGCCCGCCGCCCAGGCCGAAACCCGCAGCTACCACATCGCCGCCGGTTCTCTTGAAGACGCCCTCAACCAGTTCGGCCGCGAAAGCGGTGCACTGATCTCGTTCGGCTCGCAACTGACCCAGGGCATGAGCACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGa                                                                                                                                                                                      >  1:219973/1‑215 (MQ=255)                                       cgccCAGGCCGAAACCCGCAGCTACCACATCGCCGCCGGTTCTCTTGAAGACGCCCTCAACCAGTTCGGCCGCGAAAGCGGTGCACTGATCTCGTTCGGCTCGCAACTGACCCAGGGCATGAGCACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCa                                                                                                                                                         <  1:694411/236‑1 (MQ=255)                                                      ccGCAGCTACCACATCGCCGCCGGTTCTCTTGAAGACGCCCTCAACCAGTTCGGCCGCGAAAGCGGTGCACTGATCTCGTTCGGCTCGCAACTGACCCAGGGCATGAGCACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCa                                                                                                                                                                                              >  1:564074/1‑184 (MQ=255)                                                          ggatgcggttgagcacttcgcgggcggtggtagcaccgttgcgttcgaattcctcacggcgcactacgtcgcgggcgccggggtgctcgaatacgttgtcctgctgcgcctcggcaagccagtcgccgaccacggtggaggcgcctagctcgaCtgctgc                                                                                                                                                                                                                  >  1:334610‑M1/154‑160 (MQ=255)                                                                      cccgcccGTTCTCTTGAAGACGCCCTCAACCAGTCCGGCCGCGAAAGCGGTGCACTGATCTCGTCCGGCTCGCACCTGACCCAGGGCATGAGCACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGc                                                                                                                                                                   <  2:231145/194‑1 (MQ=255)                                                                                           gCCCTCAACCAGTTCGGCCGCGAAAGCGGTGCACTGATCTCGTTCGGCTCGCAACTGACCCAGGGCATGAGCACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGt                                                                                                     >  1:723847/1‑236 (MQ=255)                                                                                              cTCAACCAGTTCGGCCGCGAAAGCGGTGCACTGATCTCGTTCGGCTCGCAACTGACCCAGGGCATGAGCACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTc                                                                                                    <  2:479613/234‑1 (MQ=255)                                                                                                  gcgttcgaattcctcacggcgcactacgtcgcgggcgccggggtgctcgaatacgttgtcctgctgcgcctcggccagccagtcgccgaccacggtggaggcgcctagctcgaCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGc                                                                                                                                                                   <  2:314458‑M1/54‑1 (MQ=255)                                                                                                  gcgttcgaattcctcacggcgcactacgtcgcgggcgccggggtgctcgaatacgttgtcctgctgcgcctcggccagccagtcgccgaccacggtggaggcgcctagctcgaCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGc                                                                                                                                                                   >  1:314459‑M1/114‑167 (MQ=255)                                                                                                                 aaaGCGGTGCACTGATCTCGTTCGGCTCGCAACTGACCCAGGGCATGAGCACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGgcg                                                                               <  1:143365/236‑1 (MQ=255)                                                                                                                            cTGATCTCGTTCGGCTCGCAACTGACCCAGGGCATGAGCACTCAAGGCCTGGACGGTCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGa                                                                       <  1:58356/233‑1 (MQ=255)                                                                                                                                   cGTTCGGCTCGCAACTGACCCAGGGCATGAGCACTCACGGCCTGGACGGCCAGTCCGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATc                                                                                                                                                       <  2:105815/146‑1 (MQ=255)                                                                                                                                   cGTTCGGCTCGCAACTGACCCAGGGCATGAGCACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCAt                                                                                                                                                        <  2:616347/145‑1 (MQ=255)                                                                                                                                                    aCCCAGGGCATGAGCACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCg                                                                                    <  2:219973/196‑1 (MQ=255)                                                                                                                                                       cAGGGCATGAGCACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGCGGTCGGCGACTGGCTGGCCGAGGCGCCGCAGGACAACGTAGTc                                                             >  1:121176/1‑216 (MQ=255)                                                                                                                                                                  cACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCCGGACAACGTATTCGAGCAc                                                       >  1:262155/1‑211 (MQ=255)                                                                                                                                                                                   gCCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCgg                                                                                                                                                                                                              <  2:589126/43‑1 (MQ=255)                                                                                                                                                                                               cgcgCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTac  <  1:279436/235‑1 (MQ=255)                                                                                                                                                                                                     aGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGGATTCGAGCACCCCg                                                   >  1:241287/1‑180 (MQ=255)                                                                                                                                                                                                               cGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGa                       >  2:480669/1‑198 (MQ=255)                                                                                                                                                                                                                   |                                                                                                                                                                                                                       CTGCTGCTCACCACCACCCTCGGCCTGGCCATGCCCGCCGCCCAGGCCGAAACCCGCAGCTACCACATCGCCGCCGGTTCTCTTGAAGACGCCCTCAACCAGTTCGGCCGCGAAAGCGGTGCACTGATCTCGTTCGGCTCGCAACTGACCCAGGGCATGAGCACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTAC  >  AE015451/5234700‑5235124

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