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‑M2/28‑1 (MQ=255)
 cgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgagtgcaatagcaaatttcgacgcccatacgcccatgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCctgc                                                                                                    >  2:106168‑M2/111‑139 (MQ=255)
               atcacaccccatttaatatcagggaaccggacataaccccatgagtgcaatagaaaatttcgacgcccatacgcccatgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATg                                                                                      >  1:616596‑M2/97‑139 (MQ=255)
                          tttaatatcagggaaccggacataaccccatgagtgcaatagaaaatttcgacgcccatacgcccatgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTa                                                                           >  2:239717‑M2/86‑139 (MQ=255)
                              atatcagggaaccggacataaccccatgagtgcaatagaaaatttcgacgcccatacgcccatgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                                       <  2:491132‑M2/58‑1 (MQ=255)
                              atatcagggaaccggacataaccccatgagtgcaatagaaaatttcgacgcccatacgcccatgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                                       <  1:239717‑M2/58‑1 (MQ=255)
                              atatcagggaaccggacataaccccatgagtgcaatagaaaatttcgacgcccatacgcccatgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                                       <  2:807452‑M2/58‑1 (MQ=255)
                              atatcagggaaccggacataaccccatgagtgcaatagaaaatttcgacgcccatacgcccatgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                                       <  2:616596‑M2/58‑1 (MQ=255)
                               tatcagggaaccggacataaccccatgagtgcaatagaaaatttcgacgcccatacgcccatgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAAc                                                                      <  1:762566‑M2/59‑1 (MQ=255)
                                    gggaaccggacataaccccatgagtgcaatagaaaatttcgacgcccatacgcccatgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                                       <  1:39891‑M2/58‑1 (MQ=255)
                                    gggaaccggacataaccccatgagtgcaatagaaaatttcgacgcccatacgcccatgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                                       >  1:54905‑M2/76‑133 (MQ=255)
                                    gggaaccggacataaccccatgagtgcaatagaaaatttcgacgcccatacgcccatgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                                       <  2:54905‑M2/58‑1 (MQ=255)
                                    gggaaccggacataaccccatgagtgcaatagaaaatttcgacgcccatacgcccatgatgcagcagtgatagtcGCTGAAAGCCCAGCATCACGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                                       >  2:39891‑M2/76‑133 (MQ=255)
                                                      catgagtgcaatagaaaatttcgacgcccatacgcccatgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGa                                                                                 >  1:564107‑M2/58‑105 (MQ=255)
                                                      catgagtgcaatagaaaatttcgacgcccatacgcccatgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGa                                                                                 <  2:564107‑M2/48‑1 (MQ=255)
                                                      catgagtgcaatagaaaatttcgacgcccatacgcccatgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAAcgc                                               >  1:615439‑M2/58‑139 (MQ=255)
                                                          agtgcaatagaaaatttcgacgcccatacgcccatgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGt                                           <  1:106168‑M2/86‑1 (MQ=255)
                                                            tgcaatagaaaatttcgacgcccatacgcccatgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCg                                         >  2:350883‑M2/52‑139 (MQ=255)
                                                              caatagaaaatttcgacgcccatacgcccatgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCa                                       >  2:149253‑M2/50‑139 (MQ=255)
                                                                      aatttcgacgcccatacgcccatgatgcagcagtgatagtcGCTGAAAGCACAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTgg                               >  2:395227‑M2/42‑139 (MQ=255)
                                                                         ttcgacgcccatacgcccatgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGata                            <  2:560166‑M2/101‑1 (MQ=255)
                                                                              cgcccatacgcccatgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCa                       <  2:956986‑M2/106‑1 (MQ=255)
                                                                                         ccatgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCACTGACCAAAcg            <  1:350883‑M2/117‑1 (MQ=255)
                                                                                           atgatgcagcagtgatagtcGCTGAAAGCCCAGCATcc                                                                                                               >  2:331502‑M2/21‑38 (MQ=255)
                                                                                           atgatgcagcagtgatagtcGCTGAAAGCCCAGCATcc                                                                                                               <  1:331502‑M2/18‑1 (MQ=255)
                                                                                           atgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCACTGACCAAAcgcg          >  2:233351‑M2/21‑139 (MQ=255)
                                                                                            tgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAAc                                                                      >  2:54157‑M2/20‑78 (MQ=255)
                                                                                            tgatgcagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAAc                                                                      <  1:54157‑M2/59‑1 (MQ=255)
                                                                                                 cagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATggg                                                                                    >  2:437742‑M2/15‑59 (MQ=255)
                                                                                                 cagcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATggg                                                                                    <  1:437742‑M2/45‑1 (MQ=255)
                                                                                                   gcagtgatagtcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTcg  >  2:476617‑M2/13‑139 (MQ=255)
                                                                                                             tcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATg                                                                         >  2:104003‑M2/3‑58 (MQ=255)
                                                                                                             tcGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATg                                                                         <  1:104003‑M2/56‑1 (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: