Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 167,743 C→T T87I (ACC→ATC)  fhuA → ferrichrome outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913167,7430CT100.0% 87.6 / NA 28T87I (ACC→ATC) fhuAferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/14);  total (14/14)

AGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCT  >  NC_000913/167636‑167865
                                                                                                           |                                                                                                                          
aGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAg                                                                                             <  2:663182/139‑1 (MQ=255)
      cgcATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTa                                                                                       >  1:457977/1‑139 (MQ=255)
                         aCTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCg                                                                                                >  1:326072/1‑111 (MQ=255)
                         aCTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCg                                                                                                <  2:326072/111‑1 (MQ=255)
                                   cgcgACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgtgt                                                          >  1:304165/1‑139 (MQ=255)
                                   cgcgACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgtgt                                                          <  2:147634/139‑1 (MQ=255)
                                    gcgACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGtc                                                         <  2:169827/139‑1 (MQ=255)
                                     cgACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGtct                                                        >  1:581390/1‑139 (MQ=255)
                                       aCAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCg                                                               >  2:95619/1‑130 (MQ=255)
                                       aCAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCg                                                               <  1:95619/130‑1 (MQ=255)
                                                   cGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTaca                                                                    >  2:700793/1‑113 (MQ=255)
                                                   cGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTaca                                                                    <  1:700793/113‑1 (MQ=255)
                                                     gCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTgg                                        >  1:213745/1‑139 (MQ=255)
                                                              ccGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCaa                               <  2:581390/139‑1 (MQ=255)
                                                                                        caGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTcgcg     <  2:213745/139‑1 (MQ=255)
                                                                                           tCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTc                                                                                           <  1:413821/50‑1 (MQ=255)
                                                                                           tCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTc                                                                                           >  2:413821/1‑50 (MQ=255)
                                                                                           tCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCt  >  1:707993/1‑139 (MQ=255)
                                                                                             tATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGgcg                                      <  1:540930/101‑1 (MQ=255)
                                                                                             tATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGgcg                                      >  2:540930/1‑101 (MQ=255)
                                                                                               tttCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCaa                               <  1:787325/106‑1 (MQ=255)
                                                                                               tttCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCaa                               >  2:787325/1‑106 (MQ=255)
                                                                                                     ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTc                                                                                           >  2:374920/1‑40 (MQ=37)
                                                                                                     ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTc                                                                                           <  1:374920/40‑1 (MQ=37)
                                                                                                     ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgt                                                            <  2:569042/71‑1 (MQ=255)
                                                                                                     ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCgg    >  2:50716/1‑127 (MQ=255)
                                                                                                     ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCgg    <  1:50716/127‑1 (MQ=255)
                                                                                                     ttgtGATCGCCGAAGAGATGCCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgt                                                            >  1:569042/1‑71 (MQ=255)
                                                                                                           |                                                                                                                          
AGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCT  >  NC_000913/167636‑167865

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