Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	NC_002947	4,980,585	+GGC	intergenic (‑24/+44)	PP_4387 ← / ← flgE	hypothetical protein/flagellar hook protein FlgE

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 4980585	1 (0.030)	24 (0.790) +GGC	24/284	0.3	96.1%	intergenic (‑24/+44)	PP_4387/flgE	hypothetical protein/flagellar hook protein FlgE
	?	NC_002947	4980586 =	1 (0.030)					intergenic (‑25/+43)	PP_4387/flgE	hypothetical protein/flagellar hook protein FlgE

GGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGG‑‑‑CGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACC  >  NC_002947/4980439‑4980725                                                                                                                                                   |                                                                                                                                                ggCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGg                                                                                                                                                <  5:32947‑M1/148‑2 (MQ=255) ggCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGg                                                                                                                                                >  3:22439‑M1/1‑147 (MQ=255)  gCCGCACACAGGGAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGg                                                                                                                                                <  3:40092‑M1/147‑2 (MQ=255)  gCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGgg                                                                                                                                               >  7:48446‑M1/1‑146 (MQ=255)                           cgaggCTTTGCCATTGATAGAACGGGGCGGGTGCAGCGTGGGCAGGCAGAACCACCACCAGCAAGCGGAACCAGGCCAGTAACGACACCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaa                                                                                                                       <  1:33861‑M1/145‑27 (MQ=255)                            cAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaaca                                                                                                                     >  4:158679‑M1/1‑120 (MQ=255)                                 tttGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaa                                                                                                                >  3:10367‑M1/1‑115 (MQ=255)                                             aGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatc                                                                                                    >  8:3038‑M1/1‑103 (MQ=255)                                               aaCGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatca                                                                                                   >  8:160029‑M1/1‑101 (MQ=255)                                                                 ggggcagGCTGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatg                                                                                <  1:92530‑M1/148‑66 (MQ=255)                                                                             accagcaccagcAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatg                                                                    <  3:150387‑M1/148‑78 (MQ=255)                                                                                 gcaccagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCCCTATGCGGCCCCCATTACCACCACCCCAGCCCATACAAAAATCACGTCCTCTGCAAGAAGGTCTGCATGATGGTGc                                                                >  2:165799/1‑148 (MQ=255)                                                                                     cagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgcgttgcggcccccattacaaacaccaaagccaatacaatcatcaggtcagctgcatgatggtctgcagggttgtgcgttg                                                            >  7:56514‑M1/1‑63 (MQ=255)                                                                                       gcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattaaaaacaccaaagccaataaaatcatcaggtcatcttaatgatggtctgcatgatggtgctttaag                                                          >  1:56456‑M1/1‑61 (MQ=255)                                                                                         aaGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtg                                                        <  7:157675‑M1/148‑90 (MQ=255)                                                                                         aaGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagt                                                         >  4:143609‑M1/1‑59 (MQ=255)                                                                                                     gCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttg                                            <  4:10367‑M1/148‑102 (MQ=255)                                                                                                          tAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaatttcaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcatt                                       >  1:148306‑M1/1‑42 (MQ=255)                                                                                                                cAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctg                                 >  1:34112‑M1/1‑36 (MQ=255)                                                                                                                   ccGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggta                              <  8:48446‑M1/148‑116 (MQ=255)                                                                                                                      cTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagtt                           <  3:10031‑M1/148‑119 (MQ=255)                                                                                                                          aTTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctc                       >  7:12307‑M1/1‑26 (MQ=255)                                                                                                                           ttCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctct                      >  6:121454‑M1/1‑25 (MQ=255)                                                                                                                               cgTGAACCTTCACATTCTTGGggccgctttgccgccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttgagtggaaatggtcttggcattggcctggtagttgctctgcgc                  >  5:121283‑M1/1‑21 (MQ=255)                                                                                                                                      cTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttga             >  7:66582‑M1/1‑14 (MQ=255)                                                                                                                                       ttCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtgatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatca          >  2:19430‑M1/1‑13 (MQ=255)                                                                                                                                          acaTTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcagggcatctgaatgatggtctgcatgatggtgctttcagtgtaaatggtcttggcattggcctggtagttgctctgcgccttgatcagg        <  5:134707‑M1/147‑138 (MQ=255)                                                                                                                                            aTTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggttt     >  2:56754‑M1/1‑8 (MQ=255)                                                                                                                                                ttGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttacc  <  1:156676‑M1/147‑144 (MQ=255)                                                                                                                                                   |                                                                                                                                                GGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGG‑‑‑CGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACC  >  NC_002947/4980439‑4980725

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.