Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	AE015451	867,950	45 bp→TTCGGTA	9.9%	intergenic (+48/+85)	PP_0750 → / ← mqo‑I	conserved protein of unknown function/malate:quinone oxidoreductase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 867949	31 (1.200)	3 (0.120) +TTCGGTA	3/304	NT	9.9%	intergenic (+47/+130)	PP_0750/mqo‑I	conserved protein of unknown function/malate:quinone oxidoreductase
	?	AE015451	867995 =	26 (1.010)					intergenic (+93/+84)	PP_0750/mqo‑I	conserved protein of unknown function/malate:quinone oxidoreductase

TGGTCAAGGCCTACCGTGGCCTGCGCGTGGACGACTTCGAGCGCTTTGTAGGCTTCTTCGTCGAAGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGG  >  AE015451/867724‑868178                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                      tGGTCAAGGCCTACCGTGGCCTGCGCGTGGACGACTTCGAGCGCTTTGTAGGCTTCTTCGTCGAAGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCgcgc                                                                                                                                                                                                                               <  1:319168/234‑1 (MQ=255)              ccGTGGCCTGCGCGTGGACGACTTCGAGCGCTTTGTAGGCTTCTTCGTCGAAGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAgg                                                                                                                                                                                                                                       <  2:336266/213‑1 (MQ=255)                 tGGCCTGCGCGTGGACGACTTCGAGCGCTTTGTAGGCTTCTTCGTCGAAGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGcgccg                                                                                                                                                                                                                             >  1:214875/1‑220 (MQ=255)                   gCCTGCGCGTGGACGACTTCGAGCGCTTTGTAGGCTTCTTCGTCGAAGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTAc                                                                                                                                                                                                             <  2:446056/234‑1 (MQ=255)                   gCCTGCGCGTGGACGACTTCGAGCGCTTTGTAGGCTTCTTCGTCGAAGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTAc                                                                                                                                                                                                             <  2:396644/234‑1 (MQ=255)                        cgcgTGGACGACTTCGAGCGCTTTGTAGGCTTCTTCGTCGAAGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAAcc                                                                                                                                                                                                       <  1:148588/235‑1 (MQ=255)                                                 tAGGCTTCTTCGTCGAAGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAaga                                                                                                                                                                             <  1:351442/236‑1 (MQ=255)                                                  aGGCTTCTTCGTCGAAGCTGGCCTGGACCTGGACGGCAAGGAGGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCgcgc                                                                                                                                                                                                                               >  2:256538/1‑185 (MQ=255)                                                      tccttcGTCGAAGCTGGCCTGGACCTGGCCGGCAAGGACGCGCACGGCAAGCCCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAag                                                                                                                                                                              <  1:49881/228‑1 (MQ=255)                                                           cgtcgAAGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCAt                                                                                                                                                                                        <  2:379043/215‑1 (MQ=255)                                                                                         ggACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGGGGGAGCGGGCATGCCcgcg                                                                                                                                                                                 >  2:182485/1‑192 (MQ=255)                                                                                                 aCGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCACCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCg                                                                                                                                                                                                                     <  2:299763/148‑1 (MQ=255)                                                                                                                          tGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCgg                                                                                                                                                                                            <  2:312253/148‑1 (MQ=255)                                                                                                                                  aGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCGACAGGAGTTGCGCCGCACCTGCGGTCATTACCGAACCTGTgg                                                                                                                                                                                                  <  2:515860/134‑1 (MQ=255)                                                                                                                                                        cGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGtt                                                                                                                        >  1:251704/1‑186 (MQ=255)                                                                                                                                                          agagCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCg                                                                                                                                                                        <  1:506756/136‑1 (MQ=255)                                                                                                                                                           gaTCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCcgc                                                                                                                                                                                  >  1:221762/1‑125 (MQ=255)                                                                                                                                                           gaTCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGCCCTCTTCGCGGGCATGCCCGCTCCCACAGGttcggtaatgcccgcgaagagggcggtacaggctgatcaagcgcc                                                                                                                                                                                          >  2:469112‑M1/1‑72 (MQ=255)                                                                                                                                                           gaTCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGCCCTCTTCGCGGGCATGCCCGCTCCCACAGGttcggtaatgcccgcgaagagggcggtacaggctgatcaagcgcc                                                                                                                                                                                          <  1:469116‑M1/117‑46 (MQ=255)                                                                                                                                                           gaTCATCGATAAAGCCCGTGGCTCATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCcgc                                                                                                                                                                                  <  2:221758/125‑1 (MQ=255)                                                                                                                                                                        gCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATg                                                                                                                                                                                       <  1:89411/107‑1 (MQ=255)                                                                                                                                                                            gTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCcgc                                                   >  1:497448/1‑235 (MQ=255)                                                                                                                                                                             tGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAgg                                                                                                                                                 <  2:383674/140‑1 (MQ=255)                                                                                                                                                                              ggCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGcc                                                                                                                                    >  2:374936/1‑152 (MQ=255)                                                                                                                                                                                                  cGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTCCCGAACCTGTGGGAGCGGGCATGCCCGCGAaga                                                                                                                                                                             <  1:64339/91‑1 (MQ=255)                                                                                                                                                                                                  cGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAACAACGCCCCGAGGGGc                                                                                                                           >  2:198038/1‑141 (MQ=255)                                                                                                                                                                                                         ttCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGcc                                                                                                                                    <  2:251700/125‑1 (MQ=255)                                                                                                                                                                                                            gcgGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAgc                    >  2:223026/1‑234 (MQ=255)                                                                                                                                                                                                            gcgGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCG‑TTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGctctct                                                                    <  2:251267/185‑1 (MQ=255)                                                                                                                                                                                                            gcgGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCG‑TTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTa                        >  2:270908/1‑229 (MQ=255)                                                                                                                                                                                                               ggCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGctctct                                                                    <  1:5211/183‑1 (MQ=255)                                                                                                                                                                                                                   cgcCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTaa                                                                                                  <  1:81314/149‑1 (MQ=255)                                                                                                                                                                                                                      ccGCTCCCACAGGttcggtaatgcccgcgaagagggcggtacaggctgatcaagcgccaggcacaaaaaacgccccgaggggcgtttttgcgtttcagaccggcttaagcgtagtgcttggccgccgggctgctctctgccaggtccagggccgcgtcgttatcggcactgatcttgtggtacagcgcagcatcggttgccagggtcttctcgcgagccgggaagatttccttgag        >  1:452481‑M1/1‑13 (MQ=255)                                                                                                                                                                                                                            ccACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAg   <  1:182489/235‑1 (MQ=255)                                                                                                                                                                                                                             cacaGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAgg  >  1:290269/1‑235 (MQ=255)                                                                                                                                                                                                                              acaGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGctctct                                                                    <  1:256542/168‑1 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                      TGGTCAAGGCCTACCGTGGCCTGCGCGTGGACGACTTCGAGCGCTTTGTAGGCTTCTTCGTCGAAGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGG  >  AE015451/867724‑868178

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