Predicted mutation
evidence seq id position mutation annotation gene description
JC NC_000913 4,001,645 Δ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 = 40016441 (0.020)54 (1.020) 45/284 0.2 98.2% coding (219/951 nt) corA magnesium/nickel/cobalt transporter
?NC_000913 4001650 = 1 (0.020)coding (225/951 nt) corA magnesium/nickel/cobalt transporter

AGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATG  >  NC_000913/4001503‑4001776
                                                                                                                                             |                                                                                                                                    
aGTGTGGATTGAGCGTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATc                                                                                                                                     <  1:620684‑M1/143‑2 (MQ=255)
 gTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATct                                                                                                                                    >  2:772645‑M1/1‑141 (MQ=255)
   aTGGAGTGGTCTTGGCGAACGGGACGACGACGAGCGACTGCGCGTACAATCTGAACTGGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcc                                                                                                                                  <  1:447447‑M1/143‑5 (MQ=255)
     ggtctgatcGTGTCGCATCGGACGGCGACGCGCGACTGCGCGTACAATCGGAACTTGGCCGGAGCCTGGCAACCCGCCCGGCACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcctt                                                                                                                                <  2:129849‑M1/139‑7 (MQ=255)
           aTCTTGTCGAACCGGACGACGACGAGCGGCTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttctt                                                                                                                          >  1:408648‑M1/1‑131 (MQ=255)
              ttGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGCAGCATCGGCACGTTTCTTTGAAGACGACGCCGGCCTGCATccccttcttcttctt                                                                                                                       >  1:792063‑M1/1‑128 (MQ=255)
                           cgacgacggGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaa                                                                                                          <  2:400402‑M1/143‑29 (MQ=255)
                               gacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTCTGAAGACGACGACGGCCTGCATctccttcttctcctttgaagatgcggaagacc                                                                                                      >  2:153976‑M1/1‑111 (MQ=255)
                                    gCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgcc                                                                                                 <  2:174188‑M1/143‑38 (MQ=255)
                                       aCTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGTATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggt                                                                                              <  1:72341‑M1/143‑41 (MQ=255)
                                        cTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCAGCACGGAACTGGAAGACATCGAAACATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagctgcggaagatctctccggtc                                                                                             >  1:205523‑M1/1‑102 (MQ=255)
                                         ttcgcgtaCAATCTGAACTCCGCCAGCGCCTGGCAACCCGCCCGGAACTGGAAGCCCTCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaa                                                                                            <  1:108681‑M1/141‑43 (MQ=255)
                                           cgcgTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtcact                                                                                          >  2:146585‑M1/1‑99 (MQ=255)
                                                   aTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtg                                                                                  <  2:743883‑M1/143‑53 (MQ=255)
                                                   aTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtg                                                                                  <  2:937897‑M1/143‑53 (MQ=255)
                                                      tGAAGTGGGCCAGAGCCTGGCAACCCGCCCGGCACTGGAGGTCGTGGTCGCAGCGGCACGTTTCTTTGAACTCGACGACGGCCTGCATctccttctgcttctttgaagatgcggaagatcacgccggtaactccactgtggca                                                                               <  1:904600‑M1/143‑56 (MQ=255)
                                                         aCTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGAGACCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtacctccactgtggcattt                                                                            <  2:936746‑M1/143‑59 (MQ=255)
                                                         aCTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGCCGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcattt                                                                            <  1:47634‑M1/143‑59 (MQ=255)
                                                          cTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcattta                                                                           <  1:729728‑M1/143‑60 (MQ=255)
                                                               ccAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatc                                                                      >  1:863203‑M1/1‑79 (MQ=255)
                                                                    gCCTGGCAGCTCGCCCGGCACTGGAAGACATCGAAGCATCGTCACGTTTCTTTGAAGACGACGACGGCCTGCATgtcctgcttcttctttgaagatgcggaagatcacgccggtaaccccactgtggcatttaccatccgtga                                                                 <  2:242026‑M1/143‑70 (MQ=255)
                                                                     ccTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCAGCGGCACGGTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgat                                                                <  2:401963‑M1/143‑71 (MQ=255)
                                                                      cTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcctcttctttgaagatgcggaagatcacgccggtaactccactatggcatttaccatccgtgacg                                                               >  1:760183‑M1/1‑72 (MQ=255)
                                                                        ggCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggt                                                             >  2:791833‑M1/1‑70 (MQ=255)
                                                                         gCAACCCTCCCGGAACTGGAAGACATCGAAGCGTCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttctgcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtc                                                            <  2:301731‑M1/143‑75 (MQ=255)
                                                                              ccgcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagattacaatgtgatctccactgtggcatttaccatccgtgatggtcgtctg                                                       >  1:944199‑M1/1‑64 (MQ=255)
                                                                              ccgcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctg                                                       <  1:208378‑M1/143‑80 (MQ=255)
                                                                              ccgcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctg                                                       >  2:652391‑M1/1‑64 (MQ=255)
                                                                               cgcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgt                                                      <  2:841058‑M1/143‑81 (MQ=255)
                                                                                gcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtt                                                     <  2:449305‑M1/143‑82 (MQ=255)
                                                                                 cccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgttc                                                    >  2:558803‑M1/1‑61 (MQ=255)
                                                                                  tcgGAACTTGAATACATCGATGCATAGGCACGTTTCTTTGAAGACGATGACTGCTTGCAT‑‑‑‑‑CTCCTTGTTATTTTTTTAAGATGCGGAAGATCACGCCGGTAACTCTACTTTGGCATTTACTATCCGTGATGGTCGTCTGTTTa                                              <  2:463122/142‑1 (MQ=255)
                                                                                    ggAACGTGAAGACATCGAAGCATCGGCGCGTTTCTTTGAAGACGGCGACGGGCTGCATctccttcttcttctttgaagatgcggaagatcccgccggtaactccactgtggcatttaccatccgtggtggtcgtctttttact                                                 <  2:388988‑M1/143‑86 (MQ=255)
                                                                                       aCTGGAAGACATCGAAGCATCGGCACGATTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctg                                              >  2:527295‑M1/1‑55 (MQ=255)
                                                                                        cTGGAAGACAGCGAAAAACCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgc                                             <  1:253003‑M1/143‑90 (MQ=255)
                                                                                           gAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtg                                          >  1:453153‑M1/1‑51 (MQ=255)
                                                                                             aGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgag                                        >  1:356596‑M1/1‑49 (MQ=255)
                                                                                              gACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagagcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagc                                       <  1:203611‑M1/143‑96 (MQ=255)
                                                                                                 aTCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGTCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtg                                    <  1:839628‑M1/143‑99 (MQ=255)
                                                                                                  tCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtcactccactgtggcatttaccatccgtgatggtcgcctgtttactctgcgtgagcgtga                                   >  2:696004‑M1/1‑44 (MQ=255)
                                                                                                   cGACGCACCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactttggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaa                                  <  1:229819‑M1/143‑101 (MQ=255)
                                                                                                   cGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaa                                  <  1:772645‑M1/143‑101 (MQ=255)
                                                                                                   cGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaa                                  <  1:746941‑M1/143‑101 (MQ=255)
                                                                                                    gAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaac                                 <  1:676378‑M1/143‑102 (MQ=255)
                                                                                                             gCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatgggcgtctgtttactctgcgtgagcgtgaactacccgctt                        >  2:738332‑M1/1‑33 (MQ=255)
                                                                                                                 gTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccacggtggcatttaccatccgtgatggtcgtctgtttactctgcgtgggcgtggactgcccgcttttcg                    >  1:938486‑M1/1‑29 (MQ=255)
                                                                                                                 gTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagagcacgccggtaactcccctgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcg                    >  1:803242‑M1/1‑29 (MQ=255)
                                                                                                                   ttCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtc                  >  1:47221‑M1/1‑27 (MQ=255)
                                                                                                                   ttCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtc                  >  1:49394‑M1/1‑27 (MQ=255)
                                                                                                                    tCTTTGAAGGCGACGACGGCCTGCATctccttctgcttctttgatgatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtct                 <  2:234260‑M1/143‑118 (MQ=255)
                                                                                                                         gAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatc            >  1:786253‑M1/1‑21 (MQ=255)
                                                                                                                               gacgactGCTTGCATcgccttcttcttcggtggggatgcggaagatcaggccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactcggcgtgagcgtgaactgtccgcttttcgtctgtatcgtatgc      <  2:377453‑M1/143‑129 (MQ=255)
                                                                                                                               gacgacgGCCTGCGActcctgcttctcctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgc      <  1:146585‑M1/143‑129 (MQ=255)
                                                                                                                               gacgacgGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttactatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgc      <  1:746965‑M1/143‑129 (MQ=255)
                                                                                                                                   acgGCCCGCAGctccctcttcttctttgaagatgcggaagatcacgccggtaactccactgtgccattcaacctccgttagggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgc  <  1:854908‑M1/143‑133 (MQ=255)
                                                                                                                                             |                                                                                                                                    
AGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATG  >  NC_000913/4001503‑4001776

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