Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,651,191 Δ2,158 bp yhiS[yhiS] yhiS, insH11, [yhiS]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3651191 3653348 2158 33 [0] [0] 33 yhiS–[yhiS] yhiS,insH11,[yhiS]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 36511900 (0.000)33 (0.510) 30/476 0.6 100% intergenic (+528/‑101) arsC/yhiS arsenate reductase/putative uncharacterized protein YhiS
?NC_000913 3653349 = 0 (0.000)pseudogene (114/478 nt) yhiS putative uncharacterized protein YhiS

CAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAAC  >  NC_000913/3650942‑3651420
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      
cAGTGGTTTACACAACAAATTATTAAATAATNATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAac                                                                                                                                                                                                                                      <  2:231170‑M1/251‑3 (MQ=255)
  gTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCGCAAAAAAAATCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAacct                                                                                                                                                                                                                                    >  1:389066‑M1/1‑247 (MQ=255)
                aaaTTATTAAATAATTATGAGATACCTTAACCCATTTTATTATTACAAATCTAAAACATCACAAAAAAAATCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTAATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATCAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAaccttttacgagcaaa                                                                                                                                                                                                                        >  2:418740‑M1/1‑233 (MQ=255)
                     attaAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAaccttttacgagcaaaagcgga                                                                                                                                                                                                                  >  2:416282‑M1/1‑228 (MQ=255)
                      ttaAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAaccttttacaagcaaaagcggaa                                                                                                                                                                                                                 >  2:179929‑M1/1‑227 (MQ=255)
                         aataatTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAaccttttacgagcaaaagcggaaattg                                                                                                                                                                                                             <  1:156895‑M1/251‑28 (MQ=255)
                                                 gttttaTTATTACACATCTAACACTTCACCAAAAAAATCACTTCACATACTATCAGAGAACATAATAACTTTTTAAAATGCGCGCATTTTGAAAGTGGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAATTATTACAGGGTGGAGCTAACAGGGAAATCTAAGGTAATAATATCAGGAAGATAaccttttacgggcaaaagcggaaattgcttcgacaataattaaattttt                                                                                                                                                                                       >  2:395933‑M1/2‑200 (MQ=255)
                                                         attaCACATCTAAAACATCACAAAAAAAATCTCTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAaccttttacgagcaaaagcggaaattgcttcgacaataattaaatttttagaggccagg                                                                                                                                                                             <  1:356255‑M1/251‑60 (MQ=255)
                                                            acacATCTAAAACATCACAAAAAAAATCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAaccttttacgagcaaaagcggaaattgcttcgacaataattaaatttttagaggccaggat                                                                                                                                                                           >  1:258067‑M1/1‑189 (MQ=255)
                                                            acacATCTAAAACATCACAAAAAAAATCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAaccttttacgagcaaaagcggaaattgcttcgacaataattaaatttttagaggccaggat                                                                                                                                                                           <  1:43616‑M1/250‑62 (MQ=255)
                                                             cacaTCTAAAACATCACAAAAAAAATCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAaccttttacgagcaaaagcggaaattgcttcgacaataattaaatttttagaggccaggataa                                                                                                                                                                         <  1:51999‑M1/251‑64 (MQ=255)
                                                               caTCTAAAACATCACAAAAAAAATCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAaccttttacgagcaaaagcggaaattgcttcgccaataattacatttttagaggccaggatcacc                                                                                                                                                                       >  2:19544‑M1/1‑186 (MQ=255)
                                                                      aaCATCACAAAAAAAATCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAaccttttacgagcaaaagcggaaattgcttcgacaataattaaatttttagaggccaggataaccaatcttt                                                                                                                                                                <  1:152969‑M1/251‑73 (MQ=255)
                                                                                       cACTTCCCACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACAACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAaccttttacgagcaaaagcggaaattgcttcgacaagaattaaatttttagaggccaggataaccaatctttcatacaatgacatatta                                                                                                                                               <  2:372455‑M1/251‑90 (MQ=255)
                                                                                       cACTTCACACACTATCAGAGAACATAATAACTTTTTAAAAGGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAaccttttacgagcaaaagcggaaattgcttcgacaataattaaatttttagaggccaggataaccaatcttgcatacaatgacatatta                                                                                                                                               <  2:373810‑M1/251‑90 (MQ=255)
                                                                                          ttCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAaccttttacgagcaaaagcggaaattgcttcgacaataattaaatttttagaggccaggataaccaatctttcatacaatgacatattaaa                                                                                                                                             >  2:104611‑M1/1‑159 (MQ=255)
                                                                                              cacacTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAAGATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAaccttttacgagcaaaagcggaaattgcttcgacaataattaaatttttagaggccaggataatcaatctttcatacaatgacatattaaaatatc                                                                                                                                        <  2:106266‑M1/251‑97 (MQ=255)
                                                                                                 acTATCAGAGAACATAATAACTTTTTAAAATNAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAaccttttacgagcaaaagcggaaattgcttcgacaataattaaatttttagaggccaggataaccaatctttcatacaatgacatattaaaatatcagc                                                                                                                                     <  2:252418‑M1/251‑100 (MQ=255)
                                                                                                                    aCTTTTTAAAATGAGCGCATGTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAaccttttacgagcaaaagcggaaattgcttcgacaataattaaatttttagaggccaggataaccaatctttcagacaatgacatattaaaatatcagcaagaattccaaagacaatg                                                                                                                  <  2:424713‑M1/251‑119 (MQ=255)
                                                                                                                             aaTGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTGTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAtccgtttacgagcaaaagcggaaattgcttcgacaataattaaatttttagaggccaggataaccaatcgttcatacaatgacatattaaaatatcagcaagaattccaaagacaaagttataagca                                                                                                         <  1:445637‑M1/251‑128 (MQ=255)
                                                                                                                              aTGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAaccttttacgagcaaaagcggaaattgcttcgacaataattaaatttttagaggccaggataaccaatctttcatacaatgacatattaaaatatcagcaagaattccaaagacaatgttataagca                                                                                                         >  1:268765‑M1/1‑123 (MQ=255)
                                                                                                                              aTGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATAAAAGGTAATTATATCAGGAAGAAAaccttttacgcgcaaaagcggaaattgcttcgacaataattaaatttttagaggccaggataaccaatctttaatacaatgacatattaaaatatcagcaagaattccaaagacaatgttataagca                                                                                                         >  1:268738‑M1/1‑123 (MQ=255)
                                                                                                                                      aTTTTGCTAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAacattttacgagcaaaagc                                                                                                                                                                                                                     >  2:193903‑M1/1‑115 (MQ=255)
                                                                                                                                      aTTTTGATAGTTGTTAACGATTTCCTTTACGATTCTATCTCTACCCCTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAaccttttacgagcaaaagc                                                                                                                                                                                                                     <  1:193903‑M1/134‑20 (MQ=255)
                                                                                                                                                         aTTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAaccttttacgagcaaaagcggaaattgcttcgacaataattaaatttttagaggccaggataaccaatctttcatacaatgacatattaaaatatcagcaagaattccaaagacaaggttataagcaagttaaagcctttacgacactctcccga                                                                             <  1:418740‑M1/251‑156 (MQ=255)
                                                                                                                                                                                              cacaAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAaccttttacgagcaaaagcggaaattgcttcgacaataattaaatttttagaggccaggataaccaatctttcatacaatgacatattaaaatatcagcaagaattccaaagacaatgttataagcaagttaaagcctttacgacactctcccgatacaataaaattcaaacatgggcagaaatgtctgaat                                        >  2:363572‑M1/1‑59 (MQ=255)
                                                                                                                                                                                                          cGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAaccttttacgagcaaaagcggaaattgcttcgacaataattaaatttttagaggccaggataaccaatctttcatacaatgacatattaaaatatcagcaagaattccaaagacaatgttataagcaagttaaagcctttacgacactctcccgatacaataaaattcaaacatgggcagaaatgtctgaatatcaattcgaag                            <  2:189356‑M1/251‑205 (MQ=255)
                                                                                                                                                                                                                  gaTAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATAaccttttacgagcaaaagcggaaattgcttcgacaataattaaatttttagaggccaggataaccaatctttcatacaatgacatattaaaatatcagcaagaattccaaagacaatgttataagcaagttaaagcctttacgacactctcccgatacaataaaattcaaacatgggcagaaatgtctgaatatcaattcgaagtatttcaa                    <  1:380928‑M1/251‑213 (MQ=255)
                                                                                                                                                                                                                             aaTCAAAGGTAATTATATCAGGAAGATAaccttttacgagcaaaagcggaaattgcttcgacaataattaaatttttagaggccaggataaccaatctttcatacaatgacatattaaaatatcagcaagaattccaaagacaatgttataagcaagttaaagcctttacgacactctcccgatacaataaaattcaaacatgggcagaaatgtctgaatatcaattcgaagtatttcaatacgagacgct         <  1:191335‑M1/251‑224 (MQ=255)
                                                                                                                                                                                                                                   aGGTAATTATATCAGGAAGATAaccttttacgagcaaaagcggaaattgcttcgacaataattaaatttttagaggccaggataaccaatctttcatacaatgacatattaaaatatcagcaagaattccaaagacaatgttataagcaagttaaagcctttacgacactctcccgatacaataaaattcaaacatgggcagaaatgtctgaatatcaattcgaagtatttcaatacgagacgcttaaccc   <  1:457894‑M1/251‑230 (MQ=255)
                                                                                                                                                                                                                                     gTAATTATATCAGGAAGATAaccttttacgagcaaaagcggaaattgcttcgacaataattaaatttttagaggccaggataaccaatctttcatacaatgacatattaaaatatcagcaagaattccaaagacaatgttataagcaagttaaagcctttacgacactctcccgatacaataaaattcaaacatgggcagaaatgtctgaatatcaattcgaagtatttcaatacgagacgcttaaccca  >  1:227203‑M1/1‑20 (MQ=255)
                                                                                                                                                                                                                                           atatCAGGAAGATAaccttttacgagcaaaagcggaaattgcttcgacaataattaaatttttagaggccaggataaccaatctttcatacaatgacatattaaaatatcagcaagaattccaaagacaatgttataagcacgttaaa                                                                                                  >  2:72766‑M1/1‑14 (MQ=255)
                                                                                                                                                                                                                                           atatCAGGAAGATAaccttttacgagcaaaagcggaaattgcttcgacaataattaaatttttagaggccaggataaccaatctttcatacaatgacatattaaaatatcagcaagaattccaaagacaatgttataagcaagttaaa                                                                                                  <  1:72766‑M1/148‑135 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      
CAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAAC  >  NC_000913/3650942‑3651420

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