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

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 5681415	0 (0.000)	21 (0.380) +CGGG	21/250	NT	100%	intergenic (‑133/+71)	PP_4986/PP_4987	hemolysin III family channel protein/putative Chemotaxis protein
	?	NC_002947	5681416 =	0 (0.000)					intergenic (‑134/+70)	PP_4986/PP_4987	hemolysin III family channel protein/putative Chemotaxis protein

ACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCA  >  NC_002947/5681271‑5681559                                                                                                                                                 |                                                                                                                                                 aCCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcgg                                                                                                                                               <  6:133175‑M1/148‑4 (MQ=255)             gCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcag                                                                                                                                  >  7:258603‑M1/1‑133 (MQ=255)               aaaCACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCCGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagcc                                                                                                                                >  2:17328‑M1/1‑131 (MQ=255)               aaaCACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagcc                                                                                                                                <  5:321003‑M1/149‑19 (MQ=255)               aaaCACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcgggg                                                                                                                                             <  6:47336‑M1/136‑6 (MQ=255)               aaaCACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcgggg                                                                                                                                             >  5:47336‑M1/1‑131 (MQ=255)                      gACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgc                                                                                                                         <  8:258603‑M1/149‑26 (MQ=255)                         ctcgtGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgac                                                                                                                      <  5:121605‑M1/149‑29 (MQ=255)                                    ggCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCCCTGCATGAACCGCcggggctgctgcgcagccctatcgcgccacaaggcggcg                                                                                                           >  8:169675‑M1/1‑110 (MQ=255)                                               tgttgtGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccaga                                                                                                >  7:81687‑M1/1‑99 (MQ=255)                                                      aTGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccag                                                                                         >  5:362832‑M1/1‑92 (MQ=255)                                                         cTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgccgccctatcgcggcacaaggccgctcccactccagaccaccaggcc                                                                                      >  1:95340‑M1/1‑89 (MQ=255)                                                               gCTTCGTGCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTCCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgc                                                                                <  7:169675‑M1/149‑67 (MQ=255)                                                                               gAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgccgcccaatcgcgccacaaggccgctcccactccagaccaccccgccctgccccgcacatcctaacc                                                                  >  8:108164‑M1/1‑67 (MQ=255)                                                                                            tCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttg                                                      >  1:32168‑M1/1‑54 (MQ=255)                                                                                              gctaaagggGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATTAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagat                                                 <  1:101639‑M1/146‑98 (MQ=255)                                                                                                      ggCGCGTAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACGGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtca                                         <  7:108164‑M1/149‑106 (MQ=255)                                                                                                        cgcgAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacg                                       <  6:362832‑M1/149‑108 (MQ=255)                                                                                                             accggCCCCAGGCTTTAATTGCCCCTCCAAGAACCCCcggggctgctgcgcagccctatcccgacaaaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcacca                                  <  7:56057‑M1/147‑113 (MQ=255)                                                                                                               cGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagc                                <  8:81687‑M1/149‑115 (MQ=255)                                                                                                                     cAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttcc                          <  2:32168‑M1/149‑121 (MQ=255)                                                                                                                       ggattccattgcaACTGCATGACCCGCcggggctgctgcgcagccctatcgccacacacggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaa                        <  3:253055‑M1/143‑123 (MQ=255)                                                                                                                              aTTGCCACTGCATGAACCGCcggggctgctgcgcagccctagcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgacc                 <  2:95340‑M1/149‑130 (MQ=255)                                                                                                                                        cATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcttcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccgg         <  1:17328‑M1/147‑138 (MQ=255)                                                                                                                                             aCCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggtacccgc  >  2:71474‑M1/1‑5 (MQ=255)                                                                                                                                                 ccggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgc                                             <  6:80865‑M1/102‑102 (MQ=255)                                                                                                                                                 ccggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgc                                             >  5:80865‑M1/1‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                 ACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCA  >  NC_002947/5681271‑5681559

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