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

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4001645 4001649 5 32 [0] [0] 32 corA magnesium/nickel/cobalt transporter

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 40016440 (0.000)32 (0.770) 31/258 0.5 100% coding (219/951 nt) corA magnesium/nickel/cobalt transporter
?NC_000913 4001650 = 0 (0.000)coding (225/951 nt) corA magnesium/nickel/cobalt transporter

AGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAG  >  NC_000913/4001503‑4001791
                                                                                                                                             |                                                                                                                                                   
aGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttc                                                                                                                                              <  7:128053‑M1/149‑8 (MQ=255)
    tGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttct                                                                                                                                          >  5:39770‑M1/1‑138 (MQ=255)
      gATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttc                                                                                                                                        >  7:110220‑M1/1‑136 (MQ=255)
          gATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttcttt                                                                                                                                     <  6:39770‑M1/148‑17 (MQ=255)
                   gAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcgg                                                                                                                           <  8:110220‑M1/149‑27 (MQ=255)
                         gacgacgacGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctt                                                                                                                                      >  3:163402‑M1/1‑117 (MQ=255)
                         gacgacgacGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctt                                                                                                                                      <  4:163402‑M1/132‑16 (MQ=255)
                               gacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaa                                                                                                                                  <  6:191607‑M1/130‑20 (MQ=255)
                               gacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaa                                                                                                                                  >  5:191607‑M1/1‑111 (MQ=255)
                                                  aaTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaa                                                                                                                         <  6:133703‑M1/120‑29 (MQ=255)
                                                  aaTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaa                                                                                                                         >  5:133703‑M1/1‑92 (MQ=255)
                                                           ttGTCCAGAGCCTGGCAACCCGCCCGGAACTGTAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACTACGGCCTGCATctccttcttcttctttgaagatgcagaagatcacgccggtaactc                                                                                                        <  4:48331‑M1/128‑46 (MQ=255)
                                                           ttGGCCAGAGCCTGGCAACACGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcagaagatcacgccggtaacgc                                                                                                        >  3:48331‑M1/1‑83 (MQ=255)
                                                             ggCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgt                                                                                  >  6:174762‑M1/1‑81 (MQ=255)
                                                                cAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtacctccactgtggcatttaccatccgtgatg                                                                              >  5:112044‑M1/1‑78 (MQ=255)
                                                                          cAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtt                                                                    <  6:112044‑M1/149‑82 (MQ=255)
                                                                             cccgcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttga                                                                                                                                   <  7:151875‑M1/83‑19 (MQ=255)
                                                                             cccgcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttga                                                                                                                                   >  8:151875‑M1/1‑65 (MQ=255)
                                                                                           gAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgt                                                    >  8:139300‑M1/1‑51 (MQ=255)
                                                                                            aaGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtt                                                                    >  7:165315‑M1/1‑50 (MQ=255)
                                                                                            aaGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtt                                                                    <  8:165315‑M1/131‑82 (MQ=255)
                                                                                             aGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccg                                                                                                               >  7:56226‑M1/1‑49 (MQ=255)
                                                                                             aGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccg                                                                                                               <  8:56226‑M1/87‑39 (MQ=255)
                                                                                                 aTCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgc                                             >  1:78962‑M1/1‑45 (MQ=255)
                                                                                                  tCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcc                                            >  2:181237‑M1/1‑44 (MQ=255)
                                                                                                                         gAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatg                      <  7:139300‑M1/148‑128 (MQ=255)
                                                                                                                          aaGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgta                             <  7:104912‑M1/140‑121 (MQ=255)
                                                                                                                          aaGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgta                             >  8:104912‑M1/1‑20 (MQ=255)
                                                                                                                               gacgacgGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtga                                                  >  4:178004‑M1/1‑15 (MQ=255)
                                                                                                                               gacgacgGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtga                                                  <  3:178004‑M1/114‑100 (MQ=255)
                                                                                                                                  gacgGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccg              >  6:128490‑M1/1‑12 (MQ=255)
                                                                                                                                            aTctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagtcgat  <  1:181237‑M1/149‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
AGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAG  >  NC_000913/4001503‑4001791

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