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)	35 (1.090) +GGC	32/284	0.1	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

CACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCG  >  NC_002947/4980446‑4980731                                                                                                                                             |                                                                                                                                                  aggtagcggccgttgaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGc                                                                                                                                            <  4:23719‑M2/8‑1 (MQ=255) aggtagcggccggtgaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGc                                                                                                                                            >  7:158606‑M2/141‑148 (MQ=255)    tagcggccggtgaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGc                                                                                                                                         <  3:102827‑M2/11‑1 (MQ=255)               gaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACaa                                                                                                                              <  4:93457‑M2/22‑1 (MQ=255)                  caggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAaca                                                                                                                           >  6:1750‑M2/124‑148 (MQ=255)                          gccattgatagaccggggcgggggcagcgggggcaggcagaaccagcaccagcaagggtaacaaggccagtaacgacagccgcttcagtcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGc                                                                                                                    <  7:89302‑M2/32‑1 (MQ=255)                          gccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGc                                                                                                                    >  6:93921‑M2/116‑147 (MQ=255)                           ccattcatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgctacattcgcgtgaacctacacattcttggggcCCGATTGCGGCCCCAATTACAAACCACAAAGCCa                                                                                                                  >  8:87520‑M2/115‑148 (MQ=255)                            cattgatagaacggggcgggggcagcgggggcaggcagatccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCaa                                                                                                                 <  2:69041‑M2/35‑1 (MQ=255)                            cattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCa                                                                                                                  >  3:169773‑M2/114‑147 (MQ=255)                                gatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATAc                                                                                                              >  8:94892‑M2/110‑147 (MQ=255)                                  tagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACaa                                                                                                            >  4:50834‑M2/108‑147 (MQ=255)                                     aacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAatca                                                                                                         >  8:32547‑M2/105‑147 (MQ=255)                                                      gggggcaggcagtaccagccgcaggaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAatgat                                                                                       <  6:94304‑M2/61‑1 (MQ=255)                                                              gcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCa                                                                               >  4:99929‑M2/80‑148 (MQ=255)                                                              gcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCa                                                                               >  6:107145‑M2/80‑148 (MQ=255)                                                               cagaaccagcacgagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCat                                                                              <  3:130106‑M2/70‑1 (MQ=255)                                                                       gcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAAAACAATCATCAGGTAATCTGAATGATGGTCTGCATGATGGTGc                                                                      >  8:147369‑M2/71‑148 (MQ=255)                                                                         accagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCtt                                                                    <  1:29978‑M2/80‑1 (MQ=255)                                                                             gcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTACATGATGGTGCTTTCAg                                                                >  1:2394‑M2/65‑148 (MQ=255)                                                                                   ggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGaaa                                                          <  3:50834‑M2/90‑1 (MQ=255)                                                                                        aaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCGAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTc                                                     >  5:30411‑M2/54‑148 (MQ=255)                                                                                          ggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTACATGATGGTGCTTTCAGTGGAAATGGTCtt                                                   <  2:2394‑M2/97‑1 (MQ=255)                                                                                           gccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTg                                                  <  6:51766‑M2/98‑1 (MQ=255)                                                                                           gccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTg                                                  <  5:93921‑M2/98‑1 (MQ=255)                                                                                             cagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTGTGCATGATGGTGCTTTCAGTGGAAATGGTCGTGGc                                                >  6:7133‑M2/49‑148 (MQ=255)                                                                                                taacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCAtt                                             >  7:105451‑M2/46‑148 (MQ=255)                                                                                                   cgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTgg                                           <  8:158606‑M2/105‑1 (MQ=255)                                                                                                          cgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGGCTTGGCATTGGCCTGGTAg                                   <  4:179633‑M2/113‑1 (MQ=255)                                                                                                                    gcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTgcg                         >  8:115682‑M2/26‑148 (MQ=255)                                                                                                                        gaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCtt                     <  7:94892‑M2/127‑1 (MQ=255)                                                                                                                            cttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGa                   >  7:126595‑M2/18‑146 (MQ=255)                                                                                                                                acattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAg               <  5:107145‑M2/133‑1 (MQ=255)                                                                                                                                acattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAg               <  6:127030‑M2/133‑1 (MQ=255)                                                                                                                                 cattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGtt            <  8:44509‑M2/136‑1 (MQ=255)                                                                                                                                           gcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCg  >  3:94020‑M2/3‑148 (MQ=255)                                                                                                                                             |                                                                                                                                                  CACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCG  >  NC_002947/4980446‑4980731

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ 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.