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 39 [1] [0] 38 corA magnesium/nickel/cobalt transporter

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 40016441 (0.020)38 (0.940) 30/284 0.4 98.7% coding (219/951 nt) corA magnesium/nickel/cobalt transporter
?NC_000913 4001650 = 0 (0.000)coding (225/951 nt) corA magnesium/nickel/cobalt transporter

TATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCC  >  NC_000913/4001505‑4001782
                                                                                                                                           |                                                                                                                                          
tATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctc                                                                                                                                         <  1:46913‑M1/143‑4 (MQ=255)
tATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctc                                                                                                                                         >  2:537968‑M1/1‑140 (MQ=255)
 aTGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGGTTCTTTGAAGACGACGACGGCCTGCatatat                                                                                                                                        >  1:265367/1‑141 (MQ=255)
     aTTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttct                                                                                                                                    >  1:597435‑M1/1‑135 (MQ=255)
      gttatccTGTCGATCCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCTGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttctt                                                                                                                                   <  2:254910‑M1/140‑10 (MQ=255)
            ttGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctt                                                                                                                             <  2:281451‑M1/143‑16 (MQ=255)
            ttGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTCGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctt                                                                                                                             <  2:280055‑M1/143‑16 (MQ=255)
             tGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttcttt                                                                                                                            <  2:618503‑M1/143‑17 (MQ=255)
              gTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttg                                                                                                                           <  2:170553‑M1/143‑18 (MQ=255)
              gTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttg                                                                                                                           <  1:687021‑M1/143‑18 (MQ=255)
                     cGGACGACGATGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgc                                                                                                                    >  1:525956‑M1/1‑119 (MQ=255)
                     cGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGCAGACGACGACGGCCTGCATctccttcttcttctttgcagatgc                                                                                                                    >  2:91587‑M1/1‑119 (MQ=255)
                     cGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATccccttcttcttctttgaagatgc                                                                                                                    >  2:624024‑M1/1‑119 (MQ=255)
                      ggccgacgacgagcgACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCCCCCGGAACCGGAAGACATCGAAGCATCGGCACGTTCCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcg                                                                                                                   <  1:359418‑M1/140‑26 (MQ=255)
                         cgacgacggGCGACTGCGCGTACAATCTGAAGTTGGCCGGAGCCTGGCAGCCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaa                                                                                                                <  2:543339‑M1/143‑29 (MQ=255)
                             gacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGCCGACGGCCTGCATctccttcttcttctttgaagctgcggaagatc                                                                                                            >  1:608514‑M1/1‑111 (MQ=255)
                             gacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGCCGGCCTGCATctccttcttcttctttgaagatgcggaagctc                                                                                                            >  1:50861‑M1/1‑111 (MQ=255)
                             gacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatc                                                                                                            >  2:98952‑M1/1‑111 (MQ=255)
                               cgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagagcac                                                                                                          >  1:249067‑M1/1‑109 (MQ=255)
                                               cAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgcaggtaactccactg                                                                                          >  1:210691‑M1/1‑93 (MQ=255)
                                                 aTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggcaactccaccgtg                                                                                        >  2:70410‑M1/1‑91 (MQ=255)
                                                   cTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggcaactccactgtggc                                                                                      >  2:650227‑M1/1‑89 (MQ=255)
                                                     gAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccctcttcttctttgcagtcgcggaagatcacgccggtaactccactgtggcat                                                                                    <  1:503222‑M1/143‑57 (MQ=255)
                                                        cTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcattta                                                                                 <  1:101908‑M1/143‑60 (MQ=255)
                                                                  gCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtga                                                                       <  1:599273‑M1/143‑70 (MQ=255)
                                                                     tGGCCACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACCGCCTGCATctccttcttcttctgtgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatgg                                                                    <  1:479781‑M1/143‑73 (MQ=255)
                                                                           cccgcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtatctccactgtggcatttaccatccgtgatggtcgtct                                                              <  1:473190‑M1/143‑79 (MQ=255)
                                                                           cccgcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtct                                                              <  1:537968‑M1/143‑79 (MQ=255)
                                                                                ccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaacatccgtgatggtcgtctgttta                                                         >  2:71909‑M1/1‑60 (MQ=255)
                                                                                        ggaggaCATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcgtcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgt                                                 <  2:689045‑M1/143‑92 (MQ=255)
                                                                                        ggAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgt                                                 <  2:472971‑M1/143‑92 (MQ=255)
                                                                                          aaGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtga                                               <  1:596327‑M1/143‑94 (MQ=255)
                                                                                               aTCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtg                                          >  1:117238‑M1/1‑45 (MQ=255)
                                                                                                     gCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgc                                    >  1:226041‑M1/1‑39 (MQ=255)
                                                                                                          ggCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggcaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcctgagcgtgagctgcccgct                               >  2:510507‑M1/1‑34 (MQ=255)
                                                                                                                  tCTTGGGAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaaggtctcgccggtaacacaactgtgcccttttccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtct                       <  2:10034‑M1/143‑118 (MQ=255)
                                                                                                                      agaAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgtcggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtat                   <  1:77251‑M1/142‑122 (MQ=255)
                                                                                                                          gtcgaccacgccCTGCATctccttcttctgcttttgagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgta               <  2:37404‑M1/141‑126 (MQ=255)
                                                                                                                                       tGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtac  >  1:707918‑M1/1‑5 (MQ=255)
                                                                                                                                           |                                                                                                                                          
TATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCC  >  NC_000913/4001505‑4001782

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