Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	AE015451	5,681,415	+CGGG	100%	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)	18 (0.720) +CGGG	16/310	NT	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

CGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGCTCACCCACCTGCACCGCTGCCAGCTCCAACCTGCCCAGCGCCACATCCACGTAGTTCAGCTGGCTGTCGAGCTTGCACGAGCGC  >  AE015451/5681193‑5681648                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         gccacgatcagccagatggcaccaatagcagcgaggacggcaccgaccaggtgggtccaggcgttgaagcgttcaccgtaatacatgcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAg                                                                                                                                                                                                                               <  1:164921‑M2/12‑1 (MQ=255)                               gcgaggacggcaccgaccaggtgggtccaggcgttgaagcgttcaccgtaatacatgcaaacacagacctccggaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGc                                                                                                                                                                                                         <  2:219075‑M2/34‑1 (MQ=255)                                                      ggtccaggcgttgaagcgttcaccgtaatacatgcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAg                                                                                                                                                                                                              <  1:173759‑M2/29‑1 (MQ=255)                                                      ggtccaggcgttgaagcgttcaccgtaatacatgcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAg                                                                                                                                                                                                              >  2:173755‑M2/171‑199 (MQ=255)                                                                                                                          gttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTc                                                                                                                                                                                                       <  2:415535‑M2/36‑1 (MQ=255)                                                                                                                                atgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTgcgc                                                                                                                                                                             >  2:238303‑M2/97‑158 (MQ=255)                                                                                                                                  gctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCa                                                                                                                                      >  1:193444‑M2/95‑195 (MQ=255)                                                                                                                                   ctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGt                                                                                                                                        >  1:195148‑M2/94‑192 (MQ=255)                                                                                                                                      atggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGccc                                                                                                                                                                                  <  2:141895‑M2/57‑1 (MQ=255)                                                                                                                                                            cggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGc                                                                                      <  1:432289‑M2/149‑1 (MQ=255)                                                                                                                                                                    tgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCa                                                                                                                                      >  2:272805‑M2/61‑161 (MQ=255)                                                                                                                                                                                       agcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGCTCACCCACCTGCACCGCTGCCAGCTCCAACCTGCCCAg                                                <  2:281223‑M2/187‑1 (MQ=255)                                                                                                                                                                                       agcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTACCGCGACACAACGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACCTCCTCACCCAGCTCACGTTGCAg                                                                                                                                                >  2:2853‑M2/42‑132 (MQ=255)                                                                                                                                                                                                         ttgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATcg                                                                                                                                                                                                                       >  2:505305‑M2/24‑43 (MQ=255)                                                                                                                                                                                                                   atgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGc                                                                                                                                                                                                            >  2:195372‑M2/14‑44 (MQ=255)                                                                                                                                                                                                                   atgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGt                                                                                                                                                     >  2:252739‑M2/14‑99 (MQ=255)                                                                                                                                                                                                                   atgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCcgc                                                                                                                                          >  1:499742‑M2/14‑110 (MQ=255)                                                                                                                                                                                                                   atgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCcgc                                                                                                                                          <  2:499738‑M2/97‑1 (MQ=255)                                                                                                                                                                                                                      acccgccgggGCTGCGGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt                                                                                                                         <  2:20728‑M2/114‑1 (MQ=255)                                                                                                                                                                                                                         cgccgggGCTGCTGCGCAGCCCTATCGCGCCACAAGGCCGCTCCCACt                                                                                                                                                                                                  <  2:74824‑M2/41‑1 (MQ=255)                                                                                                                                                                                                                             gggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGCTCACCCACCTGCACCGCTGCCAGCTCCAACCTGCCCAGCGCCACATCCACGTAGTTCAGCTGGCTGTCGAGCTTGCACGAgcgc  <  1:423268‑M2/233‑1 (MQ=255)                                                                                                                                                                                                                               gGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCCAAcc                                                                                                                                                               >  1:64669‑M2/2‑77 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         CGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGCTCACCCACCTGCACCGCTGCCAGCTCCAACCTGCCCAGCGCCACATCCACGTAGTTCAGCTGGCTGTCGAGCTTGCACGAGCGC  >  AE015451/5681193‑5681648

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