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)	36 (0.710) +GCC	32/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

CGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCG‑‑‑CCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATC  >  NC_002947/5555461‑5555745                                                                                                                                                      |                                                                                                                                             cGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAgcg                                                                                                                                              <  3:61826/148‑1 (MQ=255)    tGTCACTGTGCAGGAAATTCAGATCGGGCTGGGCAGCGATTAGTCGGGCTGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccc                                                                                                                                          <  4:90153‑M1/149‑5 (MQ=255)        aCTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGGCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggc                                                                                                                                      <  7:239896‑M1/149‑9 (MQ=255)              cAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTAGCGAAAGGGCCCCAAAGCGgccccggccattta                                                                                                                                >  8:253568‑M1/1‑135 (MQ=255)                ggAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTTTCGCGAAAGGGCCGCATAGCGgccccggcaatttccg                                                                                                                              >  8:5827‑M1/1‑133 (MQ=255)                ggAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccg                                                                                                                              >  7:153091‑M1/1‑133 (MQ=255)                   aaTTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgt                                                                                                                           <  8:153091‑M1/149‑20 (MQ=255)                   aaTTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcg                                                                                                                            <  2:288106‑M1/148‑19 (MQ=255)                          gTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAAGGCATAGCAACACCGCTCAGGTGGGGGAGCGCCATTGTGTCGCCAAAGGGCCGCATAGCGgccccggcaatttccgggtgtat                                                                                                                       >  1:260575‑M1/1‑125 (MQ=255)                               ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaa                                                                                                               >  7:7399‑M1/1‑118 (MQ=255)                                 tGGGCTGCGATTTGGGGGGGGGACAGGAGCAACTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCGTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgc                                                                                                                     <  6:34210‑M1/141‑26 (MQ=255)                                 tGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgc                                                                                                                     >  5:34210‑M1/1‑116 (MQ=255)                                    gCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTTTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctgg                                                                                                          <  7:5827‑M1/149‑37 (MQ=255)                                     cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggg                                                                                                         >  5:330138‑M1/1‑112 (MQ=255)                                      aGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctgggg                                                                                                        <  8:7399‑M1/149‑39 (MQ=255)                                         gATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctg                                                                                                     <  1:12838‑M1/149‑42 (MQ=255)                                             aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttcc                                                                                                                               >  1:116227‑M1/1‑104 (MQ=255)                                             aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttcc                                                                                                                               <  2:116227‑M1/119‑16 (MQ=255)                                                 gggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcag                                                                                             <  4:50012‑M1/149‑50 (MQ=255)                                                  ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagc                                                                                            <  2:25400‑M1/149‑51 (MQ=255)                                                             gCAACTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcga                                                                                  <  6:330138‑M1/148‑61 (MQ=255)                                                                 cTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTTTCGCGAAAGGGCCGCATAGCGgccccgggaatttccgcctgtatggaacaaaccaggggctgcgatgcagcccgttggcgccac                                                                               >  1:34000‑M1/1‑84 (MQ=255)                                                                    cTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaag                                                                            <  2:260575‑M1/147‑67 (MQ=255)                                                                              gtgtCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgt                                                                                                                           <  8:32704‑M1/90‑20 (MQ=255)                                                                              gtgtCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgt                                                                                                                           >  7:32704‑M1/1‑71 (MQ=255)                                                                                                        cGCTCAGGTGTGGGAGCGGCCATGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatggaacaaacctgtggatgctatgcagccctttcgcgccaccaggcccacaaaggaccgcgcaggcggcaggcccggacc                                        >  8:289804‑M1/1‑47 (MQ=255)                                                                                                             aGGTGTGGGAGCGGCCTTGTTTAGCCCAAGGGCCGCAAAGCGGCCCCCGCCATTTTCCCGGGTATTCAAAAAACCATGGGCTGCTATGCAGCCCGTTAGCGACACAACGCCCACAAAAGACCCCCCAGGCGTAAAGCCCCGGACAGGcc                                   >  8:179509/1‑149 (MQ=255)                                                                                                             aGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggc                                    <  7:179509‑M1/148‑107 (MQ=255)                                                                                                              ggTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggccc                                  <  5:169562‑M1/149‑109 (MQ=255)                                                                                                                             ttGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagc                   >  5:6955‑M1/1‑26 (MQ=255)                                                                                                                               gtgtCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggccccggctggcgcgcagctg                 >  3:228721‑M1/1‑24 (MQ=255)                                                                                                                                tgtCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctga                <  7:125557‑M1/149‑127 (MQ=255)                                                                                                                                  tCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagc              >  3:115274‑M1/1‑21 (MQ=255)                                                                                                                                    gcgAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccg            >  1:59758‑M1/1‑19 (MQ=255)                                                                                                                                     cgAAAGGGCCGCCTAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccgc           <  4:115274‑M1/149‑132 (MQ=255)                                                                                                                                       ccggggccgCTATGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccgcc          >  8:167293‑M1/4‑16 (MQ=255)                                                                                                                                        cggggccGCTATGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgccgccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccgc           >  5:156045‑M1/3‑15 (MQ=255)                                                                                                                                          gggCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccgccagat      <  7:253568‑M1/149‑137 (MQ=255)                                                                                                                                             ccGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccgccagatcat   >  4:23127‑M1/1‑10 (MQ=255)                                                                                                                                              ctcatAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccgccagatcatc  <  7:250939‑M1/147‑141 (MQ=255)                                                                                                                                                      |                                                                                                                                             CGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACA‑‑GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCG‑‑‑CCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATC  >  NC_002947/5555461‑5555745

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.