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 = 40016444 (0.140)23 (0.800) 23/284 0.3 85.2% coding (219/951 nt) corA magnesium/nickel/cobalt transporter
?NC_000913 4001650 = 4 (0.140)coding (225/951 nt) corA magnesium/nickel/cobalt transporter

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

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.