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% 119.4 / NA 34Y657H (TAT→CAT) ybiWputative pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/17);  total (17/17)

ATGGTCAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCG  >  NC_000913/860510‑860763
                                                                                                                                |                                                                                                                             
aTGGTCAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCa                                                                                                                     >  1:207517/1‑139 (MQ=255)
        tACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGAcc                                                                                                             >  2:752907/1‑139 (MQ=255)
                    ggCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAgcg                                                                                                 <  1:14863/139‑1 (MQ=255)
                                cTGCCAGCGGGCTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTa                                                                                     <  2:1086317/139‑1 (MQ=255)
                                           gtgtGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAgcg                                                                                                 >  1:1661069/1‑116 (MQ=255)
                                           gtgtGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAgcg                                                                                                 <  2:1661069/116‑1 (MQ=255)
                                                   tttACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                                  <  2:207517/139‑1 (MQ=255)
                                                                            gTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAgcg                                                                                                 >  2:1484725/1‑83 (MQ=255)
                                                                            gTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAgcg                                                                                                 <  1:1484725/83‑1 (MQ=255)
                                                                                 cgcgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGcgcg                                                                                               >  1:1556607/1‑80 (MQ=255)
                                                                                 cgcgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGcgcg                                                                                               <  2:1556607/80‑1 (MQ=255)
                                                                                     ccAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGa                                <  1:212379/139‑1 (MQ=255)
                                                                                      cAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTgcc                                                                                                                         >  1:1499799/1‑49 (MQ=255)
                                                                                      cAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTgcc                                                                                                                         <  2:1499799/49‑1 (MQ=255)
                                                                                                gTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGcgcg                                                                                               <  2:1208221/65‑1 (MQ=255)
                                                                                                gTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGcgcg                                                                                               >  1:1208221/1‑65 (MQ=255)
                                                                                                     cGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGtt                >  2:508202/1‑139 (MQ=255)
                                                                                                       gagaTTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGcgcg                                                                                               >  2:1237675/1‑58 (MQ=255)
                                                                                                       gagaTTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGcgcg                                                                                               <  1:1237675/58‑1 (MQ=255)
                                                                                                       gagaTTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                                  <  1:61153/87‑1 (MQ=255)
                                                                                                       gagaTTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                                  >  2:61153/1‑87 (MQ=255)
                                                                                                       gagaTTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGc              >  1:508154/1‑139 (MQ=255)
                                                                                                          aTTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGt           >  2:1308525/1‑139 (MQ=255)
                                                                                                              atgaCGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGt                                                             <  1:606465/85‑1 (MQ=255)
                                                                                                              atgaCGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGt                                                             >  2:606465/1‑85 (MQ=255)
                                                                                                                   gTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGgcg  <  1:508202/139‑1 (MQ=255)
                                                                                                                    tACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTtcgt                                                                       >  2:1381545/1‑69 (MQ=255)
                                                                                                                    tACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTtcgt                                                                       <  1:1381545/69‑1 (MQ=255)
                                                                                                                          cgTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCcagc                                          <  2:1241627/92‑1 (MQ=255)
                                                                                                                          cgTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCcagc                                          >  1:1241627/1‑92 (MQ=255)
                                                                                                                           gTAATGGTTGCCGCCAACCGGACCACGACCGTAGcgcg                                                                                               <  1:801807/38‑1 (MQ=38)
                                                                                                                           gTAATGGTTGCCGCCAACCGGACCACGACCGTAGcgcg                                                                                               >  2:801807/1‑38 (MQ=255)
                                                                                                                             aaTGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCcagc                                          <  2:1617548/89‑1 (MQ=255)
                                                                                                                             aaTGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCcagc                                          >  1:1617548/1‑89 (MQ=255)
                                                                                                                                |                                                                                                                             
ATGGTCAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCG  >  NC_000913/860510‑860763

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