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 45 [0] [0] 45 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)45 (1.140) 40/262 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

TATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGA  >  NC_000913/4001505‑4001795
                                                                                                                                                 |                                                                                                                                                 
tgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTcc                                                                                                                                                <  5:159022‑M2/4‑1 (MQ=255)
      atggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttc                                                                                                                                          >  2:145988‑M2/140‑149 (MQ=255)
             gatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCttt                                                                                                                                    <  1:145988‑M2/16‑1 (MQ=255)
               tcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTg                                                                                                                                   >  2:61875‑M2/131‑147 (MQ=255)
                    tcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGc                                                                                                                            >  7:89604‑M2/126‑149 (MQ=255)
                    tcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGc                                                                                                                            <  5:84345‑M2/24‑1 (MQ=255)
                           ggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGAt                                                                                                                     <  7:63698‑M2/31‑1 (MQ=255)
                               gacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatc                                                                                                                                                   >  8:112206‑M2/115‑115 (MQ=255)
                               gacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatc                                                                                                                                                   <  7:112206‑M2/1‑1 (MQ=255)
                               gacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttct                                                                                                                                            <  7:151811‑M2/8‑1 (MQ=255)
                               gacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttct                                                                                                                                            >  8:151811‑M2/115‑122 (MQ=255)
                                acgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGc                                                                                                                            >  5:74206‑M2/114‑137 (MQ=255)
                                acgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGc                                                                                                                            <  6:74206‑M2/24‑1 (MQ=255)
                                     gagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGc                                                                                                                            <  4:42641‑M2/24‑1 (MQ=255)
                                     gagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGc                                                                                                                            >  3:42641‑M2/109‑132 (MQ=255)
                                               gcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgt                                                                                                 <  7:149692‑M2/51‑1 (MQ=255)
                                                 gtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttc                                                                                                                                       >  4:19269‑M2/97‑109 (MQ=255)
                                                 gtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttc                                                                                                                                       <  3:19269‑M2/13‑1 (MQ=255)
                                                     aatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGa                                                                                                                      >  3:146159‑M2/93‑122 (MQ=255)
                                                     aatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGa                                                                                                                      <  4:146159‑M2/30‑1 (MQ=255)
                                                       tctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTa                                                                                         >  1:50836‑M2/91‑149 (MQ=255)
                                                              ttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCg                                                                                  >  1:111866‑M2/84‑149 (MQ=255)
                                                                           ggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCttt                                                                                                                                    >  2:112364‑M2/71‑86 (MQ=255)
                                                                           ggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCttt                                                                                                                                    <  1:112364‑M2/16‑1 (MQ=255)
                                                                             caacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGtt                                                                   >  3:95239‑M2/69‑149 (MQ=255)
                                                                                 ccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTAct                                                               <  2:50836‑M2/85‑1 (MQ=255)
                                                                                   gcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTActct                                                             >  2:11391‑M2/63‑149 (MQ=255)
                                                                                     ccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgc                                                           <  1:11391‑M2/89‑1 (MQ=255)
                                                                                              gaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttctt                                                                                                                                        >  1:76005‑M2/52‑63 (MQ=255)
                                                                                              gaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgt                                                   <  8:89604‑M2/97‑1 (MQ=255)
                                                                                              gaagacatcgaagcatcggcacgtttatttgaagacgacgacggcctgcatCTCcttcttctt                                                                                                                                        <  2:76005‑M2/12‑1 (MQ=255)
                                                                                                    atcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgt                                                   <  2:2683‑M2/97‑1 (MQ=255)
                                                                                                    atcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgt                                                   >  1:2683‑M2/46‑142 (MQ=255)
                                                                                                      cgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGaa                                                <  4:12500‑M2/100‑1 (MQ=255)
                                                                                                      cgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGaa                                                >  3:12500‑M2/44‑143 (MQ=255)
                                                                                                        aagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGc                                        >  2:174441‑M2/42‑149 (MQ=255)
                                                                                                                     tttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTAt                           >  8:13640‑M2/29‑149 (MQ=255)
                                                                                                                        ctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATc                          <  4:81222‑M2/122‑1 (MQ=255)
                                                                                                                        ctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATc                          >  3:81222‑M2/26‑147 (MQ=255)
                                                                                                                           tgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTa                       <  7:13640‑M2/125‑1 (MQ=255)
                                                                                                                           tgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATg                     >  3:23886‑M2/23‑149 (MQ=255)
                                                                                                                                  gacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGccc              >  7:93790‑M2/16‑149 (MQ=255)
                                                                                                                                         gcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCa       <  8:93790‑M2/141‑1 (MQ=255)
                                                                                                                                          cctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAg      <  4:95239‑M2/142‑1 (MQ=255)
                                                                                                                                              catCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGGCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGa  >  4:145064‑M2/4‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGA  >  NC_000913/4001505‑4001795

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: