New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2667128 =	32 (0.600)	3 (0.060)	3/266	6.6	10.9%	coding (719/1140 nt)	hcaT	putative 3‑phenylpropionic transporter
	?	NC_000913	= 3026941	19 (0.380)					coding (1176/1320 nt)	guaD	guanine deaminase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACC  >  NC_000913/2666990‑2667270                                                                                                                                           |                                                                                                                                               gCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTCGCGCCACCACCCCCACCGAACAAAACTACCCCAc                                                                                                                                            >  2:90392/1‑143 (MQ=255)   gcccATAATGCCCCAGCGCACTACGCCGCAAAGCGCCGGGATCACCAGCACATCGCGTGCACCGCACCGGCGGAAAAGATTATTACACAGCGCACAGAAAATGACATCCGCCACCCCGCCCCCCGAGCACAAacaccacacca                                                                                                                                          >  2:817740/2‑142 (MQ=255)      cATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCCGACACGTTTACTACCCAGCGCAACGATAACGCCCTCCGCCACCACGCCCAGCCACCACAGACCCCCCAccgccc                                                                                                                                       >  2:954937/1‑142 (MQ=255)         aaTGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCgccgagg                                                                                                                                    <  1:915345/143‑1 (MQ=255)         aaTGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATAGCGGGTGCACTACAACACCGGAAAAGGTTATGAATCAGCGCACAGATAACGACCGCCGCCACCACGCCCCGCGGCCCCACACACCGCACCgccgaca                                                                                                                                    >  2:239031/1‑141 (MQ=255)               ccAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGCAACGATTATTACTCCGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAACTACCCCACCgccgggcccgagc                                                                                                                              >  2:435755/1‑142 (MQ=255)                cAGCGCACTACGCCGCAAACCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTa                                                                                                                             <  1:334424/143‑1 (MQ=255)                             cgcAAATCGCCGCGACCAACACCATCTCGCCGGCACGCCCACGGCGCACAAGGCTATGACCCAGCGCAAAGATAGTGACTTCCGCCACCCCGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGctgcctgc                                                                                                                <  1:424158/143‑1 (MQ=255)                             cgcAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCAACGCCGAGAGCGAGTAGCCAGctgcctgc                                                                                                                >  2:231909/1‑143 (MQ=255)                                      ccGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGCGTAGCCAGCTGCCTGCCAGTAAATg                                                                                                       >  1:338752/1‑143 (MQ=255)                                      ccGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATg                                                                                                       >  1:332288/1‑143 (MQ=255)                                          gaTCAACCGCAGAGCGCGCGCACTACAACCGCGGAACAGTTTATTACTCACCGCAAAGACAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGgcgc                                                                                                   <  2:116560/143‑1 (MQ=255)                                          gaTCAACAGCATAGCGCGTGCACTACCACGGCGGAAAAGGTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGgcgc                                                                                                   <  2:183112/143‑1 (MQ=255)                                           aTCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGCGTAGCCAGCTGCCTGCCAGTAAATGGCGCt                                                                                                  >  1:515512/1‑143 (MQ=255)                                             cAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTaa                                                                                                >  1:607489/1‑143 (MQ=255)                                                      atCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCCGCGAAAAGGTAAATACGTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAAt                                                                                       <  1:401388/143‑1 (MQ=255)                                                       tCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGCGTAGCCAGCTGGCCGCCAGGAAATGGCGCTCCAACCGTaagc                                                                                      >  2:57149/1‑141 (MQ=255)                                                         gcgGGCACGACAACGGCGGACAAGTTTATTTCTCAGCGCAAAGATAATTACTTCCGACACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAgg                                                                                    <  2:914079/143‑1 (MQ=255)                                                                cTACAGCGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCGCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTTAATGGCGCTAAAACCGTAAGAGGCCGCAtg                                                                             >  1:673399/1‑143 (MQ=255)                                                                    aaCGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACGTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGCCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGcc                                                                         >  1:909773/1‑143 (MQ=255)                                                                     aCGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGccc                                                                        <  1:510215/143‑1 (MQ=255)                                                                                gTTTGTCACTCCGCGCAAAGACAATGACGCCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGGAGCCAGCTGCTTGCCAGTAATTGGCGCTAAAACCGTAATAGGCCGAATGTGCCCCCTGCAATaaa                                                             <  1:227923/143‑1 (MQ=255)                                                                                       aCTCAGCGCAAAGAAAATGACTTCCGCCACGACGCCCAGCGACAACAAATACACCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGCAAATGGCGCTAAAACCGTAAGACGCCGCATGTGCGTCATGCAATaaacaaaccc                                                      >  1:915020/1‑141 (MQ=255)                                                                                         tCAGCGCAAAGATAATGACTTCCGCCACCACGCCCCGCGACAACAATTACGCTGCCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAg                                                    >  1:655525/1‑143 (MQ=255)                                                                                              gcAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAg                                               >  1:332145/1‑143 (MQ=255)                                                                                                               cGCTACCCCGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGc                              <  2:434054/143‑1 (MQ=255)                                                                                                                   accacGCTCAGCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAAc                          <  1:954937/143‑1 (MQ=255)                                                                                                                       cGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCTCCCTGCAATAAACAAACGTAGGCCAGAAAGCGCCAGTTCTGGCGAACCAgc                      <  2:502102/143‑1 (MQ=255)                                                                                                                             ccGACCACAAATCCCCCACCGCCGAGCCCGAGTAGCCAGTTGCCTGCCAGTAAATGGCGCTAACACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGTGCCAGTTCTGGCGAACCAGCGCCAAc                <  2:909773/142‑1 (MQ=255)                                                                                                                                 ccACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACg            >  1:191235/1‑143 (MQ=255)                                                                                                                                  cacaAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGc           >  1:502953/1‑143 (MQ=255)                                                                                                                                   acatacaCCCCACCGCCGAGGCCGAGTAGCCAGCTGCTTGCCAGTAAATGTCGCTAACACCCTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGACAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCa          <  2:287807/137‑1 (MQ=255)                                                                                                                                         aacccGCCGCAGAGGCCGAGGAGCCCGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCaa    <  2:762822/141‑1 (MQ=255)                                                                                                                                           cccACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCCACCACCGTGGTAGGCCGCATGGGCCGCCTGCAATGAGCAAACGCAGGCCAGAAAGCGCCAGCGCTGGCGAACCAGCGCCAACCACGCAGACCAAcc  >  2:674956/1‑143 (MQ=255)                                                                                                                                           |                                                                                                                                               GCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACC  >  NC_000913/2666990‑2667270

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.