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 31 [0] [0] 31 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)31 (0.500) 28/494 0.7 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

CCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGACGAAACGAAATACCGCGTGACAGCCATCACGCGGCAGACATTTTATTTTTTCTCGACGACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCACTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTG  >  NC_000913/4299296‑4299755
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                          
ccGTTATCGACGACCAGGTTGATTTTGCAACCTTATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCTCATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgat                                                                                                                                                                                                                   <  1:56951‑M1/251‑26 (MQ=255)
  gTTATCGACGACCAGGTTGATTTTGCATCCTGATGCCCAATAGGGGCCAACCGTGACGACTTTTTTCATCTGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATTTTACGTGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtagctccgcgatgg                                                                                                                                                                                                                 <  2:250896‑M1/251‑28 (MQ=255)
    tATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggc                                                                                                                                                                                                                <  1:135715‑M1/250‑29 (MQ=255)
    tATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggc                                                                                                                                                                                                                <  1:425583‑M1/250‑29 (MQ=255)
        gacgacCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGACCCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCAATATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGCCAGTTTTGCCGTGAggcaaggggacaggaactccgcgatggcggaat                                                                                                                                                                                                           >  1:106191‑M1/1‑218 (MQ=255)
          cgacCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaatt                                                                                                                                                                                                          >  1:469854‑M1/1‑216 (MQ=255)
              cAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacatt                                                                                                                                                                                                     <  2:32655‑M1/251‑40 (MQ=255)
                   tGATTTTGCAACCTGATGCGCAATAGGGTCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGTCGTGAggtaaggggacagtaactccgcgatggcggaattcactttttc                                                                                                                                                                                                 <  2:426056‑M1/250‑44 (MQ=255)
                    gATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttct                                                                                                                                                                                                <  1:394059‑M1/250‑45 (MQ=255)
                                                                  ttCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttt                                                                                                                                                  <  1:48587‑M1/250‑91 (MQ=255)
                                                                     aTCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttc                                                                                                                                              >  1:95293‑M1/1‑157 (MQ=255)
                                                                                  cAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggaaaggggacagtaactccgcgctggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaaccattgcgccggttttttccatattaatatt                                                                                                                                  >  2:472752‑M1/1‑144 (MQ=255)
                                                                                           aaatcaCGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttc                                                                                                                         >  1:90863‑M1/1‑135 (MQ=255)
                                                                                                  gCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCCCATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttga                                                                                                                 <  2:88391‑M1/251‑124 (MQ=255)
                                                                                                  gCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttg                                                                                                                  >  1:261570‑M1/1‑128 (MQ=255)
                                                                                                               cGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGGCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgaca                                                                                                      <  2:90863‑M1/249‑135 (MQ=255)
                                                                                                                                 tttATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTAAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtcca                                                                                   <  1:387039‑M1/250‑154 (MQ=255)
                                                                                                                                    aTGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagca                                                                                >  1:248072‑M1/1‑94 (MQ=255)
                                                                                                                                       caacattTTATGTGGGGTCATTCCCTGATATTACGGGCTCTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGATTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttctccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgaccggtaaggt                                                                            <  2:469854‑M1/249‑161 (MQ=255)
                                                                                                                                            ttttATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggtt                                                                       >  1:9978‑M1/1‑86 (MQ=255)
                                                                                                                                                     gggTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttatccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattc                                                              <  2:179593‑M1/251‑175 (MQ=255)
                                                                                                                                                                    atTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTAGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattac                                               >  1:379304‑M1/1‑62 (MQ=255)
                                                                                                                                                                        cGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagatgggcgacagagacaaaaaaagcgtccagcaaggttggttgtactattcttgccgaggcattacgacc                                           >  2:146458‑M1/1‑58 (MQ=255)
                                                                                                                                                                                   tattCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacctaataacggcg                                >  2:23639‑M1/1‑47 (MQ=255)
                                                                                                                                                                                    agtcaAAACTCGGACTAAACGCAGGCTGTCGTCAGTTTTGACGTGAggtaatgggacagtaactccgcgatggcggaattacctttttctgatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctggtcgacagagacaaaaaaagcgtccagcaagggtggttggaatattcttgccgaggcagttcgacctaataacggcgt                               <  2:170962‑M1/249‑206 (MQ=255)
                                                                                                                                                                                       cAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacctaataacggcgttt                             <  2:248072‑M1/250‑208 (MQ=255)
                                                                                                                                                                                           aCTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacctaataacggcgtttccggg                        >  1:36174‑M1/1‑39 (MQ=255)
                                                                                                                                                                                                       aaCAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacctaataacggcgtttccggggctggtgagcaa            <  2:261570‑M1/251‑225 (MQ=255)
                                                                                                                                                                                                               gtcgtcAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacctaataacggcgtttccggggct                     >  1:57714‑M1/1‑19 (MQ=255)
                                                                                                                                                                                                               gtcgtcAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacctaataacggcgtttccggggct                     <  2:57714‑M1/234‑216 (MQ=255)
                                                                                                                                                                                                                 cgtcAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacctaataacggcgtttccggggctggtgagcaaataatatgat  >  2:102271‑M1/1‑17 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                          
CCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGACGAAACGAAATACCGCGTGACAGCCATCACGCGGCAGACATTTTATTTTTTCTCGACGACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCACTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTG  >  NC_000913/4299296‑4299755

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.