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)	33 (1.310) +GGC	28/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

ACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCT  >  NC_002947/4980445‑4980729                                                                                                                                              |                                                                                                                                                caggtagcggccggtgaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTg                                                                                                                                           <  6:137495‑M2/7‑1 (MQ=255)     tagcggccggtgaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGc                                                                                                                                       <  4:88834‑M2/11‑1 (MQ=255)                 accaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACaaa                                                                                                                           <  2:148987‑M2/23‑1 (MQ=255)                  ccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACaaa                                                                                                                           <  6:80950‑M2/23‑1 (MQ=255)                   caggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAaca                                                                                                                         >  7:64831‑M2/124‑148 (MQ=255)                   caggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattagcgtgaaccttcacattcttggggcCGCTTGGCGGCCCCAATTACAAaca                                                                                                                         >  8:101482‑M2/124‑148 (MQ=255)                      gctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACAcc                                                                                                                       >  4:76563‑M2/121‑147 (MQ=255)                           gccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcccattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGc                                                                                                                  >  8:57131‑M2/116‑147 (MQ=255)                              attgatagaccggggcgggggcagcgggggcaggccgaaccagcacccgcaaggggcacaaggccagtaacgccagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCCATTACCACCACCCACGCCAAt                                                                                                              >  8:103809‑M2/113‑148 (MQ=255)                                   tagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACaa                                                                                                          <  7:114834‑M2/40‑1 (MQ=255)                                        cggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAAtcatca                                                                                                    <  1:32243‑M2/46‑1 (MQ=255)                                             cgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCa                                                                                               >  2:90378‑M2/98‑148 (MQ=255)                                                        ggggcaggcagaacccgtaccagcaaggggaacaaggccagtaacgacagccgcttcgttcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAatgatg                                                                                    <  8:79629‑M2/62‑1 (MQ=255)                                                             aggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCt                                                                                >  1:90550‑M2/82‑147 (MQ=255)                                                               gcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATTATGGTctgga                                                                             >  1:66123‑M2/80‑146 (MQ=255)                                                               gcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAAAACAATCATCAGGTAATCTGAATGATGGTCTGCa                                                                             >  8:42964‑M2/80‑148 (MQ=255)                                                                 agaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCatg                                                                           >  6:134915‑M2/78‑148 (MQ=255)                                                                    accagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCatgatg                                                                        <  3:134662‑M2/74‑1 (MQ=255)                                                                           ccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCtt                                                                  >  5:133905‑M2/68‑147 (MQ=255)                                                                             agcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCa                                                               <  1:92590‑M2/83‑1 (MQ=255)                                                                             agcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCa                                                               <  3:60089‑M2/83‑1 (MQ=255)                                                                              gcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAg                                                              >  4:78701‑M2/65‑148 (MQ=255)                                                                                  ggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTAATCTGAATGATGGTCTGGATGTTGGTGGTTTCAGTgg                                                           >  1:27290‑M2/61‑147 (MQ=255)                                                                                     gaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAAt                                                       <  3:68610‑M2/91‑1 (MQ=255)                                                                                         aaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTc                                                   >  4:17816‑M2/54‑148 (MQ=255)                                                                                         aaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCGTTCAGGGGAAATGGTc                                                   >  5:76190‑M2/54‑148 (MQ=255)                                                                                           ggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCCGCCCCAATTACAAACACCAAAGCCACTACAATCATCAGGGCCTCTGAATGATGGTGTGCATGATGGTGCTTTGAGTGGAAATTGTCtt                                                 >  5:105352‑M2/52‑148 (MQ=255)                                                                                             ccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTgg                                               <  3:78701‑M2/99‑1 (MQ=255)                                                                                             ccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTg                                                >  5:66009‑M2/50‑147 (MQ=255)                                                                                                               tcattcgcgtgacccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGtt                               <  8:29000‑M2/115‑1 (MQ=255)                                                                                                                                 acattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAg             <  3:76563‑M2/133‑1 (MQ=255)                                                                                                                                   attcttggggcCGCTTTGCGGCCCCAATTACAAACACCAACGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGttt         >  8:60677‑M2/12‑148 (MQ=255)                                                                                                                                          gggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCATGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTTGCATTGGCCTGTTAGTTGCTCTGCTCATTGATCAGGTTTACCATCt  >  8:114246‑M2/5‑148 (MQ=255)                                                                                                                                              |                                                                                                                                                ACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCT  >  NC_002947/4980445‑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.