Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 366,181 Δ93 bp coding (33‑125/3075 nt) lacZ ← beta‑galactosidase

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 366181 366273 93 18 [0] [0] 19 lacZ beta‑galactosidase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3661800 (0.000)18 (0.700) 17/486 0.5 100% coding (126/3075 nt) lacZ beta‑galactosidase
?NC_000913 366274 = 0 (0.000)coding (32/3075 nt) lacZ beta‑galactosidase

GAGTAACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGC  >  NC_000913/365932‑366403
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                               
gAGTAACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCac                                                                                                                                                                                                                               >  1:222464‑M1/1‑249 (MQ=255)
   tAACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgac                                                                                                                                                                                                                            >  1:198521‑M1/1‑246 (MQ=255)
           cGTCGGATTCTCCGGGGGAACAACCGGGGGATTGACCGTAAAGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTACCCTTGCATCGGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtg                                                                                                                                                                                                                    <  2:198521‑M1/251‑14 (MQ=255)
               ggATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagggaatc                                                                                                                                                                                                                >  1:165156‑M1/1‑234 (MQ=255)
                                                          aGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatcc                                                                                                                                                                     <  1:124669‑M1/251‑61 (MQ=255)
                                                                                   aTCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacg                                                                                                                                             >  1:249463‑M1/1‑166 (MQ=255)
                                                                                                     ccAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcctacagtgt                                                                                                                          >  2:232364‑M1/1‑148 (MQ=255)
                                                                                                            gAGGGGACGACGACAGTAGCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgcttacaattccacccaacatacgagccggaagcataaagtgtaaagcc                                                                                                                    <  1:107805‑M1/250‑110 (MQ=255)
                                                                                                             aGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTACGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCCGCGCAACTGTTGGGAAGGGCGATCAGTGCGGGCCacgacggccagtgaatcctgaatcatggtcatagctgtatcctgtgtgcaattgttatccgctcaccattccacacaaaataagagccggaagaataaagcgtaaagcctg                                                                                                                  >  1:197035‑M1/1‑140 (MQ=255)
                                                                                                                   cgacgacAGAATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCTCTTCTGGTGCCGGAAACGAGGCAAAGCGCCATTCGCCTTTCATGCGGCGCAGCTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcttagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatccgagccggaagcataaagtgtaaagcctggggtgc                                                                                                            <  1:41428‑M1/251‑118 (MQ=255)
                                                                                                                                                 aCTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaacgcctggggtgcctaatgagtgagctaactcacattaattgc                                                                              >  1:210769‑M1/1‑104 (MQ=255)
                                                                                                                                                                       cGCTTCTCGAGCCGAAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCGACTGTTGGGAAGGGCGATCGGTGCGGGCAacgacggccagtgaatccgcaatcatggtcatagctgtttcctgtgtgaaattgttatcagctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggtgtgcctaatgagtgagctaaatcacattaatttcgttgtgctcactgcccgcattc                                                        >  2:124202‑M1/1‑82 (MQ=255)
                                                                                                                                                                                       aaCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgt                                        <  1:224945‑M1/251‑186 (MQ=255)
                                                                                                                                                                                                  gctcctttcgccatcccggcggCGCACCTGTTGGGAAGGGCGATCGTTTCTGGCCacgtcggccagtgaatccgtaatcattgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccactcaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggatacctgtcgtgccagctg                             <  2:90329‑M1/245‑197 (MQ=255)
                                                                                                                                                                                                           cccattcaGGCTGCGCAACTGTTGGGAAGGGCGATCGGGGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatga                    <  2:222464‑M1/250‑206 (MQ=255)
                                                                                                                                                                                                                      ggcgcaACTGGGGGGAAGGGCGATCGGGGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaatcggccaacg         <  2:57519‑M1/250‑217 (MQ=255)
                                                                                                                                                                                                                      cgcgatactggtgGGAAGGGCGATCGGTGCGGTCCactacggccagtgaatccgtaatcattgtcatagctggttcctgtgtgaaattgtttaccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctatctcacattaattgcgttgcgcgcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaagcggccaacg         <  2:151611‑M1/245‑217 (MQ=255)
                                                                                                                                                                                                                              tGTTGGGAAGGGGGGTTGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaatcggccaacgcgcgggg  <  2:76994‑M1/250‑224 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                               
GAGTAACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGC  >  NC_000913/365932‑366403

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