Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,398,054 Δ7 bp coding (643‑649/1848 nt) mutL → methyl‑directed mismatch repair protein

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4398054 4398060 7 48 [0] [0] 48 mutL methyl‑directed mismatch repair protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 43980530 (0.000)46 (0.860) 44/292 0.4 100% coding (642/1848 nt) mutL methyl‑directed mismatch repair protein
?NC_000913 4398061 = 0 (0.000)coding (650/1848 nt) mutL methyl‑directed mismatch repair protein

AATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGATTTGCGGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGC  >  NC_000913/4397911‑4398209
                                                                                                                                                      |                                                                                                                                                    
aaaaccgaatttaaccacattgatgagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgg                                                                                                                                                      >  2:462917‑M2/151‑151 (MQ=255)
     cgaatttaaccacattgatgagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCAcc                                                                                                                                                  >  2:640739‑M2/146‑150 (MQ=255)
              ccaaattgatgagttcatccgccgcattgcgctggcgcgtttcgacgccacgatcaaactgtcgcataacggtaaaattgtgcgacagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTg                                                                                                                                        >  1:434798‑M2/137‑151 (MQ=255)
                 cattgatgagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAAc                                                                                                                                     <  2:352765‑M2/18‑1 (MQ=255)
                     gatgagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgagcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAgc                                                                                                                                 >  2:249476‑M2/130‑151 (MQ=255)
                        gagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAgcgc                                                                                                                               <  2:458872‑M2/24‑1 (MQ=255)
                          gatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTgg                                                                                                                            <  1:331825‑M2/27‑1 (MQ=255)
                              atccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGAt                                                                                                                        >  2:353926‑M2/121‑151 (MQ=255)
                                ccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGAtt                                                                                                                       <  1:635369‑M2/32‑1 (MQ=255)
                                ccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTg                                                                                                                      >  1:320382‑M2/119‑151 (MQ=255)
                                 cgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTg                                                                                                                      >  2:89214‑M2/118‑150 (MQ=255)
                                   ccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGaa                                                                                                                    >  1:643024‑M2/116‑150 (MQ=255)
                                    cgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATg                                                                                                                  >  2:104182‑M2/115‑151 (MQ=255)
                                       attgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGc                                                                                                                >  2:518397‑M2/112‑150 (MQ=255)
                                              tggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGc                                                                                                        >  1:587937‑M2/105‑151 (MQ=255)
                                                        cgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGc                                                                                              >  1:186060‑M2/95‑151 (MQ=255)
                                                           cgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTAc                                                                                           >  2:352080‑M2/92‑151 (MQ=255)
                                                           cgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTAc                                                                                           <  1:658140‑M2/60‑1 (MQ=255)
                                                              cacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTAcgcg                                                                                        <  1:89191‑M2/63‑1 (MQ=255)
                                                                cgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGc                                                                                      >  1:547840‑M2/87‑151 (MQ=255)
                                                                        ctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGc                                                                              >  2:570368‑M2/79‑151 (MQ=255)
                                                                            cgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGa                                                                           <  2:821384‑M2/76‑1 (MQ=255)
                                                                             gcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATc                                                                         >  2:319174‑M2/74‑151 (MQ=255)
                                                                                 aacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCaa                                                                      >  1:702756‑M2/70‑150 (MQ=255)
                                                                                 aacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCaa                                                                      >  1:568395‑M2/70‑150 (MQ=255)
                                                                                   cggtaaaattgtgcgtcagtaccgcgcagtgccggaaggggggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAATTc                                                                   <  1:786583‑M2/84‑1 (MQ=255)
                                                                                               gcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGccc                                                        <  1:551659‑M2/95‑1 (MQ=255)
                                                                                               gcgtcagtaccgcgcagtgccggaaggcaggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGccc                                                        <  2:737482‑M2/95‑1 (MQ=255)
                                                                                                   cagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCAc                                                    <  2:527561‑M2/99‑1 (MQ=255)
                                                                                                        ccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGCACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAg                                              <  1:138364‑M2/105‑1 (MQ=255)
                                                                                                        ccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGc                                                >  1:477488‑M2/47‑149 (MQ=255)
                                                                                                         cgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGc                                                >  1:134641‑M2/46‑148 (MQ=255)
                                                                                                         cgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCCTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGa                                             <  1:693617‑M2/106‑1 (MQ=255)
                                                                                                            gcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAAt                                          <  1:322751‑M2/109‑1 (MQ=255)
                                                                                                                   cggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTAt                                   >  2:522475‑M2/36‑151 (MQ=255)
                                                                                                                      aaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGc                                >  2:136087‑M2/33‑151 (MQ=255)
                                                                                                                       aggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCt                               >  2:500011‑M2/32‑151 (MQ=255)
                                                                                                                          cgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACg                            >  1:119886‑M2/29‑151 (MQ=255)
                                                                                                                           gggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGt                           >  1:582878‑M2/28‑151 (MQ=255)
                                                                                                                              caaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGa                         >  1:809826‑M2/25‑150 (MQ=255)
                                                                                                                                aaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACg                      <  2:320405‑M2/129‑1 (MQ=255)
                                                                                                                                  aagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTTAACGGt                    <  1:727840‑M2/131‑1 (MQ=255)
                                                                                                                                   agaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTc                   <  1:605381‑M2/132‑1 (MQ=255)
                                                                                                                                              taggcgcgGCTCCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATgcgc        >  2:151188‑M2/9‑151 (MQ=255)
                                                                                                                                               aggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGcgcg       <  2:643066‑M2/144‑1 (MQ=255)
                                                                                                                                                 gcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGAt     <  2:507570‑M2/146‑1 (MQ=255)
                                                                                                                                                    cgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGc  >  1:212708‑M2/3‑151 (MQ=255)
                                                                                                                                                    cgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGc  >  2:571659‑M2/3‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
AATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGATTTGCGGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGC  >  NC_000913/4397911‑4398209

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