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)	40 (0.640) +GGC	36/272	1.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

GCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGA‑ACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGC  >  NC_002947/4980440‑4980732                                                                                                                                                   |                                                                                                                                                    gCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGA‑ACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggc                                                                                                                                                  <  1:74400‑M1/149‑4 (MQ=255)          tggtCGCGGCCGGTGACCAGGCTTTGCCTTTGTTTGAACGGGGCGGGGGCAGCGGGGGCAGGCAGA‑ACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGg                                                                                                                                                    <  8:352222‑M1/137‑2 (MQ=255)          aGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGA‑ACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCAGCTTCATTCGCCTGAACCTTCACATTCTTGGg                                                                                                                                                    >  7:352222‑M1/1‑137 (MQ=255)            gTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGA‑ACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggc                                                                                                                                       >  8:138712‑M1/1‑135 (MQ=255)                            atgctTTGCCTTTTCTAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGA‑ACCAGCACCATCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacacc                                                                                                                       <  8:230754‑M1/147‑31 (MQ=255)                                          aTAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAGACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgcaaagcggcccc                                                                                                                                    >  1:50276‑M1/1‑106 (MQ=255)                                              aaCGGGGCGGGGGCAGCGGGGGCAGGCAGA‑ACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcat                                                                                                      >  7:132321‑M1/1‑101 (MQ=255)                                              aaCGGGGCGGGGGCAGCGGGGGCAGGCAGA‑ACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcat                                                                                                      >  4:260998‑M1/1‑101 (MQ=255)                                              aaCGGGGCGGGGGCAGCGGGGGCAGGCAGA‑ACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcat                                                                                                      >  2:266392‑M1/1‑101 (MQ=255)                                                   ggcggggGCAGCGGGGGCAGGCGGA‑ACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttggggccccaattacaaacaacaaagacaataaaaacataaggt                                                                                                 >  2:336115‑M1/1‑96 (MQ=255)                                                      gggggcagcgggggcagGCGGA‑ACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatc                                                                                             <  1:336115‑M1/149‑57 (MQ=255)                                                           cagcgggggcagGCAGA‑ACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaat                                                                                        <  7:138712‑M1/149‑62 (MQ=255)                                                                     aggcagA‑ACCAGCACCAGCAAGGGGAACAAGCCCAGTAACGACAGACGCTTCATTCCCGTGAACCTTCACATTCTTGGggccgctttgcggcccccattacaaaccccaaagccaatacaatcatgaggtcatctgaaagatggtatg                                                                               >  2:307429‑M1/1‑78 (MQ=255)                                                                       gcagA‑ACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaaatacaaacaccaaagccaatacaaacatcaggtcataagaatgatggtctgcat                                                                            >  6:314652‑M1/1‑76 (MQ=255)                                                                         agA‑ACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatga                                                                          >  1:135763‑M1/1‑74 (MQ=255)                                                                              ccagcaccagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggt                                                                      <  1:266392‑M1/149‑80 (MQ=255)                                                                                     cagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttca                                                               <  2:135763‑M1/149‑87 (MQ=255)                                                                                     cagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttca                                                               >  4:181253‑M1/1‑63 (MQ=255)                                                                                         aaGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtg                                                            >  5:120491‑M1/1‑59 (MQ=255)                                                                                           ggggAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaa                                                         <  4:321537‑M1/149‑93 (MQ=255)                                                                                           ggggAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtgga                                                          >  2:78233‑M1/1‑57 (MQ=255)                                                                                               aaCAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatgg                                                     <  1:78233‑M1/149‑97 (MQ=255)                                                                                                aCAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcgtctgaatgatggtctgcatgatggtgctttcagtggaaatggt                                                    <  8:132321‑M1/149‑98 (MQ=255)                                                                                                          tAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattg                                          <  3:260998‑M1/149‑108 (MQ=255)                                                                                                            tcgACCGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggc                                        <  6:365997‑M1/148‑110 (MQ=255)                                                                                                            aCGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaaaaaaaaacaccaaacccaaaacaatcatcaggtcatctgaattatggtctgcatgaaggtgctttcagtggaaatggtcttggcattggc                                        >  1:190726‑M1/1‑40 (MQ=255)                                                                                                                   ccGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtag                                 <  5:314652‑M1/149‑117 (MQ=255)                                                                                                                    cGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcgggcccaatttcaaaccccaaagccaaaacaatcatcagggcatctgtatcatggtctgcatgatggtgctttaagtggatatggtctaggcaatggcctggtagt                                >  6:358085‑M1/1‑32 (MQ=255)                                                                                                                              gcgTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtgg                                                           >  6:178254‑M1/1‑22 (MQ=255)                                                                                                                              gcgTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtgg                                                           <  5:178254‑M1/112‑91 (MQ=255)                                                                                                                                gTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgcct                    >  7:123888‑M1/1‑20 (MQ=255)                                                                                                                                gTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgcct                    <  2:190726‑M1/149‑130 (MQ=255)                                                                                                                                 tGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgcctt                   <  6:120491‑M1/149‑131 (MQ=255)                                                                                                                                      cTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgat                >  5:383393‑M1/1‑14 (MQ=255)                                                                                                                                      cTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttg                                                <  1:74934‑M1/115‑102 (MQ=255)                                                                                                                                      cTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttg                                                >  2:74934‑M1/1‑14 (MQ=255)                                                                                                                                         cacaTTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctggatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggt           >  1:192922‑M1/1‑11 (MQ=255)                                                                                                                                         cacaTTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcagg            <  3:181253‑M1/148‑138 (MQ=255)                                                                                                                                              aattggggccgcaaagcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggta                                  <  5:11797‑M1/119‑116 (MQ=255)                                                                                                                                              aattggggccgcaaagcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggta                                  >  6:11797‑M1/3‑6 (MQ=255)                                                                                                                                               cTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttac       >  8:115473‑M1/1‑5 (MQ=255)                                                                                                                                                  ggggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttaccagct  >  7:127095‑M1/1‑2 (MQ=255)                                                                                                                                                   gggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttaccagct  <  7:115473‑M1/148‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGA‑ACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGC  >  NC_002947/4980440‑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.