Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	NC_000913	2,722,123	Δ7 bp	coding (2171‑2177/2661 nt)	pka →	protein lysine acetyltransferase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2722122	2 (0.070)	30 (1.060)	28/248	0.2	93.8%	coding (2170/2661 nt)	pka	protein lysine acetyltransferase
	?	NC_000913	2722130 =	2 (0.070)					coding (2178/2661 nt)	pka	protein lysine acetyltransferase

TTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAACGAGAGTCGGCTGGCAGTGCGCCCTTATCCGCATCAGCTGGAAGAATGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGA  >  NC_000913/2721984‑2722263                                                                                                                                                   |                                                                                                                                      tcatcctttgctggcttctggcagtgaatttaccgcgctggatgtcacgctggatatctcgccgtttgaaggcgataacgagagtcggctggcagtgcgcccttatccgcatcagctggaagaatgggtagaattgaaaaacggtgCtt                                                                                                                                     >  8:29892‑M2/147‑149 (MQ=255) tcatcctttgctggcttctggcagtgaatttaccgcgctggatgtcacgctggatatctcgccgtttgaaggcgataacgagagtcggctggcagtgcgcccttatccgcatcagctggaagaatgggtagaattgaaaaacggtgCtt                                                                                                                                     >  2:130351‑M2/147‑149 (MQ=255)      ctttgctggcttctggcagtgaatttaccgcgctggatgtcacgctggatatctcgccgtttgaaggcgataacgagagtcggctggcagtgcgcccttatccgcatcagctggaagaatgggtagaattgaaaaacggtgCTTGTTcc                                                                                                                                >  6:52834‑M2/142‑149 (MQ=255)                     gcagtgaatttaccgcgctggatgtcacgctggatatctcgccgtttgaaggcgataacgagagtcggctggcagtgcgcccttatccgcatcagctggaagaatgggtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGcc                                                                                                                   <  1:130351‑M2/21‑1 (MQ=255)                       agtgaatttaccgcgctggatgtcacgctggatatctcgccgtttgaaggcgataacgagagtcggctggcagtgcgcccttatccgcatcagctggaagaatgggtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGCCagaa                                                                                                               <  8:18199‑M2/25‑1 (MQ=255)                            atttaccgcgctggatgtcacgctggatatctcgccgtttgaaggcgataacgagagtcggctggcagtgcgcccttatccgcatcagctggaagaatgggtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGAt                                                                                                            >  2:45354‑M2/120‑147 (MQ=255)                                   gcgctggatgtcacgctggatatctcgccgtttgaaggcgataacgagagtcggctggcagtgcgcccttatccgcatcagctggaagaatgggtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACaa                                                                                                   <  5:52834‑M2/37‑1 (MQ=255)                                         gatgtcacgctggatatctcgccgtttgaaggcgataacgagagtcggctggcagtgcgcccttatccgcatcagctggaagaatgggtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCa                                                                                              <  1:45354‑M2/42‑1 (MQ=255)                                           tgtcacgctggatatctcgccgtttgaaggcgataacgagagtcggctggcagtgcgcccttatccgcatcagctggaagaatgggtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGc                                                                                            >  2:379‑M2/105‑148 (MQ=255)                                              cacgctggatatctcgccgtttgaaggcgataacgagagtcggctggcagtgcgcccttatccgcatcagctggaagaatgggtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCa                                                                                              <  8:101957‑M2/42‑1 (MQ=255)                                              cacgctggatatctcgccgtttgaaggcgataacgagagtcggctggcagtgcgcccttatccgcatcagctggaagaatgggtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCa                                                                                              >  7:101957‑M2/102‑143 (MQ=255)                                               acgctggatatctcgccgtttgaaggcgataacgagagtcggctggcagtgcgcccttatccgcatcagctggaagaatgggtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAAtt                                                                                        >  2:56359‑M2/101‑148 (MQ=255)                                                          ctcgccgtttgaaggcgataacgagagtcggctggcagtgcgcccttatccgcatcagctggaagaatgggtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCcac                                                                                                     <  4:70060‑M2/35‑1 (MQ=255)                                                          ctcgccgtttgaaggcgataacgagagtcggctggcagtgcgcccttatccgcatcagctggaagaatgggtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCcac                                                                                                     >  3:70060‑M2/90‑124 (MQ=255)                                                                     aaggcgataacgagagtcggctggcagtgcgcccttatccgcatcagctggaagaatgtgtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAaagaa                                                                  <  7:29892‑M2/70‑1 (MQ=255)                                                                      aggcgataacgagagtcggctggcagtgcgcccttatccgcatcagctggaagaatgggtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCAttt                                                                                   >  1:36567‑M2/78‑130 (MQ=255)                                                                      aggcgataacgagagtcggctggcagtgcgcccttatccgcatcagctggaagaatgggtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCAttt                                                                                   <  2:36567‑M2/53‑1 (MQ=255)                                                                                            gcagtgcgcccttatccgcatcagctggaagaatgggtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGc                                          <  1:379‑M2/94‑1 (MQ=255)                                                                                                         tCAGCTGGAAGAATGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAAt                             >  4:30165/1‑149 (MQ=255)                                                                                                         tCAGCTGGAAGAATGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAAt                             >  6:14507/1‑149 (MQ=255)                                                                                                            cgcatcagctggaagaatgggtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTa                          >  3:18704‑M2/40‑149 (MQ=255)                                                                                                                tcagctggaagaatgggtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTAc                         <  2:89018‑M2/111‑1 (MQ=255)                                                                                                                tcagctggaagaatgggtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTAc                         >  1:89018‑M2/36‑146 (MQ=255)                                                                                                                     tggaagaatgggtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTAc                                                      >  3:27494‑M2/31‑112 (MQ=255)                                                                                                                     tggaagaatgggtagaattgaaaaacggtgCTTGTCCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTAc                                                      <  4:27494‑M2/82‑1 (MQ=255)                                                                                                                          gaatgggtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGAtt              >  2:60136‑M2/26‑147 (MQ=255)                                                                                                                          gaatgggtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGAtt              <  1:60136‑M2/122‑1 (MQ=255)                                                                                                                                gtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGc          >  2:89149‑M2/20‑145 (MQ=255)                                                                                                                                gtagaattgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGc          <  1:89149‑M2/126‑1 (MQ=255)                                                                                                                                     attgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGa  >  5:46604‑M2/15‑148 (MQ=255)                                                                                                                                      ttgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCa                      <  3:926‑M2/114‑1 (MQ=255)                                                                                                                                      ttgaaaaacggtgCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCa                      >  4:926‑M2/14‑127 (MQ=255)                                                                                                                                                   |                                                                                                                                      TTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAACGAGAGTCGGCTGGCAGTGCGCCCTTATCCGCATCAGCTGGAAGAATGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGA  >  NC_000913/2721984‑2722263

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.