Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	NC_002947	4,980,585	+GGC	100%	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)	43 (0.850) +GGC	37/260	NT	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

AACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGC  >  NC_002947/4980444‑4980732                                                                                                                                               |                                                                                                                                                   acaggtagcggccggtgaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcagtcgcgtgaaccttcacattcttggggcCGCTTTg                                                                                                                                              >  8:160097‑M2/143‑149 (MQ=255)  caggtagcggccggtgaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacaacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGc                                                                                                                                             <  3:200288‑M2/8‑1 (MQ=255)   aggtagcggccggtgaccaggcattgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaacctacacattcttgcggcCGCATTGCg                                                                                                                                            >  4:300092‑M2/141‑149 (MQ=255)          ggccggtgaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCAAAGCGGCCCCaa                                                                                                                                     <  1:238412‑M2/16‑1 (MQ=255)             cggggaccaggctatgccattgaaacaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTa                                                                                                                                  <  7:304082‑M2/19‑1 (MQ=255)               gtgacaaggctttgccattgatctaccggggcgggggcagcgggggctggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACa                                                                                                                                <  3:267539‑M2/21‑1 (MQ=255)                   ccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAaca                                                                                                                            >  5:119545‑M2/125‑149 (MQ=255)                             ccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCaaa                                                                                                                       >  1:104239‑M2/115‑144 (MQ=255)                             ccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCaaa                                                                                                                       <  2:104239‑M2/30‑1 (MQ=255)                                        acggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAatcat                                                                                                         >  3:110292‑M2/104‑147 (MQ=255)                                            ggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACaaa                                                                                                                              >  4:256356‑M2/100‑122 (MQ=255)                                            ggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACaaa                                                                                                                              <  3:256356‑M2/23‑1 (MQ=255)                                              cgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCAt                                                                                                 <  4:131995‑M2/52‑1 (MQ=255)                                                        gggggcaggcagaaccagcaccagcaaggggaacaagtccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAatgatg                                                                                       <  7:296650‑M2/62‑1 (MQ=255)                                                             caggaagaaccagcaccagcaaggggaacatggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTg                                                                                  <  3:300092‑M2/67‑1 (MQ=255)                                                                 cagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGccccaagaat                                                                                                                                 >  4:270295‑M2/79‑94 (MQ=255)                                                                 cagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGccccaagaat                                                                                                                                 <  3:270295‑M2/20‑5 (MQ=255)                                                                    aaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCa                                                                                                                                      <  4:78574‑M2/15‑1 (MQ=255)                                                                    aaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCa                                                                                                                                      >  3:78574‑M2/76‑90 (MQ=255)                                                                     accagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATgg                                                                          >  1:1091‑M2/75‑149 (MQ=255)                                                                         gcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACa                                                                                                                                <  7:281781‑M2/21‑1 (MQ=255)                                                                         gcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATAGAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCt                                                                      >  5:279464‑M2/71‑149 (MQ=255)                                                                               gcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGt                                                                >  7:97963‑M2/65‑149 (MQ=255)                                                                               gcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTCCGGCCCCAATTACAAACACCCAAGCCAATACAATCATCAGGTCAACTGAATGCTGGACTGCATGAAGGTGCTATCagg                                                                >  3:15749‑M2/65‑148 (MQ=255)                                                                               gcaagcggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattctttgggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATTGTCTGCATGATTGTGCTTTCAGt                                                                <  4:310489‑M2/85‑1 (MQ=255)                                                                                caaggggaacaaggccagtaacgccagccgcttcattcgcgtgaaccgtcacagtcttggggcCGCTTTGCGGCCCCAATTACAAACACCACAGCCAATCCACTCATCAGGTCAGCTGAATGATGGGCGGCATGATGGTGCGTTcaggg                                                               >  7:236966‑M2/64‑147 (MQ=255)                                                                                          aaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCt                                                     <  3:92731‑M2/96‑1 (MQ=255)                                                                                          aaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCt                                                     <  4:110292‑M2/96‑1 (MQ=255)                                                                                           aggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCtt                                                    >  7:36096‑M2/53‑149 (MQ=255)                                                                                           aggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCtt                                                    <  1:939‑M2/97‑1 (MQ=255)                                                                                            ggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTg                                                   <  3:158733‑M2/98‑1 (MQ=255)                                                                                            ggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTg                                                   <  2:298451‑M2/98‑1 (MQ=255)                                                                                            ggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTg                                                   <  6:279464‑M2/98‑1 (MQ=255)                                                                                                 gtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCAtt                                              <  1:141168‑M2/103‑1 (MQ=255)                                                                                                  taacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATGTGAATGATGGTCTGCATGATGGTGCGTTCAGTGGAAATGGTCTTGGCATTg                                             >  4:82639‑M2/46‑149 (MQ=255)                                                                                                    acgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGc                                           <  8:36096‑M2/106‑1 (MQ=255)                                                                                                      gacagccgcttcattcgcgtgaaaattaacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATAAAAAAATCAGGTCATAAGAATGAAGGTCTGCATGAAGGTGCGTTAAGTGGAAATGGTATTGGCATTGgccg                                         >  7:49477‑M2/42‑148 (MQ=255)                                                                                                            cgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGt                                   <  2:17046‑M2/114‑1 (MQ=255)                                                                                                             gcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGtt                                  >  1:128849‑M2/35‑149 (MQ=255)                                                                                                                     cgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTgc                           <  6:119545‑M2/122‑1 (MQ=255)                                                                                                                                  acattcttggggcCGCTTTGCGGCCCCAATTACAAACCCCCACGCCAATACAAGCCTCCGGGCATCTGAATGATGGGCTGCCTGATGGTGCTTTCAGGGGAAATGGTCTTGGCCTTGGCCGGGTAGTTGCTCTGCGCCTTGATCAGGt              <  8:236966‑M2/135‑1 (MQ=255)                                                                                                                                   cattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGttt            <  2:1091‑M2/137‑1 (MQ=255)                                                                                                                                   cattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATAAAATCAGCAGGTCATCTGAATGATGGTGTGGATGATGGTGCTTTAAGGGGAAATGGTCTTGGCAGTGGCCTGGTAGTTGCGCTGCGCCGTGATCAGGttt            >  4:262936‑M2/13‑149 (MQ=255)                                                                                                                                             gcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGc  >  2:133374‑M2/3‑149 (MQ=255)                                                                                                                                             gcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTCGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACGAGCTCGc  >  4:123556‑M2/3‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                   AACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGC  >  NC_002947/4980444‑4980732

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.