Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,335,153 Δ593 bp [ymiA]yciX [ymiA], yciX

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1335153 1335745 593 22 [0] [2] 26 [ymiA]–yciX [ymiA],yciX

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 13351520 (0.000)22 (0.590) 20/236 1.0 95.7% coding (29/165 nt) ymiA uncharacterized protein
?NC_000913 1335746 = 2 (0.050)intergenic (+288/‑85) yciX/acnA uncharacterized protein/aconitate hydratase 1; aconitase A

TGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAG  >  NC_000913/1335606‑1335872
                                                                                                                                            |                                                                                                                              
aaaagattgaattattcactggagacgattcgtaaattcgctggatttttgactaaagtttctttaggatttgtctcatcgtatgattaattacaccaatcgttcaacgtcttgctgataagtgatggcgattatatgagAAATCGTTa                                                                                                                        >  5:15855‑M2/141‑149 (MQ=255)
                 actggagacgattcgtaaattcgctggatttttgactaaagtttctttaggatttgtctcatcgtatgattaattacaccaatcgttcaacgtcttgctgataagtgatggcgattatatgagAAATCGTTAc                                                                                                                       >  3:176696‑M2/124‑133 (MQ=255)
                 actggagacgattcgtaaattcgctggatttttgactaaagtttctttaggatttgtctcatcgtatgattaattacaccaatcgttcaacgtcttgctgataagtgatggcgattatatgagAAATCGTTAc                                                                                                                       <  4:176696‑M2/10‑1 (MQ=255)
                                       gctggatttttgactaaagtttctttaggatttgtctcatcgtatgattaattacaccaatcgttcaacgtcttgctgataagtgatggcgattatatgagAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGaa                                                                                    >  6:137313‑M2/102‑146 (MQ=255)
                                                    ctaaagtttctttaggatttgtctcatcgtatgattaattacaccaatcgttcaacgtcttgctgataagtgatggcgattatatgagAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCg                                                                    >  7:116474‑M2/89‑149 (MQ=255)
                                                         ggttctttaggatttgtctcatcgtatgattaattacacaaatcgttcaacgtcttgctgataagtgattgcgattatatgagAAATCGTTACGCGATGTTTGTGATATCTTTAATATTCACCCTGAAGAGAATCAGATCTTCGCAAcc                                                               >  7:38841‑M2/84‑149 (MQ=255)
                                                                taggatttgtctcatcgtatgattaattacaccaatcgttcaacgtcttgctgataagtgatggcgattatatgagAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCa                                                         <  6:15855‑M2/72‑1 (MQ=255)
                                                                   gatttgtgtcatcgtatgattaattactccaatcgttcaacgtgtttctgatatgtgatggcgattatatgagTAATCGTTACGCTATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCTGGGCTTCGCAACGCTGTCATTaa                                                     <  8:38841‑M2/76‑1 (MQ=255)
                                                                   gatttgtctcatcgtatgattaattacaccaatcgttcaacgtcttgctgataagtgatggcgattatatgagAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATc                                                                             >  6:162828‑M2/74‑125 (MQ=255)
                                                                   gatttgtctcatcgtatgattaattacaccaatcgttcaacgtcttgctgataagtgatggcgattatatgagAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATc                                                                             <  5:162828‑M2/52‑1 (MQ=255)
                                                                         tctcatcgtatgattaattacaccaatcgttcaacgtcttgctgataagtgatggcgattatatgagAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTaa                                                     >  8:19807‑M2/68‑143 (MQ=255)
                                                                         tctcatcgtatgattaattacaccaatcgttcaacgtcttgctgataagtgatggcgattatatgagAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTaa                                                     <  7:19807‑M2/76‑1 (MQ=255)
                                                                             atcgtatgattaattacaccaatcgttcaacgtcttgctgataagtgatggcgattatatgagAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCt                                            <  8:116474‑M2/85‑1 (MQ=255)
                                                                                     attaattacaccaatcgttcaacgtcttgctgataagtgatggcgattatatgagAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATgtcgtc                                    <  2:86083‑M2/93‑1 (MQ=255)
                                                                                     attaattacaccaatcgttcaacgtcttgctgataagtgatggcgattatatgagAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATgtcgtc                                    >  1:86083‑M2/56‑148 (MQ=255)
                                                                                                     cacaAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATg                                         >  2:124534/1‑127 (MQ=255)
                                                                                                     cacaAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATg                                         <  1:124534/127‑1 (MQ=255)
                                                                                                      ttcaacgtcttgctgataagtgatggcgattatatgagAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAg                  <  5:137313‑M2/111‑1 (MQ=255)
                                                                                                                          tgatggcgattatatgagAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAg                      <  7:27876‑M2/107‑1 (MQ=255)
                                                                                                                          tgatggcgattatatgagAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAg                      >  8:27876‑M2/19‑125 (MQ=255)
                                                                                                                              ggcgattatatgagAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCaa                                  <  6:69361‑M2/95‑1 (MQ=255)
                                                                                                                              ggcgattatatgagAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCaa                                  >  5:69361‑M2/15‑109 (MQ=255)
                                                                                                                                         gagAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAggagga                                               >  7:115822‑M2/4‑85 (MQ=255)
                                                                                                                                         gagAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAggagga                                               <  8:115822‑M2/82‑1 (MQ=255)
                                                                                                                                            aaaTCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAg  <  5:34140/127‑1 (MQ=255)
                                                                                                                                            aaaTCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAg  >  6:34140/1‑127 (MQ=255)
                                                                                                                                            |                                                                                                                              
TGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAG  >  NC_000913/1335606‑1335872

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