Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,299,522 Δ12,535 bp mdtPalsB mdtP, mdtO, mdtN, ytcA, yjcS, alsK, alsE, alsC, alsA, alsB

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4299522 4312056 12535 22 [0] [0] 22 mdtP–alsB mdtP,mdtO,mdtN,ytcA,yjcS,alsK,alsE,alsC,alsA,alsB

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 42995210 (0.000)22 (0.720) 20/486 0.5 100% intergenic (‑155/+43) fdhF/mdtP formate dehydrogenase H/putative multidrug efflux pump outer membrane channel
?NC_000913 4312057 = 0 (0.000)intergenic (‑15/+44) alsB/alsR D‑allose ABC transporter periplasmic binding protein/DNA‑binding transcriptional repressor AlsR

ACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCT  >  NC_000913/4311823‑4312305
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        
cgattttgccgttatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGgacagtaac                                                                                                                                                                                                                                          <  1:303532‑M2/17‑5 (MQ=255)
        ccgttatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccatgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacatgttatgtggggtcattccctgagattacgtgcgctatgtattcaaaactctgacgaaaaacaggctgtcggcaggttttacgtgaGGTAAGGGGACAGTAACTCCGCGAt                                                                                                                                                                                                                                 <  2:301196‑M2/25‑1 (MQ=255)
                         gttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattcccggatattacgggcactattgattcaaaactcggacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGTGACAGTAACTCCGCGATGGCGGAATTGCCttttt                                                                                                                                                                                                                <  2:111260‑M2/42‑1 (MQ=255)
                                 tgcgacccgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAttt                                                                                                                                                                                                         <  2:137358‑M2/49‑1 (MQ=255)
                                          atgcgcactagggacaaaccgtgacgacttttttcatcggtctcgctccagttaatcaactcacgcctacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCt                                                                                                                                                                                               >  2:129896‑M2/193‑251 (MQ=255)
                                                   aggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgaga                                                                                                                                                                                        >  2:73571‑M2/184‑249 (MQ=255)
                                                             cgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaaca                                                                                                                                                                             >  1:149233‑M2/174‑250 (MQ=255)
                                                                gacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactattgattcaaaactctgacgaaaaacaggctggcgtcagttttgacgtgaGGTAAGGGGACAGTAACGCCGCGATGGCGGAATTACCTGTTTCTTATTTCAGCGCGTCTGGTGAGAAGAAATCAACAATTg                                                                                                                                                                         <  2:115090‑M2/81‑1 (MQ=255)
                                                                                      gctccagttaatcaaatcacgcctacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCatattaat                                                                                                                                                   >  2:284282‑M2/149‑251 (MQ=255)
                                                                                                        acgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTGGTTCAATGtt                                                                                                                                 <  1:88040‑M2/121‑1 (MQ=255)
                                                                                                         cgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCtt                                                                                                                                                                                                             <  1:195511‑M2/45‑1 (MQ=255)
                                                                                                         cgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCtt                                                                                                                                                                                                             >  2:195511‑M2/130‑174 (MQ=255)
                                                                                                          gcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCTATGGCGGAATTACCTTTGTCTGATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGt                                                                                                                                                                  <  1:276973‑M2/88‑1 (MQ=255)
                                                                                                          gcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGCATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGt                                                                                                                                                                  >  2:276973‑M2/129‑216 (MQ=255)
                                                                                                               cgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTgg                                                                                                                          <  1:163985‑M2/128‑1 (MQ=255)
                                                                                                                       cgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGAc                                                                                                                     <  2:149233‑M2/133‑1 (MQ=255)
                                                                                                                                                                 cattccctgatcttacgggcactatttattcaaaactctgacgaaacacaggctgtcgccagttttgacgtgaGGTAAGGGGCCAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTACTATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGc                                                                        >  2:276852‑M2/74‑251 (MQ=255)
                                                                                                                                                                                        atttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACg                                                 >  2:73132‑M2/51‑251 (MQ=255)
                                                                                                                                                                                            attcaaaactctgacgaaaaacaggctgtcgtcagtttagacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTGGGTTGTAATATTCGTGCCGAGGCATTACGACGTAataacttcg                                              <  1:278356‑M2/204‑5 (MQ=255)
                                                                                                                                                                                              tcaaaactctgacgaaaaacaggctgtcggcagttttgacgtgaGGCAAGGGGACAGTAACTCCGCGATGGCGGCATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGCAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGATCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTACTATTCTTGCCGAGGCATTACGACCTAATAACGGCGttt                                           >  1:142042‑M2/45‑251 (MQ=255)
                                                                                                                                                                                                                             agttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTATCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAAGAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGc            <  1:284282‑M2/238‑1 (MQ=255)
                                                                                                                                                                                                                                        gaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCCTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACTGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCt  <  1:188855‑M2/248‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        
ACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCT  >  NC_000913/4311823‑4312305

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: