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)	35 (0.590) +CGGG	30/274	1.5	97.3%	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

CCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGC‑‑‑‑GCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGC  >  NC_002947/5681272‑5681564                                                                                                                                                     |                                                                                                                                                     ccgtaatacatgcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggg                                                                                                                                                      <  2:117807‑M2/1‑1 (MQ=255) ccgtaatacatgcaaacacagacctcctcaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccatgtgacgcgaaacgggcccgatccggccccacgcgttcaatgccactccatgaaacgccgggg                                                                                                                                                      >  6:255144‑M2/149‑149 (MQ=255)  cgtaatacatgcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccggggc                                                                                                                                                     <  2:226438‑M2/2‑1 (MQ=255)      atacatgcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccggggctgc                                                                                                                                                  <  7:62033‑M2/5‑1 (MQ=255)        acatgcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCtgctgc                                                                                                                                               <  1:40485‑M2/8‑1 (MQ=255)         catgcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCtgctgc                                                                                                                                               <  2:301283‑M2/8‑1 (MQ=255)                 cacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTa                                                                                                                                      <  6:322495‑M2/17‑1 (MQ=255)                       cctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGa                                                                                                                                <  7:47664‑M2/23‑1 (MQ=255)                       cctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGa                                                                                                                                <  6:180364‑M2/23‑1 (MQ=255)                         tcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCtatcgcaga                                                                                                                               >  3:341519‑M2/125‑145 (MQ=255)                         tcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATcg                                                                                                                                   <  4:341519‑M2/20‑1 (MQ=255)                              aggtgggctgggttcctgttgtgattctgatggcgtcgttctttgtgggagcggccttgtgtcgcgtaaggggcgcgtagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACaa                                                                                                                           <  6:48190‑M2/28‑1 (MQ=255)                               ggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGc                                                                                                                        <  8:4969‑M2/31‑1 (MQ=255)                                   ggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGAGCACCCCGATCGCGACACAAGGCCGCt                                                                                                                    >  7:344692‑M2/115‑149 (MQ=255)                                    gctgggttcctgttgggatgctgatggctccgttctttgtgggagcggcctggggtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTc                                                                                                                   <  6:271374‑M2/36‑1 (MQ=255)                                                       gctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGc                                                                                                >  5:177206‑M2/95‑149 (MQ=255)                                                                          gtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCGCACATCCTaa                                                                             >  6:182046‑M2/76‑149 (MQ=255)                                                                                   gccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCCGCCCTATCGCGGCACCAGGCCGCTCCCCCTCCAGACCCCCAGGCCCTGCGCCTCCCCTCCTAACCAAGGtcc                                                                    >  1:186490‑M2/67‑148 (MQ=255)                                                                                              cgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATc                                                         >  5:136038‑M2/56‑149 (MQ=255)                                                                                              cgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATc                                                         <  6:177206‑M2/94‑1 (MQ=255)                                                                                              cgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATc                                                         >  5:106909‑M2/56‑149 (MQ=255)                                                                                                   ggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGt                                                    <  6:136038‑M2/99‑1 (MQ=255)                                                                                                             gcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAg                                          >  1:4440‑M2/41‑149 (MQ=255)                                                                                                                  cccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCt                                                                                       >  3:177241‑M2/36‑99 (MQ=255)                                                                                                                  cccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCt                                                                                       <  4:177241‑M2/64‑1 (MQ=255)                                                                                                                  cccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCaa                                                                         <  4:109238‑M2/78‑1 (MQ=255)                                                                                                                  cccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCaa                                                                         >  3:109238‑M2/36‑113 (MQ=255)                                                                                                                       gctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAAt                                <  6:106909‑M2/119‑1 (MQ=255)                                                                                                                                gCCACTGCATGAACCGCCGGGGCTGCGGCGCAGCCCAATCGCGCCACAACGCCGCGCCCCCCCCACACCACCACGCCCAGCGCCCCACAACCAAACAAAGGTCAGGTTGCAGATCCGCGGCACCCCCCCCCCGATCCAAAGccaccacg                       >  4:186570/1‑147 (MQ=255)                                                                                                                                  cactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGCCACAAGGCCGCTCCCACTCCCGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTAAGGTTTCAGATCCGCGTAACGCACCAGCCTTTCCAATTCGACAAGGTc                     >  6:302607‑M2/20‑149 (MQ=255)                                                                                                                                    ctgcatgaaccgccgggGCTGCTGCGCAGCCCTATAGCGCCACAACGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACCACCAAACCAAGGTCAGGATGCAGATCCGCGGCACCGCCCAGCCATTCCAAAGCgaccagctacg                   >  1:125531‑M2/18‑144 (MQ=255)                                                                                                                                    ctgcaggaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCg                   <  2:186490‑M2/132‑1 (MQ=255)                                                                                                                                       catgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCaa                                 <  7:272313‑M2/118‑1 (MQ=255)                                                                                                                                       catgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCaa                                 <  7:32533‑M2/118‑1 (MQ=255)                                                                                                                                       catgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCaa                                 >  8:272313‑M2/15‑132 (MQ=255)                                                                                                                                                ccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCa       >  7:248750‑M2/6‑149 (MQ=255)                                                                                                                                                     gctgctGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGc  <  3:186570/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGC‑‑‑‑GCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGC  >  NC_002947/5681272‑5681564

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.