Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 860,638 A→G Y657H (TAT→CAT)  ybiW ← putative pyruvate formate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913860,6380AG100.0% 129.2 / NA 36Y657H (TAT→CAT) ybiWputative pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/21);  total (15/21)

GATGGTCAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGTT  >  NC_000913/860509‑860770
                                                                                                                                 |                                                                                                                                    
gATGGTCAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTgccgcc                                                                                                                             <  1:26003/139‑1 (MQ=255)
     tCAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAAc                                                                                                                          <  1:120474/137‑1 (MQ=255)
     tCAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAAc                                                                                                                          >  2:120474/1‑137 (MQ=255)
          aCCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAAc                                                                                                                          >  2:358913/1‑132 (MQ=255)
          aCCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAAc                                                                                                                          <  1:358913/132‑1 (MQ=255)
              gAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGAc                                                                                                               >  1:48899/1‑139 (MQ=255)
               aGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGAcc                                                                                                              >  1:134133/1‑139 (MQ=255)
                   cGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAAc                                                                                                                          >  2:320716/1‑123 (MQ=255)
                   cGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAAc                                                                                                                          <  1:320716/123‑1 (MQ=255)
                   cGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAg                                                                                                          >  2:86470/1‑139 (MQ=255)
                            gcCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTa                                                                                                 <  2:368572/139‑1 (MQ=255)
                                 cTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATCGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTa                                                                                            <  1:5888/139‑1 (MQ=255)
                                          gggTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAg                                                                                   <  2:134133/139‑1 (MQ=255)
                                                    tttACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                                         <  1:225788/139‑1 (MQ=255)
                                                    tttACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                                         <  2:48899/139‑1 (MQ=255)
                                                      tACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTa                                                                                            <  1:6254/118‑1 (MQ=255)
                                                      tACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGACGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTa                                                                                            >  2:6254/1‑118 (MQ=255)
                                                              tCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGAACTGTTTCAGTTCGTCGATATAGGTCTGAt                                                               <  2:686288/139‑1 (MQ=255)
                                                                          aTAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCca                                                   <  1:14993/139‑1 (MQ=255)
                                                                             gTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAAc                                                                                                                          >  2:149122/1‑65 (MQ=255)
                                                                             gTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAAc                                                                                                                          <  1:149122/65‑1 (MQ=255)
                                                                                 gcgcgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGt                                            <  2:79924/139‑1 (MQ=255)
                                                                                  cgcgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTa                                           >  2:576101/1‑139 (MQ=255)
                                                                                  cgcgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTa                                           >  2:441842/1‑139 (MQ=255)
                                                                                    cgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATc                                         <  1:277234/139‑1 (MQ=255)
                                                                                      ccAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGa                                       <  1:86470/139‑1 (MQ=255)
                                                                                                  ttAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcgt                           <  2:752311/139‑1 (MQ=255)
                                                                                                      cGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcgt                           <  2:374859/135‑1 (MQ=255)
                                                                                                      cGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcgt                           >  1:374859/1‑135 (MQ=255)
                                                                                                      cGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGtt                       >  1:294589/1‑139 (MQ=255)
                                                                                                               atgaCGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTAt                                          <  2:333538/111‑1 (MQ=255)
                                                                                                               atgaCGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTAt                                          >  1:333538/1‑111 (MQ=255)
                                                                                                                   cGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGgc          >  1:274134/1‑139 (MQ=255)
                                                                                                                    gTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGgcg         <  1:450990/139‑1 (MQ=255)
                                                                                                                       cccGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCAc      <  2:565741/139‑1 (MQ=255)
                                                                                                                           cgTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGtt  >  2:70983/1‑139 (MQ=255)
                                                                                                                                 |                                                                                                                                    
GATGGTCAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGTT  >  NC_000913/860509‑860770

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