Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 257,908 Δ776 bp [crl] [crl]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 257908 258683 776 31 [1] [0] 30 [crl] [crl]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 2579071 (0.030)22 (0.820) 16/116 0.3 98.0% intergenic (+8/‑769) crl/crl pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator;regulator; Surface structures; transcriptional regulator of cryptic csgA gene for curli surface fibers/pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator;regulator; Surface structures; transcriptional regulator of cryptic csgA gene for curli surface fibers
?NC_000913 258684 = 0 (0.000)pseudogene (9/331 nt) crl pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator;regulator; Surface structures; transcriptional regulator of cryptic csgA gene for curli surface fibers

ACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGG  >  NC_000913/258619‑258749
                                                                 |                                                                 
tggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGt                                                                 <  1:270898‑M2/3‑1 (MQ=255)
tggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgcgaaggtaAGt                                                                 <  2:1048776‑M2/3‑1 (MQ=255)
 ggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTg                                                                <  2:358948‑M2/4‑1 (MQ=255)
 ggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTg                                                                <  2:853142‑M2/4‑1 (MQ=255)
      cccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAg                                                           <  1:809121‑M2/9‑1 (MQ=255)
             agcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCg                                                    >  1:913676‑M2/53‑68 (MQ=255)
               cagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGAt                                                  <  2:694588‑M2/18‑1 (MQ=255)
                  attgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGAGCAAAGATAATCGATTCt                                               <  1:1087159‑M2/21‑1 (MQ=255)
                   ttgatcaaaaaatttaccgcactaggcccgtatattcgtgatggtaAGTGCAAAGATAATCGATTCtt                                              <  1:1068898‑M2/22‑1 (MQ=255)
                   ttgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTctg                                              >  1:261098‑M2/47‑67 (MQ=255)
                   ttgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCtt                                              <  2:94427‑M2/22‑1 (MQ=255)
                        caaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCg                                         >  2:650948‑M2/42‑68 (MQ=255)
                            aaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTg                                     >  2:362720‑M2/38‑68 (MQ=255)
                             aatctaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGCTTCTTTTTCGATTGt                                    <  2:913673‑M2/32‑1 (MQ=255)
                               ttaaccgcactagtcccgtatattcgtgaaggtaAGTGCAAAGATTATCGATTCTTTTTCGTTTGTCt                                  <  1:243702‑M2/34‑1 (MQ=255)
                                   ccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCt                              <  1:825158‑M2/38‑1 (MQ=255)
                                   ccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCt                              <  1:257805‑M2/38‑1 (MQ=255)
                                       actaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTAt                          >  1:203958‑M2/27‑68 (MQ=255)
                                             cccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCa                    >  1:825400‑M2/21‑68 (MQ=255)
                                             cccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCa                    >  2:78318‑M2/21‑68 (MQ=255)
                                                  atattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTg               >  1:1045751‑M2/16‑68 (MQ=255)
                                                       cgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAAcc          >  2:1134559‑M2/11‑68 (MQ=255)
                                                          gaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCg                       >  2:33233‑M2/8‑52 (MQ=255)
                                                          gaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCg                       <  1:33233‑M2/45‑1 (MQ=255)
                                                           aaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCa      <  2:305389‑M2/62‑1 (MQ=255)
                                                           aaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCa      <  1:362720‑M2/62‑1 (MQ=255)
                                                           aaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCa      <  1:103299‑M2/62‑1 (MQ=255)
                                                            aggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCAc     >  2:770566‑M2/6‑68 (MQ=255)
                                                             ggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCAcc    >  2:866319‑M2/5‑68 (MQ=255)
                                                               taAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCgg  >  2:956734‑M2/3‑68 (MQ=255)
                                                                 |                                                                 
ACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGG  >  NC_000913/258619‑258749

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