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 28 [0] [0] 28 [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)28 (0.620) 28/466 0.5 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

AAGTAATGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGCACAGTACATACAAAGGGATTATCGCCAACATGGATATATAAGAGTAGGGGACTACTCTCTTTCTATTTTTAATCTATTATCGTTAAACAATAGATTAAATACGATAAGTTAAGCATCATCTCTCATAGTGTTTA  >  NC_000913/1586199‑1586672
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                      
aagtaaTGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccg                                                                                                                                                                                                                                 <  2:348396‑M1/251‑8 (MQ=255)
  gtaaTGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgat                                                                                                                                                                                                                               >  2:392323‑M1/1‑242 (MQ=255)
               cGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcat                                                                                                                                                                                                                  >  1:122369‑M1/1‑229 (MQ=255)
                         gCAAAGAAAGAGGCGTTTTCCTGCGTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCtctcccgatgcttttgttgcatcaggtgctgt                                                                                                                                                                                                        <  2:293488‑M1/251‑33 (MQ=255)
                                  gagGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaa                                                                                                                                                                                               <  1:221928‑M1/251‑42 (MQ=255)
                                       gTTTACCTGCTTATCCGTGACTTCCCAGGCGAGTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACTCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTTCTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacactt                                                                                                                                                                                          <  1:121943‑M1/251‑47 (MQ=255)
                                            ccTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcgg                                                                                                                                                                                     >  1:314259‑M1/1‑200 (MQ=255)
                                                 ttATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctact                                                                                                                                                                                >  2:146827‑M1/1‑195 (MQ=255)
                                                           cTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGAAAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtg                                                                                                                                                                      >  2:278912‑M1/1‑185 (MQ=255)
                                                                cAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacc                                                                                                                                                                 <  2:188255‑M1/251‑72 (MQ=255)
                                                                          cTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgtt                                                                                                                                                                                                     <  2:370035‑M1/205‑36 (MQ=255)
                                                                          cTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgtt                                                                                                                                                                                                     >  1:370035‑M1/1‑170 (MQ=255)
                                                                                                   gAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCACAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagc                                                                                                                              >  1:158595‑M1/1‑145 (MQ=255)
                                                                                                            aCTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGAACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggagctgtgttagcaaacacttgtcggctactgtcagtggtgtcacccgaaagataaaatgagagctattgttcgcttgagcaatatacgc                                                                                                                     >  1:294183‑M1/1‑136 (MQ=255)
                                                                                                                            ggCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcggttgttcgcttgagcaatatacgccaaggggtatttctg                                                                                                      <  1:424019‑M1/250‑131 (MQ=255)
                                                                                                                              ccATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgc                                                                                                   >  2:347222‑M1/1‑118 (MQ=255)
                                                                                                                                         ccGGAAAGTCAAAGCCCTCTGGCTTATTAATGTCAAACATGGCAATAACATCCTGGCGTCTATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaa                                                                                          >  2:400757‑M1/1‑107 (MQ=255)
                                                                                                                                                                 tATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagag                                                                                                                                               <  2:451503‑M1/172‑90 (MQ=255)
                                                                                                                                                                 tATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagag                                                                                                                                               >  1:451503‑M1/1‑83 (MQ=255)
                                                                                                                                                                                     aTAACATCCTGGCGTATATCCACCTGCTTACGTACTGCAGTCGCTACGGATTTAACAGCACCCactccctatgcttttgttgcatcaggtgctgtgttagcaaacactggtcgggtactgtcaggggtggctcccgaaagataaaatgggagctttggtgcgcttgagcaatatacgccaaggggtatgtctgcgcttccaggaaagtccgggagatccacggtgacagttcgaccatcaactctgcagcc                                            <  2:40657‑M1/251‑189 (MQ=255)
                                                                                                                                                                                               ggCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaaagtccgggagatccacggtgacatttcgactatcaacggtgcagccgcctgttg                                    <  1:53460‑M1/249‑197 (MQ=255)
                                                                                                                                                                                                 cGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaaagtccgggagatccacggtgacatttcgactatcaacggtgcagccgcctgtgggcat                                <  1:392323‑M1/251‑201 (MQ=255)
                                                                                                                                                                                                  gTATATTCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaaagtccgggagatccacggtgacatttcgactatcaacggtgcagccgcctgttggcat                                >  1:136135‑M1/1‑50 (MQ=255)
                                                                                                                                                                                                         cATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcgtttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccagtaaagtccgggagatccacggtgacatttcgactatcaacggtgcagccgcctgttggcataacaacat                        <  1:374006‑M1/251‑209 (MQ=255)
                                                                                                                                                                                                              tctTACGTACTGCATTCGCTACGAATTTAACAGCACCCacacccgatgcttttgtttcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagatcttttgttcgcttgagcaatttacgccaaggtgtatttctgcgcttccatgaaagtccgggagatcctctgatacatttcgtctttcaacggtgcagccgcctgttggcatatcaacatt                       <  1:278912‑M1/246‑210 (MQ=255)
                                                                                                                                                                                                              gCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgaatcaggtgctgtgcaagcaaacacttgtcagctacagtcagtggtggcacccgaaagctaaaatgagcgattttgttcgcttgagcaatatacgccaaggggaatttctgcgcttccaggaaagtccgggaaatcaacggtgacatttcgactctcaacggtgcagccgcctgttggcataacaacattgtt                    >  1:27874‑M1/1‑38 (MQ=255)
                                                                                                                                                                                                                          cATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaaagtccgggagatccacggtgacatttcgactatcaacggtgcagccgcctgttggcataacaacattgttatttgaaataata       <  2:240389‑M1/251‑226 (MQ=255)
                                                                                                                                                                                                                               gCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgacagtggtggcacccgaaagataaaatgagagcttttgttcgctttagcaatatgcgccaaggggtatttctgcgcttccaggaaagtccgggagatccacggtgacatttcgactatcaacggtgctgccgcctgttggcataacaacattgttatttgaaataatattcca  <  1:76588‑M1/251‑231 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                      
AAGTAATGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGCACAGTACATACAAAGGGATTATCGCCAACATGGATATATAAGAGTAGGGGACTACTCTCTTTCTATTTTTAATCTATTATCGTTAAACAATAGATTAAATACGATAAGTTAAGCATCATCTCTCATAGTGTTTA  >  NC_000913/1586199‑1586672

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