Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,581,699 (19-bp)2→1 intergenic (+60/+17) tktB → / ← ypfG transketolase 2/DUF1176 domain‑containing protein YpfG

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2581699 2581717 19 31 [0] [0] 31 tktB/ypfG transketolase 2/DUF1176 domain‑containing protein YpfG

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 2581698NA (NA)27 (0.430) 23/238 1.2 NA noncoding (36/36 nt) other REP179b
?NC_000913 2581718 = NA (NA)intergenic (+79/+16) tktB/ypfG transketolase 2/DUF1176 domain‑containing protein YpfG

TCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGT  >  NC_000913/2581566‑2581831
                                                                                                                                    |                                                                                                                                     
tCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAtgtcccg                                                                                                                                >  2:552753‑M1/1‑133 (MQ=255)
           ggttttaccGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAtgtcccgccattttgttt                                                                                                                     <  1:625651‑M1/137‑19 (MQ=255)
               ttACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAtgtcccgccattttgtttaacg                                                                                                                 >  2:667074‑M1/1‑118 (MQ=255)
                        agaATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTACCCGGCCTACAtgtcccgccattttgtttaacgggtgatcca                                                                                                        <  2:1097523‑M1/140‑32 (MQ=255)
                            tatTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAtgtcccgccattttgtttaacgggtgatccacaac                                                                                                    <  1:535219‑M1/140‑36 (MQ=255)
                             atTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAtgtcccgccattttgtttaacgggtgatccacaacg                                                                                                   <  1:673302‑M1/140‑37 (MQ=255)
                                  ggCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAtgtcccgccattttgtttaacgggtgatccacaacgtgggc                                                                                              >  1:644217‑M1/1‑99 (MQ=255)
                                   gCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAtgtcccgccattttgtttaacgggtgatccacaacgtgggcc                                                                                             >  2:552211‑M1/1‑98 (MQ=255)
                                    cAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAtgtcccgccattttgtttaacgggtgatccacaacgtgggcca                                                                                            >  1:1104711‑M1/1‑97 (MQ=255)
                                    cAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAtgtcccgccattttgtttaacgggtgatccacaacgtgggcca                                                                                            >  1:264940‑M1/1‑97 (MQ=255)
                                    cAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAtgtcccgccattttgtttaacgggtgatccacaacgtgggcca                                                                                            >  2:928297‑M1/1‑97 (MQ=255)
                                    cAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAtgtcccgccattttgtttaacgggtgatccacaacgtgggcca                                                                                            >  1:1209323‑M1/1‑97 (MQ=255)
                                    aaaaagCGCATAGGGTCCTGGGAGTGAGGGGTGCCTGGTGGTGATTGCCGGATGTGGATTGCCGGATGGGGCGCGGACGCGTCTTATCCGGCTTACAtggcccgccattttgtttaacgggtgatccacaacgtgggcca                                                                                            <  2:398455‑M1/139‑44 (MQ=255)
                                                    gCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAtgtcccgccattttgtttaacgggtgatccacaacgtgggccaggcgtctggcccatgc                                                                            <  2:80894‑M1/140‑60 (MQ=255)
                                                           gTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAtgtcccgccattttgtttaacgggtgatccacaacgtgggccaggcgtctggcccatgccagttat                                                                     >  2:614897‑M1/1‑74 (MQ=255)
                                                           gTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAtgtcccgccattttgtttaacgggtgatccacaacgtgggccaggcgtctggcccatgccagttat                                                                     <  2:1209323‑M1/140‑67 (MQ=255)
                                                                  gTGCCTGATGGTGATTGCGGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAtgtcccgccattttgtttaacgggtgatccacaacgtgggccaggcgtctggcccatgccagttatcgcaggt                                                              <  2:16484‑M1/140‑74 (MQ=255)
                                                                     ccTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAtgtccctccattttgtttaacgggtgatccacaacgtgggccaggcgtctggcccatgccagttatcgcaggtggg                                                           >  1:1205142‑M1/1‑64 (MQ=255)
                                                                        gATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAtgtcccgccattttgtttaacgggtgatccacaacgtgggccaggcgtctggcccatgccagttatcgcaggtgggttc                                                        >  2:702644‑M1/1‑61 (MQ=255)
                                                                         aTGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAtgtcccgccattttgtttaacgggtgatccacaacgtgggccaggcgtctggcccatgccagttatcgcaggtgggttct                                                       >  1:287239‑M1/1‑60 (MQ=255)
                                                                          tGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAtgtcccgccattttgtttaacgggtgatccacaacgtgggccaggcgtctggcccatgccagttatcgcaggtgggttctg                                                      <  1:614897‑M1/140‑82 (MQ=255)
                                                                            gTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAtttcccgccattttgtttaacgggtgatccacaacgtgggccaggcgtctggcccatgccagttatcgcaggtgggttctgcg                                                    <  2:451827‑M1/140‑84 (MQ=255)
                                                                            gTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAtgtcccgccattttgtttaacgggtgatccacaacgtgggccaggcgtctggcccatgccagttatcgcaggtgggttctgcg                                                    <  1:79094‑M1/140‑84 (MQ=255)
                                                                                                     atgcGACGCTGACGCGTCTTATCCGGCCTACAtgtcccgccattttgtttaacgggtgatccacaacgtgggccaggcgtctggcccatgccagttatcgcaggtgggttctgcggcgtaacgcaccaggcgaaaacgct                           >  1:940212‑M1/1‑32 (MQ=255)
                                                                                                       gcGACGCTGACGCGTCTTATCCGGCCTACAtgtccctccattttgtttaacgggtgatccacaacgtgggccaggcgtctggcccatgccagttatcgcaggtgggttctgcggcgtaacgcaccaggcgaaaacgctga                         <  2:1205142‑M1/140‑111 (MQ=255)
                                                                                                          aCGCTGACGCGTCTTATCCGGCCTACAtgtcccgccattttgtttaacgggtgatccacaacgtgggccaggcgtctggcccatgccagttatcgcaggtgggttctgcggcgtaacgcaccaggcgaaaacgctgaccg                      >  2:383321‑M1/1‑27 (MQ=255)
                                                                                                                aCGCGTCTTATCCGGCCTACAtgtcccgccattttgtttaacgggtgatccacaacgtgggccaggcgtctggcccatgccagttatcgcaggtgggttctgcggcgtaacgcaccaggcgaaaacgctgaccgtcaaag                >  1:185521‑M1/1‑21 (MQ=255)
                                                                                                                            cGGCCTACAtgtcccgccattttgtttaacgggtgatccacaacgtgggccaggcgtctggcccatgccagttatcgcaggtgggttctgcggcgtaacgcaccaggcgaaaacgctgaccgtcaaagcgccagcgcgcc    >  2:679378‑M1/1‑9 (MQ=255)
                                                                                                                              gCCTACAtgtcccgccattttgtttaacgggtgatccacaacgtgggccaggcgtctggcccatgccagttatcgcaggtgggttctgcggcgtaacgcaccaggcgaaaacgctgaccgtcaaagcgccagcgcgcctg  >  2:1180106‑M1/1‑7 (MQ=255)
                                                                                                                              gCCTACAtgtcccgccattttgtttaacgggtgatccacaacgtgggccaggcgtctggcccatgccagttatcgcaggtgggttctgcggcgtaacgcaccaggcgaaaacgctgaccgtcaaagcgccagcgcgcctg  >  2:568033‑M1/1‑7 (MQ=255)
                                                                                                                              gCCTACAtgtcccgccattttgtttaacgggtgatccacaacgtgggccaggcgtctggcccatgccagttatcgcaggtgggttctgcggcgtaacgcaccaggcgaaaacgctgaccgtcaaagcgccagcgcgcctg  <  1:667074‑M1/140‑134 (MQ=255)
                                                                                                                                    |                                                                                                                                     
TCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGT  >  NC_000913/2581566‑2581831

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