New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1411899 =	35 (0.880)	5 (0.160)	4/216	4.3	14.4%	coding (50/936 nt)	ttcA	tRNA s(2)C32 thioltransferase, iron sulfur cluster protein
	?	NC_000913	= 1434958	31 (0.970)					pseudogene (51/51 nt)	ttcC	pseudogene, prophage Rac integration site ttcA duplication;Phage or Prophage Related
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGT  >  NC_000913/1411765‑1412040                                                                                                                                       |                                                                                                                                              aaTGGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTAt                                                                                                                                 <  3:14582/149‑1 (MQ=255)      cgcAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGttc                                                                                                                            <  4:152778/149‑1 (MQ=255)         aGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGttcttt                                                                                                                         >  6:6823/1‑149 (MQ=255)          gAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGtt                                                                                                                             <  1:47966/144‑1 (MQ=255)          gAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGtt                                                                                                                             >  2:47966/1‑144 (MQ=255)             tctcCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCTATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTctttcttt                                                                                                                     >  3:94699/1‑149 (MQ=255)             tctcCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTctttcttt                                                                                                                     >  5:135316/1‑149 (MQ=255)             tctcCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTctttcttt                                                                                                                     >  3:23914/1‑149 (MQ=255)                                 tCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGttcttt                                                                                                                         <  2:31995/125‑1 (MQ=255)                                 tCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGttcttt                                                                                                                         >  1:31995/1‑125 (MQ=255)                                     tACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGc                                                                                                <  5:6823/146‑1 (MQ=255)                                        ccccGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCattatt                                                                                          >  1:187664/1‑149 (MQ=255)                                           cGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCa                                                                                       <  8:92495/149‑1 (MQ=255)                                                   cAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTg                                                                               <  2:187664/149‑1 (MQ=255)                                                    aaaCCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTgg                                                                              >  5:113722/1‑149 (MQ=255)                                                    aaaCCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTgg                                                                              >  1:91380/1‑149 (MQ=255)                                                      aCCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTa                                                                            <  6:135316/149‑1 (MQ=255)                                                                         cttcAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTgg                                                          <  4:23914/148‑1 (MQ=255)                                                                             aaTCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGc                                                      >  1:60315/1‑148 (MQ=255)                                                                                    ttGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTg                                              <  6:113722/149‑1 (MQ=255)                                                                                         ggcAGCATTGGCTGCGACCGCGTTACGACGCAGACGTTTTTGTAATGTGGTCAGGTTGTATGGGGCTGTCTTTGTAATTTGGTGATTTGCTTGCATTCTTTCAGTTCTCTGGTACTAAATGGGGCAAATTggg                                                         <  4:11053/131‑1 (MQ=255)                                                                                         gTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAAc                                          >  1:9829/1‑148 (MQ=255)                                                                                          tCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCATGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTtctc                                                                                 <  2:190275/108‑1 (MQ=255)                                                                                          tCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTtctc                                                                                 >  1:190275/1‑108 (MQ=255)                                                                                              cAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTgg                                                                              >  6:104571/1‑107 (MQ=255)                                                                                              cAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTgg                                                                              <  5:104571/107‑1 (MQ=255)                                                                                                    cTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACcgc                              <  2:60315/149‑1 (MQ=255)                                                                                                           cACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACttt                                               >  1:155591/1‑125 (MQ=255)                                                                                                           cACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACttt                                               <  2:155591/125‑1 (MQ=255)                                                                                                               ttACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACAt                   >  2:173365/1‑149 (MQ=255)                                                                                                                 acgacgCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATgg                 >  4:3871/1‑149 (MQ=255)                                                                                                                  cgacgCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATgg                 <  2:9829/148‑1 (MQ=255)                                                                                                                                ttGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTcgtcgt  >  6:23620/1‑149 (MQ=255)                                                                                                                                      tttGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGa                                     <  2:132551/108‑1 (MQ=255)                                                                                                                                      tttGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGa                                     >  1:132551/1‑108 (MQ=255)                                                                                                                                       |                                                                                                                                              AATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGT  >  NC_000913/1411765‑1412040

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