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

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913_3_pae_tpiA 257908–258674 258683 10–776 25 [0] [0] 25 [crl] [crl]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_tpiA = 2579070 (0.000)24 (0.810) 22/268 0.4 100% 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_3_pae_tpiA 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

TGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATT  >  NC_000913_3_pae_tpiA/257764‑258055
                                                                                                                                               |                                                                                                                                                    
tgtgTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgc                                                                                                                                                 >  1:410106‑M1/1‑144 (MQ=255)
                 aaaTTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctt                                                                                                                                >  2:18366‑M1/1‑127 (MQ=255)
                         ccAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGCCCAACCAAAGTGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgatt                                                                                                                        <  1:128417‑M1/149‑31 (MQ=255)
                                   aaaaaaGTTGCATGACAACAGGCGATAGCACTGACGTTACCGAGTGGACACCCGAAGAGCAGATTGAGCAAAAAATTTACCGCCCTAGGCCCGTATATTTGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgt                                                                                                              <  2:226845‑M1/149‑41 (MQ=255)
                                               tCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTCGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtgtggctgtttgagtcaccgt                                                                                                  >  2:376531‑M1/1‑97 (MQ=255)
                                                 acaacaGGAGATAGCAATGACGGTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtga                                                                                                <  1:280539‑M1/149‑55 (MQ=255)
                                                            tAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccg                                                                                     <  1:454181‑M1/149‑66 (MQ=255)
                                                            tAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccagcccca                                                                                     <  1:282965‑M1/149‑66 (MQ=255)
                                                                aaTGACGCTACAGAGTGGACACCCGAAGAGCAGATTGATCACAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaag                                                                                 <  2:346322‑M1/149‑70 (MQ=255)
                                                                  tGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCCCACTAGGCCCGTATTTTCGTGAAGGTAaatgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtg                                                                               >  1:422508‑M1/1‑78 (MQ=255)
                                                                     cGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgt                                                                            <  2:427561‑M1/149‑75 (MQ=255)
                                                                      gTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtg                                                                           <  1:4152‑M1/149‑76 (MQ=255)
                                                                           cGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattc                                                                      >  1:485274‑M1/1‑69 (MQ=255)
                                                                                         aaGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggat                                                        <  1:18366‑M1/149‑95 (MQ=255)
                                                                                          agagCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatg                                                       <  1:38738‑M1/149‑96 (MQ=255)
                                                                                             gCAGATTGATCAAAAAATTTACCGCCCTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggag                                                    <  1:156960‑M1/149‑99 (MQ=255)
                                                                                                        aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggcttgtggatggagctggaagcgca                                         <  2:485274‑M1/149‑110 (MQ=255)
                                                                                                         aaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcag                                        >  2:467388‑M1/1‑39 (MQ=255)
                                                                                                         aaaaaTTTACCGCACTAGGCCCCTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagagcggggattttggggcgcgtgggaggagctggaagtccag                                        >  2:161666‑M1/1‑39 (MQ=255)
                                                                                                                   cGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattcttttttttttttttttctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggttggcgctggaagcgcagggatcccgtt                              >  2:431788‑M1/1‑29 (MQ=255)
                                                                                                                   cGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataaacgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgtt                              <  1:152004‑M1/149‑121 (MQ=255)
                                                                                                                      aCTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgggaattctggggcaggtggatggagctggaagcgcaggaatcccgtttta                           >  2:433665‑M1/1‑26 (MQ=255)
                                                                                                                           gCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggatgagggggggagctggaagcgcaggaatcccgttttacctac                      >  1:13668‑M1/1‑21 (MQ=255)
                                                                                                                                       gTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagttt          >  1:339205‑M1/1‑9 (MQ=255)
                                                                                                                                               aagtgcaaagataagcgattcattttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgtt  <  2:13293‑M1/149‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATT  >  NC_000913_3_pae_tpiA/257764‑258055

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