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)	15 (0.470) +GCC	15/284	1.3	93.9%	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

CAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5555474‑5555608 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGC  >  NC_002947/5555609‑5555749                                                                                                                                        |||                                                                                                                                              CAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTAGCGAAACGGCCGCCTAGCGGCCCCGGCcaatt                                                                                                                                     >  2:64811/1‑143            GGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATG                                                                                                                          <  5:50914/148‑1            GGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTAT                                                                                                                           <  4:165669/147‑1             ggCGGCCTGGGCAGCGATTAGGCGGGCGGACATGAGCAACTCCTTTTGTCATGTGGCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGC                                                                                                                         <  1:71990/146‑1                       GCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTG                                                                                                               <  2:121776/148‑1                                                GCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGA                                                                                      <  5:136707/148‑1                                                 CAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGAC                                                                                     <  8:134388/148‑1                                                   ACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACAC                                                                                   <  2:103222/148‑1                                                     TCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTAGCGAAAGGGCCGCCTAGCGGCCCCGGCAATTTACGCGGGGATTGAAAAAAACAGGGGCTGCTATTCACCCCCTTCGCGACACAA                                                                                 >  6:165314/1‑148                                                                   GTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGGGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGCGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGA                                                                    <  7:34835/147‑1                                                                                                            CCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCCCAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGC                            >  5:165741/1‑146                                                                                                                             GCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATC         >  3:74684/1‑148                                                                                                                              CCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCA        >  2:112261/1‑148                                                                                                                                    AGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCCTCCAGC  >  1:71905/1‑148                                                                                                                                    AGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGC  <  3:128258/148‑1                                                                                                                                        |||                                                                                                                                              CAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5555474‑5555608 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGC  >  NC_002947/5555609‑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.