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

GCCCTCAACCAGTTCGGCCGCGAAAGCGGTGCACTGATCTCGTTCGGCTCGCAACTGACCCAGGGCATGAGCACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGACAACAACGGTACCGGCAGCCATGACCTGGCGCTGAACTTCGGCATCCGC  >  AE015451/5234790‑5235221                                                                                                                                                                                                                    |                                                                                                                                                                                                                             gCCCTCAACCAGTTCGGCCGCGAAAGCGGTGCACTGATCTCGTTCGGCTCGCAACTGACCCAGGGCATGAGCACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGt                                                                                                                                                                                                      >  1:723847/1‑236 (MQ=255)    cTCAACCAGTTCGGCCGCGAAAGCGGTGCACTGATCTCGTTCGGCTCGCAACTGACCCAGGGCATGAGCACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTc                                                                                                                                                                                                     <  2:479613/234‑1 (MQ=255)                       aaaGCGGTGCACTGATCTCGTTCGGCTCGCAACTGACCCAGGGCATGAGCACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGgcg                                                                                                                                                                                <  1:143365/236‑1 (MQ=255)                                  cTGATCTCGTTCGGCTCGCAACTGACCCAGGGCATGAGCACTCAAGGCCTGGACGGTCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGa                                                                                                                                                                        <  1:58356/233‑1 (MQ=255)                                                          aCCCAGGGCATGAGCACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCg                                                                                                                                                                                     <  2:219973/196‑1 (MQ=255)                                                             cAGGGCATGAGCACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGCGGTCGGCGACTGGCTGGCCGAGGCGCCGCAGGACAACGTAGTc                                                                                                                                                              >  1:121176/1‑216 (MQ=255)                                                                        cACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCCGGACAACGTATTCGAGCAc                                                                                                                                                        >  1:262155/1‑211 (MQ=255)                                                                                                     cgcgCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTac                                                                                                   <  1:279436/235‑1 (MQ=255)                                                                                                           aGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGGATTCGAGCACCCCg                                                                                                                                                    >  1:241287/1‑180 (MQ=255)                                                                                                                     cGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGa                                                                                                                        >  2:480669/1‑198 (MQ=255)                                                                                                                            CTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGa                                                                                                                                                                                                                           <  2:8183‑M1/92‑1 (MQ=255)                                                                                                                                  ggCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGgcgc                                                                                                                                                >  2:209599/1‑161 (MQ=255)                                                                                                                                  ggCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGgcgc                                                                                                                                                <  1:209599/161‑1 (MQ=255)                                                                                                                                      gCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGGAGGACAACGTATTCGAGCACCCCGGCGCCCGCGCCggc                                                                                                                                      >  2:166172/1‑165 (MQ=255)                                                                                                                                       cGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTTAGCGCCCCCCCCGTTGTCGGCGACTCGCCGGCCGAGGCGCAGCAGGCCAACg                                                                                                                                                                   >  2:65644/1‑137 (MQ=255)                                                                                                                                         gCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGActggct                                                                                                                                                                                          <  2:269546/112‑1 (MQ=255)                                                                                                                                               aaGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAAc                                                           <  2:723846/233‑1 (MQ=255)                                                                                                                                                              gctgaacgcgttgtcggtttcctggcgtgcttgcaggccgctgccgcgcagcagTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGc                                                                                                              >  1:314459‑M2/55‑167 (MQ=255)                                                                                                                                                              gctgaacgcgttgtcggtttcctggcgtgcttgcaggccgctgccgcgcagcagTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGc                                                                                                              <  2:314458‑M2/113‑1 (MQ=255)                                                                                                                                                                  cAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGcc                                                                                                                               <  2:605846/146‑1 (MQ=255)                                                                                                                                                                        gTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGacaa                                                >  2:727938/1‑219 (MQ=255)                                                                                                                                                                                    aCCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAAc                                                                                                                <  2:537179/143‑1 (MQ=255)                                                                                                                                                                                     cccATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGaca                                                 <  2:269882/205‑1 (MQ=255)                                                                                                                                                                                     cccATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGACAACAACGGTACCGGCAGCCATGACCTGGCGCt                  >  1:267199/1‑236 (MQ=255)                                                                                                                                                                                           acaGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGac                                                  >  1:665304/1‑198 (MQ=255)                                                                                                                                                                                           acaGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGac                                                  <  2:665303/198‑1 (MQ=255)                                                                                                                                                                                           acaGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGGGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCggggc                                                               >  1:232488/1‑184 (MQ=255)                                                                                                                                                                                            ttgcaggccgctgccgcgcagcagTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGACAACAACGGTACCGGCAGCCATGACCTGGCGCTGAACt             >  1:475361‑M2/25‑234 (MQ=255)                                                                                                                                                                                             aGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTg                                                                                   <  1:166172/163‑1 (MQ=255)                                                                                                                                                                                                 ccTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAAc                                                                                                                <  1:455531/130‑1 (MQ=255)                                                                                                                                                                                                  ctgctgctgGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGc                                                                                                              <  2:365431/131‑1 (MQ=255)                                                                                                                                                                                                  ctgctgctgGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGc                                                                                                              >  2:293033/1‑131 (MQ=255)                                                                                                                                                                                                  ctgctgctgGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGc                                                                                                              <  1:293034/131‑1 (MQ=255)                                                                                                                                                                                                      tgctgGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGACAACAACGGTACCGGCAGCCATGACCTGGCGCTGAACTTCGGCATCcgc  >  1:275953/1‑235 (MQ=255)                                                                                                                                                                                                             gcagcagTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTTGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATcc                                                                      >  1:334610‑M2/8‑160 (MQ=255)                                                                                                                                                                                                                    |                                                                                                                                                                                                                             GCCCTCAACCAGTTCGGCCGCGAAAGCGGTGCACTGATCTCGTTCGGCTCGCAACTGACCCAGGGCATGAGCACTCAAGGCCTGGACGGCCAGTACGACACGCGCCAGGGGCTCGACGCCCTGCTGCGCGGCAGCGGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGACAACAACGGTACCGGCAGCCATGACCTGGCGCTGAACTTCGGCATCCGC  >  AE015451/5234790‑5235221

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