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% 94.2 / NA 30*355* (TAG→TAAyqiIfimbrial protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/16);  total (14/16)

CTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAGGTTCTATGTTCTTCAT  >  NC_000913/3191569‑3191827
                                                                                                                               |                                                                                                                                   
cttTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACa                                                                                                                          <  1:1410703/139‑1 (MQ=255)
             tAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACgg                                                                                                             <  2:816323/139‑1 (MQ=255)
                  cgcCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCAt                                                                                                        <  1:684302/139‑1 (MQ=255)
                                             cAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACa                                                                                                                          <  2:1323935/94‑1 (MQ=255)
                                             cAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACa                                                                                                                          >  1:1323935/1‑94 (MQ=255)
                                                aaCTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTaaaa                                                                                                                <  2:693507/101‑1 (MQ=255)
                                                aaCTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTaaaa                                                                                                                >  1:693507/1‑101 (MQ=255)
                                                aaCTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGAATCCTGAAAATCACCTATAACTAATATCTAATACa                                                                                                                          >  1:1534732/1‑90 (MQ=255)
                                                aaCTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGAATCCTGAAAATCACCTATAACTAATATCTAATACa                                                                                                                          <  2:1534732/90‑1 (MQ=255)
                                                     aGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGATATATCAGTGCTCTAACTCCAGccct                                                                     >  2:1201569/1‑139 (MQ=255)
                                                                 gCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGcc                                                                                                  <  1:1442441/98‑1 (MQ=255)
                                                                 gCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGcc                                                                                                  >  2:1442441/1‑98 (MQ=255)
                                                                            tgttgAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCa                                                                                        >  1:455324/1‑97 (MQ=255)
                                                                            tgttgAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCa                                                                                        <  2:455324/97‑1 (MQ=255)
                                                                                gAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATAc                                          <  1:10526/139‑1 (MQ=255)
                                                                                gAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTACCTGCATCATCTGTTGACAATAc                                          >  2:1332083/1‑139 (MQ=255)
                                                                                   ggACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAgg                                                                                                    >  1:991164/1‑78 (MQ=255)
                                                                                   ggACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAgg                                                                                                    <  2:991164/78‑1 (MQ=255)
                                                                                        gTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATa                                  >  1:1189701/1‑139 (MQ=255)
                                                                                                   aCGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGtt                                                  <  2:367902/112‑1 (MQ=255)
                                                                                                   aCGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGtt                                                  >  1:367902/1‑112 (MQ=255)
                                                                                                   aCGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCaaa                       >  1:280397/1‑139 (MQ=255)
                                                                                                     gTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTACCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAAc                     <  1:1332083/139‑1 (MQ=255)
                                                                                                        aTCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAg                  >  1:771756/1‑139 (MQ=255)
                                                                                                        aTCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAg                  <  2:771756/139‑1 (MQ=255)
                                                                                                                   cACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAgg                                                                                                    <  1:93263/46‑1 (MQ=255)
                                                                                                                   cACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAgg                                                                                                    >  2:93263/1‑46 (MQ=255)
                                                                                                                        ataACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGctctc                                                                    <  1:506422/73‑1 (MQ=255)
                                                                                                                        ataACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGctctc                                                                    >  2:506422/1‑73 (MQ=255)
                                                                                                                        ataACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAGGTTCTATGTTCTTCAt  <  2:280397/139‑1 (MQ=255)
                                                                                                                               |                                                                                                                                   
CTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAGGTTCTATGTTCTTCAT  >  NC_000913/3191569‑3191827

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