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)	33 (0.540) +GGC	28/338	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

GGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAG  >  AE015451/4980381‑4980801                                                                                                                                                                                                             |                                                                                                                                                                                                                         ggCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcg                                                                                                                                                                                                              >  1:391873‑M1/1‑205 (MQ=255)                 taagagacaggaggggcgcgcccAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcat                                                                                                                                                                           <  1:329241‑M1/226‑48 (MQ=255)                          agagTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctt                                                                                                                                                                                                                  >  1:315654‑M1/1‑180 (MQ=255)                              tGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaa                                                                                                                                                              <  1:1129974‑M1/236‑61 (MQ=255)                                    aCCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacacca                                                                                                                                                                                           >  2:1031356‑M1/1‑170 (MQ=255)                                          ccTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaac                                                                                                                                                                                               <  2:1128768‑M1/191‑28 (MQ=255)                                                    ggTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatc                                                                                                                                                                             <  1:122853‑M1/199‑46 (MQ=255)                                                         tgGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagcca                                                                                                                                                                                     <  2:529327‑M1/186‑38 (MQ=255)                                                                  cacaGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatgg                                                                                                                          <  1:430735‑M1/236‑97 (MQ=255)                                                                  cacaGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatgg                                                                                                                                                        >  2:427218‑M1/1‑140 (MQ=255)                                                                     aGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaat                                                                                                                                                                              <  1:31322‑M1/181‑45 (MQ=255)                                                                          gCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgga                                                                                                                                              >  2:953193‑M1/1‑132 (MQ=255)                                                                               cGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggc                                                                                                             >  1:952561‑M1/1‑127 (MQ=255)                                                                               cGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtct                                                                                                                                                     >  1:516280‑M1/1‑127 (MQ=255)                                                                                                    gATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgcct                                                                                         <  1:616716‑M1/235‑130 (MQ=255)                                                                                                            ggggcggggGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACAACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgc                                                                                                                                                   <  2:466114‑M1/169‑72 (MQ=255)                                                                                                                           ggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcatt                                                                                                                <  2:190986‑M1/189‑107 (MQ=255)                                                                                                                                               cagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgcgttg                                                                                                                                                                                                                >  2:971425‑M1/1‑63 (MQ=255)                                                                                                                                                          aCAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggacgc                                                                                                                                                                                                                    >  2:21679‑M1/1‑52 (MQ=255)                                                                                                                                                                ccAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaac                                                                                                                                                                                               >  2:963690‑M1/1‑46 (MQ=255)                                                                                                                                                                ccAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaac                                                                                                                                                                                               <  1:963691‑M1/73‑28 (MQ=255)                                                                                                                                                                   gTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttaccagctcgccggtcaggtcgacgttggagtcctccagggc                                     >  1:510262‑M1/1‑43 (MQ=255)                                                                                                                                                                    tAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttaccagctcgccggtcaggtcgacgttggagtcctccagggcaccaccggaaatc                        >  1:651539‑M1/1‑42 (MQ=255)                                                                                                                                                                      aCGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgc                                                                                             >  2:895549‑M1/1‑40 (MQ=255)                                                                                                                                                                      aCGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttg                                                                                                   >  2:1074491‑M1/1‑40 (MQ=255)                                                                                                                                                                      aCGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaag                                                                                                                                                                                        <  1:542587‑M1/74‑35 (MQ=255)                                                                                                                                                                      aCGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaag                                                                                                                                                                                        >  2:542586‑M1/1‑40 (MQ=255)                                                                                                                                                                         aCAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttaccagctcgccggtcaggtcgacgttggagtcctccagggcaccaccggaaatctg                      <  2:651538‑M1/233‑197 (MQ=255)                                                                                                                                                                            gCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgct                                                                                                                                        <  2:1016392‑M1/116‑83 (MQ=255)                                                                                                                                                                            gCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgct                                                                                                                                        >  1:1016393‑M1/1‑34 (MQ=255)                                                                                                                                                                                cTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtg                                                                                                                                 <  2:913382‑M1/119‑90 (MQ=255)                                                                                                                                                                                cTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtg                                                                                                                                 >  1:913383‑M1/1‑30 (MQ=255)                                                                                                                                                                                                  tCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcatggtagttgctctgcgccttgatcaggtttaccagctcgccggtcaggtcgacgttggagtcctccagggcaccaccggaaatctggccgagggtaccggtgtccg  >  1:1031658‑M1/1‑12 (MQ=255)                                                                                                                                                                                                     caTTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcggcttgatcaggttta                                                                             >  1:901956‑M1/1‑9 (MQ=255)                                                                                                                                                                                                     caTTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttaccagctagccggtcaggtcgacgttggagtcc                                             >  2:678047‑M1/1‑9 (MQ=255)                                                                                                                                                                                                           tGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgcctggatcagggttaccagctcggcggtcaggtcgacgttggagtcctccagggcaccaccggcaatctggccgagggtaccgggg      >  2:123745‑M1/1‑3 (MQ=255)                                                                                                                                                                                                             |                                                                                                                                                                                                                         GGCAGGCGGCGTTGTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGGAGTCCTCCAGGGCACCACCGGAAATCTGGCCGAGGGTACCGGTGTCCGGAG  >  AE015451/4980381‑4980801

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.