Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	NC_000913_3_hsa_tpiA	1,726,433	Δ11 bp	33.3%	coding (411‑421/600 nt)	nemR →	transcriptional repressor for the nemRA‑gloA operon, quinone‑, glyoxal‑, and HOCl‑activated

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_hsa_tpiA	= 1726432	16 (0.460)	8 (0.230)	8/276	NT	33.3%	coding (410/600 nt)	nemR	transcriptional repressor for the nemRA‑gloA operon, quinone‑, glyoxal‑, and HOCl‑activated
	?	NC_000913_3_hsa_tpiA	1726444 =	16 (0.460)					coding (422/600 nt)	nemR	transcriptional repressor for the nemRA‑gloA operon, quinone‑, glyoxal‑, and HOCl‑activated

CCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGC  >  NC_000913_3_hsa_tpiA/1726292‑1726579                                                                                                                                             |                                                                                                                                                    ccAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCACGGCGTGATCGCCCTGCTCTCTCAgg                                                                                                                                             >  1:496300/1‑149 (MQ=255)          acTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTtggaaaatggccgtgag                                                                                                                                    >  1:553109‑M1/1‑132 (MQ=255)                 cAGTTTTGCCAACATTGAACCCTTAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTtggaaaatggccgtgagaaccatt                                                                                                                             <  2:85759‑M1/149‑25 (MQ=255)                              aTGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTAAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGACCCATGGt                                                                                                                >  2:42843/1‑149 (MQ=255)                                      aTCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGACCCATTCTTTAACCtt                                                                                                        >  2:337791/1‑149 (MQ=255)                                           tGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTAAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTtgtg                                                                                                   >  1:191840/1‑149 (MQ=255)                                                       aGTAAAATTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGGAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCGCCCAGGCGCTGGAAAAGGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGc                                                                                       <  1:490658/149‑1 (MQ=255)                                                       aGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTtggaaaatggccgtgagaaccattgtttaaccttttgtggcgaaccgctgcaacaggcacaa                                                                                       >  1:636481‑M1/1‑87 (MQ=255)                                                        gTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAAAAGATATGCGGAGGGCGAGCGATAAAGGTGCTCGCGGCGTGATCGGCCTGCTCTCTCAGGCGGTGGAAAATGGCCTTGCGAACCCTTGTATAACCTTTTGCGGCCAACCATTGCa                                                                                      >  2:323031/1‑149 (MQ=255)                                                         tAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTtggaacatggccgtgagaaccattgtttaaccttttttggcgaacctctgcaacaggcacaggt                                                                                     >  1:448227‑M1/1‑85 (MQ=255)                                                            aaCTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCCGGCGACGGGAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAg                                                                                  <  1:517810/149‑1 (MQ=255)                                                                 tctGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGcaca                                                                             >  1:219105/1‑149 (MQ=255)                                                                 tctGCCGAAGTGTGCGAACTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGcaca                                                                             >  1:162055/1‑149 (MQ=255)                                                                    gCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTtggaaaatggccgtgagaaccattgtttaaccttttgtggcgaaccgctgcaacaggcacaagtgctttacgcac                                                                          <  1:610621‑M1/149‑76 (MQ=255)                                                                     ccGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGATGCTCGCGGCGTGAGCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACAGCTGCAGCAGGCACAAggg                                                                         >  1:580481/1‑147 (MQ=255)                                                                                      tcattagaTATTCGCAGCGCGATGGATAAAGGTGCGCGCAGCGTGGTCGCCCTGCTCTCGCAGGCGCTGGAAAATGGCCGTGGGAACCATTGTTTAACCTTCAGTAGCCAACCGCGGCAACAGGCACAAGTGCTTTACGCCCTGTGGCt                                                        >  2:300454/6‑149 (MQ=255)                                                                                         gaagaTATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTggg                                                     >  1:46147/1‑149 (MQ=255)                                                                                                  cgcAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTtggaaaatggccgtgagaaccattgtttaaccttttgtggcgaaccgctgcaacaggcacaagtgctttacgcactgtggcttggcgcgaatctgcaggccaaaa                                            <  1:84395‑M1/149‑106 (MQ=255)                                                                                                            tGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCCCGAATCTGCAGGCCaaa                                  >  1:213188/1‑149 (MQ=255)                                                                                                                  aaGGTGTTCGCGGCGTGATCGCCCTGCTtggaaaatggccgtgagaaccattgtttaaccttttgtggcgaaccgctgcaacaggcacaagtgctttacgcactgtggcttggcgcgaatctgcaggccaaaatttcgcgcagtttcga                            <  1:143320‑M1/149‑122 (MQ=255)                                                                                                                  aaGGTGCTCGCGGCGTGGTCGCCCTGCTCTCTCAGCTGCTGCAAAACGGCCGAGAGCAGGATTGTTTAACCTTTTGTGGCGAACCGCAGCAACCGGCGCAAGTCCTTTAATCACTTAGGATTGGCGCGAATCTGCAGGCCAACATTTcg                            >  2:436574/1‑149 (MQ=255)                                                                                                                  aaGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTcg                            >  1:255992/1‑149 (MQ=255)                                                                                                                         tCGCGGCGTGATCGCCCTGCTtggaaaatggccgtgagaaccattgtttaaccttttgtggcgaaccgctgcaacaggcacaagtgctttacgcactgtggcttggcgcgaatctgcaggccaaaatttcgcgcagtttcgagccactg                     >  1:498757‑M1/1‑21 (MQ=255)                                                                                                                             ggcggcaTCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGa                 <  1:2749/143‑1 (MQ=255)                                                                                                                                            ctctctCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCAAATCTGCAGGCCAAAATTTTGCGCAGTTTCGAGCCACTGGAAAAcac  >  2:62798/1‑147 (MQ=255)                                                                                                                                             |                                                                                                                                                    CCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGC  >  NC_000913_3_hsa_tpiA/1726292‑1726579

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