New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4635603	15 (0.410)	14 (0.380) +G	14/480	0.9	51.0%	coding (83/474 nt)	creA	PF05981 family protein CreA
	?	NC_000913	4638790 =	12 (0.330)					coding (613/1353 nt)	creD	putative inner membrane protein CreD

CGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGC  >  NC_000913/4638540‑4639034                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     cGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCg                                                                                                                                                                                                                                                      >  1:320643/1‑251 (MQ=255)       cGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCAGTGAACCCGGCACCGGGTTAGAGCAAGGCGTGCAGGGCTTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCg                                                                                                                                                                                                                                                <  1:185980/251‑1 (MQ=255)                tttttgctgttttagcattcagtgtttttttcttagtagattatgttttagggcttgacggaaataaaagtattgagattttgttcttaatcaatatgttatttaccgtgacgatctaattgctcgtgtaatagataaaaatggtaacaatatgaaatacaagcatttgatcctgtctttaagccttataatggtggggcgattggctcatgcagaagagtttggttcggtcgagGGTACCGGCGATCtt                                                                                                                                                                                                                                        <  2:230729‑M2/15‑1 (MQ=255)                 cgTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCtttt                                                                                                                                                                                                                                      >  1:87417/1‑250 (MQ=255)                                attcagggtttttttcttagtagagtatgttttagggcttgacggaaataaaagtattgagattttgttcttaatcaatatgttatttaccgtgacgaactaattgctcgtgtaatagataaaaatggtaacaatatgaaatacaagcatttgatcctgtctttaagcctgataatgctggggccattggctcatgcagaagagattggttcggtcgagGGTACCGGCGATCTTTCTGTGGTGCCTggcgg                                                                                                                                                                                                                       <  2:160034‑M2/32‑1 (MQ=255)                                                                    gtgATTAGCGGTGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGcc                                                                                                                                                                                   <  1:89768/251‑1 (MQ=255)                                                                     ggatgaGCGTCTGGGATGCGCGTGTTATTGGTGTTGTGAAAGCGGCTGAAGTTAACGGAACGGCTCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGGGCATATCCCTTGACCTGAAGTGGAGTGGCGGAAGCAGAACGTGAAGCTGAATATGGCCGTTAATTTAAGCGGTACCGGCGATCCTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGATCTTTACCAGCAACTGGCCg                                                                                                                                                                                  <  2:251936/246‑1 (MQ=255)                                                                          aGCGTCGGGGAGGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATcc                                                                                                                                                                             <  2:320643/251‑1 (MQ=255)                                                                                         tgagattttgttcttaatcaatatgttatttaccgtgacgaactaattgctcgtgtaatagaaaaaaatggtaacaatatgaaatacaagcatttgatcctgtctttaagcctgataatgctggggccattggctcatgcagaagagattggttcggtcgagGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATAACCTTAACCAGCAACTGGCCGCCTCCCAGTTTTTTAGGTg                                                                                                                                                               >  2:126190‑M2/163‑250 (MQ=255)                                                                                           agattttgttcttaatcaatatgttatttaccgtgacgacccaattgcccgtgtaatagataaaaatggtaacaatatgaaatacaagcatttgatcctgtctttaagcctgataatgctggagccattggctcatgcagaagagattggttcggtcgagGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGAtt                                                                                                                                                            >  1:231847‑M2/161‑251 (MQ=255)                                                                                                              tatgttatttaccgtgacgaaataattgctcgggtaatagatcaaactggtaacaatatgaactacaagcatttgatcctgtctttaagcctgataatgctggggccattggcaaatacagaagagattgattcggtcgagGGTACCGGCGATCTTTCTGTGGGGCCTGGAGGGCGGAATACCGAAATGACCTTAACCAGCAAATGGCCGCATACCAGTTTTTTAGGGGATTTTATACCAGCCAAACggg                                                                                                                                          >  2:327963‑M2/142‑250 (MQ=255)                                                                                                                                                       taaaaatggtaacaatatgaaatacaagcatttgatcctgtctttaagcctgataatgctggggccattggctcatgcagaagagattggttcggtcgagGGTACCGGCGATCTTTCTGTGGTGCCTGGCGAGAGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGGTTTTCTACCAGCCAAACGGGGAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCt                                                                                                >  1:167628‑M2/101‑251 (MQ=255)                                                                                                                                                                    aatatgaaatacaagcatttgatcctgtctttacgcctgatactgctgggaccaatggctcatgcagaagagattggttcggtcgaaGGTACCGGCGATCTTTCTGTGGTGACTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGCGATTTTCTACCAGCCAAACGGGAAGCTAGCGAGACAGGATTTCAGGCGCACCGGCAAAGCAGCTGGTTTGCtaata                                                                                    >  2:293366‑M2/88‑250 (MQ=255)                                                                                                                                                                          gggCAGGGCGTGCATATCCCTTGACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTACCCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCATGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAAGCTCgg                                                                             <  2:83479/251‑1 (MQ=255)                                                                                                                                                                          aaatacaagcatttgatcctgtctttaagcctgataatgctggggccattggctcatgcagaagagattggttcggtcgagGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACGGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCgg                                                                             <  1:251567‑M2/170‑1 (MQ=255)                                                                                                                                                                                     gCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTg                                                                  <  2:87417/251‑1 (MQ=255)                                                                                                                                                                                                tACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCaa                                                        <  1:368242/250‑1 (MQ=255)                                                                                                                                                                                                  ttaagcctgataatgctggggccattggctcatgcagaagagattggttcggtcgagGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGGCCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGGGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATg                                                      >  2:122676‑M2/58‑250 (MQ=255)                                                                                                                                                                                                                 ctggggccattggctcatgcagaagagattggttcggtcgagGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGaaa                                       >  1:103406‑M2/43‑250 (MQ=255)                                                                                                                                                                                                                  ggAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTACTAATTTCGGTGAGCGTTTTTCTTCAGGCCATGATACCCGCTGGGAAAAc                                     >  2:219285/1‑251 (MQ=255)                                                                                                                                                                                                                           aCCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGc                                                                                                                                  <  1:338039/149‑1 (MQ=255)                                                                                                                                                                                                                           aCCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGc                                                                                                                                  >  2:338039/1‑149 (MQ=255)                                                                                                                                                                                                                                        agagattggttcggtcgagGGTAACGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCCAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGc                 >  1:262900‑M2/20‑249 (MQ=255)                                                                                                                                                                                                                                             ttggttcggtcgagGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTttagcgtcgggatcggg          >  2:119832‑M2/15‑243 (MQ=255)                                                                                                                                                                                                                                             ttggttcggtcgagGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTGGCGTCg                  <  1:119832‑M2/229‑1 (MQ=255)                                                                                                                                                                                                                                                     gtcgagGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGc  <  1:9129‑M2/245‑1 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     CGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGC  >  NC_000913/4638540‑4639034

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