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

ATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAG  >  NC_000913/4001510‑4001787
                                                                                                                                            |                                                                                                                                         
tatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatctc                                                                                                                                         >  2:537968‑M2/141‑143 (MQ=255)
tatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatctc                                                                                                                                         <  1:46913‑M2/3‑1 (MQ=255)
     attgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttct                                                                                                                                    >  1:597435‑M2/136‑143 (MQ=255)
      gttatcctgtcgatccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatctgcacgtttctttgaagacgacgacggcctgcatCTCcttctt                                                                                                                                   <  2:254910‑M2/9‑1 (MQ=255)
            ttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttctt                                                                                                                             <  2:281451‑M2/15‑1 (MQ=255)
            ttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctcgcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttctt                                                                                                                             <  2:280055‑M2/15‑1 (MQ=255)
             tgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCttt                                                                                                                            <  2:618503‑M2/16‑1 (MQ=255)
              gtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTg                                                                                                                           <  2:170553‑M2/17‑1 (MQ=255)
              gtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTg                                                                                                                           <  1:687021‑M2/17‑1 (MQ=255)
                     cggacgacgatgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGc                                                                                                                    >  1:525956‑M2/120‑143 (MQ=255)
                     cggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgcagacgacgacggcctgcatCTCCTTCTTCTTCTTTGCAGATGc                                                                                                                    >  2:91587‑M2/120‑143 (MQ=255)
                     cggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCCCCTTCTTCTTCTTTGAAGATGc                                                                                                                    >  2:624024‑M2/120‑143 (MQ=255)
                      ggccgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccccccggaaccggaagacatcgaagcatcggcacgttcctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCg                                                                                                                   <  1:359418‑M2/25‑1 (MQ=255)
                         cgacgacgggcgactgcgcgtacaatctgaagttggccggagcctggcagcccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGaa                                                                                                                <  2:543339‑M2/28‑1 (MQ=255)
                             gacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgccgacggcctgcatCTCCTTCTTCTTCTTTGAAGCTGCGGAAGATc                                                                                                            >  1:608514‑M2/112‑143 (MQ=255)
                             gacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgccggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGCTc                                                                                                            >  1:50861‑M2/112‑143 (MQ=255)
                             gacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATc                                                                                                            >  2:98952‑M2/112‑143 (MQ=255)
                               cgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGAGCAc                                                                                                          >  1:249067‑M2/110‑143 (MQ=255)
                                               caatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCAGGTAACTCCACtg                                                                                          >  1:210691‑M2/94‑143 (MQ=255)
                                                 atctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGCAACTCCACcgtg                                                                                        >  2:70410‑M2/92‑143 (MQ=255)
                                                   ctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGCAACTCCACTGTGGc                                                                                      >  2:650227‑M2/90‑143 (MQ=255)
                                                     gaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCCTCTTCTTCTTTGCAGTCGCGGAAGATCACGCCGGTAACTCCACTGTGGCAt                                                                                    <  1:503222‑M2/56‑1 (MQ=255)
                                                        cttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTa                                                                                 <  1:101908‑M2/59‑1 (MQ=255)
                                                                  gcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGa                                                                       <  1:599273‑M2/69‑1 (MQ=255)
                                                                     tggccacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgaccgcctgcatCTCCTTCTTCTTCTGTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATgg                                                                    <  1:479781‑M2/72‑1 (MQ=255)
                                                                           cccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTATCTCCACTGTGGCATTTACCATCCGTGATGGTCGTCt                                                              <  1:473190‑M2/78‑1 (MQ=255)
                                                                           cccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCt                                                              <  1:537968‑M2/78‑1 (MQ=255)
                                                                                ccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTAACATCCGTGATGGTCGTCTGTTTa                                                         >  2:71909‑M2/61‑143 (MQ=255)
                                                                                        ggaggacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCGTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgt                                                 <  2:689045‑M2/91‑1 (MQ=255)
                                                                                        ggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgt                                                 <  2:472971‑M2/91‑1 (MQ=255)
                                                                                          aagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtga                                               <  1:596327‑M2/93‑1 (MQ=255)
                                                                                               atcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgtg                                          >  1:117238‑M2/46‑143 (MQ=255)
                                                                                                     gcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGc                                    >  1:226041‑M2/40‑143 (MQ=255)
                                                                                                          ggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGCAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCCTGAGCGTGAGCTGCCCGCt                               >  2:510507‑M2/35‑143 (MQ=255)
                                                                                                                  tcttgggagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGGTCTCGCCGGTAACACAACTGTGCCCTTTTCCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCt                       <  2:10034‑M2/117‑1 (MQ=255)
                                                                                                                      agaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGTCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTAt                   <  1:77251‑M2/121‑1 (MQ=255)
                                                                                                                          gtcgaccacgccctgcatCTCCTTCTTCTGCTTTTGAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTa               <  2:37404‑M2/125‑1 (MQ=255)
                                                                                                                                       tgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGtac  >  1:707918‑M2/6‑142 (MQ=255)
                                                                                                                                            |                                                                                                                                         
ATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAG  >  NC_000913/4001510‑4001787

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: