Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	NC_002947	5,555,608	+GCC	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.030)	18 (0.510) +GCC	18/282	1.2	94.8%	intergenic (‑84/+87)	PP_4887/PP_4888	hypothetical protein/methyl‑accepting chemotaxis transducer
	?	NC_002947	5555609 =	1 (0.030)					intergenic (‑85/+86)	PP_4887/PP_4888	hypothetical protein/methyl‑accepting chemotaxis transducer

GGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5555476‑5555608 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTG  >  NC_002947/5555609‑5555730                                                                                                                                      |||                                                                                                                           GGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCC                                                                                                                >  4:47764/1‑148   AAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCTAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGC                                                                                                              <  2:194647/148‑1     ATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGT                                                                                                            <  3:47764/148‑1                                  GGGCGGACATGAGCACCTCCTTTTGTCATGGGGCGCCGACAGCATACCAACCCCGCTCAGGTGTGGGCGCGGCCTTGTGGCGCTAACGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCA                                                                               <  1:213326/148‑1                                     CGGACATGCGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGGCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCC                                                                            <  1:73815/148‑1                                            GAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGC                                                                     <  8:4777/148‑1                                                 ACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACAC                                                                <  6:32483/148‑1                                                   TCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTAGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTTGGGCTGCTATGCAGCCCGTTAGCGACAAAA                                                              >  6:98816/1‑148                                                                 GTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGG                                                                                          >  6:52929/1‑106                                                                  TCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGG                                                                                         <  5:52929/106‑1                                                                   CGCCGACAGCATAGCAACACCGCTCTGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCTACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACC                                               <  8:106322/147‑1                                                                         CAGCATAGCAACACCGCGCAGGTGTGGGAGCGGCCTTGTGGCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGG                                        <  1:48897/148‑1                                                                          AGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAAT                                                                                                                    >  2:10339/1‑71                                                                           GCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATT                                                                                                                   <  1:10339/71‑1                                                                              TAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAA                                                    <  7:19699/131‑1                                                                                                GTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAG                 <  2:10851/148‑1                                                                                                           CTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCGGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCA      <  1:16398/148‑1                                                                                                               TGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGAAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTACGGACCGGAACTGCCCCAGCCTCGCGCGCCGCTc  >  2:108438/1‑147                                                                                                                                      |||                                                                                                                           GGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5555476‑5555608 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTG  >  NC_002947/5555609‑5555730

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.