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	3 (0.050)	23 (0.430) +CGGG	21/268	2.5	88.8%	intergenic (‑133/+71)	PP_4986/PP_4987	hemolysin III family channel protein/putative Chemotaxis protein
	?	NC_002947	5681416 =	3 (0.050)					intergenic (‑134/+70)	PP_4986/PP_4987	hemolysin III family channel protein/putative Chemotaxis protein

GCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAA‑‑‑CC‑GC‑‑‑‑GCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGC  >  NC_002947/5681283‑5681564                                                                                                                                              |                                                                                                                                                     acatgcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccggggct                                                                                                                                                    <  1:159002‑M2/3‑1 (MQ=255) acatgcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccggggct                                                                                                                                                    >  2:159002‑M2/142‑144 (MQ=255)   atgcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTgcgc                                                                                                                                             >  6:308376‑M2/140‑149 (MQ=255)     gcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccaatgcattaaccgccgggGCTGATGCGCAg                                                                                                                                           >  6:276104‑M2/138‑149 (MQ=255)       aaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGcc                                                                                                                                         <  4:186771‑M2/14‑1 (MQ=255)          cacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTa                                                                                                                                      <  2:96999‑M2/17‑1 (MQ=255)                cctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGa                                                                                                                                <  5:308376‑M2/23‑1 (MQ=255)                    ctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttgatttccacttcatgaaccgcccggGCTGCTGCGCAGCCCTATCGccacaa                                                                                                                             >  5:187373‑M2/123‑147 (MQ=255)                       aggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtggcgaaaggggcgcgaagcggccccaggctttaatttgcactgcctgaaccgccgggGCGTCTGCGCCGCCCGATagcgccaccacg                                                                                                                         >  1:248134‑M2/120‑145 (MQ=255)                                          ttgttGATGCCTTGTTTCTTTGTGGTTGGGGCGTGGTGTGGGTAACGGGGCGGTAAGCGGCCCCAGGCTTTCATTCCCACTCATTCAACGCCCGGG‑‑‑‑GCTGCTGCGCAGCCATATCGCCAAAAAAGGCCGCTCCCACTCCAGAccacca                                                                                                   <  6:187373/144‑1 (MQ=255)                                              atgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGa                                                                                                                                <  7:185470‑M2/23‑1 (MQ=255)                                                gctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGc                                                                                                >  7:359283‑M2/95‑149 (MQ=255)                                                  tgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccctaggctttcattgccactgcatgaaccgccggggc                                                                                                                                                     <  2:333788‑M2/2‑1 (MQ=255)                                                  tgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccctaggctttcattgccactgcatgaaccgccggggc                                                                                                                                                     >  1:333788‑M2/93‑94 (MQ=255)                                                         ttCTTTGTGGGAGCGGCCTTGGGTCGCGAAAGGGGCGCGACGCGGCCCCAGGCTTTCATTCCCACTCCATCAC‑‑‑CC‑GCCGGGGCTGCTGCGCAGCCCTAGCGCGACACAAGGCCGCTCCCACTcc                                                                                                            <  5:364173/124‑1 (MQ=17)                                                                      ggcgcggccttggggcgcgaaaggggcgcgaaccggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCa                                                                          <  5:276104‑M2/77‑1 (MQ=255)                                                                                gcgctAATGGGGCGCTAAGGGGCCCCAGGCTTTAATTCCCACTCATTCAC‑‑‑CC‑GCCGGGGGTGCTGCGAAGCCCTATCGCTAAAAATGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCa                                                             <  7:20198/147‑1 (MQ=255)                                                                                      gcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCACGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTg                                                               >  4:277967‑M2/57‑144 (MQ=255)                                                                                             gggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCATGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGAt                                                          >  1:23382‑M2/50‑142 (MQ=255)                                                                                             gggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGAt                                                          <  2:23382‑M2/93‑1 (MQ=255)                                                                                                 gcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGcgtaa                                                  >  6:71910‑M2/46‑144 (MQ=255)                                                                                                 gcgaagcggccccaggctttcattgccactgcatgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCa                                                  <  5:71910‑M2/101‑1 (MQ=255)                                                                                                                                catgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAAt                                <  6:113713‑M2/119‑1 (MQ=255)                                                                                                                                catgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAAt                                >  5:113713‑M2/15‑133 (MQ=255)                                                                                                                                catgaaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAAt                                >  7:336793‑M2/15‑133 (MQ=255)                                                                                                                                    aaccgccgggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCa                                                                    <  7:362897‑M2/83‑1 (MQ=255)                                                                                                                                           gggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCAc      >  4:108750‑M2/4‑148 (MQ=255)                                                                                                                                            ggGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCt    >  3:203817‑M2/3‑149 (MQ=255)                                                                                                                                              gctgctGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGc  >  6:219568/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                     GCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAA‑‑‑CC‑GC‑‑‑‑GCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGC  >  NC_002947/5681283‑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.