Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 568,449 IS2 (–) +5 bp coding (135‑139/333 nt) emrE → multidrug/betaine/choline efflux transporter EmrE

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 568449 =0 (0.000)11 (1.090) 8/298 0.2 100% coding (135/333 nt) emrE multidrug/betaine/choline efflux transporter EmrE
?NC_000913 = 1469240 NA (NA)noncoding (1/1331 nt) IS2 repeat region
* ? NC_000913 = 5684530 (0.000)5 (0.500) 3/296 0.8 100% coding (139/333 nt) emrE multidrug/betaine/choline efflux transporter EmrE
?NC_000913 1467910 = NA (NA)noncoding (1331/1331 nt) IS2 repeat region

ATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAATTATTTGTTATTGTGCATCATTCTGGTTATTAGCTCAGACGCTGGCTTATATTCCTACAGGGATTGCTTATGCTATCTGGTCAGGAGTCGGTATTGTCCTGATTAGCTTACTGTCATGGGGATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTT  >  NC_000913/568233‑568658
                                                                                                                                                                                                                        |                                                                                                                                                                                                                 
gctacaccatgttgccgggcaacgagggagacagtcatccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaGTTATTAGCTCAGACGCTGGCTTATATTCCTAc                                                                                                                                                                                   <  1:91260‑M1/33‑1 (MQ=255)
             gatgaAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAATTATTTGTTATTGTGCATCATTCTGGTTATtggatttgcccctatatttccagacatctgttatcacttaan                                                                                                                                                                     <  1:160114‑M1/250‑43 (MQ=255)
                                       cccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtataaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaGTTATTAGCTCAGACGCTGGCTTATATTCCTACAGGGATTGCTTATGCTATCTGGTCAGGAGTCGGTATTGTnc                                                                                                                                          >  1:199576‑M1/178‑251 (MQ=255)
                                                           gTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGCGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGACGGTTTTACACGGTTATGGCCATCTGTTGGTACAATTATTTGTTATTGTGCATCATTCTGGTTATtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacc                                                                                                                      >  1:222145‑M1/1‑162 (MQ=255)
                                                                               ntgtggtacgccgtctgcgtttctccggccctaagacctcagtcatctgttctccactgaatagtctaaaaactagaattaagactatcacttatttaagtgatattggttgtctggagcttcagggggccagtctaGTTATTAGCTAAGACGCTGGCTTATATTCCTACAGGGATTGCTTATGCTATCTGGTc                                                                                                                                                           <  1:32358‑M1/57‑1 (MQ=255)
                                                                                                                aTACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTTTTGGTACAATTATTTGTTATTGTGCATCATTCTGGTTATtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggt                                                                 <  1:227155‑M1/251‑143 (MQ=255)
                                                                                                                                           ntagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaGTTATTAGCTCAGACGCTGGCTTATATTCCTACAGGGATTGCTTATGCTATCTGGTCAGGAGTCGGTATTGTCCTGATTAGCTTACTGTCATGGGGATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGGGTGTTGATTATTAAttt                                      <  1:152499‑M1/174‑1 (MQ=255)
                                                                                                                                                                   ctatcacttatttaagtgatattggttgtctggagatcaagggggccagtataGTTATTAGCTCAGACGCTGGCTTATATTCCTACAGGGATTGCTTATGCTCTCTGGTCAGGAGTCGGAATTGTCCTGATTAGCTTACTGTCATGGGGATTTTTCGGc                                                                                                          <  1:96937‑M1/106‑1 (MQ=255)
                                                                                                                                                                     atcacttatttaagtgatattggttgtctggagattcagggggccagtctaGTTATTAGCTCa                                                                                                                                                                                                        <  1:92512‑M1/12‑1 (MQ=255)
                                                                                                                                                                                aagtgatattggttgtctggagattcagggggccagtctaGTTATTAGCTCAGACGCTGGCTTATATTCCTACAGGGATTTCTTATGCTATCTAGTCAGGAGTCGGTATTGTCCTGATTAGCTTACTGGCATGCGGATGTGTCGGCCAACGGCTTGACCTGACAGCCATTATATCCATGATCTTGATTTGTGCCGGTGGGTTGATTATTAATGTATTTTCACGAAGCACACCACATTAAAATAAgttgtt  >  1:205922‑M1/41‑250 (MQ=255)
                                                                                                                                                                                 aTCTGTTGGTACAATTATTTGTTGTTGTGCATCATTCTGGTTATtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccn                                                                                                                                                       <  1:156732‑M1/100‑57 (MQ=255)
                                                                                                                                                                                 aTCTGTTGGTACAATTATTTGTTGTTGTGCATCATTCTGGTTATtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccg                                                                                                                                                       >  2:156732‑M1/1‑44 (MQ=255)
                                                                                                                                                                                  ntgatattggttgtctggagattcagggggccagtctaGTTATTAGCTCAGACGCTGGCTTATATTCCTACAGGGATTGCt                                                                                                                                                                         <  1:152743‑M1/43‑1 (MQ=255)
                                                                                                                                                                                  gtgatattggttgtctggagattcagggggccagtctaGTTATTAGCTCAGACGCTGGCTTATATTCCTACAGGGATTGCt                                                                                                                                                                         >  2:152743‑M1/39‑81 (MQ=255)
                                                                                                                                                                                         tggttgtcgtgagattcagggggccagtctaGTTATTAGCTCAGACGCTg                                                                                                                                                                                                 >  1:162709‑M1/32‑50 (MQ=255)
                                                                                                                                                                                                                ccagtctaGTTATTAGCTCAGACGCTGGCTTATATTCCTACAGGGATTGCTTATGCTATCTGGTCAGGAGTCGGTATTGTCCTGTTTAGCTTACTGTCATGGGGATTTTTCGGCCAACGGCTGGACCTg                                                                                           >  1:141318‑M1/9‑129 (MQ=255)
                                                                                                                                                                                                                        |                                                                                                                                                                                                                 
ATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAATTATTTGTTATTGTGCATCATTCTGGTTATTAGCTCAGACGCTGGCTTATATTCCTACAGGGATTGCTTATGCTATCTGGTCAGGAGTCGGTATTGTCCTGATTAGCTTACTGTCATGGGGATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTT  >  NC_000913/568233‑568658

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