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	1 (0.020)	37 (0.630) +GCC	30/260	NT	97.4%	intergenic (‑84/+87)	PP_4887/PP_4888	hypothetical protein/methyl‑accepting chemotaxis transducer
	?	NC_002947	5555609 =	1 (0.020)					intergenic (‑85/+86)	PP_4887/PP_4888	hypothetical protein/methyl‑accepting chemotaxis transducer

ACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCG‑‑‑CCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGC  >  NC_002947/5555468‑5555749                                                                                                                                             |                                                                                                                                                 aCTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccgg                                                                                                                                           <  5:132990‑M1/148‑8 (MQ=255)   tgtgCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaa                                                                                                                                        >  2:229413‑M1/1‑139 (MQ=255)        aGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttcc                                                                                                                                   >  8:158460‑M1/1‑134 (MQ=255)        aGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccgacaatttcc                                                                                                                                   >  3:44989‑M1/1‑134 (MQ=255)            aaTTCAGGTCGGACTGGGCTGCGATTAGTCGGGCGGAAATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcg                                                                                                                                <  8:79755‑M1/148‑19 (MQ=255)                 aGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtat                                                                                                                           <  7:188148‑M1/148‑24 (MQ=255)                   gTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgc                                                                                                                         >  4:86565‑M1/1‑123 (MQ=255)                    tCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaa                                                                                                                       <  1:229413‑M1/149‑28 (MQ=255)                    tCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaa                                                                                                                       <  6:74791‑M1/149‑28 (MQ=255)                        ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaa                                                                                                                   <  4:80392‑M1/149‑32 (MQ=255)                           gggCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacct                                                                                                                >  1:30669‑M1/1‑115 (MQ=255)                            ggCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctg                                                                                                               <  6:16586‑M1/149‑36 (MQ=255)                             gCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctgg                                                                                                              <  7:158460‑M1/149‑37 (MQ=255)                                  gATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGCAAGGGCCGCATAGCGgccccggcaatttgcgcgggtatggcacaaacctggggctg                                                                                                         >  7:32164‑M1/1‑108 (MQ=255)                                   aTTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgc                                                                                                        <  3:86565‑M1/149‑43 (MQ=255)                                               gACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgcccccgtcaatttccgcgtgtatgcaacaaacctggggctgcatagcagcccc                                                                                            <  1:279606‑M1/149‑55 (MQ=255)                                               gACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgcccccgtcaatttccgcgtgtatgcaacaaacctggggctgcatagcagcccc                                                                                            <  1:284322‑M1/149‑55 (MQ=255)                                                 cATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgcatagcagccccagg                                                                                          <  6:22320‑M1/149‑57 (MQ=255)                                                   tGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgc                                                                                        <  2:30669‑M1/149‑59 (MQ=255)                                                    gAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTAGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcg                                                                                       >  7:339884‑M1/1‑90 (MQ=255)                                                    gAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTAGCGAAAGGGCCGCAAAGCGgccccgggcatttgcgcgggtttggaaaaaaccaggggcggcgatgcagcccgttcgcg                                                                                       >  3:264120‑M1/1‑90 (MQ=255)                                                             cTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTTTAGCCAAATGGCCCCATAGCGGCCCCGGGAATTTTCGCGGGTTTTCCACCAACCAGGGGCGGGTATGCCGCCCGTTCGCGGCACACggg                                                                              >  5:202133/1‑148 (MQ=255)                                                                ttGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggccc                                                                            <  7:70982‑M1/148‑71 (MQ=255)                                                                    cctgtgtcGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGGCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaa                                                                       <  3:40526‑M1/147‑76 (MQ=255)                                                                               gACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcag                                                            <  5:210647‑M1/149‑87 (MQ=255)                                                                                aCAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcagg                                                           <  7:241768‑M1/149‑88 (MQ=255)                                                                                    cATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggctt                                                        >  2:45222‑M1/1‑58 (MQ=255)                                                                                          aaCACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggccc                                                 <  4:347180‑M1/149‑98 (MQ=255)                                                                                           acacCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttc                                                                                          >  6:31993‑M1/1‑51 (MQ=255)                                                                                           acacCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttc                                                                                          <  5:31993‑M1/107‑57 (MQ=255)                                                                                                            gAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctgg                               <  3:309942‑M1/149‑116 (MQ=255)                                                                                                              gCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgtaggcttcaggcccggaactggcccaggctggcg                             <  1:45222‑M1/149‑118 (MQ=255)                                                                                                               cGGCCTTGTGCCGCTAACGGGCCGCATCGCGcccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgc                            <  5:250424‑M1/149‑119 (MQ=255)                                                                                                                   cTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcag                        <  4:264120‑M1/149‑123 (MQ=255)                                                                                                                           gcgAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccg                >  6:21730‑M1/1‑19 (MQ=255)                                                                                                                                    ccGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccgccagatcat       >  6:254020‑M1/1‑10 (MQ=255)                                                                                                                                     cGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccgccagatcatc      <  7:322852‑M1/149‑141 (MQ=255)                                                                                                                                     cGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccgccagatcatc      <  5:216504‑M1/149‑141 (MQ=255)                                                                                                                                        aTAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccgccagatcatccag   <  6:202133‑M1/149‑144 (MQ=255)                                                                                                                                         tAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccgccagatcatccagc  <  8:339884‑M1/149‑145 (MQ=255)                                                                                                                                             |                                                                                                                                                 ACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCG‑‑‑CCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGC  >  NC_002947/5555468‑5555749

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.