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 34 [0] [0] 34 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)34 (0.490) 32/492 0.9 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

GCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGACGAAACGAAATACCGCGTGACAGCCATCACGCGGCAGACATTTTATTTTTTCTCGACGACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCACTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGG  >  NC_000913/4299281‑4299757
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                            
gCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaagggg                                                                                                                                                                                                                                     >  1:6429‑M1/1‑241 (MQ=255)
   cGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacag                                                                                                                                                                                                                                 <  1:253769‑M1/251‑14 (MQ=255)
       cGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaac                                                                                                                                                                                                                             >  1:95967‑M1/1‑234 (MQ=255)
        gATTTAGCAGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGTGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtcaggggacagtaact                                                                                                                                                                                                                            >  2:129196‑M1/1‑233 (MQ=255)
            ttGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgc                                                                                                                                                                                                                        >  1:433245‑M1/1‑229 (MQ=255)
            ttGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTGTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccg                                                                                                                                                                                                                         <  1:533899‑M1/250‑22 (MQ=255)
                           acCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacc                                                                                                                                                                                                         <  1:405881‑M1/251‑38 (MQ=255)
                                 ttGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgctggcggaattacctttttc                                                                                                                                                                                                   >  2:393278‑M1/1‑208 (MQ=255)
                                      tttGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatt                                                                                                                                                                                              >  1:573989‑M1/1‑203 (MQ=255)
                                        tGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttattt                                                                                                                                                                                             <  2:276180‑M1/250‑50 (MQ=255)
                                          cAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaatttcctttttcttatttcag                                                                                                                                                                                          <  1:107543‑M1/251‑53 (MQ=255)
                                           aaCCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagc                                                                                                                                                                                         <  2:4984‑M1/251‑54 (MQ=255)
                                            aCCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcg                                                                                                                                                                                        <  2:551486‑M1/251‑55 (MQ=255)
                                                gATGCGCAATATGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACTCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcc                                                                                                                                                                                    <  2:564282‑M1/251‑59 (MQ=255)
                                                   gcgcAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGCggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctgg                                                                                                                                                                                 >  1:303741‑M1/1‑190 (MQ=255)
                                                   gcgcAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctgg                                                                                                                                                                                 >  2:89801‑M1/1‑190 (MQ=255)
                                                                   ccGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaaca                                                                                                                                                                 <  2:356357‑M1/251‑78 (MQ=255)
                                                                          gacTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACCTTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAACCAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggagagaagaaatcaacaattgcgc                                                                                                                                                          >  1:421563‑M1/1‑167 (MQ=255)
                                                                            cTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccg                                                                                                                                                        >  1:302146‑M1/1‑165 (MQ=255)
                                                                                            cGCTCCAGTTAATCAAATCACGCATACGCGCTCGCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaa                                                                                                                                        <  2:101179‑M1/251‑103 (MQ=255)
                                                                                                              caCGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGTggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgt                                                                                                                      >  2:152189‑M1/1‑131 (MQ=255)
                                                                                                                   aTACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGAGTCATTCCCTGATATTACTGGCTCTATTTAGTCAAAACTCTGACGAAAAACAGGCGGTCGTCAGTTTTTACGTGAggtaaggggacagtaactccgcgatggcggaattacttttttcttattgcagcgcgggtggtgagaagaaatcaacaattgcgctggttttttccatattaatattagcttgttcaatgttgagc                                                                                                                 <  2:516995‑M1/251‑126 (MQ=255)
                                                                                                                                                          aTTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggt                                                                          >  2:537251‑M1/1‑87 (MQ=255)
                                                                                                                                                           ttttATGTGGGGTCCTTCCCTGATATTACGGGCACTATTTATTCACAACTCTGCCGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaatgggacagtaactccgcgatggcggaattacctttttcttatttcagcgagcctggtgagaagaaatcaaccattgcgccggttttttccatattaagattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggtt                                                                         >  2:166302‑M1/1‑86 (MQ=255)
                                                                                                                                                                       tCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttg                                                             >  2:171240‑M1/1‑74 (MQ=255)
                                                                                                                                                                        cATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgc                                                            <  1:89801‑M1/251‑179 (MQ=255)
                                                                                                                                                                                    ttACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacg                                                >  2:111305‑M1/1‑61 (MQ=255)
                                                                                                                                                                                               atttattCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacctaataacg                                     <  2:295170‑M1/251‑202 (MQ=255)
                                                                                                                                                                                                           ctctGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacctaataacggcgtttccgggg                         >  1:285417‑M1/1‑38 (MQ=255)
                                                                                                                                                                                                               gACGAAAAACAGGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattaccttttgcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacctaataacggcgtttccggggctg                      <  1:430851‑M1/250‑217 (MQ=255)
                                                                                                                                                                                                                      aaCACGCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatcttagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacctaataacggcgtttccggggctggtgagct               >  2:544250‑M1/1‑27 (MQ=255)
                                                                                                                                                                                                                          ggCTGTCGTCAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacctaataacggcgtttccggggctggtgagcaaataa          <  2:6429‑M1/251‑229 (MQ=255)
                                                                                                                                                                                                                              gtcgtcAGTTTTGACGTGAggtaagggtacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgatgcattacgacctaataacggcgtttccggtgctggtgagcaaataatatg      <  2:549969‑M1/251‑233 (MQ=255)
                                                                                                                                                                                                                                  tcAGTTTTGACGTGAggtaaggggacagtaactccgcgatggcggaattacctttttcttatttcagcgcgcctggtgagaagaaatcaacaattgcgccggttttttccatattaatattagcttgttcaatgttgagctgggcgacagagacaaaaaaagcgtccagcaaggttggttgtaatattcttgccgaggcattacgacctaataacggcgtttccggggctggtgagcaaataatatgatcg  >  1:185317‑M1/1‑15 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                            
GCCCGGACGATTTTGCCGTTATCGACGACCAGGTTGATTTTGCAACCTGATGCGCAATAGGGGCAAACCGTGACGACTTTTTTCATCGGTCTCGCTCCAGTTAATCAAATCACGCATACGCGCTCTCGACTACAGTATGCATCTTTTATGCCACATTTTATGTGGGGTCATTCCCTGATATTACGGGCACTATTTATTCAAAACTCTGACGAAAAACAGGCTGTCGTCAGTTTTGACGTGACGAAACGAAATACCGCGTGACAGCCATCACGCGGCAGACATTTTATTTTTTCTCGACGACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCACTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGG  >  NC_000913/4299281‑4299757

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