Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 612,345 A→G V50V (GTT→GTC fepA ← ferric enterobactin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913612,3450AG100.0% 103.9 / NA 30V50V (GTT→GTCfepAferric enterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/13);  total (17/13)

CATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGG  >  NC_000913/612273‑612422
                                                                        |                                                                             
cATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAGACGCCAgg                                                                       >  1:120525/1‑81 (MQ=255)
   ggTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAGACGCCAGgcgc                                                                    >  1:120540/1‑81 (MQ=255)
       cGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAGACGCCAGGCGCCTg                                                                 >  1:120535/1‑80 (MQ=255)
                  cGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAGACGCCAGGCGCCTGTAAGTTc                                                          <  1:120498/76‑1 (MQ=255)
                  cGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAGACGCCAGGCGCCTGTAAGTTCt                                                         <  2:120499/77‑1 (MQ=255)
                  cGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAGACGCCAGGCGCCTGTAAGTTCTGCTc                                                     >  2:120532/1‑81 (MQ=255)
                   gACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAGACGCCAGGCGCCTGTAAGTTCTGCTcg                                                    >  2:120546/1‑81 (MQ=255)
                   gACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAGACGCCAGGCGCCTGTAAGTTCTGCTc                                                     <  1:120500/80‑1 (MQ=255)
                    acacATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAGACGCCAGGCGCCTGTAAGTTCTGCTcgg                                                   <  2:120501/81‑1 (MQ=255)
                    acacATCGCGGGCAACCGGGCTTTTGCGGATTTCATCTGCGGTGATGGTCGAGACGCCAGGCGCCTGTAAGTTCTGCTcgg                                                   >  1:120539/1‑81 (MQ=255)
                            cgGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAGACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGt                                            >  2:120541/1‑80 (MQ=255)
                              ggCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAGACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAAc                                         <  1:120502/81‑1 (MQ=255)
                                 aaCCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAGACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGa                                       >  1:120552/1‑80 (MQ=255)
                                 aaCCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAGACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGAc                                      <  1:120503/81‑1 (MQ=255)
                                   ccGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAGACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACaa                                    >  1:120533/1‑81 (MQ=255)
                                      ggTTTTTGCGGATTTCATCTGCGGTGATTGTCGAGACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAg                                 >  2:120542/1‑81 (MQ=255)
                                          tttGCGGATTTCATCTGCGGTGATGGTCGAGACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTAt                              <  1:120504/80‑1 (MQ=255)
                                          tttGCGGATTTCATCTGCGGTGATGGTCGAGACGCCAGGCGCCTGTAAGCTCTGCTCGGCGGCGGTAACGACAATAGTATc                             >  2:120547/1‑81 (MQ=255)
                                           ttGCGGATTTCATCTGCGGTGATGGTCGAGACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCg                            >  1:120549/1‑81 (MQ=255)
                                           ttGCGGATTTCATCTGCGGTGATGGTCGAGACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCg                            >  1:120548/1‑81 (MQ=255)
                                                 aTTTCATCTGCGGTGATGGTCGAGACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCAtgt                      <  2:120505/81‑1 (MQ=255)
                                                  tttCATCTGCGGTGATGGTCGAGACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCAtgt                      >  2:120538/1‑80 (MQ=255)
                                                   ttCATCTGCGGTGATGGTCGAGACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCAtgtg                     >  2:120545/1‑80 (MQ=255)
                                                      aTCTGCGGTGATGGTCGAGACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGaaa                  <  2:120506/80‑1 (MQ=255)
                                                      aTCTGCGGTGATGGTCGAGACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGaaa                  <  1:120507/80‑1 (MQ=255)
                                                      aTCTGCGGTGATGGTCGAGACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGaaa                  >  1:120553/1‑80 (MQ=255)
                                                            ggTGATGGTCGAGACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGt           <  1:120508/81‑1 (MQ=255)
                                                                   gTCGAGACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTAtcggtc    <  1:120509/81‑1 (MQ=255)
                                                                    tCGAGACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTAtcggtcg   <  2:120510/81‑1 (MQ=255)
                                                                     cGAGACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTAtcggtcgg  >  2:120543/1‑81 (MQ=255)
                                                                        |                                                                             
CATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGG  >  NC_000913/612273‑612422

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