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

GCCTTGCAGGCCCGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGCTCACCCACCTGCACCGCTGCCAGCTCCAACCTGCCCAGCGCCACATCCACGTAGTTCAGCT  >  AE015451/5681181‑5681625                                                                                                                                                                                                                                           |                                                                                                                                                                                                                   gCCTTGCAGGCCCGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcg                                                                                                                                                                                                                  >  1:297177‑M1/1‑235 (MQ=255)             cgcgGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGGCACTGCATGAACCGCcggggctgcggcg                                                                                                                                                                                                       >  1:140291‑M1/1‑223 (MQ=255)                     cGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctat                                                                                                                                                                                              <  1:219102‑M1/237‑23 (MQ=255)                                     aCCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcg                                                                                                                                                                                            <  2:276505‑M1/223‑25 (MQ=255)                                                     aCGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacccaaggccgctcccactccagaccac                                                                                                                                                              >  1:441655‑M1/1‑183 (MQ=255)                                                                  ggTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgc                                                                                                                                                  <  2:297175‑M1/236‑67 (MQ=255)                                                                           aGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggtgctgctgcgcagccctatcgcgacacaaggccgctcccac                                                                                                                                                                        >  1:139336‑M1/1‑161 (MQ=255)                                                                                   aaGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcc                                                                                                                                                 >  1:113630‑M1/1‑153 (MQ=255)                                                                                    aGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaagg                                                                                                                                <  1:272192‑M1/236‑85 (MQ=255)                                                                                                                               gCTGGGTTCCTGTTGTGATGCTGGTGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATGGCCACTGCATGAACCGCcggggctgctgcgcagcccta                                                                                                                                                                                               <  2:82933‑M1/130‑22 (MQ=255)                                                                                                                                  gggTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCATTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaa                                                                                                                                                                                    <  1:134616‑M1/138‑33 (MQ=255)                                                                                                                                           ttgtGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttcc                                                                                            >  1:98904‑M1/1‑97 (MQ=255)                                                                                                                                                aTGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagcccaccaggccctgcgcctcacctccta                                                                                                                                       >  2:264709‑M1/1‑92 (MQ=255)                                                                                                                                                aTGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcacc                                                                                                     >  2:358027‑M1/1‑92 (MQ=255)                                                                                                                                                  gCTGATGGCTTCGTTCTTTGTGGGAGCGGCATTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCCTTCATTGCCACTGCATGAACCGCcggggctgc                                                                                                                                                                                                           >  2:93949‑M1/1‑90 (MQ=255)                                                                                                                                                           ttCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGTAGCGGCCCCAGGCCTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggctgcagatccgcgtcacgcaccagcctttccaatgc                                                                                       >  1:73901‑M1/1‑81 (MQ=255)                                                                                                                                                                   ttGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACcgc                                                                                                                                                                                                                    <  2:170562/73‑1 (MQ=255)                                                                                                                                                                   ttGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACcgc                                                                                                                                                                                                                    >  1:170563/1‑73 (MQ=255)                                                                                                                                                                               ccTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgct                                                                                                                                                                                                          <  1:20473‑M1/71‑11 (MQ=255)                                                                                                                                                                               ccTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgct                                                                                                                                                                                                          >  2:20473‑M1/1‑61 (MQ=255)                                                                                                                                                                                 ttGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCAATGCATGAACCGCcggggatgct                                                                                                                                                                                                          <  2:98885‑M1/69‑11 (MQ=255)                                                                                                                                                                                            aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggtacccgcgccacctgctcacccccccgcacccctgccagctccaacctgccc                       >  1:139011‑M1/1‑48 (MQ=255)                                                                                                                                                                                            aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggtacccgcgccacctgctcacccacctgcaccgctgccagctccaacctgccc                       >  1:432288‑M1/1‑48 (MQ=255)                                                                                                                                                                                                     gAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggtacccg                                                                     <  2:176499‑M1/182‑144 (MQ=255)                                                                                                                                                                                                        gCGGCCCCAGCCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccact                                                                                                                                                                       <  1:36845‑M1/81‑46 (MQ=255)                                                                                                                                                                                                                             cACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggtacccgcgccacctgctcacccacctgcaccgctgccagctccaacctgcccagcgccacatccacgtagttc  <  1:358029‑M1/225‑211 (MQ=255)                                                                                                                                                                                                                                   aTGAACCGCcggggctgctgcgcagccctatcgcg                                                                                                                                                                                          <  1:175206‑M1/35‑27 (MQ=255)                                                                                                                                                                                                                                   aTGAACCGCcggggctgctgcgcagccctatcgcg                                                                                                                                                                                          >  2:175205‑M1/1‑9 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                   GCCTTGCAGGCCCGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGCTCACCCACCTGCACCGCTGCCAGCTCCAACCTGCCCAGCGCCACATCCACGTAGTTCAGCT  >  AE015451/5681181‑5681625

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.