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 (1.160) 33/222 0.1 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

AGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGA  >  NC_000913/4001538‑4001795
                                                                                                                |                                                                                                                                                 
cgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttc                                                                                                                                          <  2:50941‑M2/10‑1 (MQ=255)
cgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttc                                                                                                                                          >  1:50941‑M2/113‑122 (MQ=255)
             cgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtg                                                                                                  <  4:141241‑M2/50‑1 (MQ=255)
                 tacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTgg                                                                                               >  2:89426‑M2/96‑148 (MQ=255)
                      tctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttc                                                                                                                                             >  1:114440‑M2/91‑97 (MQ=255)
                      tctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttc                                                                                                                                             <  2:114440‑M2/7‑1 (MQ=255)
                               ggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACt                                                                                                   <  7:38324‑M2/49‑1 (MQ=255)
                               ggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACt                                                                                                   >  8:38324‑M2/82‑130 (MQ=255)
                                           gcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCt                                                                      <  1:89426‑M2/78‑1 (MQ=255)
                                               cccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCACGGatggtaaa                                                                        <  4:112311‑M2/76‑4 (MQ=255)
                                               cccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCACGGatggtaaa                                                                        >  3:112311‑M2/66‑138 (MQ=255)
                                                ccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGc                                                                                                                            <  1:6261‑M2/24‑1 (MQ=255)
                                                ccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGc                                                                                                                            >  2:6261‑M2/65‑88 (MQ=255)
                                                    ccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGa                                                                               <  2:42438‑M2/69‑1 (MQ=255)
                                                    ccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGa                                                                               >  1:42438‑M2/61‑129 (MQ=255)
                                                         actggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgc                                                           >  2:183546‑M2/56‑144 (MQ=255)
                                                         actggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgc                                                           <  1:183546‑M2/89‑1 (MQ=255)
                                                             gaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGc                                                                                              <  6:73282‑M2/54‑1 (MQ=255)
                                                             gaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGc                                                                                              >  5:73282‑M2/52‑105 (MQ=255)
                                                               agacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCttt                                                                                                                                    >  3:29002‑M2/50‑65 (MQ=255)
                                                               agacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCttt                                                                                                                                    <  4:29002‑M2/16‑1 (MQ=255)
                                                                   atcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGAt                                                                                                                     <  8:41642‑M2/31‑1 (MQ=255)
                                                                   atcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgtg                                                                                                <  7:49216‑M2/52‑1 (MQ=255)
                                                                   atcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgtg                                                                                                >  8:49216‑M2/46‑97 (MQ=255)
                                                                                    tttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGa                                                                               >  4:123145‑M2/29‑97 (MQ=255)
                                                                                       ctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGa                                                                                                                                  <  7:144343‑M2/18‑1 (MQ=255)
                                                                                       ctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGa                                                                                                                                  >  8:144343‑M2/26‑43 (MQ=255)
                                                                                        tttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCgg                                                                                                             >  2:174002‑M2/25‑63 (MQ=255)
                                                                                        tttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCgg                                                                                                             <  1:174002‑M2/39‑1 (MQ=255)
                                                                                        tttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAGCTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTa                       >  1:143070‑M2/25‑149 (MQ=255)
                                                                                             agacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGgt                                                                           >  4:12982‑M2/20‑92 (MQ=255)
                                                                                             agacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGgt                                                                           <  3:12982‑M2/73‑1 (MQ=255)
                                                                                               acgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGc                >  2:131172‑M2/18‑149 (MQ=255)
                                                                                                        gcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAGCTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCa       <  2:143070‑M2/141‑1 (MQ=255)
                                                                                                             catCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGa  <  1:131172‑M2/146‑1 (MQ=255)
                                                                                                                |                                                                                                                                                 
AGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGA  >  NC_000913/4001538‑4001795

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