New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	3663643 =	22 (0.540)	4 (0.110)	4/304	NT	16.4%	coding (267/891 nt)	PP_3227	Transcriptional regulator, LysR family
	?	AE015451	3663659 =	21 (0.570)					coding (251/891 nt)	PP_3227	Transcriptional regulator, LysR family

CCTTCGGGGAAAATCTTGCCAGGCGGGAACAGGGCCGCGGCAGCGTCCAGCTCACCGTTTTCCATGCGCGCGATCAGGCTGCTGCCCCAGCCGCTGGTGACCTGGGTGCGCAGTTGGGGGTAGGCCTCACGGATCTGTGCCAGGGCGTCGAGCAGCACAACGTCGCCCAATGTCTGCGGCACACCCAGGCGCAAGGTGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCTTGCACTGCTCATAGACCCGCAAGCCCATGGGGGTGGGCTTGAGCGGCTTGGTGTTGCGGTCGAGGAGGGTGGCGCCGAGGTCTTCCTCGAAATTCTGCACCCGGCGGGTAATGGCCGGCTGCGTCAGTTGCAGGGCCTCGGCGGCAAGGTTGGTGGACTGGCAACGAATGACGGCGA  >  AE015451/3663443‑3663881                                                                                                                                                                                                                         |                                                                                                                                                                                                                               ccTTCGGGGAAAATCTTGCCAGGCGGGAACAGGGCCGCGGCAGCGTCCAGCTCACCGTTTTCCATGCGCGCGATCAGGCTGCTGCCCCAGCCGCTGGTGACCTGGGTGCGCAGTTGGGGGTAGGCCTCACGGATCTGTGCCAGGGCGTCGAGCAGCACAACGTCGCCCAATGTCTGCGGCACACCCAGGCGCAAGGTGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACg                                                                                                                                                                                                                <  2:534274/233‑1 (MQ=255)        ggAAAATCTTGCCAGGCGGGAACAGGGCCGCGGCAGCGTCCAGCTCACCGTTTTCCATGCGCGCGATCAGGCTGCTGCCCCAGCCGCTGGTGACCTGGGTGCGCAGTTGGGGGTAGGCCTCACGGATCTGTGCCAGGGCGTCGAGCAGCACAACGTCGCCCAATGTCTGCGGCACACCCAGGCGCAAGGTGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCt                                                                                                                                                                                                      >  1:275173/1‑236 (MQ=255)                                                                        gATCAGGCTGCGGCCCCAGCCGCTGGTGGCCTGGGTGCGCAGTTGGGGGTAGGCCTCACGGATCTGTGCCAGGGCGTCGAGCAGCACAACGTCGCCCAATGTCTGCGGCACACCCAGGCGCAAGGTGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACg                                                                                                                                                                                                                <  2:238118/162‑1 (MQ=255)                                                                          tCAGGCTGCTGCCCCAGCCGCTGGTGACCTGGGTGCGCAGTTGGGGGTAGGCCTCACGGATCTGTGCCAGGGCGTCGAGCAGCACAACGTCGCCCAATGTCTGCGGCACACCCAGGCGCAAGGTGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACg                                                                                                                                                                                                                <  2:585928/160‑1 (MQ=255)                                                                             tcggcgccaccctcctcgaccgcaacaccaagccgcccaagcccccccccatgggcttgcgggtctatgagcagtgcaaggcgatcctgcgcgagatagactcgctgcgtgaactggtagccaaCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAg                                                                                                                                                                                                             >  2:382454‑M1/125‑160 (MQ=255)                                                                                tgctgcCCCAGCCGCTGGTGACCTGGGTGCGCAGTTGGGGGTAGGCCTCACGGATCTGTGCCAGGGCGTCGAGCAGCACAACGTCGCCCAATGTCTGCGGCACACCCAGGCGCAAGGTGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCTTGCACTGCTCATAGACCCGCAAGCCCATGGGGGTGGGCTTGAGCGGCTTGgt                                                                                                                               >  1:168469/1‑235 (MQ=255)                                                                                         aGCCGCTGGTGACCTGGGTGCGCAGTTGGGGGTAGGCCTCACGGATCTGTGCCAGGGCGTCGAGCAGCACAACGTCGCCCAATGTCTGCGGCACACCCAGGCGCAAGGTGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCa                                                                                                                                                                                            <  2:305163/165‑1 (MQ=255)                                                                                                      cTGGGTGCGCAGTTGGGGGTAGGCCTCACGGATCTGTGCCAGGGCGTCGAGCAGCACAACGTCGCCCAATGTCTGCGGCACACCCAGGCGCAAGGTGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCtt                                                                                                                                                                                  <  2:153418/162‑1 (MQ=255)                                                                                                            gcgcAGTTGGGGGTAGGCCTCACGGATCTGTGCCAGGGCGTCGAGCAGCACAACGTCGCCCAATGTCTGCGGCACACCCAGGCGCAAGGTGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCTTGCACTGCTCATAGACCCGCAAGCCCATGGGGGTGGGCTTGAGCGGCTTGGTGTTGCGGTCGAGGAGGGTGGCGCCGAg                                                                                                    >  1:386499/1‑234 (MQ=255)                                                                                                                     gcccacccccatgggcttgcgggtctatgagcagtgcaaggcgagcctgcgcgagatagactcgctgcgtgaactggtagccaacgacggcgcgccgtcgTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCTTg                                                                                                                                                                                 <  2:738717‑M2/48‑1 (MQ=255)                                                                                                                               tCACGGATCTGTGCCAGGGCGTCGAGCAGCACAACGTCGCCCAATGTCTGCGGCACACCCAGGCGCAAGGTGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCa                                                                                                                                                                                            <  2:647571/127‑1 (MQ=255)                                                                                                                               tCACGGATCTGTGCCAGGGCGTCGAGCAGCACAACGTCGCCCAATGTCTGCGGCACACCCAGGCGCAAGGTGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCa                                                                                                                                                                                            >  1:647572/1‑127 (MQ=255)                                                                                                                               tCACGGATCTGTGCCAGGGCGTCGAGCAGCACAACGTCGCCCAATGTCTGCGGCACACCCAGGCGCAAGGTGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCTTGCACTGCTCATAGACCCGCAAGCCCATgggg                                                                                                                                                   >  2:501709/1‑168 (MQ=255)                                                                                                                                          cgtctatgagccgtgcaagtcgatcctgcgcgagatagactcgccgcgtgaactggcagccaaAGACGGCGCGCCGTCGATGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCa                                                                                                                                                                                            >  1:680711‑M1/64‑116 (MQ=255)                                                                                                                                             cAGGGCGTCGAGCAGCACAACGTCGCCCAATGTCTGCGGCACACCCAGGCGCAAGGTGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCTTGCACTGCTCATAGACCCGCAAGCCCATgggg                                                                                                                                                   >  1:130907/1‑154 (MQ=255)                                                                                                                                             cAGGGCGTCGAGCAGCACAACGTCGCCCAATGTCTGCGGCACACCCAGGCGCAAGGTGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCTTGCACTGCTCATAGACCCGCAAGCCCATGGGGGTGGGCTTGAGCGGCTTGGTGTTGCGGTCGAGGAGGGTGGCGCCGAGGTCTTCCTCGAAATTCTGCACCCGGCGGGTAAt                                                                   >  2:402757/1‑234 (MQ=255)                                                                                                                                             cAGGGCGTCGAGCAGCACAACGTCGCCCAATGTCTGCGGCACACCCAGGCGCAAGGTGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCTTGCACTGCTCATAGACCCGCAAGCCCATGGGGGCGGGCTTGAGCGGCTTGGTGTTGCGGTCGAGGAGGGGGGCGCCGAGGTCTTCCTCg                                                                                          >  1:69305/1‑211 (MQ=255)                                                                                                                                                       tgcaaggcgatcctgcgcgagatagcctcgctgcgtgaactggtagccaacgacggcgcgccgtcgTTGGCTACCAGTTCACGCAGCGAGTCTAtct                                                                                                                                                                                                  <  2:529666‑M2/31‑1 (MQ=255)                                                                                                                                                       agcagcACAACGTCGCCCAATGTCTGCGGCACACCCAGGCGCAAGGTGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCTTGCACTGCTCATAGACCCGCAAGCCCACGGGGGTGGGCTTGAGCGGCTTGGTGTTGCGGTCGAGGAGGGTGGCGCCGAGGTCTTc                                                                                              >  2:446408/1‑197 (MQ=255)                                                                                                                                                         caaggcgatcctgcgcgagatagactcgctgcgtgaactggtagccaaCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGCTCGCCTTGCACTGCTCATAGACCCGCAAGCCCATGGGGGTGGGc                                                                                                                                             <  2:738717‑M1/100‑1 (MQ=255)                                                                                                                                                          ggcaCAACGTGGCGCAATGTCTGCGGCACACCCAGGCGCAAGGGGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCTTGCACTGCTCATAGACCCg                                                                                                                                                                <  2:23397/127‑1 (MQ=255)                                                                                                                                                               accgccGCCCAATGTCTGCGGCACACCCAGGCGCAAGGTGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCTTGCACTGCTCATAGACCCGCAAGCCCATGGGGGTGGGCtt                                                                                                                                           <  2:263192/139‑1 (MQ=255)                                                                                                                                                                    tgcgcgagatagactcgctgcgtgaactggtagccatcgacggcgcgccgtctTTGGCTGCCAGTTCACGCGGCGAGTCTATCTCGCGCAGGATCGACTTGCACGGCTCATAGacg                                                                                                                                                                  >  1:680711‑M2/54‑115 (MQ=255)                                                                                                                                                                    cGCCCAATGTCTGCGGCACACCCAGGCGCAAGGTGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCTTGCACTGCTCATAGACCCGCAAGCCCATGGGGGTGGGCtt                                                                                                                                           >  2:393913/1‑139 (MQ=255)                                                                                                                                                                         aaTGTCTGCGGCACACCCAGGCGCAAGGTGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCTTGCACTGCTCCTAGACCCGCAAGCCCATg                                                                                                                                                      >  2:52504/1‑123 (MQ=255)                                                                                                                                                                          agatagactcgctgcgtgaactggtagccaaCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGGCTATCTCGCGCAGGATCGCCTTGCa                                                                                                                                                                               <  2:529666‑M1/66‑1 (MQ=255)                                                                                                                                                                                     ctgcgtgaactggtagccaacgacggcgcgccgtcgTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCTTGCACTGCTCATAGACCCGCAAGCCCATGGGGGGGGGCTTGGGCGGCTTGGTGTTGCGGTCGAGGAGGGTGGCGCCGa                                                                                                     >  2:382454‑M2/37‑160 (MQ=255)                                                                                                                                                                                               gcAAGGTGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCTTGCACTGCTCATAGACCCGCAAGCCCATGGGGGTGGGCTTGAGCGGCCTGGTGTTGCGGTCGAGGAGGGTGGCGCCGAGGTCTTCCTCGAAATTCTGCACCCGGCGGGTAATGGCCGGCTg                                                          >  1:511880/1‑193 (MQ=255)                                                                                                                                                                                                   ggTGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCTTGCACTGCTCATAGACCCGCAAGCCCATGGGGGTGGGCTTGAGCGGCTTGGTGTTGCGGTCGAGGAGGGTGGCGCCGAGGTCTTCCTCGAAATTCTGCACCCGGCGGGTAATGGCCGGCTGCGTCAGTTGCAGGGCCTCGGCGGCAAGGTTGGTGGACTGGCAACGa            >  1:746818/1‑235 (MQ=255)                                                                                                                                                                                                      gCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCTTGCACTGCTCATAGACCCGCAAGCCCATGGGGGTGGGCTTGAGCGGCTTGGTGTTGCGGTCGAGGAGGGTGGCGCCGAGGTCTTCCTCGAAATTCTGCACCCGGCGGGTAATGGCCGGCTGCgg                                                       >  1:108406/1‑188 (MQ=255)                                                                                                                                                                                                              gcgcgccgtcgTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCTTGCACTGCTCATAGACCCGCAAGCCCATGGGGGTGGGCTTGAGCGGCTTGGTGTTGCGGTCGAGGAGGGTGGCGCCGAGGTCTTCCTCGAAATTCTGCACCCGGCGGGTAATGGCCGGCTGCGTCAGTTGCAGGGCCTCGGCGGCAAGGTTGGTGGACTGGCAACGAATGACGGcga  >  2:544090‑M2/12‑234 (MQ=255)                                                                                                                                                                                                                gcgcCGTCGATTGCTACCAGATCACGCAGCGAGTCTATCTCGCGCa                                                                                                                                                                                            <  1:279876/46‑1 (MQ=255)                                                                                                                                                                                                                  gccgtcgTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCTTg                                                                                                                                                                                 >  2:672055‑M2/8‑55 (MQ=255)                                                                                                                                                                                                                  GCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCtt                                                                                                                                                                                  <  2:232971‑M1/54‑1 (MQ=255)                                                                                                                                                                                                                  GCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCtt                                                                                                                                                                                  >  1:232972‑M1/1‑54 (MQ=255)                                                                                                                                                                                                                  GCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCTTg                                                                                                                                                                                 >  1:672056‑M1/1‑55 (MQ=255)                                                                                                                                                                                                                   CCGGCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCTTGCACTGCTCATAGACCCGCAAGCCCATGGGGGTGGGCTTGAGCGGCTTGGTGTTGCGGTCGAGGAGGGTGGCGCCGAGGTCTTCCTCGAAATTCTGCACCCGGCGGGTAATGGCCGGCTGCGTCAGTTGCAGGGCCTCGGCGGCAAGGTtggt                         >  2:168468‑M1/1‑206 (MQ=255)                                                                                                                                                                                                                         |                                                                                                                                                                                                                               CCTTCGGGGAAAATCTTGCCAGGCGGGAACAGGGCCGCGGCAGCGTCCAGCTCACCGTTTTCCATGCGCGCGATCAGGCTGCTGCCCCAGCCGCTGGTGACCTGGGTGCGCAGTTGGGGGTAGGCCTCACGGATCTGTGCCAGGGCGTCGAGCAGCACAACGTCGCCCAATGTCTGCGGCACACCCAGGCGCAAGGTGCCCGACGGCGCGCCGTCGTTGGCTACCAGTTCACGCAGCGAGTCTATCTCGCGCAGGATCGCCTTGCACTGCTCATAGACCCGCAAGCCCATGGGGGTGGGCTTGAGCGGCTTGGTGTTGCGGTCGAGGAGGGTGGCGCCGAGGTCTTCCTCGAAATTCTGCACCCGGCGGGTAATGGCCGGCTGCGTCAGTTGCAGGGCCTCGGCGGCAAGGTTGGTGGACTGGCAACGAATGACGGCGA  >  AE015451/3663443‑3663881

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