Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	NC_002947	5,555,608	+GCC	100%	intergenic (‑84/+87)	PP_4887 ← / ← PP_4888	hypothetical protein/methyl‑accepting chemotaxis transducer

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 5555608	0 (0.000)	37 (0.730) +GCC	29/254	NT	100%	intergenic (‑84/+87)	PP_4887/PP_4888	hypothetical protein/methyl‑accepting chemotaxis transducer
	?	NC_002947	5555609 =	0 (0.000)					intergenic (‑85/+86)	PP_4887/PP_4888	hypothetical protein/methyl‑accepting chemotaxis transducer

ACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCT  >  NC_002947/5555468‑5555753                                                                                                                                             |                                                                                                                                                  aCTGTGCCGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggc                                                                                                                                           <  6:104396‑M1/149‑9 (MQ=255)     tgCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaat                                                                                                                                        <  4:38628‑M1/148‑12 (MQ=255)       cAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaattt                                                                                                                                      <  3:285728‑M1/148‑14 (MQ=255)        aGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttcc                                                                                                                                    >  3:83342‑M1/1‑134 (MQ=255)          gAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgc                                                                                                                                  >  5:132516‑M1/1‑132 (MQ=255)                             gCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctgg                                                                                                               <  6:267813‑M1/149‑37 (MQ=255)                                   aTTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTTTGTCGCGAAAGGGCCGCATATCGgccccggcaatttccgagtggatgcaaacaaacaggggctg                                                                                                          >  5:148473‑M1/1‑107 (MQ=255)                                       gTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCGGACAGCATACCAACACCGCTCAGGTGTGGGAGCGGCTTTGTGTCGCGATAGGGCCTCATAGCGgcctcggcaatttccgcgtgtatgcaacaaacctgg                                                                                                               <  8:204325‑M1/139‑37 (MQ=255)                                       gTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatg                                                                                                     >  2:70084‑M1/1‑103 (MQ=255)                                       gTCGGGCGGACATGAGCAACTCCATTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctgg                                                                                                               >  7:204325‑M1/1‑103 (MQ=255)                                           ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGTGAAAGGGCCGCCAAGCGgccccggcaaattacgcgggtatggaacgaaccaggggctgggatgaagc                                                                                                 >  8:113121‑M1/1‑99 (MQ=255)                                            gCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTTTAGCGCAAGGGCCCCATATCGgccccggcaatttgcgcgtgtatgcaacaaaccaggggctgctatgcagcc                                                                                                >  3:38838‑M1/1‑98 (MQ=255)                                                   tGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgc                                                                                         <  6:132516‑M1/149‑59 (MQ=255)                                                    gAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgc                                                                                         >  4:266397‑M1/1‑90 (MQ=255)                                                    gAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgc                                                                                         <  3:266397‑M1/148‑59 (MQ=255)                                                      gCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagc                                                                                                 <  8:190785‑M1/138‑51 (MQ=255)                                                      gCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagc                                                                                                 >  7:190785‑M1/1‑88 (MQ=255)                                                       cAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaac                                                                                                                       >  4:313195‑M1/1‑87 (MQ=255)                                                       cAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaac                                                                                                                       <  3:313195‑M1/115‑29 (MQ=255)                                                           tCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTAGCGAAAGGGCCCCATAGCGgccccggcaatttacgcgtgtttgcaacaaaactggggctgctattcagcccgttcgcgacaca                                                                                   >  3:289522‑M1/1‑83 (MQ=255)                                                                     aTGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaag                                                                       <  1:70084‑M1/149‑77 (MQ=255)                                                                               gACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctgggg                                                                                                             <  1:124552‑M1/101‑39 (MQ=255)                                                                               gACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctgggg                                                                                                             >  2:124552‑M1/1‑63 (MQ=255)                                                                                          aaCACCGCTCAGGTGTGGGAGCGGCCTTGTGTAGCGAAAGGACCGCCTAGCGgccccggcaatttgcgcgggtatgcaacaaaccaggggctgcgatgcagccctttcgcgccacaaggcccccaaaggaccccgcaggcgtgaggccc                                                  >  4:75231‑M1/1‑52 (MQ=255)                                                                                                    atgtgtggGAGCGGCCTTGTGTCGCGAAAGGGCCGCATATCGgccccttaaatttccgcgtgtattcaacaaccctgtggctgatatgaagccctttcgcgacacaaggcccacaaatgaccgcgcaggcttcaggccaggaactggac                                        <  6:148473‑M1/147‑108 (MQ=255)                                                                                                    aGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggc                                         <  4:83342‑M1/148‑107 (MQ=255)                                                                                                    aGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactgg                                          <  7:113121‑M1/147‑106 (MQ=255)                                                                                                         tgGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggc                                   <  4:38838‑M1/149‑113 (MQ=255)                                                                                                                 gCCTTGTGTCGCGAAAGAGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactgg                                          <  5:205632‑M1/134‑106 (MQ=255)                                                                                                                 gCCTTGTGTCGCGAAAGAGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactgg                                          >  6:205632‑M1/1‑29 (MQ=255)                                                                                                                          ctcgaaAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagcc                  <  7:61452‑M1/147‑130 (MQ=255)                                                                                                                           gcgAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagcc                  <  5:108143‑M1/148‑130 (MQ=255)                                                                                                                            cgAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccgc                >  2:233070‑M1/1‑18 (MQ=255)                                                                                                                             gAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgccgctgagcc                  >  8:61452‑M1/1‑17 (MQ=255)                                                                                                                               cggggccGCTATGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccg                 >  7:41081‑M1/3‑15 (MQ=255)                                                                                                                               acgggccGCATAGCGgccccggcaatttccgcgtgtatgcaacaaccctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccgcc               <  5:48920‑M1/145‑133 (MQ=255)                                                                                                                                 gggCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcggctcagctgcgcgccagcc           <  4:289522‑M1/149‑137 (MQ=255)                                                                                                                                      gCATCGCGcccccgtcaattaccgcgtttatccaacaaacctggggcggctatgcagccctttcgctacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccgccagatcatcc      <  3:3699‑M1/149‑142 (MQ=255)                                                                                                                                          aGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaacgaccgcgcaggcttcaggccaggaactggcccaggctggcgcgaagcagagccgccagatcatccagca  >  4:183179‑M1/1‑4 (MQ=255)                                                                                                                                             |                                                                                                                                                  ACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCT  >  NC_002947/5555468‑5555753

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.