Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 168,350 Δ4 bp coding (867‑870/2244 nt) fhuA → ferrichrome outer membrane transporter

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 168350 168353 4 20 [0] [0] 20 fhuA ferrichrome outer membrane transporter

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 1683490 (0.000)20 (0.640) 19/266 0.8 100% coding (866/2244 nt) fhuA ferrichrome outer membrane transporter
?NC_000913 168354 = 0 (0.000)coding (871/2244 nt) fhuA ferrichrome outer membrane transporter

GCGTTATGCTATTGCACCGGCGTTCACCTGGCGTCCGGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCCGGAAACCGGTTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCCGCTGCCGAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTAATGAGAAGATGGTCGGCTACAGCTTCGATCACGAATTTAACGACACCTTTACTGTGCGTCAGAA  >  NC_000913/168215‑168478
                                                                                                                                      |                                                                                                                                 
gcGTTATGCTATTGCACCGGCGTTCACCTGGCGTCCGGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCCGGAAACCGGTTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCCGCTaacggtaagcgtct                                                                                                                     <  1:87495‑M1/149‑15 (MQ=255)
  gTTATGCTATTGCACCGGCGTTCACCTGGCGTCCGGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCCGGAAACCGGTTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCCGCTaacggtaagcgtctgc                                                                                                                   <  4:101738‑M1/149‑17 (MQ=255)
                              gCGTCCGGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCCGGAAACCGGTTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCCGCTaacggtaa                                                                                                                           >  2:79301‑M1/1‑105 (MQ=255)
                                    ggATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCCGGAAACCGGTTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCCGCTaacggtaagcgtctgccgacagactttaatgaaggggcgaagaacaacac                                                                                 >  1:119340‑M1/1‑99 (MQ=255)
                                          tAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCCGGAAACCGGTTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCCGCTaacggtaagcgtctgccgacagactttaatgaa                                                                                                  <  8:10076‑M1/126‑34 (MQ=255)
                                          tAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCCGGAAACCGGTTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCCGCTaacggtaagcgtctgccgacagactttaatgaa                                                                                                  >  7:10076‑M1/1‑93 (MQ=255)
                                                     ttACCTTCCTTTCTTACTTCCAGAACGAGCCGGAAACCGGTTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCCGCTaacggtaagcgtctgccgacagactttaatgaaggggcgaagaacaacacctattctcgtaatgag                                                                 <  4:52740‑M1/148‑67 (MQ=255)
                                                          ttccttTCTTACTTCCAGAACGAGCCGGAAACCGGTTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCCGCTaacggtaagcgtctgccgacagactttaatgaaggggcgaagaacaacacctattctcgtaatgagaagatg                                                           >  6:1890‑M1/1‑77 (MQ=255)
                                                                     cTTCCAGAACGAGCCGGAAACCGGTTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCCGCTaacggtaagcgtctgccgacagactttaatgaaggggcgaagaacaacacctattctcgtaatgaga                                                                <  2:134472‑M1/133‑68 (MQ=255)
                                                                     cTTCCAGAACGAGCCGGAAACCGGTTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCCGCTaacggtaagcgtctgccgacagactttaatgaaggggcgaagaacaacacctattctcgtaatgaga                                                                >  1:134472‑M1/1‑66 (MQ=255)
                                                                         cAGAACGAGCCGGAAACCGGTTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCCGCTaacggtaagcgtctgccgacagactttaatgaaggggcgaagaacaacacctattctcgtaa                                                                     <  5:586‑M1/124‑63 (MQ=255)
                                                                         cAGAACGAGCCGGAAACCGGTTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCCGCTaacggtaagcgtctgccgacagactttaatgaaggggcgaagaacaacacctattctcgtaa                                                                     >  6:586‑M1/1‑62 (MQ=255)
                                                                          aGAACGAGCCGGAAACCGGTTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCCGCTaacggtaagcgtctgccgacagactttaatgaaggggcgaagaacaacacctattctcgtaatgagaagatggtcggctacagcttc                                            >  6:101021‑M1/1‑61 (MQ=255)
                                                                               gAGCCGGAAACCGGTTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCCGCTaacggtaagcgtctgccgacagactttaatgaaggggcgaagaacaacacctattctcgtaatgagaagatggtcggctacagcttcgatcac                                      <  2:119340‑M1/149‑94 (MQ=255)
                                                                                 gCCGGAAACCGGTTATTACGGCTGGTTGCCGGAAGAGGGAACCGTTGAGCCGCTaacggtaagcgtctgccgacagactttaatgaaggggcgaagaacaacacctattctcgtaatgagaagatggtcggctacagcttcgatcacg                                     >  4:96161‑M1/1‑54 (MQ=255)
                                                                                 gCCGGAAACCGGTTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCCGCTaacggtaagcgtctgccgacagactttaatgaaggggcgaagaacaacacctattctcgtaatgagaagatggtcggctacagcttcgatcacg                                     >  3:26836‑M1/1‑54 (MQ=255)
                                                                                          cGGTTATTACGGCTGGTTGCCGGAAGAGGGAACCGTTGAGCCGCTaacggtaagcgtctgccgacagactttaatgaaggggcgaagaacaacacctattctcgtaatgagaagatggtcggctacagcttcgatcacgaatttaacga                           <  3:96161‑M1/149‑105 (MQ=255)
                                                                                                    ggCTGGTTGCCGAAAGAGGGAACCGTTGAGCCGCTaacggtaagcgtctgccgacagactttaatgaaggggcgaagaacaacacctattctcgtaatgagaagatggtcggctacagcttcgatcacgaatttaacgacacctttact                 <  5:101021‑M1/149‑115 (MQ=255)
                                                                                                                 aaGAGGGAACCGTTGAGCCGCTaacggtaagcgtctgccgacagactttaatgaaggggcgaagaacaacacctattctcgtaatgagaagatggtcggctacagcttcgatcacgaatttaacgacacctttactgtgcgtcagaac     >  2:35393‑M1/1‑22 (MQ=255)
                                                                                                                   gagGGAACCGTTGAGCCGCTaacggtaagcgtctgccgacagactttaatgaaggggcgaagaacaacacctattctcgtaatgagaagatggtcggctacagcttcgatcacgaatttaacgacacctttactgtgcgtcagaacctg  <  5:1890‑M1/149‑130 (MQ=255)
                                                                                                                                      |                                                                                                                                 
GCGTTATGCTATTGCACCGGCGTTCACCTGGCGTCCGGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCCGGAAACCGGTTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCCGCTGCCGAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTAATGAGAAGATGGTCGGCTACAGCTTCGATCACGAATTTAACGACACCTTTACTGTGCGTCAGAA  >  NC_000913/168215‑168478

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: