Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,191,696 G→A *355* (TAG→TAA yqiI → fimbrial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,191,6960GA100.0% 81.3 / NA 27*355* (TAG→TAAyqiIfimbrial protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/12);  total (15/12)

CTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAGGTTCTATGTTCTTCATCC  >  NC_000913/3191590‑3191829
                                                                                                          |                                                                                                                                     
cTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAg                                                                                                       >  2:786194/1‑139 (MQ=255)
  tCGTTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAACACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGc                                                                                                     >  2:539345/1‑139 (MQ=255)
            cACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATc                                                                                           >  1:201465/1‑139 (MQ=255)
                                aGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAATGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAg                                                                                         <  2:601501/121‑1 (MQ=255)
                                aGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAg                                                                                         >  1:601501/1‑121 (MQ=255)
                                     gAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAgg                                                                                                      >  1:206150/1‑103 (MQ=255)
                                     gAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAgg                                                                                                      <  2:206150/103‑1 (MQ=255)
                                          cGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAgg                                                                                                      <  1:226880/98‑1 (MQ=255)
                                          cGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAgg                                                                                                      >  2:226880/1‑98 (MQ=255)
                                              aGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGcatcat                                                         <  1:786194/139‑1 (MQ=255)
                                                           gAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATAc                                            >  1:349646/1‑139 (MQ=255)
                                                              ggACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAg                                         >  1:693264/1‑139 (MQ=255)
                                                                aCCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGc                                       >  2:1717/1‑139 (MQ=255)
                                                                              aCGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATAc                                            <  1:96079/120‑1 (MQ=255)
                                                                              aCGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATAc                                            >  2:96079/1‑120 (MQ=255)
                                                                                gTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACgg                                                                                                               <  1:568130/51‑1 (MQ=255)
                                                                                gTAATCCTGAAAATCACCTATAACCAATATCTAATACAAACACTAAAACgg                                                                                                               >  2:568130/1‑51 (MQ=255)
                                                                                  aaTCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACa                                                >  2:448822/1‑112 (MQ=255)
                                                                                  aaTCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACa                                                <  1:448822/112‑1 (MQ=255)
                                                                                   aTCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGctct                                                                       >  1:683906/1‑88 (MQ=255)
                                                                                   aTCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGctct                                                                       <  2:683906/88‑1 (MQ=255)
                                                                                       tGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCt                                                       >  1:708307/1‑100 (MQ=255)
                                                                                       tGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCt                                                       <  2:708307/100‑1 (MQ=255)
                                                                                         aaaaTCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAg                                                                                         >  2:543808/1‑64 (MQ=255)
                                                                                         aaaaTCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAg                                                                                         <  1:543808/64‑1 (MQ=255)
                                                                                                  tataACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAGGTTCTATGTTCTTCa     <  2:693264/139‑1 (MQ=255)
                                                                                                     aaCTAATATCTAATACAAACACTAAAACGGGCCATCAGGACCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAGGTTCTATGTTCTTCATcc  <  2:201465/139‑1 (MQ=255)
                                                                                                          |                                                                                                                                     
CTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAGGTTCTATGTTCTTCATCC  >  NC_000913/3191590‑3191829

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