New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3555446	20 (0.580)	4 (0.130)	4/226	3.9	15.6%	coding (2362/2706 nt)	malT	mal regulon transcriptional activator
	?	NC_000913	= 3555468	25 (0.790)					coding (2384/2706 nt)	malT	mal regulon transcriptional activator
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGGCAGGCCGGACGTAAAAGTGACGCCCAGCGCGTGTTGCTGGACGCATTAAAACTGGCGAATCGCACCGGATTTATCAGCCATTTTGTCATCGAAGGCGAAGCGATGGCGCAACAACTGCGTCAGCTGATTCAGCTTAATACGCTGCCGGAACTGGAACAGCATCGCGCGCAGCGTATTCTGCGAGAAATCAATCAACATCATCGGCATAAATTCGCCCATTTCGATGAGAATTTCGTTGAACGTCTGCTAAATCATCCTGAAGTACCTGAACTGATCCGCACCAG  >  NC_000913/3555299‑3555585                                                                                                                                                    |                                                                                                                                            tGGCAGGCCGGACGTAAAAGTGACGCCCAGCGCGTGTTGCTGGACGCATTAAAACTGGCGAATCGCACCGGATTTATCAGCCATTTTGTCATCGAAGGCGAAGCGATGGCGCAACAACTGCGTCAGCTGATTCAGCTTAATACGCTGcc                                                                                                                                            >  1:174771/1‑149 (MQ=255)                aaaaGTGACGCCCAGCGCGTGTTGCTGGACGCATTAAAACTGGCGAATCGCACCGGATTTATCAGCCATTTTGTCATCGAAGGCGAAGCGATGGCGCAACAACTGCGTCAGCTGATTCAGCTTAATACGCTGCCGGAACTGGAACAGCa                                                                                                                             >  5:162653/1‑149 (MQ=255)                                           gACGCATTAAAACTGGCGAATCGCACCGGATTTATCAGCCATTTTGTCATCGAAGGCGAAGCGATGGCGCAACAACTGCGTCAGCTGATTCAGCTTAATACGCTGCCGGAACTGGAACAGCATCGCGCGCAGCGTATTCTGCGAGaa                                                                                                    <  6:162653/147‑1 (MQ=255)                                            aCGCATTAAAACTGGCGAATCGCACCGGATTTATCAGCCATTTTGTCATCGAAGGCGAAGCGATGGCGCAACAACTGCGTCAGCTGATTCAGCTTAATACGCTGCCGGAACTGGAACAGCATCGCGCGCAGCGTATTCTGCGAGAaatc                                                                                                 <  2:174771/149‑1 (MQ=255)                                                     aaCTGGCGAATCGCACCGGATTTATCAGCCATTTTGTCATCGAAGGCGAAGCGATGGCGCAACAACTGCGTCAGCTGATTCAGCTTAATACGCTGCCGGAACTGGAACAGCATCGCGCGCAGCGTATTCTGCGAGAAATCAATCAAcat                                                                                        >  1:120541/1‑149 (MQ=255)                                                                  cACCGGATTTATCAGCCATTTTGTCATCGAAGGCGAAGCGATGGCGCAACAACTGCGTCAGCTGATTCAGCTTAATACGCTGCCGGAACTGGAACAGCATCGCGCGCAGc                                                                                                                  <  8:17338/110‑1 (MQ=255)                                                                  cACCGGATTTATCAGCCATTTTGTCATCGAAGGCGAAGCGATGGCGCAACAACTGCGTCAGCTGATTCAGCTTAATACGCTGCCGGAACTGGAACAGCATCGCGCGCAGc                                                                                                                  >  7:17338/1‑110 (MQ=255)                                                                            aTCAGCCATTTTGTCATCGAAGGCGAAGCGATGGCGCAACAACTGCGTCAGCTGATTCAGCTTAATACGCTGCCGGAACTGGAACAGCATCGCGCGCAGCGTATTCTGCGAGAaatcaat                                                                                              >  1:17003/1‑120 (MQ=255)                                                                            aTCAGCCATTTTGTCATCGAAGGCGAAGCGATGGCGCAACAACTGCGTCAGCTGATTCAGCTTAATACGCTGCCGGAACTGGAACAGCATCGCGCGCAGCGTATTCTGCGAGAaatcaat                                                                                              <  2:17003/120‑1 (MQ=255)                                                                               aGCCATTTTGTCATCGAAGGCGAAGCGATGGCGCAACAACTGCGTCAGCTGATTCAGCTTAATACGCTGCCGGAACTGGAACAGCATCGCGCGCAGCGTATTCTGCGAGAAAtcaatcaa                                                                                           <  3:121192/120‑1 (MQ=255)                                                                               aGCCATTTTGTCATCGAAGGCGAAGCGATGGCGCAACAACTGCGTCAGCTGATTCAGCTTAATACGCTGCCGGAACTGGAACAGCATCGCGCGCAGCGTATTCTGCGAGAAAtcaatcaa                                                                                           >  4:121192/1‑120 (MQ=255)                                                                                                    gAAGCGATGGCGCAACAACTGCGTCAGCTGATTCAGCTTAATACGCTGCCGGAACTGGAACAGCATCGCGCGCAGCGTATTCTGCGAGAAATCAATCAACATCATCGGCATAAATTCGCCCATTTCGATGAGAATTTCGTTGAACGTCt                                         <  2:120541/149‑1 (MQ=255)                                                                                                             gcgcAACAACTGCGTCAGCTGATTCAGCTTAATACGCTGCCGGAACTGGAAc                                                                                                                                 <  6:104160/52‑1 (MQ=255)                                                                                                                    aacTGCGTCAGCTGATTCAGCTTAATACGCTGCCGGAACTGGAACAGCATCGCGCGCAGCGTATTCTGCGAGAAATCAATCAACATCATCGGCATAAATTCGCCCATTTCGATGAGAATTTCGTTGAACGTCt                                         <  1:38251/133‑1 (MQ=255)                                                                                                                    aacTGCGTCAGCTGATTCAGCTTAATACGCTGCCGGAACTGGAACAGCATCGCGCGCAGCGTATTCTGCGAGAAATCAATCAACATCATCGGCATAAATTCGCCCATTTCGATGAGAATTTCGTTGAACGTCt                                         >  2:38251/1‑133 (MQ=255)                                                                                                                                 gATTCAGCTTAATACGCTGCCGGAACTGGAACAGCATCGCGCGCAGCGTATTCTGCGAGAAATCAATCAACATCATCGGCATAAATTCGCCCATTTCGATGAGAATTTCGTTGAAc                                             <  2:188762/116‑1 (MQ=255)                                                                                                                                 gATTCAGCTTAATACGCTGCCGGAACTGGAACAGCATCGCGCGCAGCGTATTCTGCGAGAAATCAATCAACATCATCGGCATAAATTCGCCCATTTCGATGAGAATTTCGTTGAAc                                             >  1:188762/1‑116 (MQ=255)                                                                                                                                           aaTACGCTGCCGGAACTGGAACAGCATCGCGCGCAGCGTATTCTGCGAGAAATCAATCAACATCATCGGCATAAATTCGCCCATTTCGATGAGAATTTCGTTGAACGTCTGCTAAATCATCCTGAAGTACCTGAACTGATCCGCACCAg  >  3:44163/1‑149 (MQ=255)                                                                                                                                                gCTGCCGGAACTGGAACAGCATCGCGCGCAGCGTATTCTg                                                                                                          >  2:142873/1‑40 (MQ=255)                                                                                                                                                gCTGCCGGAACTGGAACAGCATCGCGCGCAGCGTATTCTg                                                                                                          <  1:142873/40‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                            TGGCAGGCCGGACGTAAAAGTGACGCCCAGCGCGTGTTGCTGGACGCATTAAAACTGGCGAATCGCACCGGATTTATCAGCCATTTTGTCATCGAAGGCGAAGCGATGGCGCAACAACTGCGTCAGCTGATTCAGCTTAATACGCTGCCGGAACTGGAACAGCATCGCGCGCAGCGTATTCTGCGAGAAATCAATCAACATCATCGGCATAAATTCGCCCATTTCGATGAGAATTTCGTTGAACGTCTGCTAAATCATCCTGAAGTACCTGAACTGATCCGCACCAG  >  NC_000913/3555299‑3555585

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