Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,857,142 7 bp→GATAGTC coding (50‑56/2562 nt) mutS → methyl‑directed mismatch repair protein

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2857142 2857148 7 33 [0] [0] 33 mutS methyl‑directed mismatch repair protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 28571410 (0.000)31 (0.690)
+GATAGTC
22/202 0.8 100% coding (49/2562 nt) mutS methyl‑directed mismatch repair protein
?NC_000913 2857149 = 0 (0.000)coding (57/2562 nt) mutS methyl‑directed mismatch repair protein

ACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCTCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTCG  >  NC_000913/2857038‑2857275
                                                                                                       |                                                                                                                                      
aCGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctg                                                                                                     <  2:915426‑M1/139‑36 (MQ=255)
 cGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGCAAATTTCGACGCCCATACGCCCATGATGCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgc                                                                                                    >  2:106168‑M1/1‑103 (MQ=255)
               aTCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatg                                                                                      >  1:616596‑M1/1‑89 (MQ=255)
                          tttAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatgggtgattttta                                                                           >  2:239717‑M1/1‑78 (MQ=255)
                              atatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatgggtgatttttatgaa                                                                       <  1:239717‑M1/139‑66 (MQ=255)
                              atatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatgggtgatttttatgaa                                                                       <  2:807452‑M1/139‑66 (MQ=255)
                              atatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatgggtgatttttatgaa                                                                       <  2:616596‑M1/139‑66 (MQ=255)
                              atatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatgggtgatttttatgaa                                                                       <  2:491132‑M1/139‑66 (MQ=255)
                               tatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatgggtgatttttatgaac                                                                      <  1:762566‑M1/139‑67 (MQ=255)
                                    gggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatgggtgatttttatgaa                                                                       <  2:54905‑M1/133‑66 (MQ=255)
                                    gggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatgggtgatttttatgaa                                                                       >  1:54905‑M1/1‑68 (MQ=255)
                                    gggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatgggtgatttttatgaa                                                                       <  1:39891‑M1/133‑66 (MQ=255)
                                    gggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTgatagtcgctgaaagcccagcatcacgagatcctgctgttttaccggatgggtgatttttatgaa                                                                       >  2:39891‑M1/1‑68 (MQ=255)
                                                      cATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatgggtgatttttatgaactgttttatgacgacgcaaaacgc                                               >  1:615439‑M1/1‑50 (MQ=255)
                                                      cATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatgggtga                                                                                 <  2:564107‑M1/105‑56 (MQ=255)
                                                      cATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatgggtga                                                                                 >  1:564107‑M1/1‑50 (MQ=255)
                                                          aGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatgggtgatttttatgaactgttttatgacgacgcaaaacgcgcgt                                           <  1:106168‑M1/139‑94 (MQ=255)
                                                            tGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatgggtgatttttatgaactgttttatgacgacgcaaaacgcgcgtcg                                         >  2:350883‑M1/1‑44 (MQ=255)
                                                              cAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatgggtgatttttatgaactgttttatgacgacgcaaaacgcgcgtcgca                                       >  2:149253‑M1/1‑42 (MQ=255)
                                                                      aaTTTCGACGCCCATACGCCCATGATGCAGCAGTgatagtcgctgaaagcacagcatcccgagatcctgctgttttaccggatgggtgatttttatgaactgttttatgacgacgcaaaacgcgcgtcgcaactgctgg                               >  2:395227‑M1/1‑34 (MQ=255)
                                                                         ttCGACGCCCATACGCCCATGATGCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatgggtgatttttatgaactgttttatgacgacgcaaaacgcgcgtcgcaactgctggata                            <  2:560166‑M1/139‑109 (MQ=255)
                                                                         ttCGACGCCCATACGCCCATGATGCAGCAGTgatagtc                                                                                                                                 <  2:831467‑M1/38‑8 (MQ=255)
                                                                         ttCGACGCCCATACGCCCATGATGCAGCAGTgatagtc                                                                                                                                 >  1:831467‑M1/1‑31 (MQ=255)
                                                                              cgcccatacgcccatGATGCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatgggtgatttttatgaactgttttatgacgacgcaaaacgcgcgtcgcaactgctggatatttca                       <  2:956986‑M1/139‑114 (MQ=255)
                                                                                         ccatGATGCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatgggtgatttttatgaactgttttatgacgacgcaaaacgcgcgtcgcaactgctggatatttcactgaccaaacg            <  1:350883‑M1/139‑125 (MQ=255)
                                                                                           atgatgCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatgggtgatttttatgaactgttttatgacgacgcaaaacgcgcgtcgcaactgctggatatttcactgaccaaacgcg          >  2:233351‑M1/1‑13 (MQ=255)
                                                                                           atgatgCAGCAGTgatagtcgctgaaagcccagcatcc                                                                                                               <  1:331502‑M1/38‑26 (MQ=255)
                                                                                           atgatgCAGCAGTgatagtcgctgaaagcccagcatcc                                                                                                               >  2:331502‑M1/1‑13 (MQ=255)
                                                                                            tgatgCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatgggtgatttttatgaac                                                                      >  2:54157‑M1/1‑12 (MQ=255)
                                                                                            tgatgCAGCAGTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatgggtgatttttatgaac                                                                      <  1:54157‑M1/78‑67 (MQ=255)
                                                                                                 cagcagTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatggg                                                                                    >  2:437742‑M1/1‑7 (MQ=255)
                                                                                                 cagcagTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatggg                                                                                    <  1:437742‑M1/59‑53 (MQ=255)
                                                                                                   gcagTgatagtcgctgaaagcccagcatcccgagatcctgctgttttaccggatgggtgatttttatgaactgttttatgacgacgcaaaacgcgcgtcgcaactgctggatatttcactgaccaaacgcggtgcttcg  >  2:476617‑M1/1‑5 (MQ=255)
                                                                                                       |                                                                                                                                      
ACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCTCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTCG  >  NC_000913/2857038‑2857275

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