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 30 [0] [0] 30 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)30 (0.870) 28/246 0.4 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

TCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCA  >  NC_000913/4001515‑4001790
                                                                                                                                       |                                                                                                                                            
attgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatc                                                                                                                                              <  7:53806‑M2/1‑1 (MQ=255)
attgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatc                                                                                                                                              >  8:53806‑M2/136‑136 (MQ=255)
         gtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATg                                                                                                                        >  2:132926‑M2/127‑149 (MQ=255)
          tcgaaccggactacgacgagcgactgcgcgtacaatctgaacttggccagagactggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGc                                                                                                                       <  3:14748‑M2/24‑1 (MQ=255)
               ccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAg                                                                                                                  <  1:132926‑M2/29‑1 (MQ=255)
                cggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCgg                                                                                                                     >  2:190117‑M2/120‑145 (MQ=255)
                                actgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTaa                                                                                                     >  8:35461‑M2/104‑145 (MQ=255)
                                actgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTaa                                                                                                     <  7:35461‑M2/42‑1 (MQ=255)
                                    cgcgtacaatctgatcttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgactacggcctgcatCTCCTTCTTCTTCttt                                                                                                                               <  5:142957‑M2/16‑1 (MQ=255)
                                    agagtacaatctgaacttggacagagccttgcaacccgaccggaactgtaatacatcgaagcatcggcacgtttctttgaagacgaagacggccggcatCTCCTTCTTCTTCttt                                                                                                                               >  6:142957‑M2/100‑115 (MQ=255)
                                          caatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTgg                                                                                          >  3:4525‑M2/94‑146 (MQ=255)
                                          caatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTgg                                                                                          <  4:4525‑M2/53‑1 (MQ=255)
                                             tctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttct                                                                                                                                 >  6:49525‑M2/91‑104 (MQ=255)
                                             tctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttct                                                                                                                                 <  5:49525‑M2/14‑1 (MQ=255)
                                             tctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTa                                                                                    >  2:181546‑M2/91‑149 (MQ=255)
                                                         cagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGc                                                                                                           >  2:91407‑M2/79‑114 (MQ=255)
                                                         cagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGc                                                                                                           <  1:91407‑M2/36‑1 (MQ=255)
                                                          agagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATgg                                                                       <  1:190117‑M2/72‑1 (MQ=255)
                                                                   caacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGtt                                                              <  1:181546‑M2/81‑1 (MQ=255)
                                                                      cccgcccggaactggaagacatcgaagcatcggcacgtttctgtgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGCTCACGCCg                                                                                                         <  7:120054‑M2/38‑1 (MQ=255)
                                                                      cccgcccggaactggaagacatcgaagcatcggcaagtttctttgaagacgacgacggcctgcctCTCCTTCTTCTTCTTTGAAGATGAGTAAGATAACGCCg                                                                                                         >  8:120054‑M2/66‑103 (MQ=255)
                                                                       ccgcccggaactggaagacatcgaagcctcggcacgtttctttgaagacgacgccggcctgcatCTCCTACTTCTTCTTGGAAGATGCGGAAGATCACGCCGGTCACTCCACGGTGGCATTTACCATCCGTGATGGTCGTCTGTTTAct                                                          >  7:125078‑M2/65‑149 (MQ=255)
                                                                            cggcactggaagacatcgaagcatcggcaagtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcg                                                     <  8:125078‑M2/90‑1 (MQ=255)
                                                                                          atcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGc                                       >  6:48935‑M2/46‑149 (MQ=255)
                                                                                                   tcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATg                                                                        <  4:124057‑M2/71‑1 (MQ=255)
                                                                                                   tcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATg                                                                        >  3:124057‑M2/37‑107 (MQ=255)
                                                                                                             tctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCg                    >  4:97481‑M2/27‑149 (MQ=255)
                                                                                                                      acgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGc           >  5:88399‑M2/18‑149 (MQ=255)
                                                                                                                        gacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGccc         >  3:109306‑M2/16‑149 (MQ=255)
                                                                                                                               gcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCa  <  4:109306‑M2/141‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                            
TCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCA  >  NC_000913/4001515‑4001790

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