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 = 40016434 (0.130)23 (0.770) 23/284 0.3 85.2% coding (219/951 nt) corA magnesium/nickel/cobalt transporter
?NC_000913 4001649 = 4 (0.130)coding (225/951 nt) corA magnesium/nickel/cobalt transporter

ATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTAT  >  NC_000913/4001513‑4001774
                                                                                                                                  |                                                                                                                                   
acgttggcGAACGGGACGACGACGAGCGACTGCGCGTACAATCGGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTCTGAAGACGACGACGGCCTGCAGctccttcttctt                                                                                                                         <  2:347332‑M1/140‑13 (MQ=255)
aTCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttctt                                                                                                                         >  1:132372‑M1/1‑131 (MQ=255)
 tCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAACACATCGAAGCATCGGCACGTTTCTTTGCAGACGACGACGCCCTGCATctccttcttcttc                                                                                                                        >  1:65342‑M1/1‑130 (MQ=255)
            cGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgc                                                                                                             <  2:19627‑M1/143‑25 (MQ=255)
                agacgacgagCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTTGCAACCCGCCCGCAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaa                                                                                                         <  2:317139‑M1/142‑29 (MQ=255)
                   caaggagcgaCTGCGAGGACAGTCTGAACTAGACCAGAGCCGGGCAACCCGCCCGGAACTGGAAGACATCGAAGCACCGGCACGTTTCTTTGCAGACGACGACGGCCTGCATctccttcttcttcttgcaagatgcggaagat                                                                                                      <  1:281694‑M1/139‑32 (MQ=255)
                                ggcgtACACTCTGACCTTGGCCAGAGCCTGGCAACCCTCCTGGATCCGGAAGATATCGAAGCATCTGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaact                                                                                         <  1:281253‑M1/142‑45 (MQ=255)
                                         tCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGCACTAGCAGGCCAAGACGCACCGGCACGTTTCATAAATTTATCTGTTTACCTGCAT‑‑‑‑‑CTCCTTCTTCTCCTTGGACGGTGCGGAAGATCACGCCGGTAACTCCACTGTgg                                                                           >  2:525118/1‑143 (MQ=255)
                                           tGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactcccccggggca                                                                              >  2:43385‑M1/1‑88 (MQ=255)
                                           tGAACTGGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggca                                                                              <  2:251604‑M1/143‑56 (MQ=255)
                                            gAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACTAC‑‑‑‑GCAAGAATATGCAAGTATCGTCTTCTTTGAAGTTGCGGAAGGTCACGCCGGTACCTCGACTGTGGcgt                                                                        >  2:382187/1‑141 (MQ=255)
                                            gAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaacgcaactggggcat                                                                             >  2:38816‑M1/1‑87 (MQ=255)
                                                     cAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcccgccggtaactccgctgtggcatttaccatcc                                                                    >  2:114059‑M1/1‑78 (MQ=255)
                                                        agCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGAAGACGGCCTGCATctccttcttctttttttaagatgcggaagatcacgccggtaactccactgtgatatctatcatccgtg                                                                 >  1:55926‑M1/1‑75 (MQ=255)
                                                          ccTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGAAGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagataacgccggtaactccgccggggcatttgacaagcgtgat                                                               >  2:276937‑M1/1‑73 (MQ=255)
                                                                  cccgcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtct                                                       <  2:132372‑M1/143‑79 (MQ=255)
                                                                       ccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgttta                                                  <  2:410931‑M1/143‑84 (MQ=255)
                                                                                gAAGACATCGAAGCATCGGCACGTTGCTTTGAAGACGACGACGGCCCGC‑‑ATT‑‑‑TCCTTCTTCTTCTGTGAAGCTGCGGAAGATCACGCCGGTAACACCACTGGGGCATTTGGGATCCGCGATGGTCGGCTGTGTAGTCGgagtg                                    >  2:178003/1‑143 (MQ=255)
                                                                                          aaGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaact                               >  2:332837‑M1/1‑41 (MQ=255)
                                                                                            gCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATAACTCTACCTGGTGGTGCATAT‑‑GGTTGTTGGCACCCCAAACCCCACTGTGGCATTTACCATACGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGc                        >  2:525096/1‑143 (MQ=255)
                                                                                             cATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcc                            >  2:211242‑M1/1‑38 (MQ=255)
                                                                                                      gTTTCTTTGAAGACGACGACGGCCTGCATcgccttcttcttcttagaagaagcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcg                   >  2:184655‑M1/1‑29 (MQ=255)
                                                                                                        ttCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtc                 <  2:525933‑M1/143‑117 (MQ=255)
                                                                                                           cttgAAGACGACGACGGCCTGCATctccctcttctgctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgt              <  2:482952‑M1/142‑120 (MQ=255)
                                                                                                            ttGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgttagcgtgaactgcccgctttccgtctgta             >  2:288509‑M1/1‑23 (MQ=255)
                                                                                                             tGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgacctgcccgcgtttcgtctgtat            >  2:236662‑M1/1‑22 (MQ=255)
                                                                                                                       gacgGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtg  >  1:312423‑M1/1‑12 (MQ=255)
                                                                                                                                  |                                                                                                                                   
ATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTAT  >  NC_000913/4001513‑4001774

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.