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 29 [0] [0] 29 [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)29 (0.550) 27/486 0.6 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

ATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGCCATGTTACGGGTGCATGAAGGGCTAAATACATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACGCCTGGGCGCGT  >  NC_000913/404574‑405057
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
tatttattatgcgtatgcttcaaaaacaaaattatctgcaacatgttgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCg                                                                                                                                                                                                                                           >  1:123472‑M2/243‑251 (MQ=255)
                gcttcaaaaacaaaattatctgcaacatgttgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgagaagggagggatcaacttagtacgtcccggaggagggatgacaaTGACGTCCGGCGAAAGTCATCAgg                                                                                                                                                                                                                            <  1:281639‑M2/24‑1 (MQ=255)
                      aaaacaaaattctctgcaacatgatgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCt                                                                                                                                                                                                                     >  1:133558‑M2/221‑251 (MQ=255)
                                            gccgatctactgtgcttgtatgaagttgatgcgtggtgttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataagatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCTAAAGTCAGCAGGATTTGCGGAGTTTTTCTGATTTTTCCATa                                                                                                                                                                                               <  2:191086‑M2/53‑1 (MQ=255)
                                               gaaatactgtgcttttatgaatttgatgcgtgtttttctccataacttatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGCTTTGCTGCGTTTTTCTGATTTTTCCATacac                                                                                                                                                                                            >  2:282508‑M2/196‑251 (MQ=255)
                                                    actgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagactaattgacccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctggggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTtttccataagttcgg                                                                                                                                                                                        >  1:191086‑M2/191‑243 (MQ=255)
                                                                                 tttctccataaattatatatgtccacattcggacttaggggacagaatacttgaaccattcctctgtaacgcagcataatcgttagcacgaaacataccatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGCCGCCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGc                                                                                                                                                          >  1:442331‑M2/162‑251 (MQ=255)
                                                                                               atatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGt                                                                                                                                            >  2:380377‑M2/148‑251 (MQ=255)
                                                                                               atatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCAGCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGt                                                                                                                                            <  1:472601‑M2/104‑1 (MQ=255)
                                                                                                       ccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTAc                                                                                                                                    <  1:130113‑M2/112‑1 (MQ=255)
                                                                                                                           aagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgtgttctatgcattgataattgatggatcaacttattacgtccctgaggatggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATGTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAAGGGTGGTGCCTc                                                                                                                <  2:280647‑M2/132‑1 (MQ=255)
                                                                                                                             gaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTccc                                                                                                              >  1:157482‑M2/118‑251 (MQ=255)
                                                                                                                             gaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTccc                                                                                                              >  1:21617‑M2/118‑251 (MQ=255)
                                                                                                                                     gaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAGGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGc                                                                                                       >  2:263051‑M2/110‑250 (MQ=255)
                                                                                                                                                         agcataatcgttagcgcgaaacctaatatgttttctatgcattgataattgatggctcaacttattaagtccctgaggagggatgacaaTGACGTCCGGCGACAGTCATCAGGCTTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTAGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTAc                                                                                                                                    >  2:210358‑M2/90‑201 (MQ=255)
                                                                                                                                                         agcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTGTTCTGATTTGTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTAc                                                                                                                                    <  1:210358‑M2/112‑1 (MQ=255)
                                                                                                                                                               atcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGAt                                                                            <  1:306663‑M2/168‑1 (MQ=255)
                                                                                                                                                                              cataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCa                                                             <  2:133558‑M2/183‑1 (MQ=255)
                                                                                                                                                                               ataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTcc                                                              >  1:369004‑M2/68‑249 (MQ=255)
                                                                                                                                                                                  atatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCAc                                                         >  1:116329‑M2/65‑251 (MQ=255)
                                                                                                                                                                                                 attgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCTAAAGTCATCAGGAGTTGCTGAGTTTTTCTGAGTTTTCCATACACTGCGTCAt                                                                                                                                                                                    <  2:51906‑M2/64‑1 (MQ=255)
                                                                                                                                                                                                 attgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCAt                                                                                                                                                                                    >  1:51906‑M2/50‑113 (MQ=255)
                                                                                                                                                                                                  ttgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAAccc                                         >  1:260866‑M2/49‑251 (MQ=255)
                                                                                                                                                                                                        attgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTAcc                                   >  2:386040‑M2/43‑251 (MQ=255)
                                                                                                                                                                                                                      ttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCaa                      >  1:254685‑M2/29‑250 (MQ=255)
                                                                                                                                                                                                                        attacgtccctgcggagggctaacaaTGACGTCCGGCGACAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCCGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCggcaactt                   >  2:246158‑M2/27‑248 (MQ=255)
                                                                                                                                                                                                                                  tgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACGCCTg         >  2:210835‑M2/17‑251 (MQ=255)
                                                                                                                                                                                                                                    aggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACTCCTggg       <  1:386040‑M2/237‑1 (MQ=255)
                                                                                                                                                                                                                                         ggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACGCCTGGGCGCGt  >  1:356580‑M2/10‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
ATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGCCATGTTACGGGTGCATGAAGGGCTAAATACATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACGCCTGGGCGCGT  >  NC_000913/404574‑405057

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