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

TATGCTATTGCACCGGCGTTCACCTGGCGTCCGGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCCGGAAACCGGTTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCCGCTGCCGAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTAATGAGAAGATGGTCGGCTACAGCTTCGATCACGAATTTAACGACACCTTTACTGTGCGTCAGAACCTG  >  NC_000913/168219‑168482
                                                                                                                                       |                                                                                                                                
gcgttatgctattgcaccggcgttcacctggcgtccggatgataaaaccaattttaccttcctttcttacttccagaacgagccggaaaccggttattacggctggttgccgaaagagggaaccgttgagccgctAACGGTAAGCGTCt                                                                                                                     <  1:87495‑M2/14‑1 (MQ=255)
  gttatgctattgcaccggcgttcacctggcgtccggatgataaaaccaattttaccttcctttcttacttccagaacgagccggaaaccggttattacggctggttgccgaaagagggaaccgttgagccgctAACGGTAAGCGTCTGc                                                                                                                   <  4:101738‑M2/16‑1 (MQ=255)
                              gcgtccggatgataaaaccaattttaccttcctttcttacttccagaacgagccggaaaccggttattacggctggttgccgaaagagggaaccgttgagccgctAACGGTaa                                                                                                                           >  2:79301‑M2/106‑113 (MQ=255)
                                    ggatgataaaaccaattttaccttcctttcttacttccagaacgagccggaaaccggttattacggctggttgccgaaagagggaaccgttgagccgctAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAacac                                                                                 >  1:119340‑M2/100‑149 (MQ=255)
                                          taaaaccaattttaccttcctttcttacttccagaacgagccggaaaccggttattacggctggttgccgaaagagggaaccgttgagccgctAACGGTAAGCGTCTGCCGACAGACTTTAATGaa                                                                                                  <  8:10076‑M2/33‑1 (MQ=255)
                                          taaaaccaattttaccttcctttcttacttccagaacgagccggaaaccggttattacggctggttgccgaaagagggaaccgttgagccgctAACGGTAAGCGTCTGCCGACAGACTTTAATGaa                                                                                                  >  7:10076‑M2/94‑126 (MQ=255)
                                                     ttaccttcctttcttacttccagaacgagccggaaaccggttattacggctggttgccgaaagagggaaccgttgagccgctAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTAATgag                                                                 <  4:52740‑M2/66‑1 (MQ=255)
                                                          ttcctttcttacttccagaacgagccggaaaccggttattacggctggttgccgaaagagggaaccgttgagccgctAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTAATGAGAAGATg                                                           >  6:1890‑M2/78‑149 (MQ=255)
                                                                     cttccagaacgagccggaaaccggttattacggctggttgccgaaagagggaaccgttgagccgctAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTAATgaga                                                                <  2:134472‑M2/67‑1 (MQ=255)
                                                                     cttccagaacgagccggaaaccggttattacggctggttgccgaaagagggaaccgttgagccgctAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTAATgaga                                                                >  1:134472‑M2/67‑133 (MQ=255)
                                                                         cagaacgagccggaaaccggttattacggctggttgccgaaagagggaaccgttgagccgctAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTaa                                                                     <  5:586‑M2/62‑1 (MQ=255)
                                                                         cagaacgagccggaaaccggttattacggctggttgccgaaagagggaaccgttgagccgctAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTaa                                                                     >  6:586‑M2/63‑124 (MQ=255)
                                                                          agaacgagccggaaaccggttattacggctggttgccgaaagagggaaccgttgagccgctAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTAATGAGAAGATGGTCGGCTACAGCTTc                                            >  6:101021‑M2/62‑148 (MQ=255)
                                                                               gagccggaaaccggttattacggctggttgccgaaagagggaaccgttgagccgctAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTAATGAGAAGATGGTCGGCTACAGCTTCGATCAc                                      <  2:119340‑M2/93‑1 (MQ=255)
                                                                                 gccggaaaccggttattacggctggttgccggaagagggaaccgttgagccgctAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTAATGAGAAGATGGTCGGCTACAGCTTCGATCACg                                     >  4:96161‑M2/55‑148 (MQ=255)
                                                                                 gccggaaaccggttattacggctggttgccgaaagagggaaccgttgagccgctAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTAATGAGAAGATGGTCGGCTACAGCTTCGATCACg                                     >  3:26836‑M2/55‑148 (MQ=255)
                                                                                          cggttattacggctggttgccggaagagggaaccgttgagccgctAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTAATGAGAAGATGGTCGGCTACAGCTTCGATCACGAATTTAACGa                           <  3:96161‑M2/104‑1 (MQ=255)
                                                                                                    ggctggttgccgaaagagggaaccgttgagccgctAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTAATGAGAAGATGGTCGGCTACAGCTTCGATCACGAATTTAACGACACCTTTACt                 <  5:101021‑M2/114‑1 (MQ=255)
                                                                                                                 aagagggaaccgttgagccgctAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTAATGAGAAGATGGTCGGCTACAGCTTCGATCACGAATTTAACGACACCTTTACTGTGCGTCAGAAc     >  2:35393‑M2/23‑148 (MQ=255)
                                                                                                                   gagggaaccgttgagccgctAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTAATGAGAAGATGGTCGGCTACAGCTTCGATCACGAATTTAACGACACCTTTACTGTGCGTCAGAACCTg  <  5:1890‑M2/129‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                
TATGCTATTGCACCGGCGTTCACCTGGCGTCCGGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCCGGAAACCGGTTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCCGCTGCCGAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTAATGAGAAGATGGTCGGCTACAGCTTCGATCACGAATTTAACGACACCTTTACTGTGCGTCAGAACCTG  >  NC_000913/168219‑168482

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: