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)	25 (0.510) +GGC	24/332	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

AATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGC  >  AE015451/4980372‑4980819                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           aaTGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaatt                                                                                                                                                                                                                       <  1:295437‑M1/235‑22 (MQ=255)      gcgcGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaa                                                                                                                                                                                                                  >  2:742864‑M1/1‑209 (MQ=255)           gCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaa                                                                                                                                                                                                            <  1:742867‑M1/236‑33 (MQ=255)               gcggcgTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaatta                                                                                                                                                                                                                      <  1:144732‑M1/222‑23 (MQ=255)                                        ggCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaa                                                                                                                                                                                >  1:589091‑M1/1‑175 (MQ=255)                                                                tgtgGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttc                                                                                                                                                       >  1:355681‑M1/1‑151 (MQ=255)                                                                            acaGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaat                                                                                                                                                                                                     <  2:627156‑M1/178‑40 (MQ=255)                                                                              aGGTAGCGGACGGTGACCAGGCTTTGCCATTGATAGACCGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatct                                                                                                                                                                                   <  2:226391‑M1/194‑58 (MQ=255)                                                                                 tAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgat                                                                                                                                                                            <  2:622922‑M1/198‑65 (MQ=255)                                                                                                    tttGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatc                                                                                                                                                                                               >  1:901366‑M1/1‑115 (MQ=255)                                                                                                              aTAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgc                                                                                                                    <  1:299422‑M1/225‑121 (MQ=255)                                                                                                              aTAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctgg                                                                                                                           <  2:589088‑M1/218‑114 (MQ=255)                                                                                                                          gggggcagcgggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggcc                                                                                                                                                                                                                                        >  2:138950‑M1/1‑93 (MQ=255)                                                                                                                            gggcagcgggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcat                                                                                                                                                                                             <  2:920621‑M1/138‑48 (MQ=255)                                                                                                                                          ggcagAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatg                                                                                                                                                               >  2:401941‑M1/1‑77 (MQ=255)                                                                                                                                             agAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctga                                                                                                                                                                                 <  2:626786‑M1/133‑60 (MQ=255)                                                                                                                                             agAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctga                                                                                                                                                                                 >  1:626789‑M1/1‑74 (MQ=255)                                                                                                                                                       cgaggaaggggaACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtgga                                                                                                                                                 <  2:838836‑M1/150‑92 (MQ=255)                                                                                                                                                          gcaAGGGTAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttaccagctcgccggtcag                                                                               <  2:730884‑M1/218‑158 (MQ=255)                                                                                                                                                                  aaCAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgcctt                                                                                                          >  1:312266‑M1/1‑53 (MQ=255)                                                                                                                                                                  aaCAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctc                                                                                                                  >  1:701‑M1/1‑53 (MQ=255)                                                                                                                                                                      aGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgcctggatcaggtttaccagctcgccggtcaggtcgacgttggagtcctccagggcaccac                                                  >  1:888068‑M1/1‑49 (MQ=255)                                                                                                                                                                          cAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttaccagctcgccggtcaggtcgacgttggagtc                                                                >  1:700083‑M1/1‑45 (MQ=255)                                                                                                                                                                                               tCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttaccagctcgccggtca                                                                                <  2:849959‑M1/180‑157 (MQ=255)                                                                                                                                                                                                       aCCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggggctttcagtggaaatggtcttggcattggcctt                                                                                                                            >  2:309827‑M1/1‑16 (MQ=255)                                                                                                                                                                                                                ttCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcgg                                                                                                    >  1:229627‑M1/1‑7 (MQ=255)                                                                                                                                                                                                                      gggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttaccagctcgccggtcaggtcgacgttggagtcctccagggcaccaccggaaatctggccgagggtaccggtgtccggagcgccgatcaccggca  >  1:647830‑M1/1‑1 (MQ=255)                                                                                                                                                                                                                      gggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttaccagctcgccggtcaggtcgacgttggagtcctccagggcaccac                                                  >  2:267807‑M1/1‑1 (MQ=255)                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           AATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCACGC  >  AE015451/4980372‑4980819

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