Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	NC_000913	2,720,868	(CGCTGG)2→3	coding (916/2661 nt)	pka →	protein lysine acetyltransferase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2720855	1 (0.030)	24 (0.790) +GCGCTG	21/232	0.4	96.0%	coding (903/2661 nt)	pka	protein lysine acetyltransferase
	?	NC_000913	2720856 =	1 (0.030)					coding (904/2661 nt)	pka	protein lysine acetyltransferase

TATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCACGAGCTGTTTTCGGCGGTGGAAACCCTTAGCCATATGCGCCCGCTACGTGGCGACCGGCTGATGATTATCAGCAACGGTGCTGCGCCTGCCGCGCTGGCGCTGGATGCCTTATGGTCACGCAATGGCAAGCTGGCAACGCTAAGCGAAGAAACCTGCCAGAAACTGCGCGATGCACTGCCAGAACATGTGGCAATATCTAACCCGCTCGATCTACGCGATGACGCCAGCAGT  >  NC_000913/2720726‑2720996                                                                                                                                   |                                                                                                                                             gcgtgccggtttgttgcgggtacaggacacccacgagctgttttcggcggtggaaacccttagccatatgcgcccgctacgtggcgaccggctgatgattatcagcaacggtgctgcgcctgccgcgctgGCGCTGGCGCTGGATGCCt                                                                                                                            <  2:7304‑M2/19‑1 (MQ=255)               tgcgggtacaggacacccacgagctgttttcggcggtggaaacccttagccatatgcgcccgctacgtggcgaccggctgatgattatcagcaacggtgctgcgcctgccgcgctgGCGCTGGCGCTGGATGCCTTATGGTCACGCAAt                                                                                                              <  1:128539‑M2/33‑1 (MQ=255)                gcgggtacaggacacccacgagctgttttcggcggtggaaacccttagccatatgcgcccgccacgtggcgaccggctgatgattatcagcaacggtgctgcgcctgccgcgctgGCGCTGGCGCTGGATGCCTTATGGTCACGCAATg                                                                                                             >  2:76542‑M2/116‑149 (MQ=255)                 cgggtacaggacacccacgagctgttttcggcggtggaaacccttagccatatgcgcccgcgacgtggcgaccggctgatgattatcagcaacggtgctgcgcctgccgcgctgGCGCTGGCGCTGGATGCCTTATGGTCACGCAATg                                                                                                             <  5:107034‑M2/34‑1 (MQ=255)                        aggacacccacgagctgttttcggcggtggaaacccttagccatatgcgcccgctacgtggcgaccggctgatgattatcagcaacggtgctgcgcctgccgcgctgGCGCTGGCGCTGGATGCCTTATGGTCACGCAATGGCAAGCTg                                                                                                     <  7:42128‑M2/42‑1 (MQ=255)                        aggacacccacgagctgttttcggcggtggaaacccttagccatatgcgcccgctacgtggcgaccggctgatgattatcagcaacggtgctgcgcctgccgcgctgGCGCTGGCGCTGGATGCCTTATGGTCACGCAATGGCAAGCTg                                                                                                     <  1:139208‑M2/42‑1 (MQ=255)                        aggacacccacgagctgttttcggcggtggaaacccttagccatatgcgcccgctacgtggcgaccggctgatgattatcagcaacggtgctgcgcctgccgcgctgGCGCTGGCGCTGGATGCCTTATGGTCACGCAATGGCAAGCTg                                                                                                     >  8:36177‑M2/108‑149 (MQ=255)                              cccacgagctgttttcggcggtggaaacccttagccatatgcgcccgctacgtggcgaccggctgatgattatcagcaacggtgctgcgcctgccgcgctgGCGCTGGCGCTGGATGCCTTATGGTCACGCAATGGCAAGCTGGCAACg                                                                                               <  7:36177‑M2/48‑1 (MQ=255)                                 acgagctgttttcggcggtggaaacccttagccatatgcgcccgctacgtggcgaccggctgatgattatcagcaacggtgctgcgcctgccgcgctggcgctggcgc                                                                                                                                     >  3:57861‑M2/99‑108 (MQ=255)                                 acgagctgttttcggcggtggaaacccttagccatatgcgcccgctacgtggcgaccggctgatgattatcagcaacggtgctgcgcctgccgcgctggcgctggcgc                                                                                                                                     <  4:57861‑M2/10‑1 (MQ=255)                                 acgagctgttttcggcggtggaaacccttagccatatgcgcccgctacgtggcgaccggctgatgattatcagcaacggtgctgcgcctgccgcgctgGCGCTGGCGCTGGATGCCTTATGGTCACGCAATGGCAAGCTGGCAACGCt                                                                                             >  3:100150‑M2/99‑148 (MQ=255)                                    agctgttttcggcggtggaaacccttagccatatgcgcccgctacgtggcgaccggctgatgattatcagcaacggtgctgcgcctgccgcgctgGCGCTGGCGCTGGATGCCTTATGGTCACGCAATGGCAAGCTGGCAACGCTAAGc                                                                                         >  5:10983‑M2/96‑149 (MQ=255)                                              ggcggtggaaacccttagccatatgcgcccgctacgtggcgaccggctgatgattatcagcaacggtgctgcgcctgccgcgctgGCGCTGGCGCTGGATGCCTTATGGTCACGCAATGGCAAGCTGGCAACGCTAAGCGAAGAAAcc                                                                                >  2:69659‑M2/86‑148 (MQ=255)                                              ggcggtggaaacccttagccatatgcgcccgctacgtggcgaccggctgatgattatcagcaacggtgctgcgcctgccgcgctgGCGCTGGCGCTGGATGCCTTATGGTCACGCAATGGCAAGCTGGCAACGCTAAGCGAAGAAAcc                                                                                <  1:69659‑M2/63‑1 (MQ=255)                                                cggtggaaacccttagccatatgcgcccgccacgtggcgaccggctgatgattatcagcaacggtgctgcgcctgccgcgctgGCGCTGGCGCTGGATGCCTTATGGTCACGCAATGGCAAGCTGGCAACGCTAAGCGAAGAAACCTGc                                                                             <  1:76542‑M2/66‑1 (MQ=255)                                                  gtggaaacccttagccatatgcgcccgctacgtggcgaccggctgatgattatcagcaacggtgctgcgcctgccgcgctgGCGCTGGCGCTGGATGCCTTATGGTCACGCAATGGCAAGCTGGCAACGCTAAGCGAAGAAACCTGCCa                                                                           <  6:10983‑M2/68‑1 (MQ=255)                                                         cccttagccatatgcgcccgctacgtggcgaccggctgatgattatcagcaacggtgctgcgcctgccgcgctgGCGCTGGCGCTGGATGCCTTATGGTCACGCAATGGCAAGCTGGCAACGCTAAGCGAAGAAACCTGCCAGAAACTg                                                                    <  4:100150‑M2/75‑1 (MQ=255)                                                            ttagccatatgcgcccgctacgtggcgaccggctgatgattatcagcaacggtgctgcgcctgccgcgctgGCGCTGGCGCTGGATGCCTTATGGTCACGCAATGGCAAGCTGGCAACGCTAAGCGAAGAAACCTGCCCGAAACTgcgc                                                                 >  7:76085‑M2/72‑149 (MQ=255)                                                              agccatatgcgcccgctacgtggcgaccggctgatgattatcagcaacggtgctgcgcctgccgcgctgGCGCTGGCGCTGGATGCCTTATGGTCACGCAATGGCAAGCTGGCAACGCTAAGCGAAGAAACCTGCCAGAAACTgcgc                                                                 <  2:123571‑M2/78‑1 (MQ=255)                                                                            gctacgtggcgaccggctgatgattatcagcaacggtgctgcgcctgccgcgctgGCGCTGGCGCTGGATGCCTTATGGTCACGCAATGGCAAGCTGGCAACGCTAAGCGAAGAAACCTGCCAGAAACTGCGCGATGCACTGCCAGaa                                                  >  1:1376‑M2/56‑148 (MQ=255)                                                                                                   tgatgaTTATCAGCAACGGTGCTGCGCCTGCCGCGCTGGCGCTGGATGCCTTATGGTCACGCAATGGCAAGCTGGCAACGCTAAGCGAAGAAACCTGCCAGAAACTGCGCGATGCACTGCCAGAACATGTGGCAATATCTAACCCGCTc                          >  7:28277/1‑149 (MQ=255)                                                                                                        agcaacggtgctgcgcctgccgcgctgGCGCTGGCGCTGGATGCCTTATGGTCACGCAATGGCAAGCTGGCAACGCTAAGCGAAGAAACCTGCCAGAAACTGCGCGATGCACTGCCAGAACATGTGGCAATATCTAACCCGCTCGATc                      <  2:1376‑M2/121‑1 (MQ=255)                                                                                                          caacggtgctgcgcctgccgcgctgGCGCTGGCGCTGGATGCCTTATAGTCACGCAATGGCAAGCTGGCAACGCTAAGCGAAGAAACCTGCCAGAAACTGCGCGATGCACTGCCAGAACATGTGGCAATATCTAACCCGCTCGATCTAc                   >  8:10524‑M2/26‑149 (MQ=255)                                                                                                                      gcctgccgcgctgGCGCTGGCGCTGGATGCCTTATGGTCACGCAATGGCAAGCTGGCAACGCTAAGCGAAGAAACCTGCCAGAAACTGCGCGATGCACTGCCAGAACATGTGGCAATATCTAACCCGCTCGATCTACGCGATGACGCca       >  4:62866‑M2/14‑149 (MQ=255)                                                                                                                             gcgctgGCGCTGGCGCTGGATGCCTTATGGTCACGCAATGGCAAGCTGGCAACGCTAAGCGAAGAAACCTGCCAGAAACTGCGCGATGCACTGCCAGAACATGTGGCAATATCTAACCCGCTCGATCTACGCGATGACGCCAGCAGt  >  4:60129‑M2/7‑147 (MQ=255)                                                                                                                             gcgctgGCGCTGGCGCTGGATGCCTTATGGTCACGCAATGGCAAGCTGGCAACGCTAAGCGAAGAAACCTGCCAGAAACTGCGCGATGCACTGCCAGAACATGTGGCAATATCTAACCCGCTCGATCTACGCGATGACGCCAGCAGt  <  3:60129‑M2/141‑1 (MQ=255)                                                                                                                                   |                                                                                                                                             TATTCAGCGTGCCGGTTTGTTGCGGGTACAGGACACCCACGAGCTGTTTTCGGCGGTGGAAACCCTTAGCCATATGCGCCCGCTACGTGGCGACCGGCTGATGATTATCAGCAACGGTGCTGCGCCTGCCGCGCTGGCGCTGGATGCCTTATGGTCACGCAATGGCAAGCTGGCAACGCTAAGCGAAGAAACCTGCCAGAAACTGCGCGATGCACTGCCAGAACATGTGGCAATATCTAACCCGCTCGATCTACGCGATGACGCCAGCAGT  >  NC_000913/2720726‑2720996

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