New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3467138	33 (1.100)	4 (0.140)	4/206	2.9	12.2%	intergenic (‑147/+22)	bfd/chiA	bacterioferritin‑associated ferredoxin/periplasmic endochitinase
	?	NC_000913	= 3467144	27 (0.970)					intergenic (‑153/+16)	bfd/chiA	bacterioferritin‑associated ferredoxin/periplasmic endochitinase
Rejected: Frequency below/above cutoff threshold.

ATACTTTCTTCAAAAATAAATGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGTTCTTTAGTCAGGCCGTGGAATTCCGGAGACAGGTTAGTGCCGCCCGGCAGGTCGCGCGCGATTGAC  >  NC_000913/3467004‑3467280                                                                                                                                             |                                                                                                                                         ataCTTTCTTCAAAAATAAATGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGgtatcgc                                                                                                                                   <  2:115868‑M2/148‑8 (MQ=255)            aaaaaTAAATGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTa                                                                                                                        <  1:115868/148‑1 (MQ=255)             aaaaTAAATGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTAtt                                                                                                                      >  8:58722/1‑149 (MQ=255)                                 ttATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTa                                                                                                                 >  6:91801/1‑134 (MQ=255)                                 ttATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTa                                                                                                                 <  5:91801/134‑1 (MQ=255)                                             cGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATAAGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGaa                                                                                                   <  1:165608/136‑1 (MQ=255)                                             cGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATAAGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGaa                                                                                                   >  2:165608/1‑136 (MQ=255)                                              ggAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAg                                                                                                                                          <  8:108134/96‑1 (MQ=255)                                                aaaGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCg                                                                                   >  4:42655/1‑149 (MQ=255)                                                        ttttttAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGc                                                                                    >  4:49371/1‑140 (MQ=255)                                                        ttttttAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGc                                                                                    <  3:49371/140‑1 (MQ=255)                                                         tttttAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTa                                                                                 <  6:111274/142‑1 (MQ=255)                                                         tttttAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTa                                                                                 >  5:111274/1‑142 (MQ=255)                                                                gCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAaca                                                                                                                              >  6:43166/1‑90 (MQ=255)                                                                gCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAaca                                                                                                                              <  5:43166/90‑1 (MQ=255)                                                                                tGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGTTCTTTAGTCAGGCCGt                                                    <  3:42655/148‑1 (MQ=255)                                                                                  aaTAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTaaa                                                                                                            <  2:70421/90‑1 (MQ=255)                                                                                                  tACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGTTCTTTAGTCAGGCCGTGGAATTCCGGAGACAg                                    >  3:84202/1‑146 (MQ=255)                                                                                                  tACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGTTCTTTAGTCAGGCCGTGGAATTCCGGAGACAg                                    <  4:84202/146‑1 (MQ=255)                                                                                                                 aTTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGAttt                                                                                              >  1:123830/1‑73 (MQ=255)                                                                                                                 aTTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGAttt                                                                                              <  2:123830/73‑1 (MQ=255)                                                                                                                 aTTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCg                                                                            >  7:59636/1‑91 (MQ=255)                                                                                                                 aTTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCg                                                                            <  8:59636/91‑1 (MQ=255)                                                                                                                     aCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCg                                                                            <  4:147907/87‑1 (MQ=255)                                                                                                                     aCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCg                                                                            >  3:147907/1‑87 (MQ=255)                                                                                                                     aCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGTTCTTTAGTCAGGCCGTGGAATTCCGGAGACAGGTTAGTGCCGCCCGGCAGGTc               <  5:165471/148‑1 (MQ=255)                                                                                                                       aaaaaaTGCGGCGATACTGGAAGgtatcgccgcattttttgtcaatgggctacttttggtagcccattatttttattcatatttctatt                                                                        <  5:149342‑M2/89‑67 (MQ=255)                                                                                                                       aaaaaaTGCGGCGATACTGGAAGgtatcgccgcattttttgtcaatgggctacttttggtagcccattatttttattcatatttctatt                                                                        >  6:149342‑M2/1‑23 (MQ=255)                                                                                                                                ggcgATACTGGAAGgtatcgccgcattttttgtcaatgggctacttttggtagcccattattt                                                                                         <  5:66679‑M2/63‑50 (MQ=255)                                                                                                                                ggcgATACTGGAAGgtatcgccgcattttttgtcaatgggctacttttggtagcccattattt                                                                                         >  6:66679‑M2/1‑14 (MQ=255)                                                                                                                                ggcgATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGTTCTTTAGTCAGGCCGTGGAATTCCGGAGACAGGTTAGTGCCGCCCGGCAGGTCGCGCGCGa       >  2:169057/1‑145 (MQ=255)                                                                                                                                ggcgATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGTTCTTTAGTCAGGCCGTGGAATTCCGGAGACAGGTTAGTGCCGCCCGGCAGGTCGCGCGCGa       <  1:169057/145‑1 (MQ=255)                                                                                                                                   gATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGTTCTTTAGTCAGGCCGTGGAATTCCGGAGACAGGTTAGTGCCGCCCGGCAGGTCGCGCGCGATTGAc  >  1:139381/1‑147 (MQ=255)                                                                                                                                             |                                                                                                                                         ATACTTTCTTCAAAAATAAATGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGTTCTTTAGTCAGGCCGTGGAATTCCGGAGACAGGTTAGTGCCGCCCGGCAGGTCGCGCGCGATTGAC  >  NC_000913/3467004‑3467280

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