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 35 [0] [0] 35 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)35 (0.960) 33/258 0.2 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

GATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAG  >  NC_000913/4001513‑4001791
                                                                                                                                   |                                                                                                                                                   
gATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttcttt                                                                                                                                     <  2:16926‑M1/148‑17 (MQ=255)
 aTCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttg                                                                                                                                    <  5:9793‑M1/148‑18 (MQ=255)
  tCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcctt                                                                                                                                               >  6:64541‑M1/1‑130 (MQ=255)
  tCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcctt                                                                                                                                               <  5:64541‑M1/136‑7 (MQ=255)
         gAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcgg                                                                                                                           <  3:126935‑M1/149‑27 (MQ=255)
             cGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttc                                                                                                                                              >  8:49322‑M1/1‑119 (MQ=255)
             cGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttc                                                                                                                                              <  7:49322‑M1/126‑8 (MQ=255)
                      acgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccgg                                                                                                              <  6:64351‑M1/149‑40 (MQ=255)
                             aCTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactc                                                                                                        >  4:46322‑M1/1‑103 (MQ=255)
                                   cggacaaTCTGACCTTGGCCAGAGCCTGGCAACCCGCCCGGAACAGGAAGACAGCGAAGCATCGGCACGGTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtg                                                                                                 <  7:139509‑M1/146‑53 (MQ=255)
                                          tCTTAACTTGGACAGATCCTGGAAACCCGACCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGAAGACGGCCTGCATctccttcttcgtctttgaagatgcggaaga                                                                                                                       <  8:34613‑M1/120‑31 (MQ=255)
                                          tCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGCCATCGAAGCATCGGCACGTTTCTGTGAAGCCGACGACGGCCGGCAGctccttcgtcttctgtgcagatgcggaaga                                                                                                                       >  7:34613‑M1/1‑90 (MQ=255)
                                               aCTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgc                                                                                                                 <  5:168041‑M1/121‑37 (MQ=255)
                                               aCTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgc                                                                                                                 >  6:168041‑M1/1‑85 (MQ=255)
                                                cTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatcc                                                                                    >  5:132676‑M1/1‑84 (MQ=255)
                                                     ccAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGCCGGCCTGCATcgccttcttcttcgttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgat                                                                               >  7:45094‑M1/1‑79 (MQ=255)
                                                                            aaCTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtga                                                        <  3:46322‑M1/149‑94 (MQ=255)
                                                                                              cATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccgtcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccctccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttt                                      <  8:45094‑M1/149‑112 (MQ=255)
                                                                                                       gTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgta                             >  5:83505‑M1/1‑29 (MQ=255)
                                                                                                          tCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatc                           <  6:132676‑M1/148‑123 (MQ=255)
                                                                                                            tttGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgta                        >  5:12303‑M1/1‑24 (MQ=255)
                                                                                                             ttGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaa                                                 >  4:146137‑M1/1‑23 (MQ=255)
                                                                                                             ttGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaa                                                 <  3:146137‑M1/123‑101 (MQ=255)
                                                                                                              tGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgta                        <  6:12303‑M1/147‑126 (MQ=255)
                                                                                                               gAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgc                     >  6:164566‑M1/1‑21 (MQ=255)
                                                                                                                aaGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgc                     <  6:83505‑M1/148‑129 (MQ=255)
                                                                                                                aaGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgctt                                       >  2:74890‑M1/1‑20 (MQ=255)
                                                                                                                aaGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgctt                                       <  1:74890‑M1/130‑111 (MQ=255)
                                                                                                                           ggCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagc          >  1:16969‑M1/1‑9 (MQ=255)
                                                                                                                           ggCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagc          >  1:116862‑M1/1‑9 (MQ=255)
                                                                                                                              cTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagt      >  6:84767‑M1/1‑6 (MQ=255)
                                                                                                                                gCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagtc     >  4:138142‑M1/1‑4 (MQ=255)
                                                                                                                                  aTctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagtcgat  <  5:164566‑M1/149‑148 (MQ=255)
                                                                                                                                   tctccttcttcttctttgaagatgcggaagatcacgccggtaactccgctgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaa                                                 >  7:99891‑M1/1‑1 (MQ=255)
                                                                                                                                   tctccttcttcttctttgaagatgcggaagatcacgccggtaactccgctgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaa                                                 <  8:99891‑M1/101‑101 (MQ=255)
                                                                                                                                   |                                                                                                                                                   
GATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAG  >  NC_000913/4001513‑4001791

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: