Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 1,871,057 IS1 (+) +9 bp coding (673‑681/1476 nt) cdgI → putative c‑di‑GMP binding protein CdgI

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1871057 =0 (0.000)15 (0.570) 14/490 0.7 100% coding (673/1476 nt) cdgI putative c‑di‑GMP binding protein CdgI
?NC_000913 1978503 = NA (NA)noncoding (768/768 nt) IS1 repeat region
* ? NC_000913 = 18710650 (0.000)20 (0.760) 18/490 0.4 100% coding (681/1476 nt) cdgI putative c‑di‑GMP binding protein CdgI
?NC_000913 = 1979270 NA (NA)noncoding (1/768 nt) IS1 repeat region

TGATCCAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAACTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAGTATAATTATACTATCTGGTATATCAGTCGCGGGATTGAAGTTTCCAGTAAACT  >  NC_000913/1870836‑1871304
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
cggatggctggccgctgtatgaatcccgcctgaagggaaagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccGTCGTCGATTACTGTCCGGACAACCATACt                                                                                                                                                                                                                            >  1:46196‑M1/222‑251 (MQ=255)
                    gaatcccgcctgaagggaaagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACt                                                                                                                                                                                                         >  1:40842‑M1/202‑250 (MQ=255)
                        cccgcctgaagggaaagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATAc                                                                                                                                                                                                     >  1:238356‑M1/198‑250 (MQ=255)
                                 agggaaagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTg                                                                                                                                                                                           >  2:171647‑M1/189‑251 (MQ=255)
                                           gcccgtaatcagcaagcgatctacgcagcgccttgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccGTCGTCGATTACTGTCCGGCCAACCATACTGCGACCTGGGGAATCCACTATACAAAAATCTTGGTTtgtcggt                                                                                                                                                                                 <  2:189809‑M1/73‑4 (MQ=255)
                                                    cAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATggtgatgctgccaacttactgatttagtgtatgat                                                                                                                                                                                                              >  1:208514‑M1/1‑178 (MQ=255)
                                                    cAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATggtgatgctgccaacttactgatttagtgtatgat                                                                                                                                                                                                              <  2:208514‑M1/213‑36 (MQ=255)
                                                        aagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATCGGTTTGTCTGTGGGCATTTTTACt                                                                                                                                                                    >  2:144206‑M1/166‑251 (MQ=255)
                                                               gTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctc                                                                                                                                                             >  1:247359‑M1/1‑167 (MQ=255)
                                                                  tAACCTCGCTGGCGCTCTTTTGTTAAGGCAAAGCCAACATCCTTGCCACCAATAAGAAACCAACGGGACTCATGTTACTGACGATGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGAAGTTATAACGCTGACTCTTCTTTGTATGTCGTCGATggtgatgctgccaacttactgatttcgtgtatgatggcgtttttgaggtgctccagtggcttctgtttctatcagctgtcactactg                                                                                                                                                          >  2:265944‑M1/1‑164 (MQ=255)
                                                                                       acctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGcc                                                                                                                                     <  2:224510‑M1/117‑1 (MQ=255)
                                                                                                tgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACt                                                                                                                             <  2:58119‑M1/125‑1 (MQ=255)
                                                                                                cagACAACATCCTTGCCAACAATTAGAAACAACCGGGACTCCTGTTGCTGGCGCTGCTCCCTTTTTTAGTTTTTCCCCTTCTGGCACCCAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATggttatgctgccaacttactgatttagtgtatgatggtgtttttgtggtgctcccgtggcttctgtttctattagctgtccctcctgttcagctactgacggggtggtgcttaacgg                                                                                                                            >  2:15779‑M1/2‑134 (MQ=255)
                                                                                                         acctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTa                                                                                                                   >  1:58119‑M1/117‑251 (MQ=255)
                                                                                                                     tgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGACTCAACTATACAAAAATATTGGTTTGTCTGTGAGCATTTTTACTGTTATTTATAATCATGCGATCACGTTAACCCAGCGAACTATGGCCGTTATTAGcattattg                                                                                                       <  2:54705‑M1/147‑2 (MQ=255)
                                                                                                                          aaaaaaaCGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATggtgattctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgct                                                                                                  <  1:129941‑M1/251‑144 (MQ=255)
                                                                                                                           aaaaaaCGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgct                                                                                                  >  2:75997‑M1/1‑107 (MQ=255)
                                                                                                                                  gctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGc                                                                                          <  1:210666‑M1/160‑1 (MQ=255)
                                                                                                                                      aTCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACACTCTGAGCAGTTATCACGCTGACTATTCTTTGTATGTCGTCGATggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctc                                                                                      >  1:236095‑M1/1‑96 (MQ=255)
                                                                                                                                               cTGGAGCTGATACCTTTTTTAGTTTTTCACCTTCTGGCACACCCTCTGAGCAGTTATAACGCTGCCTATTCTTTGTATGTCGTCGCTggtgatgctgcacacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaacagcaccgccggccatcagcgctatctctgctctcactgccgt                                                                              >  2:169101‑M1/1‑87 (MQ=255)
                                                                                                                                                      tGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatg                                                                      >  1:9078‑M1/1‑80 (MQ=255)
                                                                                                                                                             ttttttAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATggcgatgctgccaacttactgatttagtgtatgatggtgtttttgaggttctccagtggcttctgtttctatcagctgtccctcctgttaagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacctggcaactg                                                               >  1:229265‑M1/1‑73 (MQ=255)
                                                                                                                                                                tttAGTTTTTCCCGTTCTGGCAAACAACCTGAGCAGTTATAACGCTGTCTTTTCTTTGTATGTAGTCGATggtgatgctgccaacttactgattaagtgtatgatggagtttttgaggtgctccagtggcttctgtttctatcagctgtcactcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgc                                                              >  2:117535‑M1/1‑70 (MQ=255)
                                                                                                                                                                tttAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgc                                                              >  1:225682‑M1/1‑70 (MQ=255)
                                                                                                                                                                           tcgggcattatctgaacataaaacactatcaataagttggagtcattaccGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAGTATAAt                                                 <  1:217377‑M1/201‑1 (MQ=255)
                                                                                                                                                                                   gCACACAATCGGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctc                                         <  2:9078‑M1/251‑201 (MQ=255)
                                                                                                                                                                                             tGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaacccggtac                               >  2:41625‑M1/1‑41 (MQ=255)
                                                                                                                                                                                              gAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaacccggtacg                              <  1:36400‑M1/251‑212 (MQ=255)
                                                                                                                                                                                                    ttATAACGCTGACTATTCTTTGTATGTCGTCGATggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaacccggtacgcacc                          >  2:222244‑M1/1‑34 (MQ=255)
                                                                                                                                                                                                        caataagttggagtcattaccGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAGTATAATTATACTATCTGGTATATCAGTCGCGGGAt                    <  1:224513‑M1/230‑1 (MQ=255)
                                                                                                                                                                                                        caataagttggagtcattaccGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACGCTGGATGAGTATAATTATACTATCTGGTATATCAGTCGCGGGa                     >  1:62811‑M1/22‑250 (MQ=255)
                                                                                                                                                                                                              gACTATTCTTTGTATGTCGTCGATggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctagcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaacccggtacgcaccagaaaatcattg              <  2:247359‑M1/251‑228 (MQ=255)
                                                                                                                                                                                                                      cattaccGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAGTATAATTATACTATCTGGTATATCAGTCGCGGGATTGAAGTTTCCAGt       >  1:250963‑M1/8‑250 (MQ=255)
                                                                                                                                                                                                                       ttGTATGTCGTCGATggtgatgctgccaacttattgatttagtgtatgatggtgtttttgaggggctccagtggcttcgttttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaagccgatacgcaccagaaaagcattgatatggc       <  2:236095‑M1/249‑235 (MQ=255)
                                                                                                                                                                                                                           aTGTCGTCGATggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaccatggcaactgcagttcacttacaccgcttctcaacccggtacgcaccagaaaatcattgatatggccatga  >  2:258029‑M1/1‑11 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
TGATCCAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAACTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAGTATAATTATACTATCTGGTATATCAGTCGCGGGATTGAAGTTTCCAGTAAACT  >  NC_000913/1870836‑1871304

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