New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 5201567	9 (0.400)	15 (0.680)	15/310	0.5	59.0%	intergenic (‑53/‑127)	PP_4582/PP_4583	conserved membrane protein of unknown function/putative Peptidase
	?	AE015451	5201666 =	12 (0.550)					intergenic (‑152/‑28)	PP_4582/PP_4583	conserved membrane protein of unknown function/putative Peptidase

GACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGAC  >  AE015451/5201356‑5201744                                                                                                                                                                                                                    |                                                                                                                                                                                  gacgaTGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCacaggcaactccaatcaa                                                                                                                                                                 <  1:271534‑M1/230‑19 (MQ=255)                          cccGCCAGCACCACGCGTAACGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGgcgc                                                                                                                                                                                 <  1:126403/188‑1 (MQ=255)                                                                                                       cccTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTAccc                                                                                               <  1:268906/194‑1 (MQ=255)                                                                                                        ccTAACGCGGCTACTGCGTACATGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAGACACTTTTCCGGGCacaggcaactccaatcaagagaccttcaatgcca                                                                                                                                                 <  2:250314‑M1/143‑35 (MQ=255)                                                                                                         cTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTGCTCGGTAATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCGCTTCGCGGGTAAACCCGCTCCTACCTCGGt                                                                                                                          <  2:209347/165‑1 (MQ=255)                                                                                                                     cTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCacaggcaactccaatcaagagaccttcaatgccaaccaagccccaacccccgattgcccaagccgaccaggccaccgac                                                                                                    >  2:403268‑M1/1‑96 (MQ=255)                                                                                                                     cTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCgggca                                                                                                                                             >  2:42377/1‑132 (MQ=255)                                                                                                                     cTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGtt                                                                                                              >  1:212879/1‑165 (MQ=255)                                                                                                                     cTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGtt                                                                                                              <  2:212878/165‑1 (MQ=255)                                                                                                                     cTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTgg                                                                                                            >  2:268904/1‑167 (MQ=255)                                                                                                                           aCAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCacaggcaactccaatcaagagaccttcaatgccaaccaagccccaacccccgattgcccaagccgaccaggccaccgacccgtatgcctggctgcaacagcgcgacacccccgaggtactcgcttacctgcaagctgaaaacgccta                                <  1:403270‑M1/237‑148 (MQ=255)                                                                                                                              aGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCacaggcaactccaatcaagagaccttcaatgccaaccaagccccaacccccgattgcccaagccgaccaggccaccgacccgtatgcctggctgcaacagcgcgacacccccgaggtactcgcttacctgcaagctga                                         >  2:185235‑M1/1‑87 (MQ=255)                                                                                                                                   cAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCacaggcaactccaatcaagagaccttcaatgccaaccaagccccaacccccgattgcccaagccgaccaggccaccgacccgtatgcctggctgcaacagcgcgacacccccgaggtactcgcttacctgcaagctgaaaacgcctaccaggaa                         >  1:243091‑M1/1‑82 (MQ=255)                                                                                                                                                        tcGCTCATTGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATCAACACCTTTCCGGGCacaggcaactccaatcaagagaccttcaatgccaaccaagccccaacccccgattgcccaagccgaccaggcc                                                                                                          <  2:436420‑M1/134‑74 (MQ=255)                                                                                                                                                         cGGTCATGGTGTTCGTGCAATCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCacaggcaactccaatcaagagaccttcaatgccaaccaagccccaacccccgattgcccaagccgaccaggccaccgacccgtttgcctggctgcaacagcgcgacacccccgaggtactcgcttacctgcaagct                                           <  1:174494‑M1/196‑137 (MQ=255)                                                                                                                                                                         gCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCacaggcaactccaatcaagagaccttcaatgccaaccaagccccaacccccgattgcccaagccgaccaggccacc                                                                                                       <  2:203181‑M1/120‑77 (MQ=255)                                                                                                                                                                         gCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCacaggcaactccaatcaagagaccttcaatgccaaccaagccccaacccccgattgcccaagccgaccaggccacc                                                                                                       >  1:203182‑M1/1‑44 (MQ=255)                                                                                                                                                                               cAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCctc                                                                                                                                                       <  1:293539/66‑1 (MQ=255)                                                                                                                                                                               cAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCctc                                                                                                                                                       >  2:293537/1‑66 (MQ=255)                                                                                                                                                                                              cTATGATAAACACTTTTCCGGGCacaggca                                                                                                                                                                            <  2:140961‑M1/30‑8 (MQ=255)                                                                                                                                                                                               tATGATAAACACTTTTCCGGGCacaggcaactccaatcaagagaccttcaatgcaaaccaagccccaaaccccgattgccaaaaccgacc                                                                                                               >  2:406082‑M1/1‑22 (MQ=255)                                                                                                                                                                                                aTGATAAGCGCTTTGCCGGGCgcaggcaactccaatcaagagcccttcaatgccaaccaagccccaacccccgattgcccaagccgaccaggccaccgacccgtatgcctggctgcaacagcg                                                                             <  2:56363‑M1/123‑103 (MQ=255)                                                                                                                                                                                                              tCCGGGCacaggcaactccaatcaagagaccttcaatgccaaccaagccccagcccccgattgcccaagccgaccaggccaccgacccgtatgcctggctgcaacagcgcgacacccccgaggtactcgcttacctgcaagctgaaaacg                                    <  1:42377‑M1/150‑144 (MQ=255)                                                                                                                                                                                                               ccGGGCacaggcaactccaaccaagagaccttcaatgccaaccaagccccaacccccgattgcccaagccgaccaggcccccgacccgtatgcctggctgcaacagcgcgacgcccccgaggtactcgcttacctgcaagctgaaaacgcctaccaggaagcctgcctggccgaccaggcacc  <  1:137748‑M1/183‑178 (MQ=255)                                                                                                                                                                                                                    |                                                                                                                                                                                  GACGATGGTCGGGGTCAGCAGCAGCACCGCCAGCACCACGCGTAGCGGTTCGCGCAGCCAGCGCCACATCCAGCCGGTGAGCTTGAAGCCCACCAGCAGGCACCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGAC  >  AE015451/5201356‑5201744

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