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

ATGGTCAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCG  >  NC_000913/860510‑860768
                                                                                                                                |                                                                                                                                  
aTGGTCAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCa                                                                                                                          >  1:1018361/1‑139 (MQ=255)
  ggTCAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAAc                                                                                                                        <  2:1178983/139‑1 (MQ=255)
     cAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCgg                                                                                                                     <  2:971492/139‑1 (MQ=255)
       gTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGAc                                                                                                                   <  2:73896/139‑1 (MQ=255)
        tACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGAcc                                                                                                                  >  2:9853/1‑139 (MQ=255)
        tACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGAcc                                                                                                                  >  2:240597/1‑139 (MQ=255)
             gAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGAc                                                                                                             >  2:839410/1‑139 (MQ=255)
                  cGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAg                                                                                                        >  2:535961/1‑139 (MQ=255)
                    ggCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCTGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAgcg                                                                                                      <  2:740010/139‑1 (MQ=255)
                      cTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGcgcg                                                                                                    <  1:8712/139‑1 (MQ=255)
                                              tgGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTtcgt                                                                            <  1:839410/139‑1 (MQ=255)
                                                     tACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGAtata                                                                     >  2:393649/1‑139 (MQ=255)
                                                      aCGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAg                                                                    >  1:210620/1‑139 (MQ=255)
                                                                      gCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAgcg                                                        <  2:165765/135‑1 (MQ=255)
                                                                      gCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGATTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAgcg                                                        >  1:165765/1‑135 (MQ=255)
                                                                           aGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCcagc                                               <  1:747466/139‑1 (MQ=255)
                                                                            gTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTa                                                                                               >  1:595620/1‑90 (MQ=255)
                                                                            gTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTa                                                                                               <  2:595620/90‑1 (MQ=255)
                                                                            gTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                                       <  2:1191988/114‑1 (MQ=255)
                                                                            gTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                                       >  1:1191988/1‑114 (MQ=255)
                                                                              cTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCcagcag                                             <  1:755558/138‑1 (MQ=255)
                                                                              cTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCcagcag                                             >  2:755558/1‑138 (MQ=255)
                                                                                 cgcgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATgg                                                                                            <  1:435428/88‑1 (MQ=255)
                                                                                 cgcgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATgg                                                                                            >  2:435428/1‑88 (MQ=255)
                                                                                 cgcgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTa                                         >  1:647284/1‑139 (MQ=255)
                                                                                    gcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCg                                      <  1:535961/139‑1 (MQ=255)
                                                                                     ccAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGa                                     <  2:210620/139‑1 (MQ=255)
                                                                                                   aGCGGAGATTGATGACGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcgtcg                       <  2:343049/139‑1 (MQ=255)
                                                                                                              atgaCGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                                       <  2:1076225/80‑1 (MQ=255)
                                                                                                              atgaCGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                                       >  1:1076225/1‑80 (MQ=255)
                                                                                                              atgaCGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcgtcg                       >  2:660026/1‑128 (MQ=255)
                                                                                                              atgaCGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcgtcg                       <  1:660026/128‑1 (MQ=255)
                                                                                                              atgaCGTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACtt            <  1:393649/139‑1 (MQ=255)
                                                                                                                   gTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGt                                          >  1:936840/1‑104 (MQ=255)
                                                                                                                   gTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGt                                          <  2:936840/104‑1 (MQ=255)
                                                                                                                   gTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGgcg       <  2:1018361/139‑1 (MQ=255)
                                                                                                                   gTACCCGCGTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGgcg       <  1:240597/139‑1 (MQ=255)
                                                                                                                        cgcgTAATGGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCg  >  1:62666/1‑139 (MQ=255)
                                                                                                                                |                                                                                                                                  
ATGGTCAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCG  >  NC_000913/860510‑860768

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