New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3394456 =	31 (0.810)	3 (0.080)	3/270	5.3	9.2%	coding (1803/3801 nt)	yhdP	DUF3971‑AsmA2 domains protein
	?	NC_000913	3394495 =	30 (0.820)					coding (1764/3801 nt)	yhdP	DUF3971‑AsmA2 domains protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TCTTTCAGCGGTGTCTCATCAAAGTAAGGGCCAACGGCTTTACCCGGACCTTTAATGTCAGCGTCAATCAGCAGTTTTTCTTTTGAGTAGTCAGGGATCACTGCGGTAAGATTACTCGCGCGCACGCCGCCCAGATTAACGCCATCGGTTTTCATCCATAAACCGTCGTTAATAAAGTCCAGTTCAATATCA‑AGGTTAGTTAATGCAGGCCAGTCCGGCTGGAAGGCAAACTTCGCGTTGCGCAGCGGCACCAGCACTTCAAACTGACCTTCGTT  >  NC_000913/3394322‑3394596                                                                                                                                       |                                                                                                                                              tCTTTCAGCGGTGTCTCATCAAAGTAAGGGCCAACGGCTTTACCCGGACCTTTAATGTCAGCGTCAATCAGCAGTTTTTCTTTTGAGTAGTCAGGGATCCCTGCGGTAAGATTACTCGCCCGCACCCCGCTCCGCTTAACGcg                                                                                                                                       >  1:518033/1‑142 (MQ=255)              gctgaccaAAGTATTGGCCAACGGCTTTACCCGGACCTTTAAGGTCAGCGTCAATCAGCAGTTTTTGTTTTGAGTAGTCAGGGATCACTGCGGTAAGATTACTCGCGCGCACGCCGCCCAGATTAACGCCGTCGGTTTTCATc                                                                                                                          <  2:319955/137‑1 (MQ=255)                       aGTAAGGGCCAACGGCTTTACCCGGACCTTTAATGTCAGCGTCAATCAGCAGTTTTTCTTTTGAGTAGTCAGGGATCACTGCGGTAAGATTACTCGCGCGCACGCCGCCCAGATTAACGCCATCGGTTTTCATCCATAAACcg                                                                                                                 <  1:167074/143‑1 (MQ=255)                         tAATGGCCAACGGTTTTACCTGTTCCTTTAATGTCAGCGTCAATCAGCAGTTTTTCTTGTGAGGAGTCAGGGATCACTGCGGTAAGATTACTCGCGCGCACGCCGCCCAGATTAACGCCATCGGTTTTCATCCATAAACcgtc                                                                                                               <  1:469786/143‑1 (MQ=255)                         tAAGGGCCAACGGCTTAACGCGGACCATTAAAGACAGCGTCAATCAGCAGTTTTTCTTTTGAGTAGTCAGGGATCACTGCGGTAAGATTACTCGCGCGCACGCCGCCCAGATTAACGCCATCGGTTTTCATCCATAAACcgtc                                                                                                               <  1:120946/143‑1 (MQ=255)                                 aaCGGCTTTACCCGGACCTTTAATGTCAGCGTCAATCAGCAGTTTTTCTTTTGAGTAGTCAGGGATCACTGCGGTAAGATTACTCGCGCGCACGCCGCCCAGATTAACGCCATCGGTTTTCATCCATAAACCGTCGTtaataa                                                                                                       >  1:307013/1‑143 (MQ=255)                                              ggttcTTTAATGTCAGTGTTTATCTTCCGTTTTTCTTTTGAGTAGGCAGGGAGCACGGCGGTAAGACTACGCGCGCGCACGCCGCCCAGATTAACGCCATCGGTTTTCATCCATAAACCGTCGTTAATAAAGTCCAGTTCAat                                                                                          <  1:80477/139‑1 (MQ=255)                                                aCCTTTAATGTCAGCGTCAATCAGCAGTTTTTCTTTTGAGTAGTCAGGGATCACTGCGGTAAGATTACTCGCGCGCACGCCGCCCAGATTAACGCCATCGGTATTCATCCAGAAACCGTCGTTCATAAAGCCCAGTTCAatat                                                                                        >  2:374571/1‑143 (MQ=255)                                                   tttAATGTCAGCGTCAATCAGCAGTTTTTCTTTTGAGTAGTCAGTGATCACTGCGGTAAGATTACTCGCGCGCAAGCCGCCCAGATTAACGCCATCGGTTTTCATCCATAAAACGTCGTGAATAAAATTCAACTAAATATCa‑a                                                                                    >  2:191709/1‑143 (MQ=255)                                                   tttAATGTCAGCGTCAATCAGCAGTTTTTCTTTTGAGTAGTCAGGGATCACTGCGGTAAGATTACTCGCGCGCACGCCGCCCAGATTAACGCCATCGGTTTTCATCCATAAACCGTCGTTCATATCGTCCAGTTCAATATCa‑a                                                                                    >  1:419865/1‑143 (MQ=255)                                                        tGTCAGCGTCAATCAGCAGTTTTTCTTTTGAGTAGTCAGGGATCACTGCGGTAAGATTACTCGCGCGCACGCCGCCCAGATTAACGCCATCGGTTTTCATCCATCAACCGTCGTTAATACAGTCGAGTTCCATATCA‑CGgcta                                                                               >  2:655788/1‑143 (MQ=255)                                                        tGTCAGCGTCAATCAGCAGTTTTTCTTTTGAGTAGTCAGGGATCACTGCGGTAAGATTACTCGCGCGCACGCCGCCCAGATTAACGCCATCGGTTTTCATCCATAACCCGTCGTTAATACACTCCCGTTCAATATCA‑CGgttg                                                                               >  1:120993/1‑142 (MQ=255)                                                         gTCAGCGTCAATCAGCAGTTTTACTTTTGAGTAGTCAGGGATCACGGCGGTAAGATTACTCGCGCGCACGCCGCCCAGATTAACGCCATCGGTTTTCATCCATAAACCGTCGTTAATAAAGTCCAGTTCAATATCA‑AGgttag                                                                              <  1:355946/143‑1 (MQ=255)                                                              cGTCAATCAGAAGAATTTCTTTTGAGTAGTCAGGGATCACTGCGGTCAGATTACTCGCGCGCACGCCGCCCAGATTAACGCCATCGGTTTTCATCCATAAACCGTCGCTAATAAAATCCCGTTCAATATCC‑ACgttagttagt                                                                         >  2:88505/1‑141 (MQ=255)                                                               gTCTATCGGCAGTTTTTCTTGTGAGTAGTCAGGGATCACTGCGGTAAGATTACTCGCGCGCACGCCGCCCAGATTAACGCCATCGGTTTTCATCCATAAACCGTCGTTAATAAAGTCCAGTTCAATATCA‑AGGTTAGTTAATg                                                                        <  2:423700/143‑1 (MQ=255)                                                                           tttttCTTTTGAGTAGTCAGGGATCACTGCGGTAAGATTACTCGCGCGCACGCCGCCCAGATTAACGCCATCGGTTTTCATCCATAAACCGTCGTTACTAAAGTCCAGTTCAATATCA‑AGGTTAGCTAATGCAGGCCAGTCCg                                                            >  2:388517/1‑143 (MQ=255)                                                                              ttcttcTGAGTAGTCAGGGATCACTGCGGTAAGATTACTCGCGCGCACGCCGCCCAGATTAACGCCATCGGTTTTCATCCATAAACCGTCGTTCATAAAGTCCACTTCAATATCA‑AGGATAGTTAATGCAGGCCAGTCCGGCt                                                         >  1:434641/1‑143 (MQ=255)                                                                                             aGGGATCACTGAGGTAAGAATACTCGCGCGCACGCCGCACAGATTAATGCCATGGGGTTTCATCAATAAACCGTCGTTAATAAAGTCCAGTTCAATATAA‑AGGTTAGTTAATGCAGGGGAGTCCGGATGGAAGGCAAAATTcg                                          >  2:690715/1‑143 (MQ=255)                                                                                               ggttCACTGCTGTAAGATTACTCTCGCGCACGCCGCCTAGATTAACGCCAGCGGGTATCATCCATAAACCGTCGTTAATAAAGTCCAGTTCAATATCA‑AGGTTAGTTAATGCAGGCCAGTCCGGCTGGAAGGCAAACTTcgcg                                        <  1:426212/140‑1 (MQ=255)                                                                                                      tGCGGTAAGATTACTCGCGCGCACGCCGCCCAGATTAACGCCATCGGTTTTCATCCATAAACCGTCGTTAATAAAGGCCAGTTCAATATCA‑AGGTTAGTTAATGCACGCCAGTCCGGCTGGAAGGCAAAATTCGCGGTgcgcg                                 >  2:458856/1‑142 (MQ=255)                                                                                                        cGGTAAGATTACTCGCGCGCACGCCGCCCAGATTAACGCCATCGGTTTTCATCCATAAACCGTCGTTAATAAAGTCCAGTTCAATATCA‑TGGTTAGTTAATGCGTGCCAGTCCTGCTGGAAGGCAAACTTCGCGTTGCGCAGc                               >  1:30770/1‑143 (MQ=255)                                                                                                             aGATTACTCGCGCGCACGCCGCCCAGATTAACGCCATCGGTTTTCATCCATAAACCGTCGTTAATAAAGTCCAGTTCAATATCA‑AGGTTAGTTAATGCAGGCAAGTCCGGCTGGAAAGCAAACATCGCGTTGCGCAGCGGcgc                          >  2:385248/1‑141 (MQ=255)                                                                                                               aTTACTCGCGCGCACGCCACCCAGATTAACGCCATCGGGTTTCATCCATAAACCGTCGTTAATAAAGTCCAGTTCAATATGC‑AGGATAGTTAATGCACGCCAGTCCGGCTGGAAGGCACACTGCGCGTTACGCAGCGGCagcc                        >  2:346603/1‑140 (MQ=255)                                                                                                                 tACTCGCGCGCACGCCGCCCAGATTAACGCCATCGGTTTTCATCCATAAACCGTCGTTAATAAAGTCCAGTTCAATATCA‑AGGTTTGTTAATGCAGGCCAGTCCGGCTGGAAGGCAAACTTCGCGTTGCGCAGCGGCACCAGc                      >  1:378717/1‑143 (MQ=255)                                                                                                                          gcACGCCGCCCAGATTAACGCCATCGGTTTTCATCCATAAACCGTCGTTAATAAAGACCAGATCAATATCA‑AGGGTAGCTAATGCAGGCCAGTCAGGCTGGAAGGCAAACTTGCCGTTGCGCAGCGGCACCCGCCCTACCAAc             >  2:665843/1‑143 (MQ=255)                                                                                                                           cACGCCGCCTAGATTAACGCCATCGGTTTTCATCTATAAATCGTCGTGAATAAAGTCCAGTACAATATCA‑AGGTTAGTTAATGCAGGCCAGTCCGGCTGGAAGGCAAACTTCGCGTTGCGCAGCGGCACCAGCACTTCAAACt            <  1:407969/143‑1 (MQ=255)                                                                                                                            aCGCCGCCCAGATTAACGCCATCGGTGTTCATCCATAAACCGTCGTTAATAAAGTCCAGTTCAATATCA‑CGGTTAGTTAATGCAGGCCAGTCCGGCTGGAAGGCAAACTTCGCGTCGCGCAGCGGCACCAGCCCTTCAAACTg           >  2:728883/1‑143 (MQ=255)                                                                                                                             cgccgcCCAGATTAACGCCATCGGTTTTCATCCATAAACCGTCGTTAATAAAGTCCAGTTCAATATCA‑AGGTTAGTTAATGCAGGCCAGTCCGGCTGGAAGGAAAAGTAGACGTTTCTCAGCGGCACCAGCACTTCAAACTGa          <  2:398033/143‑1 (MQ=255)                                                                                                                                   ccAGATTAACGCCATCGGTTTTCATCCATATACCGTCGTTAATAAAGTCCAGTTCAATATCA‑AGGTTAGTTAATGCAGGCCAGTCCGGCTGGAAGACAAACTTCGCGTTGCGCAGCGGCACCAGCACTTCAAACTGACCTTCg    <  1:33020/143‑1 (MQ=255)                                                                                                                                    cataTTAACGCCATCGGTTTTCATCCATAAACCGTGGCTAATAAAGTCCAGTTCAAGATCA‑AGGTTAGTTAATGCAGGCCAGTCCGGCTGGAAGGCAAACTTCGCGTTGCGCAGCGGCACCAGCACTTCAAACTGACCTTCGt   <  1:645435/140‑1 (MQ=255)                                                                                                                                    cAGATTAACGCCATCGGTTTTCATCCATAAACCGTCGTTAATAAAGTCCAGTTCAATATCA‑AGGTTAGTTAATGCAGGCCAGTCCGGCTGGAAGGCAAACTTCGCGTTGCGCAGCGGCACCAGCACTTCAAACTGACCTTCGt   <  1:710469/143‑1 (MQ=255)                                                                                                                                     aGATTAACGCCATCGGTTTTCATCCATAAACCGTCGTTAATAAAGCACTGAT‑AAGACAACATGTTAGTTAATGCAGGCCAGTCCGGCTGGAAGGCAAACTTCGCGTTGCGCAGCGGCACCAGCACTTCAAACTGACCTTCGtt  <  1:727015/143‑1 (MQ=255)                                                                                                                                       |                                                                                                                                              TCTTTCAGCGGTGTCTCATCAAAGTAAGGGCCAACGGCTTTACCCGGACCTTTAATGTCAGCGTCAATCAGCAGTTTTTCTTTTGAGTAGTCAGGGATCACTGCGGTAAGATTACTCGCGCGCACGCCGCCCAGATTAACGCCATCGGTTTTCATCCATAAACCGTCGTTAATAAAGTCCAGTTCAATATCA‑AGGTTAGTTAATGCAGGCCAGTCCGGCTGGAAGGCAAACTTCGCGTTGCGCAGCGGCACCAGCACTTCAAACTGACCTTCGTT  >  NC_000913/3394322‑3394596

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.