Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 390,251 Δ14,565 bp [yaiT][dgcC] 16 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 390251 404815 14565 33 [0] [0] 32 [yaiT]–[dgcC] [yaiT],ytiB,insF2,insE2,yaiT,iprA,ampH,sbmA,yaiW,yaiY,yaiZ,ddlA,iraP,phoA,psiF,[dgcC]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3902500 (0.000)32 (0.540) 29/494 0.7 100% intergenic (‑523/‑1) hemB/yaiT porphobilinogen synthase/putative autotransporter YaiT
?NC_000913 404816 = 0 (0.000)coding (1114/1116 nt) dgcC diguanylate cyclase DgcC

GATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGCCATGTTACGGGTGCATGAAGGGCTAAATACATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACGCCTGGGCGCGTGGCAT  >  NC_000913/404573‑405062
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
atatttattatgcgtatgcttcaaaaacaaaattatctgcaacatgttgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTcc                                                                                                                                                                                                                                                 <  1:455543‑M2/8‑1 (MQ=255)
               atgcttcaaaaacaaaattatctgcaacatgttgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGtcatca                                                                                                                                                                                                                                   <  1:228793‑M2/22‑1 (MQ=255)
               atgcttcaaaaacaaaattatctgcaacatgttgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataagtgaaccattcgtctgtaacgcagcgtaatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGtcatca                                                                                                                                                                                                                                   <  1:392047‑M2/22‑1 (MQ=255)
                             aaattatctgcaacatgttgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGttt                                                                                                                                                                                                                    <  1:313879‑M2/37‑1 (MQ=255)
                                      gcaacatgttgaaatactttgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaGGACGTCCGGCGAAAGTCAGCAGGATTTGCTGAGTTTTTCTGAt                                                                                                                                                                                                             <  2:171491‑M2/44‑1 (MQ=255)
                                                    tactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCg                                                                                                                                                                                             >  1:349963‑M2/192‑251 (MQ=255)
                                                        gtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGCTTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCAt                                                                                                                                                                                         >  2:211729‑M2/188‑251 (MQ=255)
                                                          gcttttatgaatttgatgcgtgtttttatccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCg                                                                                                                                                                                       >  1:464199‑M2/186‑251 (MQ=255)
                                                                   aatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATc                                                                                                                                                                               >  2:308396‑M2/177‑250 (MQ=255)
                                                                                         ataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTtctc                                                                                                                                                        <  2:315024‑M2/97‑1 (MQ=255)
                                                                                             attatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTcc                                                                                                                                                    >  2:53742‑M2/151‑251 (MQ=255)
                                                                                             attatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTcc                                                                                                                                                    >  1:239469‑M2/151‑251 (MQ=255)
                                                                                               tatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCa                                                                                                                                                   <  1:308396‑M2/102‑1 (MQ=255)
                                                                                                                ggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCa                                                                                                                                  >  1:306104‑M2/132‑250 (MQ=255)
                                                                                                                      aggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAAtggtgg                                                                                                                           <  2:457105‑M2/126‑1 (MQ=255)
                                                                                                                      aggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAAtggtgg                                                                                                                           <  1:53742‑M2/126‑1 (MQ=255)
                                                                                                                                 taattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGt                                                                                                                >  1:371298‑M2/115‑251 (MQ=255)
                                                                                                                                         ccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGcagcc                                                                                                        >  2:362018‑M2/107‑249 (MQ=255)
                                                                                                                                                                   gttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGc                                                                                                            >  1:70261‑M2/81‑221 (MQ=255)
                                                                                                                                                                   gttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGc                                                                                                            <  2:70261‑M2/141‑1 (MQ=255)
                                                                                                                                                                        cgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTccc                                                                         >  1:198066‑M2/76‑251 (MQ=255)
                                                                                                                                                                           gaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggaggggtgacaaTGACGTCCGGCGAAAGTCCTCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGCGAGCAGACAATATCTTTCAGTGCCGCCGGATGTCTTCCCGtt                                                                      >  2:364850‑M2/73‑251 (MQ=255)
                                                                                                                                                                                   atatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGAt                                                                                                                                                                                                                               >  2:434864‑M2/65‑90 (MQ=255)
                                                                                                                                                                                   atatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGAt                                                                                                                                                                                                                               <  1:434864‑M2/26‑1 (MQ=255)
                                                                                                                                                                                      tgttttctatgcattgataattgatggatcaacttattacgtccctgagcagggatgacaaTGACGTCCGGCGAAAGTAATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCAt                                                           >  1:365238‑M2/62‑251 (MQ=255)
                                                                                                                                                                                                gcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACCCTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGaa                                                 <  1:419404‑M2/200‑1 (MQ=255)
                                                                                                                                                                                                     gataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCAt                                            >  1:190795‑M2/47‑251 (MQ=255)
                                                                                                                                                                                                               ggatcaaattattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCtgg                                  >  2:161043‑M2/37‑250 (MQ=255)
                                                                                                                                                                                                                                         gggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTACACGCCTGGGcgcg        >  2:40805‑M2/11‑251 (MQ=255)
                                                                                                                                                                                                                                           gatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACGCCTGGGCGCGTg      >  2:86582‑M2/9‑251 (MQ=255)
                                                                                                                                                                                                                                            atgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCAAACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACGCCTGGGCGCGTgg     <  1:101801‑M2/244‑1 (MQ=255)
                                                                                                                                                                                                                                               acaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACGCCTGGGCGCGTGGCAt  >  2:463916‑M2/5‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
GATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGCCATGTTACGGGTGCATGAAGGGCTAAATACATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACGCCTGGGCGCGTGGCAT  >  NC_000913/404573‑405062

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