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	0 (0.000)	28 (0.590) +GGC	26/278	0.4	100%	intergenic (‑24/+44)	PP_4387/flgE	hypothetical protein/flagellar hook protein FlgE
	?	NC_002947	4980586 =	0 (0.000)					intergenic (‑25/+43)	PP_4387/flgE	hypothetical protein/flagellar hook protein FlgE

GCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCT  >  NC_002947/4980440‑4980729                                                                                                                                                   |                                                                                                                                                gaacacaggtagcggccggtgaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcc                                                                                                                                                 >  1:36179‑M2/147‑147 (MQ=255)        ggtagcggccggtgaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGc                                                                                                                                          >  1:262212‑M2/140‑147 (MQ=255)         gatgcgcccggtgacctggctgcctctttgatagaccgggccggggtcagcgggggcaggcagaaccagcaccagcaaggggaactaggccagtcacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCg                                                                                                                                         <  7:234201‑M2/9‑1 (MQ=255)               gccggtgaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCa                                                                                                                                   >  5:69332‑M2/133‑147 (MQ=255)                    tgaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgcaccttcacattcttggggcCGCTTTGCGGCCCCAATTa                                                                                                                               >  5:134513‑M2/128‑146 (MQ=255)                      accaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACaa                                                                                                                            <  1:255431‑M2/22‑1 (MQ=255)                         aggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAaca                                                                                                                         >  2:235865‑M2/123‑147 (MQ=255)                          gcctttgccattgaaagcacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAacac                                                                                                                        <  5:259418‑M2/26‑1 (MQ=255)                             tttgccattgatacaacgggtggggggcatcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCaa                                                                                                                     <  3:231584‑M2/29‑1 (MQ=255)                                 ccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGc                                                                                                                  >  5:214123‑M2/115‑146 (MQ=255)                                    ttgatagaaccgcgcggcggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAAt                                                                                                              <  8:193723‑M2/36‑1 (MQ=255)                                                               ggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATgg                                                                                   <  2:262212‑M2/63‑1 (MQ=255)                                                                   ggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCa                                                                                               <  5:179281‑M2/51‑1 (MQ=255)                                                                                   gcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttgggggCGCGTTGCGGCCCCAATTACAAACACCACAGCCAATACAATCATCAGGTGATCTGAATGATGGGCTGCATGATGGTGCGTTCa                                                               >  6:178402‑M2/65‑147 (MQ=255)                                                                                      aggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTg                                                            <  6:214123‑M2/86‑1 (MQ=255)                                                                                        gggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGa                                                          <  6:69332‑M2/88‑1 (MQ=255)                                                                                               aggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTc                                                   >  5:234395‑M2/53‑147 (MQ=255)                                                                                                cgccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCt                                                  <  4:260819‑M2/96‑1 (MQ=255)                                                                                                  ccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTg                                                >  3:62961‑M2/50‑147 (MQ=255)                                                                                                  ccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATAATCAGGTCATATGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTg                                                >  6:20546‑M2/50‑147 (MQ=255)                                                                                                         ctacagccgctgcattcgagtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACGAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTgg                                         <  6:134513‑M2/105‑1 (MQ=255)                                                                                                          gacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATGAGGTCATCTGAATGATGGTCGGCATGATGGTGGTTTGAGTGGATATGGTGTTGGCTTtggg                                        >  4:48591‑M2/42‑146 (MQ=255)                                                                                                                cgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCAAGTGAATGATGGTCTGCATGATGGTGCGTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTa                                  >  8:163652‑M2/36‑147 (MQ=255)                                                                                                                cgcttcattcgcgagaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAAAGAAGGTCTGCATGAAGGTGCTTTCAGTGGAATTGGTCTTGGCATTGGCCAGGTa                                  >  4:56568‑M2/36‑147 (MQ=255)                                                                                                                 gcttcattcgcgtgaacgttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACATTCATCAGGTCATCAGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAg                                 <  2:134246‑M2/113‑1 (MQ=255)                                                                                                                               aaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTAAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCt                    >  3:101372‑M2/21‑146 (MQ=255)                                                                                                                                accttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCtt                   <  3:128828‑M2/127‑1 (MQ=255)                                                                                                                                         ttcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGtt          >  3:108586‑M2/11‑146 (MQ=255)                                                                                                                                                ggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCt  >  2:247945‑M2/4‑147 (MQ=255)                                                                                                                                                   |                                                                                                                                                GCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCT  >  NC_002947/4980440‑4980729

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.