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

CTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCA  >  AE015451/4980369‑4980816                                                                                                                                                                                                                          |                                                                                                                                                                                                                                       aatgggcgcggcaggcggcgttgttgtaggggccagagtggcgcacccagcctttgcccgggttgtgggccgaacacaggtagcggccggtgaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAAtt                                                                                                                                                                                                                       <  1:295437‑M2/18‑1 (MQ=255)      gcgcggcaggcggcgttgttgtaggggccagagtggcgcacccagcctttgcccgggttgtgggccgaacacaggtagcggccggtgaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACaaa                                                                                                                                                                                                                  >  2:742864‑M2/213‑235 (MQ=255)           gcaggcggcgttgttgtaggggccagagtggcgcacccagcctttgcccgggttgtgggccgaacacaggtagcggccggtgaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCaa                                                                                                                                                                                                            <  1:742867‑M2/29‑1 (MQ=255)               gcggcgttgttgtaggggccagagtggcgcacccagcctttgcccgggttgtgggccgaacacaggtagcggccggtgaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTa                                                                                                                                                                                                                      <  1:144732‑M2/19‑1 (MQ=255)                                        ggcgcacccagcctttgcccgggttgtgggccgaacacaggtagcggccggtgaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGaa                                                                                                                                                                                >  1:589091‑M2/179‑235 (MQ=255)                                                                tgtgggccgaacacaggtagcggccggtgaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTc                                                                                                                                                       >  1:355681‑M2/155‑236 (MQ=255)                                                                            acaggtagcggccggtgaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAAt                                                                                                                                                                                                     <  2:627156‑M2/36‑1 (MQ=255)                                                                              aggtagcggacggtgaccaggctttgccattgatagaccggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCt                                                                                                                                                                                   <  2:226391‑M2/54‑1 (MQ=255)                                                                                 tagcggccggtgaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAatgat                                                                                                                                                                            <  2:622922‑M2/61‑1 (MQ=255)                                                                                                    tttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAatc                                                                                                                                                                                               >  1:901366‑M2/119‑160 (MQ=255)                                                                                                              atagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTgg                                                                                                                           <  2:589088‑M2/110‑1 (MQ=255)                                                                                                              atagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGc                                                                                                                    <  1:299422‑M2/117‑1 (MQ=255)                                                                                                                          gggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcc                                                                                                                                                                                                                                        >  2:138950‑M2/97‑97 (MQ=255)                                                                                                                            gggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAatcat                                                                                                                                                                                             <  2:920621‑M2/44‑1 (MQ=255)                                                                                                                                          ggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCatgatg                                                                                                                                                               >  2:401941‑M2/81‑154 (MQ=255)                                                                                                                                             agaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGa                                                                                                                                                                                 <  2:626786‑M2/56‑1 (MQ=255)                                                                                                                                             agaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGa                                                                                                                                                                                 >  1:626789‑M2/78‑133 (MQ=255)                                                                                                                                                       cgaggaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGa                                                                                                                                                 <  2:838836‑M2/88‑1 (MQ=255)                                                                                                                                                          gcaagggtaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAg                                                                               <  2:730884‑M2/154‑1 (MQ=255)                                                                                                                                                                  aacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGctc                                                                                                                  >  1:701‑M2/57‑175 (MQ=255)                                                                                                                                                                  aacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCtt                                                                                                          >  1:312266‑M2/57‑183 (MQ=255)                                                                                                                                                                      aggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTGGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGcaccac                                                  >  1:888068‑M2/53‑235 (MQ=255)                                                                                                                                                                          cagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGtc                                                                >  1:700083‑M2/49‑217 (MQ=255)                                                                                                                                                                                               tcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCa                                                                                <  2:849959‑M2/153‑1 (MQ=255)                                                                                                                                                                                                       accttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGGGCTTTCAGTGGAAATGGTCTTGGCATTGGCctt                                                                                                                            >  2:309827‑M2/20‑127 (MQ=255)                                                                                                                                                                                                                ttcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGatcgg                                                                                                    >  1:229627‑M2/11‑141 (MQ=255)                                                                                                                                                                                                                      gggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGcaccac                                                  >  2:267807‑M2/5‑187 (MQ=255)                                                                                                                                                                                                                      gggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCa  >  1:647830‑M2/5‑235 (MQ=255)                                                                                                                                                                                                                        gcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCa  >  2:551770‑M2/3‑233 (MQ=255)                                                                                                                                                                                                                          |                                                                                                                                                                                                                                       CTCAATGGGCGCGGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAGCGCCGATCACCGGCA  >  AE015451/4980369‑4980816

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.