Predicted mutation
evidence seq id position mutation annotation gene description
JC NC_000913 4,001,644 Δ5 bp coding (220‑224/951 nt) corA → magnesium/nickel/cobalt transporter

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 40016431 (0.020)51 (1.100) 43/284 0.1 98.1% coding (219/951 nt) corA magnesium/nickel/cobalt transporter
?NC_000913 4001649 = 1 (0.020)coding (225/951 nt) corA magnesium/nickel/cobalt transporter

CAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCG  >  NC_000913/4001501‑4001783
                                                                                                                                              |                                                                                                                                            
cAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGACAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGCTTCTCTGAAGACGACGACGGCCTGCat                                                                                                                                              >  2:742663/1‑143 (MQ=255)
    aTGGATTGATCTTGTCGCACCGGACGACGCCGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcc                                                                                                                                          <  2:658365‑M1/143‑5 (MQ=255)
       gcttgaTCTTGTCGAACCGGACCACGACGAGCCACTTCGCGTACAATCTGAACTTGGCCAGCGCATGGCAACCCGCCCGGAACTGGAAGACATCGACGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttc                                                                                                                                       <  2:647652‑M1/141‑8 (MQ=255)
          tggtgctggtgacgcgggcgatgacgagcGAGTGCGCGTAGAATCTGAACTTGGCCAGAGCGTGGCAACCCGCCCGGAACTGGAAGACAGCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttc                                                                                                                                    <  2:393790‑M1/129‑11 (MQ=255)
           gATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttct                                                                                                                                   >  1:398885‑M1/1‑132 (MQ=255)
                 gTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttg                                                                                                                             <  1:357216‑M1/143‑18 (MQ=255)
                 gTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGCCACGTTGCTTTGCAGACGACGACGGCCTGCCGcgcctgcttctccttcg                                                                                                                             >  2:597948‑M1/1‑126 (MQ=255)
                           tcgacgacgagCTGCTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGTCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcgga                                                                                                                   <  2:454644‑M1/142‑28 (MQ=255)
                            cgacgacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGTGCCTGGCAACCCGCCCGGAACTCGAAGACATCGAAGCATCGGCACGTTTCTTTGACGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaa                                                                                                                  <  1:35582‑M1/143‑29 (MQ=255)
                            cgacgacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaa                                                                                                                  <  2:237818‑M1/143‑29 (MQ=255)
                            cgacgacgaCCGACTGCCCGTACCATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaa                                                                                                                  <  1:387439‑M1/143‑29 (MQ=255)
                              acgacgaGCTACTGCGCGTACGAGCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaaga                                                                                                                <  1:24158‑M1/143‑31 (MQ=255)
                                gacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCCTCGGCACGTGTCTTTGAAGACGACGCCGGCCTGCATccccttcttcctccttgacgatgcggaagatc                                                                                                              >  2:430770‑M1/1‑111 (MQ=255)
                                 acgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCCGCATctccttcttcttctttgaagatgcggaagatca                                                                                                             <  2:582390‑M1/143‑34 (MQ=255)
                                   gaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGCAGACATCGAAGCACCGGCACGTTTCTTTGAAGCCGACGACGGCCTGCATctcctccttcttctttgaagatgcggaagatcacg                                                                                                           >  2:558683‑M1/1‑108 (MQ=255)
                                    aGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgc                                                                                                          <  1:751932‑M1/143‑37 (MQ=255)
                                     gcgcgtgcgcgTACGATCTGAACGTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcatgcc                                                                                                         <  2:733492‑M1/138‑38 (MQ=255)
                                     gCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgcc                                                                                                         >  1:131968‑M1/1‑106 (MQ=255)
                                      cGACTGCGCGTACAAGCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccg                                                                                                        <  1:554471‑M1/143‑39 (MQ=255)
                                          tGCGCGTACAATCTGAGCTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaa                                                                                                    <  2:736215‑M1/143‑43 (MQ=255)
                                                 tggatctGATCTTCGACAGGGCCTGGCAACGCGCGCGGAACTGGAAGACATGGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgtggaagatcacgccggtaactccact                                                                                             <  2:4913‑M1/140‑50 (MQ=255)
                                                   caTCTGAACTCGGCCAGAGCCTGGCAACCCGCCCGGAACGGGGAGACATCGAAGCCTCGGTCCGTTTCATTGTATTG‑‑GATGGCATCATTTGTCTCCTTCTTCTTCTTTGAAGACGCGGTATATCACGCCGGTAACTCCACtgt                                                                                      <  2:784378/142‑1 (MQ=255)
                                                     tCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactcgactgtgg                                                                                         >  1:487185‑M1/1‑90 (MQ=255)
                                                      cgaaactcGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcgtctgtgaaaatcatcaacatcacgccggtaactccactgtggc                                                                                        <  2:99277‑M1/140‑55 (MQ=255)
                                                               gCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagaggcggaagctcacgccggcaactccactgtggcagttaccat                                                                               >  2:614384‑M1/1‑80 (MQ=255)
                                                                   gagCCGTGCTACCCGCCGGGAACTGGAATGCGTCTGAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgt                                                                           <  2:127738‑M1/143‑68 (MQ=255)
                                                                     gCCTGGCAAGCCGGTCTGAATTGGAAGACATCGAAGCATGGGCACGTTTCTTTGGAGACGACGACGGCCTGCATctccgtcttcttctttgaagatgcggaagatgacgccggtaactccactgtggcatttaccatccgtga                                                                         <  2:81320‑M1/143‑70 (MQ=255)
                                                                       gtgGCACGCCGCCCGGAGCTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcctgcttcttctttgaagatgcggaagatcacgccggtaactccacggtggcatttaccatccgtgatg                                                                       <  2:359291‑M1/142‑72 (MQ=255)
                                                                       cTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatg                                                                       >  1:680134‑M1/1‑72 (MQ=255)
                                                                       cTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagacgcggaatatcacgccggtaactccactgtggcatttaccatccgtgatg                                                                       <  1:765287‑M1/143‑72 (MQ=255)
                                                                          gCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtc                                                                    <  1:598323‑M1/143‑75 (MQ=255)
                                                                           cAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatcacggatggtaa                                                                   <  2:595458‑M1/143‑76 (MQ=255)
                                                                                cgcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgt                                                              <  1:337657‑M1/143‑81 (MQ=255)
                                                                                       aaCTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactct                                                       <  1:17456‑M1/143‑88 (MQ=255)
                                                                                          tGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttctgcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcg                                                    <  2:435981‑M1/143‑91 (MQ=255)
                                                                                           ggaggaGATCGAAGCGTCGGTATGCTTGTTCGAAGACGACGAGGGCCTGCACctccttcttcttctttgaagatgcggaagatcatgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgt                                                   <  2:343187‑M1/143‑92 (MQ=255)
                                                                                           ggAAGACATCGAAGCATCTGCACGTTTCTTTGAAGACGACGCCGGCCTGCATctccttctgcttcttggaagatgcggaatatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgt                                                   <  2:550685‑M1/143‑92 (MQ=255)
                                                                                           ggAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcattcaccatccgtgatggtcgtggatttactctgcgt                                                   >  2:378686‑M1/1‑52 (MQ=255)
                                                                                            gAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtg                                                  >  1:103849‑M1/1‑51 (MQ=255)
                                                                                            gAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtg                                                  >  2:292404‑M1/1‑51 (MQ=255)
                                                                                             cagACACTGAAGCATCGGCACTTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtga                                                 <  2:752020‑M1/142‑94 (MQ=255)
                                                                                             aaGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtga                                                 >  1:138083‑M1/1‑50 (MQ=255)
                                                                                                  ttggacgggtcggcacgttTCTTTGATGACGACGCCGGCCTGCATctccttcttcttctttgaagatgcggaagagcacgccggtaactccactgtggcatctaccatccgccatggtcgtctgtttactctgcgtgagcgtg                                            <  2:703671‑M1/134‑99 (MQ=255)
                                                                                                     gAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcattcaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaac                                         >  2:582987‑M1/1‑42 (MQ=255)
                                                                                                        gCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcctgcgcgtgagctgc                                      >  2:415033‑M1/1‑39 (MQ=255)
                                                                                                            cGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggatgttcccgcccgtccctccactgtggcatttcccattcttgatggtagactgtttactctgcgtgagtgtgaactgaccga                                  >  1:64814‑M1/1‑35 (MQ=255)
                                                                                                                cacgtcctgggaagacgacgaCGGCCTGCATctccgtcttcttcgttgaagatgcggaagatcacgccggtaactccgctgtggcatttaccatccgtgatggtcggctgtttactctgcgtgagcgtgaactgcccgctttt                              <  1:620219‑M1/133‑113 (MQ=255)
                                                                                                                aCGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgctttt                              <  1:262808‑M1/143‑113 (MQ=255)
                                                                                                                   tttCTTTGAAGACGACGCCGCCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgt                           <  1:154534‑M1/143‑116 (MQ=255)
                                                                                                                   tttCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtgactcccctgtggcatttcccctccgtgatggtcgtctgtttactctgcgtgagcgtgacctgcccgcttttcgt                           >  1:203558‑M1/1‑28 (MQ=255)
                                                                                                                       tttGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcctttaccatccgtgatggtcgtctgtctactctgcccgagcgtgacctgcccgcttttcgtccgt                       >  1:742275‑M1/1‑24 (MQ=255)
                                                                                                                        ttGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgcctgtttactctgcgtgagcgtgaactgcccgcttttcctctgtg                      >  2:653858‑M1/1‑23 (MQ=255)
                                                                                                                                            cATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgagatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgcagcc  >  2:494159‑M1/1‑3 (MQ=255)
                                                                                                                                              |                                                                                                                                            
CAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCG  >  NC_000913/4001501‑4001783

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