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 = 40016432 (0.030)43 (0.650) 37/284 1.0 95.6% coding (219/951 nt) corA magnesium/nickel/cobalt transporter
?NC_000913 4001649 = 2 (0.030)coding (225/951 nt) corA magnesium/nickel/cobalt transporter

GGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGT  >  NC_000913/4001507‑4001791
                                                                                                                                              |                                                                                                                                              
agtatgggttgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgccgacggcctgcatc                                                                                                                                                <  2:410063‑M2/1‑1 (MQ=255)
   atggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTcc                                                                                                                                             <  2:878070‑M2/4‑1 (MQ=255)
      gattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggcccgatgctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttc                                                                                                                                          >  1:17898‑M2/137‑143 (MQ=255)
      gattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttc                                                                                                                                          >  2:532202‑M2/137‑143 (MQ=255)
            tcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttc                                                                                                                                    >  2:215246‑M2/131‑143 (MQ=255)
            tctggtggaaccggacgacgacgcgcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttc                                                                                                                                    <  2:222698‑M2/13‑1 (MQ=255)
              tggtcgaaccgtacgacgccgagcgcctgcgcgtacaagctgaactttgccagagcctggcaccccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttctt                                                                                                                                  <  1:485760‑M2/15‑1 (MQ=255)
                gtcgaatcacacgacgacgagcgactgctcgtacaatctgaacttgtccatagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTg                                                                                                                                <  2:336622‑M2/17‑1 (MQ=255)
                 tcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCTTTCTTCTTCTTTGa                                                                                                                               >  2:217558‑M2/126‑143 (MQ=255)
                 tcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGa                                                                                                                               >  2:418098‑M2/126‑143 (MQ=255)
                 gcgagccggacgacgacgggcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaacgggaagacagcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGa                                                                                                                               <  1:1037688‑M2/18‑1 (MQ=255)
                 gcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGa                                                                                                                               <  1:650826‑M2/18‑1 (MQ=255)
                   gaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTCCTTCCTCGAAg                                                                                                                             >  1:613102‑M2/124‑143 (MQ=255)
                                   ggcgactgggcgtacaaactgaccttgggcagagccggggaagccgcgcggaactggaagacatcgaagcatcggcaggtttctttgaagacgacgacggcctgcatCTCCGTCTTCTTCTTTGAAGATGCGGAAGATCACGc                                                                                                             <  1:158480‑M2/36‑1 (MQ=255)
                                     ctgcggcgcgtacaatctgaacttggccagggcctggcaacccgcccggaactggaagacatcgaagcttcggcacgtttctttgaagacgacgacggcctgcatCTCCTGCTTCTTCTTTGAAGATGCGGAAGATCACGCCg                                                                                                           <  1:440231‑M2/38‑1 (MQ=255)
                                      gactgcgcgtacaatctgaacctggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCgg                                                                                                          <  1:983517‑M2/39‑1 (MQ=255)
                                             cgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgcagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTCACTcc                                                                                                   >  2:218523‑M2/98‑143 (MQ=255)
                                                acaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACt                                                                                                <  2:772332‑M2/49‑1 (MQ=255)
                                                  aatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGCCGATGCGGAAGATCACGCCGGGAACCCCACtgt                                                                                              >  2:839782‑M2/93‑143 (MQ=255)
                                                   agctgaacgtggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCCTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgtg                                                                                             <  2:457799‑M2/52‑1 (MQ=255)
                                                     ctgaacttggccagggcctggcaacccgcccggaactggaagacagcgaagcatcggcacgtttctttgcagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGc                                                                                           <  1:46975‑M2/54‑1 (MQ=255)
                                                      tgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTCCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCa                                                                                          <  1:356471‑M2/55‑1 (MQ=255)
                                                          tctcgcgcagggcctggcacccggcccggaaccggaagacatcgaagcatcggcacgttctttggagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTa                                                                                      <  1:796373‑M2/59‑1 (MQ=255)
                                                           ttggctagagtctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTAc                                                                                     <  1:301462‑M2/60‑1 (MQ=255)
                                                               ccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATc                                                                                 >  2:441725‑M2/80‑143 (MQ=255)
                                                                ccTGGCAACCCGCCCGGAACTGGAAGACATAGAAGAATCGGCACTTTTCTTTGAAGACGACGACGGCCTGCATTTT‑‑‑‑‑CGTCGTCTTGGTTGAAGTTGCGGAAGATCACGCCGTGAACTCGACTGTGGCAGTTACTATCCGTGAt                                                                           >  1:21175/1‑143 (MQ=255)
                                                                    gactcgcaacccgcccggaacgggaagacatcgcagcatcggcacgtttctttgaagacgacgacggcctgcatCACCTTCTGCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACGCCACTGTGGCATTTACCATCCGTGa                                                                            <  2:728178‑M2/69‑1 (MQ=255)
                                                                      ctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACAGTGCCATTTACCATCCGTGATg                                                                          <  1:414951‑M2/71‑1 (MQ=255)
                                                                         gcaacccgcccggaactggaagacatcgcagcatcggcacgtttctttgaagacgacgacggcctgcatGTCCTTCTTCTTCTTCGAAGATGCGGAAGATCACGCCGGTAACTCCACCGTGGCATGCACCAACCGCGatggcc                                                                       >  1:811401‑M2/70‑141 (MQ=255)
                                                                             gccggcgggaacgggaatgcatcgaagcctcggcacgtttctttgacgacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCTACTGTGGCATTTACCATCCGTGATGGTCGTCt                                                                   <  1:44713‑M2/78‑1 (MQ=255)
                                                                                gcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCgtcttct                                                                >  2:1110882‑M2/63‑140 (MQ=255)
                                                                                     gcaaGCCAGCGACGCGTCTGGACGTTTCGTGGAGGACGACGACGGCCTGCGT‑‑‑‑‑CTCCTTCTTCTTCTTTGAAGATGCGGATGCTCACGCCGGTAACTCCACTGTGGCATTTACCAACCGTGATGGTCGTCTGTTGACTCTgcgt                                                      <  1:368956/141‑1 (MQ=255)
                                                                                            aaggcatgcaagcatcggaacggttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGACGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtga                                                    <  2:377058‑M2/93‑1 (MQ=255)
                                                                                                     aagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGACGATCACGCCGGTAACTCCACTGTGGCATTTACCATTCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACt                                           >  2:800359‑M2/42‑143 (MQ=255)
                                                                                                     aagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACt                                           <  1:418098‑M2/102‑1 (MQ=255)
                                                                                                           cggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGc                                     >  2:313910‑M2/36‑143 (MQ=255)
                                                                                                            ggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTGTTTGAAGATGCGGAAGATGACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCt                                    >  2:825901‑M2/35‑143 (MQ=255)
                                                                                                              cacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCttt                                  >  1:396987‑M2/33‑143 (MQ=255)
                                                                                                              cacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACAGTGGCATTTAACATCCGTGATGGTCGTCCGTTTACTCTGCGTGAGCGTGACCTGCCCactgt                                  >  1:801893‑M2/33‑141 (MQ=255)
                                                                                                                    tctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCt                            <  2:625891‑M2/117‑1 (MQ=255)
                                                                                                                        tgaagacgacgccggcctgcctCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGtan                        <  1:592754‑M2/121‑2 (MQ=255)
                                                                                                                        tgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGtag                        >  2:686397‑M2/23‑142 (MQ=255)
                                                                                                                             acgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTACCTCCACTGTGGCATCTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTAt                   <  2:901749‑M2/126‑1 (MQ=255)
                                                                                                                             acgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGtac                   >  2:542577‑M2/18‑142 (MQ=255)
                                                                                                                             acgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGGGAACTGCCCGCTTTTCGTCCGTCTGGtac                   >  1:961511‑M2/18‑142 (MQ=255)
                                                                                                                                              cTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCCTATGCGTGACCGTAGCCAgg  >  2:949066/1‑142 (MQ=255)
                                                                                                                                              |                                                                                                                                              
GGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGT  >  NC_000913/4001507‑4001791

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.