Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	AE015451	4,980,585	+GGC	100%	intergenic (‑24/+44)	PP_4387 ← / ← flgE	conserved exported protein of unknown function/Flagellar hook protein FlgE

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 4980585	0 (0.000)	29 (0.340) +GGC	29/350	NT	100%	intergenic (‑24/+44)	PP_4387/flgE	conserved exported protein of unknown function/Flagellar hook protein FlgE
	?	AE015451	4980586 =	0 (0.000)					intergenic (‑25/+43)	PP_4387/flgE	conserved exported protein of unknown function/Flagellar hook protein FlgE

AATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCC  >  AE015451/4980372‑4980820                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            aaTGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgct                                                                                                                                                                                                                                      <  1:155837‑M1/221‑8 (MQ=255)          ggCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacacca                                                                                                                                                                                                              >  1:843743‑M1/1‑205 (MQ=255)                                       tGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctga                                                                                                                                                                                  <  1:335994‑M1/235‑60 (MQ=255)                                            cACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatg                                                                                                                                                                            >  1:741737‑M1/1‑171 (MQ=255)                                              cccAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacac                                                                                                                                                                                                                <  1:118080‑M1/198‑30 (MQ=255)                                                                   gggCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatg                                                                                                                                                                               <  2:826326‑M1/210‑63 (MQ=255)                                                                           cacaGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgc                                                                                                                                                                                                                                  <  2:486210‑M1/151‑12 (MQ=255)                                                                                       ccGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcatt                                                                                                                                   <  2:771883‑M1/234‑107 (MQ=255)                                                                                        cGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattaccaacaccaaagccaatacaatcatcaggtcatc                                                                                                                                                                                     >  2:1104896‑M1/1‑127 (MQ=255)                                                                                            gACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctgg                                                                                                                            <  1:390907‑M1/236‑114 (MQ=255)                                                                                                                            gggcagcgggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCAGATTCTTGGggccgctttccggcccccattacaaacaccaaaggcaatacaatcatc                                                                                                                                                                                             >  1:617991‑M1/1‑91 (MQ=255)                                                                                                                                    ggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggc                                                                                                                                      >  1:1403498‑M1/1‑83 (MQ=255)                                                                                                                                                    gcaccagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttaccagctcgccggtcaggtcgacgttgg                                                                     >  1:1254720‑M1/1‑67 (MQ=255)                                                                                                                                                     caccagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttaccagctcgccggtcaggtcgacgttg                                                                      <  2:1254717‑M1/233‑168 (MQ=255)                                                                                                                                                             aGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagt                                                                                                                        >  2:1135202‑M1/1‑58 (MQ=255)                                                                                                                                                              ggggAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcagggttaa                                                                                               >  2:157570‑M1/1‑57 (MQ=255)                                                                                                                                                                  aaCAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctctgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggttctttcagtggaaatggtcttggcattggcctgg                                                                                                                            >  2:67660‑M1/1‑53 (MQ=255)                                                                                                                                                                   aCAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgcc                                                                                                             >  2:664712‑M1/1‑52 (MQ=255)                                                                                                                                                                     aaGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttaccagctcgccggtcaggtcgacgttggagtcctccagggcac                                                      >  2:1310573‑M1/1‑50 (MQ=255)                                                                                                                                                                        gCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttaccagctcgccggtcaggtcggcggtggaggc                                                                 >  2:212894‑M1/1‑47 (MQ=255)                                                                                                                                                                           aGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttg                                                                                                                                        <  1:836114‑M1/145‑102 (MQ=255)                                                                                                                                                                                     gCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttaccagctcgccggtcaggtcgacgttggagtcctccagggccccaccggaaatctggccgag                                   >  1:1125056‑M1/1‑34 (MQ=255)                                                                                                                                                                                      ccGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatgggctgcatgattgtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgctcaggtttaccagctcgccgg                                                                                    >  1:485799‑M1/1‑33 (MQ=255)                                                                                                                                                                                       cGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttaccagctcgccggtcaggtcgacgttggagtcctccagggcaccaccggaaatctggccgagg                                  >  1:1478968‑M1/1‑32 (MQ=255)                                                                                                                                                                                          ttCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttaccagctcgccggtcaggtcgacgttggagtcctccagggcaccaccggaaatctggccgagggtac                              >  1:210395‑M1/1‑29 (MQ=255)                                                                                                                                                                                            cATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaattatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttaccagctcgccggtcaggtcgacgttggagtcctccagggcaccaccggaaatctggccgag                                   >  2:865325‑M1/1‑27 (MQ=255)                                                                                                                                                                                                  cgTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtgg                                                                                                                                                   <  1:680425‑M1/111‑91 (MQ=255)                                                                                                                                                                                                  cgTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtgg                                                                                                                                                   >  2:680423‑M1/1‑21 (MQ=255)                                                                                                                                                                                                       aCCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttaccagctcgccggtcaggtcgacgttggagtcctccagggcaccaccggaactctggc                                       >  1:1377433‑M1/1‑16 (MQ=255)                                                                                                                                                                                                                 tCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttaccagctcgccggtcaggtcgacgttggagtcctccagggcaccaccggaaatctggccgagggtaccggtgtccggagcgccgatcacc       <  1:1135205‑M1/236‑231 (MQ=255)                                                                                                                                                                                                                      gggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttaccagctcgccggtcaggtcgacgttggagtcctccagggcaccaccggaaatctggccgagggtaccggtgtccggagcgccgatcaccggcac  <  1:1310576‑M1/236‑236 (MQ=255)                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            AATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGCC  >  AE015451/4980372‑4980820

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