Predicted mutation
evidence seq id position mutation freq annotation gene description
JC AE015451 1,531,732 Δ107 bp 23.3% coding (38‑144/630 nt) PP_1344 ← conserved protein of unknown function

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AE015451 = 15317318 (0.140)8 (0.140) 7/328 NT 23.3% coding (145/630 nt) PP_1344 conserved protein of unknown function
?AE015451 1531839 = 45 (0.810)coding (37/630 nt) PP_1344 conserved protein of unknown function

ACAAAGTACCATCGCGCCAGGACGCCACATGGCAGTGCTCGCGGGCAGCCGGTTGCAACTGGCTTTCCAGCAGGCGCTGCAAGTGCTCCAGGCGCTCGGCCTGGTTCAGCAGCAGGCGCAGCGGGCGCACCTGGCGCAGCAGGGCGGAAGGTGGCTGGGCGGGGGAGGGTTTGTAGGCCATGGGGGTACGCGTGAAGTGACGTGGGCACAGTCTAACATCTCCGTTGCCTGTGCCGGCCTCTTCGCGGGGGAACCCGCTCCTGCCAGGGTTATCGGTATCTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAGCCAAAGACAAGACTGGTTCTCGATAGCTTTCATGTTGCCAAGCCTTGAACTCAGTGAAAATGCCCCTATCTAAAACAGGCCCCTGCCAAAACGCTGTA  >  AE015451/1531521‑1531940
                                                                                                                                                                                                                  |                                                                                                                                                                                                                 
acaaAGTACCATCGCGCCAGGACGCCACATGGCAGTGCTCGCGGGCAGCCGGTTGCAACTGGCTTTCCAGCAGGCGCTGCAAGTGCTCCAGGCGCTCGGCCTGGTTCAGCAGCAGGCGCAGCGGGCGCACCTGGCGCAGCAGGGCGGAAGGTGGCTGGGCGGGGGAGGGTTTGTAGGCCATGGGGGTACGCGTGAAGTGACGTGGGCACAGTCTAACATCTCCGTTGCCTGTGcc                                                                                                                                                                                           <  1:955058/235‑1 (MQ=255)
                             tGGCAGTGCTCGCGGGCAGCCGGTTGCAACTGGCTTTCCAGCAGGCGCTGCAAGTGCTCCAGGCGCTCGGCCTGGTTCAGCAGCAGGCGCAGCGGGCGCACCTGGCGCAGCAGGGCGGAAGGTGGCTGGGCGGGGGAGGGTTTGTAGGCCATGGGGGTACGCGTGAAGTGACGTGGGCACAGgcagccaaagacaagactggttctcgatagctttcatgttgccaagccttga                                                                                                                                                               <  2:366209‑M1/234‑53 (MQ=255)
                                             cAGCCGGTTGCAACTGGCTTTCCAGCAGGCGCTGCAAGTGCTCCAGGCGCTCGGCCTGGTTCAGCAGCAGGCGCAGCGGGCGCACCTGGCGCAGCAGGGCGGAAGGTGGCTGGGCGGGGGAGGGTTTGTAGGCCATGGGGGTACGCGTGAAGTGACGTGGGCACAGTCTAACATCTCCGTTGCCTGTGCCGGCCTCTTCGCGGGGGAACCCGCTCCTGCCAGGGTTAt                                                                                                                                                     <  2:1082749/228‑1 (MQ=255)
                                                                                                                     gcAGCGGGCGCACCTGGCGCAGCAGGGCGGAAGGTGGCTGGGCGGGGGAGGGTTTGTAGGCCATGGGGGTACGCGTGAAGTGACGTGGGCACAGgcagccaaagacaagactggttctcgatagctttcatgttgccaagccttgaactcagtgaaaatg                                                                                                                                                 <  1:358775‑M1/160‑67 (MQ=255)
                                                                                                                     gcAGCGGGCGCACCTGGCGCAGCAGGGCGGAAGGTGGCTGGGCGGGGGAGGGTTTGTAGGCCATGGGGGTACGCGTGAAGTGACGTGGGCACAGgcagccaaagacaagactggttctcgatagctttcatgttgccaagccttgaactcagtgaaaatg                                                                                                                                                 >  2:358772‑M1/1‑94 (MQ=255)
                                                                                                                     gcAGCGGGCGCACCTGGCGCAGCAGGGCGGAAGGTGGCTGGGCGGGGGAGGGTTTGTAGGCCATGGGGGTACGCGTGAAGTGACGTGGGCACAGTCTAACATCTCCGTTGCCTGTGCCGGCCTCTTCGCGGGGGAACCCGCTCCTGCCAGGGTTATCGGTATCTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAGCCAAAGACAAGACTGGTTCTCGATAGCttt                                                                      >  1:183670/1‑235 (MQ=255)
                                                                                                                               atccTGGCGCAGCAGGGCGGAAGGTGGCTGGGCGGGGGAGGGTTTGTAGGCCATGGGGGTACGCGTGAAGTGACGTGGGCACAGTCTAACATCTCCGTTGCCTGTGCCGGCCTCTTCGCGGGGGAACCCGCTCCTGCCAg                                                                                                                                                           <  1:275017/138‑1 (MQ=255)
                                                                                                                                           cgggggggAGGGTGGCGGGGCGGGGGAGGGTTTGTAGGCCATGGGGGTACGCGTGAAGTGACGTGGGCACAGTCTAACATCTCCGTTGCCTGTGCCGGCCTCTTCGCGGGGGAACCCGCTCCTGCCAGGGTTATCGGTATCTTGTAGGAGCgg                                                                                                                                  <  1:162192/151‑1 (MQ=255)
                                                                                                                                               gCGGAAGGTGGCTGGGCGGGGGAGGGTTTGTAGGCCATGGGGGTACGCGTGAAGTGACGTGGGCACAGgcagccaaagacaagactggttctcgatagctttcatgttgccaagccttgaactcagtgaaaatg                                                                                                                                                 <  1:117441‑M1/134‑67 (MQ=255)
                                                                                                                                               gCGGAAGGTGGCTGGGCGGGGGAGGGTTTGTAGGCCATGGGGGTACGCGTGAAGTGACGTGGGCACAGgcagccaaagacaagactggttctcgatagctttcatgttgccaagccttgaactcagtgaaaatg                                                                                                                                                 >  2:117439‑M1/1‑68 (MQ=255)
                                                                                                                                                 ggAAGGTGGCTGGGCGGGGGAGGGTTTGTAGGCCATGGGGGTACGCGTGAAGTGACGTGGGCACAGTCTAACATCTCCGTTGCCTGTGCCGGCctctgc                                                                                                                                                                                  >  2:162190/1‑97 (MQ=255)
                                                                                                                                                   agggtgGCTGGGCGGGGGAGGGTTTGTAGGCCATGGGGGTGCGCGTGAAGTGACGTGGGCACAGgcagccaaagacaagactggttctcgatagctttcatgttgccaagccttgaactcagtgaaaatgcccctatctaaaacaggcccctgccaaaacgctgtatcagcatcgtggaaatccaccactttcctcaccatcgtttccggg                                                                <  2:211335‑M1/209‑148 (MQ=255)
                                                                                                                                                                                     gggggTACGCGTGAAGTGACGTGGGCACAGTCTAACa                                                                                                                                                                                                            >  1:645254/1‑37 (MQ=255)
                                                                                                                                                                                         gTACGCGTGAAGTGACGTGGGCACAGgcagccaaagacaagactggttctcgagagctttcatgttgccaagccttgaactcagtgaaaatgccactatctaaaacaggcccctgccaaaacgctgtagcagcatcgtggaaatccaccactttcctcaccatcgtatccgggtagaatgctcgttcg                                                 >  2:126063‑M1/1‑26 (MQ=255)
                                                                                                                                                                                             gcgTGAAGTGACGTGGGCACAGTCTAACATCTCCGTTGCCTGTGCCGGCCTCTTCGCGGGGGAACCCGCTCCTGCCAGGGTTATCGGTATCTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAGCCAAAGACAAGACTGGTTCTCGATAGCTTTCATGTTGCCAAGCCTTGAAc                                                  >  1:1050735/1‑183 (MQ=255)
                                                                                                                                                                                                  aaGTGACGTGGGCACAGgcagccaaagacaagactggttctcgatagctttcatgttgccaagccttgaactcagtgaaaatgcccctatctaaaacaggcccctgccaaaacgctgtatcagcatcgtggaaatccaccactttcctcaccatcgtttccgggtacaatgctcgttcgcatgcggcctatgggctgcacgggcgactcatggggccgccctccat  <  2:947285‑M1/226‑210 (MQ=255)
                                                                                                                                                                                                                caGgcagccaaagacaagactggttctcgatagctttcatgttgccaagccttgaactcagtgaaaatgcccctatctaaaacaggcccctgccaaaacgctgtatcagcatcgtggaaatccaccactttcctcaccatcgtttccgggtagaatgctcgttcgcatgcggcct                                       <  2:843612‑M1/175‑173 (MQ=255)
                                                                                                                                                                                                                  |                                                                                                                                                                                                                 
ACAAAGTACCATCGCGCCAGGACGCCACATGGCAGTGCTCGCGGGCAGCCGGTTGCAACTGGCTTTCCAGCAGGCGCTGCAAGTGCTCCAGGCGCTCGGCCTGGTTCAGCAGCAGGCGCAGCGGGCGCACCTGGCGCAGCAGGGCGGAAGGTGGCTGGGCGGGGGAGGGTTTGTAGGCCATGGGGGTACGCGTGAAGTGACGTGGGCACAGTCTAACATCTCCGTTGCCTGTGCCGGCCTCTTCGCGGGGGAACCCGCTCCTGCCAGGGTTATCGGTATCTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAGCCAAAGACAAGACTGGTTCTCGATAGCTTTCATGTTGCCAAGCCTTGAACTCAGTGAAAATGCCCCTATCTAAAACAGGCCCCTGCCAAAACGCTGTA  >  AE015451/1531521‑1531940

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