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.940) +GGC	28/284	0.2	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

GCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCC  >  NC_002947/4980457‑4980733                                                                                                                                  |                                                                                                                                                    ggtgaccaggctttgccattgatagcacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggacagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTa                                                                                                                                   <  3:121318‑M2/19‑1 (MQ=255)     accaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACaaa                                                                                                                               <  1:159195‑M2/23‑1 (MQ=255)            tttgccattgatagaccggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCaaa                                                                                                                        <  2:175051‑M2/30‑1 (MQ=255)                 cattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcccattcttggggcCGCTTTGCGGCCCCAATTACAAACACAAAAGCCa                                                                                                                    >  3:108545‑M2/114‑147 (MQ=255)                  attgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCCATTACAACCACCAAAGCCAAt                                                                                                                  >  4:122591‑M2/113‑148 (MQ=255)                   ttgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGGGGGCCCCATTACAAACCCCCAAGCCAAt                                                                                                                  >  1:166818‑M2/112‑147 (MQ=255)                           acggggcgggggcagcgggggcaggcagaaccagcaccagcaaggtgaccaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAatcatc                                                                                                         <  3:114706‑M2/45‑1 (MQ=255)                           acggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAatca                                                                                                           >  3:65569‑M2/104‑146 (MQ=255)                                        agcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaacctacacattcttggggcCGCTTTGCGGCCCCAAATACAAACACCAAAGCCAACACAATAATCAGGTCATCAGaa                                                                                             >  3:35416‑M2/91‑147 (MQ=255)                                                    cagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaacctgcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCat                                                                                <  8:130942‑M2/70‑1 (MQ=255)                                                         ccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATgg                                                                           >  7:124597‑M2/74‑148 (MQ=255)                                                              accagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCtt                                                                      >  7:154240‑M2/69‑148 (MQ=255)                                                                        tgaacaccgccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAACGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGaaa                                                            <  6:125075‑M2/90‑1 (MQ=255)                                                                              aggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGGATGATGGTGCTTTCAGTGGAAATGGTCt                                                      >  8:68804‑M2/53‑148 (MQ=255)                                                                                gccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACCCCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTg                                                    >  8:16637‑M2/51‑148 (MQ=255)                                                                                            agccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATAAAATCATCAGGTCATCTGAAAGATGGTCTGCATTCTGGTGCTTTAAGTGGAAATGGTAGTGGCATTTTCCTgg                                        >  5:110632‑M2/39‑148 (MQ=255)                                                                                                         gcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTgcg                           >  7:127428‑M2/26‑148 (MQ=255)                                                                                                          cgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTgcgc                          <  2:155785‑M2/124‑1 (MQ=255)                                                                                                           gtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCCAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGcc                         >  3:136982‑M2/24‑148 (MQ=255)                                                                                                             gaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCtt                       <  4:3421‑M2/127‑1 (MQ=255)                                                                                                               accttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTg                      <  1:23304‑M2/128‑1 (MQ=255)                                                                                                                 cttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATAATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTTGCCTGGTAGTTGCTCTGCGCCTtgc                     >  6:89199‑M2/18‑145 (MQ=255)                                                                                                                  ttcacattcttggggcCGCTGTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGAGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCa                  <  3:60481‑M2/132‑1 (MQ=255)                                                                                                                     agattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGGCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAg                 <  4:57477‑M2/133‑1 (MQ=255)                                                                                                                     acattcttggggcCGCTTTGCGGACCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAgg                <  7:132488‑M2/134‑1 (MQ=255)                                                                                                                       attcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGttt             >  2:144814‑M2/12‑148 (MQ=255)                                                                                                                           ttggggcCGCTTTGCGGCCCCAATTATAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCa         >  8:99316‑M2/8‑148 (MQ=255)                                                                                                                              gggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCt      >  2:116465‑M2/5‑148 (MQ=255)                                                                                                                               ggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTAcc          >  5:95114‑M2/4‑143 (MQ=255)                                                                                                                                gcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCa         <  6:95114‑M2/141‑1 (MQ=255)                                                                                                                                 cCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGc   >  4:106940‑M2/2‑148 (MQ=255)                                                                                                                                  cgcTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGGTCTGCGCCTTGATCAGGTTTACCAGCTCGcc  >  6:144453/1‑148 (MQ=255)                                                                                                                                  |                                                                                                                                                    GCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCC  >  NC_002947/4980457‑4980733

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.