Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	NC_002947	5,681,415	+CGGG	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	1 (0.020)	23 (0.490) +CGGG	23/276	0.5	96.0%	intergenic (‑133/+71)	PP_4986/PP_4987	hemolysin III family channel protein/putative Chemotaxis protein
	?	NC_002947	5681416 =	1 (0.020)					intergenic (‑134/+70)	PP_4986/PP_4987	hemolysin III family channel protein/putative Chemotaxis protein

CGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGC‑‑‑‑GCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTA  >  NC_002947/5681273‑5681550                                                                                                                                               |                                                                                                                                            cGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggg                                                                                                                                         <  4:218218‑M1/147‑5 (MQ=255)            cAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGTGAAGCGGCCCCAGGCTTTCATTTCCACTGCATGAACCGCcggggctgctgcgca                                                                                                                              >  5:88061‑M1/1‑132 (MQ=255)            cAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgca                                                                                                                              <  8:43070‑M1/147‑16 (MQ=255)             aaaCACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCccgggctgctgggcag                                                                                                                             >  5:105567‑M1/1‑131 (MQ=255)                      cctcctGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgc                                                                                                                    <  8:218885‑M1/147‑26 (MQ=255)                           tGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacac                                                                                                               <  3:137269‑M1/147‑31 (MQ=255)                               gTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCACcggggctgctgcgcagccctatcgcgacacaagg                                                                                                           <  2:156381‑M1/147‑35 (MQ=255)                                 gggCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggcc                                                                                                         <  5:251209‑M1/147‑37 (MQ=255)                                  ggCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCATGTGTCGCGAAAGGGGCGCGAAGCGGCCCCACGCTTTAATTGCCACTGCATGAACCGCccgggcagctgggcagcccaatagcgccacaagcccc                                                                                                        >  1:16525‑M1/1‑110 (MQ=255)                                        gTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAACGGGCGCGCACCCGCCCCCCGCTTTACATGCCACTGCAAGAACCGCcggggctgcagcgccgccctatcgcgacccaacgccgctccc                                                                                                   >  6:80383‑M1/1‑104 (MQ=255)                                                   gATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGGGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccac                                                                                       <  5:65220‑M1/147‑55 (MQ=255)                                                                                         gtCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATCACCCGCcggggcggctgcgcagccctagcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggtt                                                  <  2:180961‑M1/146‑92 (MQ=255)                                                                                           cgcgAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgccacacaaggccgctaccactccagaccaccacgccctgcgccacacatccaaaccaaggtcacggt                                                  >  8:126565‑M1/1‑53 (MQ=255)                                                                                           cgcgAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTTCGCCGGCCAATCGCGACACAAGGACCCTACCACTCCAGACCACCAGCCCCTGCGCCAGACCTACAAAACAAAGTCAGGTTGc                                                >  8:243070/1‑146 (MQ=19)                                                                                              gAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccaga                                                                                           >  8:132605‑M1/1‑50 (MQ=255)                                                                                              gAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccaga                                                                                           <  7:132605‑M1/100‑51 (MQ=255)                                                                                                aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctccccctccagaccaccaggcccagcgcctcacatcctaaccaaggtcaggttgcagatcc                                          >  2:70369‑M1/1‑48 (MQ=255)                                                                                                                  ccAGGCTTTCATTGCCACTGCATGAACCGCcggggctgcttcgcagccctatcgcgacacaaggccgctcccacaccagaccaccaggcccagcgccacacatcctaaccaaggtaaggttgcagctacgcgtcacgcaccagcct                         >  3:193555‑M1/1‑30 (MQ=255)                                                                                                                       cTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaa                   >  1:932‑M1/1‑25 (MQ=255)                                                                                                                        tttCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaat                  <  7:126565‑M1/147‑124 (MQ=255)                                                                                                                             ttGCCACTGCATGAACCGCgggggctgcggcgctgccctatcgctacaaatggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgac             <  2:1344‑M1/147‑129 (MQ=255)                                                                                                                                 cACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccag          >  8:71289‑M1/1‑15 (MQ=255)                                                                                                                                 cACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccag          <  7:71289‑M1/146‑132 (MQ=255)                                                                                                                                       aTGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggt   >  3:59820‑M1/1‑9 (MQ=255)                                                                                                                                        tGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggta  >  1:24507‑M1/1‑8 (MQ=255)                                                                                                                                               |                                                                                                                                            CGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGC‑‑‑‑GCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTA  >  NC_002947/5681273‑5681550

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.