New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3683802 =	29 (0.740)	5 (0.130)	5/272	4.7	15.7%	coding (1817/1989 nt)	yhjK	cyclic‑di‑GMP phosphodiesterase
	?	NC_000913	3683818 =	26 (0.690)					coding (1801/1989 nt)	yhjK	cyclic‑di‑GMP phosphodiesterase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAGCGCCACCCGAACTCCGGCATTGCGCAGC  >  NC_000913/3683692‑3683960                                                                                                                               |                                                                                                                                               cgCCCTCCGCAATCCCACCGCCCGCGTTTGCCCGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCCGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCAGCAAGCCt                                                                                                                                <  2:386048/142‑1 (MQ=255)     cTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTcaa                                                                                                                            <  1:118084/143‑1 (MQ=255)           aaTACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCATCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCaacaaaca                                                                                                                      >  1:307998/1‑143 (MQ=255)                         cTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCCGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTc                                                                                                        <  1:12263/143‑1 (MQ=255)                          ttttGCCAGACAGTCATCTTGTGCTTCAGTCTCAATGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGTAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCa                                                                                                       <  2:399874/143‑1 (MQ=255)                                cagccagTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTCTCTCCCGGCCAGCCTTCGACAACCATTTTGTCGACTTCCAGTACGt                                                                                                 >  1:335191/1‑143 (MQ=255)                                     attggcgTTTGGCTGCAGTCTCCACGCCCTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGAtt                                                                                            <  2:268499/139‑1 (MQ=255)                                       tCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTCGTCGACTTTCAGTACGTCGACTgg                                                                                          >  2:556782/1‑143 (MQ=255)                                         gcgTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCa                                                                                        >  1:677029/1‑143 (MQ=255)                                               tgCTTCAGTCTCCACGCCTTCGGCAATCATGTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGAtt                                                                                  <  1:604950/143‑1 (MQ=255)                                                  ttCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGCTTTTCACTACGGCGATTGGCAACCATCTCa                                                                               >  2:140211/1‑143 (MQ=255)                                                            tcGCCTTCTGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGTAACGATTTCATATGCTGCAg                                                                     <  2:143203/142‑1 (MQ=255)                                                                   cGGCAATCATGTGTAAGTGCAGGCCCTGCGCCAGCCTGACAATCCAACGGGGGGGGCGTTTGTCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTGTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGc                                                              <  1:133606/143‑1 (MQ=255)                                                                           aTTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGc                                                      <  2:307998/143‑1 (MQ=255)                                                                            tttGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATCTGCTGCCGCTGCCGCAGCCCTGCg                                                     >  1:530920/1‑143 (MQ=255)                                                                              tGTAAGTTCAGGCTCTGCGCGAGCATGATAATTGTAGCAATCATGCTGCTATCTCCCGGCAAGACTTCAACAAACATTTGGTCGATTTTCAGTACGTCGATTGGCAACGATTTAATATGCTGCAGCTGACGCAGCCCTGCGTa                                                   <  2:174551/143‑1 (MQ=255)                                                                               gTACGTTCAGACGCTTCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTCCGTCGATTGGCAACGCTTTCATATGCTGCAGCTGACGCCGCCCTCCGCAg                                                  >  2:258749/1‑143 (MQ=255)                                                                                 aatttCAGGCTCCCCGCCTCCTTGAGAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGcc                                                <  1:188909/140‑1 (MQ=255)                                                                                 aaGTTCAGGCTCCGCGCCAGCATGATAATTGCAGCAATCATGCTGCCATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGcc                                                <  2:333243/143‑1 (MQ=255)                                                                                                    gCAGGATAATTGGAGCAAGCGTGCTGCTATCTCCCGGCAAGCCTTCAACAAACAGTTTGTCGATTTTCAGGACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAgc                             <  1:219227/143‑1 (MQ=255)                                                                                                         attagtGCAGCAGTCTTTCTGCTGTCTCCCGGCAAGCCTTTAACGAGCATGTTGTCGATTTTCAGTACGTCGCTTGGCTACGATTTCAGCTGCTGTAGCTCACGCAGCCCTGCGTATCCCATGCCGAAATCATCCAGCGCCAc                        <  2:192605/138‑1 (MQ=255)                                                                                                         aTAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAACGCCAc                        >  1:593807/1‑143 (MQ=255)                                                                                                           aaTTGCAGCAATCATGCTGCTATCTCCCGGCTAGCCTTCAACAAACATTTTGTCGATTTTCAGGACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAGCGCCAccc                      <  1:506409/143‑1 (MQ=255)                                                                                                               gcagcaATCATGCTGCTATCTCCCGGCAAGCCTTCCACACACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCCTATGCTGCAGCTGACGCAGCCCTGCGTAGCCCCTGCCGAAATCATCCCCCGCCACCCGAAc                  >  1:33198/1‑143 (MQ=255)                                                                                                                  gcaATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAGCGCCACCCGAACTcc               >  2:607896/1‑143 (MQ=255)                                                                                                                    actCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAGCGCCACCCGAACTCCgg             <  1:40034/141‑1 (MQ=255)                                                                                                                         tgctgcTTTCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTTGGTGCGCCCACCGCCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAGCGCCACCCGAACTCCGGCATTg        <  1:710876/143‑1 (MQ=255)                                                                                                                               tATCTCCCGGCGAGCCGTCAACGATCATTTTGTCGATTTTCGGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAGCGCCACCCGAACTCCGGCATTGCGCAGc  <  2:174285/143‑1 (MQ=255)                                                                                                                               |                                                                                                                                               AGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAGCGCCACCCGAACTCCGGCATTGCGCAGC  >  NC_000913/3683692‑3683960

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