Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	AE015451	5,681,415	+CGGG	intergenic (‑133/+71)	PP_4986 ← / ← PP_4987	Channel protein, hemolysin III family/putative Chemotaxis protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 5681415	0 (0.000)	25 (1.140) +CGGG	22/308	0.1	100%	intergenic (‑133/+71)	PP_4986/PP_4987	Channel protein, hemolysin III family/putative Chemotaxis protein
	?	AE015451	5681416 =	0 (0.000)					intergenic (‑134/+70)	PP_4986/PP_4987	Channel protein, hemolysin III family/putative Chemotaxis protein

GGCCCGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGCTCACCCACCTGCACCGCTGCCAGCTCCAACCTGCCCAGCGCCACATCCACGTAGTTCAGCTGGCTGTCGAGCTTGCACGAGCGCG  >  AE015451/5681189‑5681649                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                          cgcggccacgatcagccagatggcaccaatagcagcgaggacggcaccgaccaggtgggtccaggcgttgaagcgttcaccgtaatacatgcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattggcactgcatgaaccgccgggGCTgcggcg                                                                                                                                                                                                                                   >  1:140291‑M2/228‑236 (MQ=255)         cgatcagccagatggcaccaatagcagcgaggacggcaccgaccaggtgggtccaggcgttgaagcgttcaccgtaatacatgcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTAt                                                                                                                                                                                                                          <  1:219102‑M2/18‑1 (MQ=255)                         accaatagcagcgaggacggcaccgaccaggtgggtccaggcgttgaagcgttcaccgtaatacatgcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATcg                                                                                                                                                                                                                        <  2:276505‑M2/20‑1 (MQ=255)                                         acggcaccgaccaggtgggtccaggcgttgaagcgttcaccgtaatacatgcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACCCAAGGCCGCTCCCACTCCAGaccac                                                                                                                                                                                          >  1:441655‑M2/188‑237 (MQ=255)                                                      ggtgggtccaggcgttgaagcgttcaccgtaatacatgcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTgcgc                                                                                                                                                                              <  2:297175‑M2/62‑1 (MQ=255)                                                               aggcgttgaagcgttcaccgtaatacatgcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccggtGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCAc                                                                                                                                                                                                    >  1:139336‑M2/166‑205 (MQ=255)                                                                       aagcgttcaccgtaatacatgcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGcc                                                                                                                                                                             >  1:113630‑M2/158‑220 (MQ=255)                                                                        agcgttcaccgtaatacatgcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAgg                                                                                                                                                            <  1:272192‑M2/80‑1 (MQ=255)                                                                                                                   gctgggttcctgttgtgatgctggtggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcatggccactgcatgaaccgccgggGCTGCTGCGCAGCCCTa                                                                                                                                                                                                                           <  2:82933‑M2/17‑1 (MQ=255)                                                                                                                      gggttcctgttgtgatgctgatggcttcgttctttgtgggagcggcattgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACaa                                                                                                                                                                                                                <  1:134616‑M2/28‑1 (MQ=255)                                                                                                                               ttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTcc                                                                                                                        >  1:98904‑M2/102‑217 (MQ=255)                                                                                                                                    atgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGCCCACCAGGCCCTGCGCCTCACCTCCTa                                                                                                                                                                   >  2:264709‑M2/97‑169 (MQ=255)                                                                                                                                    atgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCAcc                                                                                                                                 >  2:358027‑M2/97‑203 (MQ=255)                                                                                                                                      gctgatggcttcgttctttgtgggagcggcattgtgtcgcgaaaggggcgcgaagcggccccaggccttcattgccactgcatgaaccgccggggctgc                                                                                                                                                                                                                                       >  2:93949‑M2/95‑99 (MQ=255)                                                                                                                                               ttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgtagcggccccaggccttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGCTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGc                                                                                                                   >  1:73901‑M2/86‑206 (MQ=255)                                                                                                                                                                   ccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccggggctgct                                                                                                                                                                                                                                      >  2:20473‑M2/66‑71 (MQ=255)                                                                                                                                                                   ccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccggggctgct                                                                                                                                                                                                                                      <  1:20473‑M2/6‑1 (MQ=255)                                                                                                                                                                     ttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccaatgcatgaaccgccggggatgct                                                                                                                                                                                                                                      <  2:98885‑M2/6‑1 (MQ=255)                                                                                                                                                                                aaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGCTCACCCCCCCGCACCCCTGCCAGCTCCAACCTGccc                                                   >  1:139011‑M2/53‑237 (MQ=255)                                                                                                                                                                                aaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGCTCACCCACCTGCACCGCTGCCAGCTCCAACCTGccc                                                   >  1:432288‑M2/53‑237 (MQ=255)                                                                                                                                                                                         gaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCcg                                                                                                 <  2:176499‑M2/139‑1 (MQ=255)                                                                                                                                                                                            gcggccccagcctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACt                                                                                                                                                                                                   <  1:36845‑M2/41‑1 (MQ=255)                                                                                                                                                                                                                 cactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGCTCACCCACCTGCACCGCTGCCAGCTCCAACCTGCCCAGCGCCACATCCACGTAGTTc                              <  1:358029‑M2/206‑1 (MQ=255)                                                                                                                                                                                                                       atgaaccgccgggGCTGCTGCGCAGCCCTATcgcg                                                                                                                                                                                                                      >  2:175205‑M2/14‑35 (MQ=255)                                                                                                                                                                                                                       atgaaccgccgggGCTGCTGCGCAGCCCTATcgcg                                                                                                                                                                                                                      <  1:175206‑M2/22‑1 (MQ=255)                                                                                                                                                                                                                                 gggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGCTCACCCACCTGCACCGCTGCCAGCTCCAACCTGCCCAGCGCCACATCCACGTAGTTCAGCTGGCTGTCGAGCTTGCACGAGcgcg  >  1:343547‑M2/4‑237 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                          GGCCCGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGCTCACCCACCTGCACCGCTGCCAGCTCCAACCTGCCCAGCGCCACATCCACGTAGTTCAGCTGGCTGTCGAGCTTGCACGAGCGCG  >  AE015451/5681189‑5681649

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