Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,586,443 Δ555 bp [ydeP][ydeQ] [ydeP], [ydeQ]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1586443 1586997 555 25 [0] [0] 25 [ydeP]–[ydeQ] [ydeP],[ydeQ]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 15864420 (0.000)25 (0.440) 24/472 0.9 100% coding (45/2280 nt) ydeP putative oxidoreductase YdeP
?NC_000913 1586998 = 0 (0.000)coding (737/915 nt) ydeQ putative fimbrial adhesin protein YdeQ

AGGTAAGTAATGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGCACAGTACATACAAAGGGATTATCGCCAACATGGATATATAAGAGTAGGGGACTACTCTCTTTCTATTTTTAATCTATTATCGTTAAACAATAGATTAAATACGATAAGTTAAGCATCATCTCTCATAGTGTTTAAATTAA  >  NC_000913/1586195‑1586678
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            
aGGTAAGTAATGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCact                                                                                                                                                                                                                                           <  2:154900‑M1/251‑4 (MQ=255)
 ggTAAGTAATGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactc                                                                                                                                                                                                                                          >  1:536159‑M1/1‑247 (MQ=255)
      gtaaTGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgat                                                                                                                                                                                                                                     <  2:407618‑M1/251‑10 (MQ=255)
      gtaaTGATTGAACCGTATCCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgat                                                                                                                                                                                                                                     <  2:220075‑M1/251‑10 (MQ=255)
                   cGTATTCGCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTGCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcat                                                                                                                                                                                                                        <  2:365349‑M1/251‑23 (MQ=255)
                      aTTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatca                                                                                                                                                                                                                      <  2:453271‑M1/250‑25 (MQ=255)
                         ctcAGCAAAGAAAGAGGCGTTTACCTGCTTATCAGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtg                                                                                                                                                                                                                  >  2:400385‑M1/1‑223 (MQ=255)
                            aGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTGCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGGCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgctgttgttgcatcaggtgctg                                                                                                                                                                                                               <  2:575000‑M1/251‑32 (MQ=255)
                                                          cGTGACTTCCCAGTCTATTGCTTGTGCGCCGTTTTCACAAATGGCGAATGTCGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGGCAAACATGGCAATAATATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcgtttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcag                                                                                                                                                                                 <  1:590586‑M1/251‑62 (MQ=255)
                                                                    cAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAAGAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttggtgcatcaggtgctgttttagcaatcacttgtcggctactgtcagtggtggcacc                                                                                                                                                                       <  1:486241‑M1/251‑72 (MQ=255)
                                                                             gCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagat                                                                                                                                                               >  1:170360‑M1/1‑171 (MQ=255)
                                                                                    cgcCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgag                                                                                                                                                       <  1:597265‑M1/251‑88 (MQ=255)
                                                                                        gTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagct                                                                                                                                                   <  1:251854‑M1/251‑92 (MQ=255)
                                                                                                    gTCGAAAGACGCACTGTGCTTAGGATTGGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCTGCGGCTTACGTGCTGCATGCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacatgtgtcggctactgtcagtggttgcacccgaaagataaaatgagagcttttggtcgctt                                                                                                                                        <  2:519690‑M1/250‑103 (MQ=255)
                                                                                                       gAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgcaagataaaatgagagcttttgttcgcttgagc                                                                                                                                    >  2:22031‑M1/1‑145 (MQ=255)
                                                                                                                   gtgCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtgttggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgcca                                                                                                                         <  2:508330‑M1/250‑118 (MQ=255)
                                                                                                                                ggCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgc                                                                                                           <  1:424808‑M1/251‑132 (MQ=255)
                                                                                                                                                         aGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaaagtccgggagatcc                                                                                  >  1:338087‑M1/1‑95 (MQ=255)
                                                                                                                                                                              cAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaaagtccgggagatccacggtgacatttcgactatca                                                             <  1:376428‑M1/251‑178 (MQ=255)
                                                                                                                                                                                     ggCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaaagtccgggagatccacggtgacatttcgactatcaacggtgc                                                      <  1:400385‑M1/251‑185 (MQ=255)
                                                                                                                                                                                          tAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaaagtccgggagatccacggtgacatttcgactatcaacggtgcagccg                                                 >  1:261135‑M1/1‑62 (MQ=255)
                                                                                                                                                                                                             cATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaaagtccgggagatccacggtgacatttcgactatcaacggtgcagccgcctggtggcataacaacat                              <  2:536159‑M1/251‑209 (MQ=255)
                                                                                                                                                                                                                             gCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaaagtccgggagatccacggtgacatttcgactatcaacggtgcagccgcctgttggcataacaacattgttatttgaaata                <  1:22031‑M1/249‑223 (MQ=255)
                                                                                                                                                                                                                              cATTCGCTCCGGATTTAACAGCCCCCactcccgatgcttttgttgcatcaggtgctgttttagcaaacacttgtcggctactgtcagtgggggcacccgaaagataaaatgagagcttttgttcgctttagcaatatacgccaaggtgtatgtctgtgcttccaggaaagtcctggagagccatggtgacatgtcgactatcaacggtgcagtcgcctgttggcataacaacattgttatgtgaaataata             <  2:350517‑M1/251‑226 (MQ=255)
                                                                                                                                                                                                                                          aTTTAACAGCACCCactcccgatgcttttgttgcatcaggagctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaaagtccgggagatccacggtgacatttcgactatcaacggtgcagccgcctgttggcataacaacattgttatttgaaatattattccacgtaaa  >  1:322922‑M1/1‑14 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            
AGGTAAGTAATGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGCACAGTACATACAAAGGGATTATCGCCAACATGGATATATAAGAGTAGGGGACTACTCTCTTTCTATTTTTAATCTATTATCGTTAAACAATAGATTAAATACGATAAGTTAAGCATCATCTCTCATAGTGTTTAAATTAA  >  NC_000913/1586195‑1586678

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