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)	26 (0.520) +CGGG	25/252	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

ACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTG  >  NC_002947/5681271‑5681563                                                                                                                                                 |                                                                                                                                                     aCCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTTATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcgg                                                                                                                                                   <  5:166032‑M1/148‑4 (MQ=255)  ccGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcgggg                                                                                                                                                 <  1:2023‑M1/149‑6 (MQ=255)           aTGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgc                                                                                                                                        >  4:292483‑M1/1‑135 (MQ=255)             gCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcag                                                                                                                                      >  5:234048‑M1/1‑133 (MQ=255)             gCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgca                                                                                                                                       <  8:51010‑M1/148‑16 (MQ=255)             gCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcg                                                                                                                                         >  7:51010‑M1/1‑133 (MQ=255)              cAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagc                                                                                                                                     <  6:23070‑M1/149‑18 (MQ=255)               aaaCACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcg                                                                                                                                                    >  6:255049‑M1/1‑131 (MQ=255)               aaaCACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcg                                                                                                                                                    <  5:255049‑M1/133‑3 (MQ=255)                   acaGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcgggg                                                                                                                                                 >  1:181011‑M1/1‑127 (MQ=255)                   acaGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcgggg                                                                                                                                                 <  2:181011‑M1/132‑6 (MQ=255)                    caGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatc                                                                                                                               <  3:292483‑M1/149‑24 (MQ=255)                                   gggCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgac                                                                                                                          >  6:280169‑M1/1‑111 (MQ=255)                                   gggCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgac                                                                                                                          <  5:280169‑M1/139‑29 (MQ=255)                                    ggCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgct                                                                                                               <  6:234048‑M1/149‑40 (MQ=255)                                     gcggggttcCTGTTGTGTTGCTGAGGGGTTCGTTCTTGGGGGGAGCGGCCTTGTGTCGCTAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTCCATCAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctc                                                                                                              <  8:294020‑M1/146‑41 (MQ=255)                                         ggTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGAcggggctgctgcgcagccctatagcgccacaaggcctctacca                                                                                                           >  5:240774‑M1/1‑105 (MQ=255)                                              cttttgtgatGCTGATGTCTTCGTTCTTTGTCGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGATTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggcctctcccactcca                                                                                                      <  6:240774‑M1/145‑49 (MQ=255)                                                           gATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCAATGCATTAAACGGcggggctgctgcgcaggcctatcgcgacacaagtccgctcccactccagaccaccaggccct                                                                                        >  7:75724‑M1/1‑87 (MQ=255)                                                                 ttCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcg                                                                                                                                                    <  3:132547‑M1/83‑3 (MQ=255)                                                                 ttCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcg                                                                                                                                                    >  4:132547‑M1/1‑81 (MQ=255)                                                                                    gCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggt                                                                 >  7:243015‑M1/1‑62 (MQ=255)                                                                                                 aaaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccacca                                                                                              <  5:22854‑M1/105‑57 (MQ=255)                                                                                                   aGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcagg                                                             <  5:70497‑M1/136‑90 (MQ=255)                                                                                                   aGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcagg                                                             >  6:70497‑M1/1‑47 (MQ=255)                                                                                                    ggggCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgccacaaggccgctcccactccagaccacccggccctgcgcctcacctcctaaccaaggtcaggttgcagatccgcgt                                               >  6:78447‑M1/1‑46 (MQ=255)                                                                                                                  ccccAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctt                                 >  4:153842‑M1/1‑32 (MQ=255)                                                                                                                      aGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttcc                              >  1:217147‑M1/1‑28 (MQ=255)                                                                                                                       ggCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaa                            <  3:153842‑M1/149‑123 (MQ=255)                                                                                                                        gCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaat                           <  8:243015‑M1/149‑124 (MQ=255)                                                                                                                            tCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcga                       >  3:125207‑M1/1‑22 (MQ=255)                                                                                                                              aTTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgac                      >  3:206095‑M1/1‑20 (MQ=255)                                                                                                                                   cACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcccgcaccagccgttccaatgcgacccggtc                >  7:284658‑M1/1‑15 (MQ=255)                                                                                                                                          tGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcagggtgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggtacc         <  4:206095‑M1/149‑142 (MQ=255)                                                                                                                                                 ccggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccacaaggccctgcgccacacatcctaaacaaggtcaggatgcagatccgcgtcacgcaccagcctttccaattcgaccacgtcacgtacccccgcca  >  2:263273‑M1/1‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     ACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTG  >  NC_002947/5681271‑5681563

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.