New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3683801 =	30 (0.730)	5 (0.130)	5/272	4.8	15.2%	coding (1817/1989 nt)	yhjK	cyclic‑di‑GMP phosphodiesterase
	?	NC_000913	3683817 =	27 (0.690)					coding (1801/1989 nt)	yhjK	cyclic‑di‑GMP phosphodiesterase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAGCGCCACCCGAAC  >  NC_000913/3683661‑3683943                                                                                                                                             |                                                                                                                                               tttCAATAGGGAGTGGGCGAGCAAAGAGGTAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTgca                                                                                                                                              <  1:299585/143‑1 (MQ=255) tttCAACAGGGAGTGGGCTAGCAAAAAGGAAGCCCGGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTGGTGCTTCAGTCTCCACGCCTTCGGCAATCATTCGTAAGTTCAGGCTCTGCGCCAGCATGATAATTgca                                                                                                                                              <  1:523896/143‑1 (MQ=255)    cAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTgcagca                                                                                                                                           <  1:81710/143‑1 (MQ=255)           gattGGGCGAGCAAGAAGTAAGGCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCAtgc                                                                                                                                    <  1:305022/140‑1 (MQ=255)             gTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCGTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCAtgctg                                                                                                                                  <  1:663655/143‑1 (MQ=255)              tGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATCGCAGCAATCAtgctgc                                                                                                                                 >  2:456043/1‑143 (MQ=255)                               cgCCCTCCGCAATCCCACCGCCCGCGTTTGCCCGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCCGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCAGCAAGCCt                                                                                                                <  1:386048/142‑1 (MQ=255)                                   cTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTcaa                                                                                                            <  2:118084/143‑1 (MQ=255)                                         aaTACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCATCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCaacaaaca                                                                                                      >  2:307998/1‑143 (MQ=255)                                                       cTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCCGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTc                                                                                        <  2:12263/143‑1 (MQ=255)                                                        ttttGCCAGACAGTCATCTTGTGCTTCAGTCTCAATGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGTAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCa                                                                                       <  1:399874/143‑1 (MQ=255)                                                              cagccagTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTCTCTCCCGGCCAGCCTTCGACAACCATTTTGTCGACTTCCAGTACGt                                                                                 >  2:335191/1‑143 (MQ=255)                                                                   attggcgTTTGGCTGCAGTCTCCACGCCCTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGAtt                                                                            <  1:268499/139‑1 (MQ=255)                                                                     tCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTCGTCGACTTTCAGTACGTCGACTgg                                                                          >  1:556782/1‑143 (MQ=255)                                                                       gcgTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCa                                                                        >  2:677029/1‑143 (MQ=255)                                                                             tgCTTCAGTCTCCACGCCTTCGGCAATCATGTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGAtt                                                                  <  2:604950/143‑1 (MQ=255)                                                                                ttCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGCTTTTCACTACGGCGATTGGCAACCATCTCa                                                               >  1:140211/1‑143 (MQ=255)                                                                                          tcGCCTTCTGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGTAACGATTTCATATGCTGCAg                                                     <  1:143203/142‑1 (MQ=255)                                                                                           ctccagctgCAATCATGCTTGAGTTCATGCTCTGCGCCAGCATCATCATTGCTGCAATCATGTTGCTTTCTCCCGGCATGCCTTCAACGATCATTTTGTCGATTTTCAGTACGTTGATTGGTAACGATTTCATATGCTGCAGc                                                    <  1:231171/135‑1 (MQ=255)                                                                                                 cGGCAATCATGTGTAAGTGCAGGCCCTGCGCCAGCCTGACAATCCAACGGGGGGGGCGTTTGTCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTGTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGc                                              <  2:133606/143‑1 (MQ=255)                                                                                                         aTTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGc                                      <  1:307998/143‑1 (MQ=255)                                                                                                          tttGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATCTGCTGCCGCTGCCGCAGCCCTGCg                                     >  2:530920/1‑143 (MQ=255)                                                                                                            tGTAAGTTCAGGCTCTGCGCGAGCATGATAATTGTAGCAATCATGCTGCTATCTCCCGGCAAGACTTCAACAAACATTTGGTCGATTTTCAGTACGTCGATTGGCAACGATTTAATATGCTGCAGCTGACGCAGCCCTGCGTa                                   <  1:174551/143‑1 (MQ=255)                                                                                                             gTACGTTCAGACGCTTCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTCCGTCGATTGGCAACGCTTTCATATGCTGCAGCTGACGCCGCCCTCCGCAg                                  >  1:258749/1‑143 (MQ=255)                                                                                                               aatttCAGGCTCCCCGCCTCCTTGAGAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGcc                                <  2:188909/140‑1 (MQ=255)                                                                                                               aaGTTCAGGCTCCGCGCCAGCATGATAATTGCAGCAATCATGCTGCCATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGcc                                <  1:333243/143‑1 (MQ=255)                                                                                                                                  gCAGGATAATTGGAGCAAGCGTGCTGCTATCTCCCGGCAAGCCTTCAACAAACAGTTTGTCGATTTTCAGGACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAgc             <  2:219227/143‑1 (MQ=255)                                                                                                                                       attagtGCAGCAGTCTTTCTGCTGTCTCCCGGCAAGCCTTTAACGAGCATGTTGTCGATTTTCAGTACGTCGCTTGGCTACGATTTCAGCTGCTGTAGCTCACGCAGCCCTGCGTATCCCATGCCGAAATCATCCAGCGCCAc        <  1:192605/138‑1 (MQ=255)                                                                                                                                       aTAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAACGCCAc        >  2:593807/1‑143 (MQ=255)                                                                                                                                         aaTTGCAGCAATCATGCTGCTATCTCCCGGCTAGCCTTCAACAAACATTTTGTCGATTTTCAGGACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAGCGCCAccc      <  2:506409/143‑1 (MQ=255)                                                                                                                                             gcagcaATCATGCTGCTATCTCCCGGCAAGCCTTCCACACACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCCTATGCTGCAGCTGACGCAGCCCTGCGTAGCCCCTGCCGAAATCATCCCCCGCCACCCGAAc  >  2:33198/1‑143 (MQ=255)                                                                                                                                             |                                                                                                                                               TTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAATCATTTGTAAGTTCAGGCTCTGCGCCAGCATGATAATTGCAGCAATCATGCTGCTATCTCCCGGCAAGCCTTCAACAAACATTTTGTCGATTTTCAGTACGTCGATTGGCAACGATTTCATATGCTGCAGCTGACGCAGCCCTGCGTAGCCCATGCCGAAATCATCCAGCGCCACCCGAAC  >  NC_000913/3683661‑3683943

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.