Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 807,445 C→A A155A (GCG→GCT ybhB ← putative kinase inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913807,4450CA100.0% 112.0 / NA 31A155A (GCG→GCTybhBputative kinase inhibitor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (23/8);  total (23/8)

AGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATCGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCTTCATCGA  >  NC_000913/807367‑807521
                                                                              |                                                                            
aGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATAGc                                                                            <  1:165016/81‑1 (MQ=255)
  aCAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATAGCAg                                                                          <  1:165017/81‑1 (MQ=255)
        aGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATAGCAGTAATc                                                                     >  2:165045/1‑80 (MQ=255)
          ttgttTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATAGCAGTAATCGAg                                                                  >  2:165039/1‑81 (MQ=255)
                ttAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATAGCAGTAATCGAGGCGCtt                                                            >  1:165054/1‑81 (MQ=255)
                    tACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATAGCAGTAATCGAGGCGCTTGc                                                          >  2:165043/1‑79 (MQ=255)
                    tACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATAGCAGTAATCGAGGCGCTTGCCa                                                        <  2:165018/81‑1 (MQ=255)
                       cTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATAGCAGTAATCGAGGCGCTTGCCagag                                                     <  1:165019/81‑1 (MQ=255)
                       cTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATAGCAGTAATCGAGGCGCTTGCCa                                                        >  2:165044/1‑78 (MQ=255)
                          tAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATAGCAGTAATCGAGGCGCTTGCCagagag                                                   <  2:165020/80‑1 (MQ=255)
                          tAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATAGCAGTAATCGAGGCGCTTGCCAGAGAGt                                                  >  2:165059/1‑81 (MQ=255)
                              tGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATAGCAGTAATCGAGGCGCTTGCCAGAGAGTGGa                                               >  1:165041/1‑80 (MQ=255)
                                  aCCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATAGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATg                                           >  1:165093/1‑80 (MQ=255)
                                   ccAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATAGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGaa                                         >  2:165052/1‑81 (MQ=255)
                                    cAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATAGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAAc                                        >  1:165049/1‑81 (MQ=255)
                                         ggCATTGCGCCATCTGGCAGAGTGATTAACTAAACATAGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTa                                    >  1:165063/1‑80 (MQ=255)
                                         ggCATTGCGCCATCTGGCAGAGTGATTAACTAAACATAGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTa                                    >  2:165056/1‑80 (MQ=255)
                                          gCATTGCGCCATCTGGCAGAGTGATTAACTAAACATAGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTaaa                                  >  1:165050/1‑81 (MQ=255)
                                                cgcCATCTGGCAGAGTGATTAACTAAACATAGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGAc                            >  2:165048/1‑81 (MQ=255)
                                                  ccATCTGGCAGAGTGATTAACTAAACATAGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCa                          >  2:165047/1‑81 (MQ=255)
                                                   cATCTGGCAGAGTGATTAACTAAACATAGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCAt                         <  1:833238/81‑1 (MQ=255)
                                                    aTCTGGCAGAGTGATTAACTAAACATAGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATc                        >  1:165087/1‑81 (MQ=255)
                                                       tGGCAGAGTGATTAACTAAACATAGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATcgcg                     >  2:165046/1‑81 (MQ=255)
                                                         gCAGAGTGATTAACTAAACATAGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGcgcc                   >  1:165051/1‑81 (MQ=255)
                                                          cAGAGTGATTAACTAAACATAGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGcgccg                  >  1:165055/1‑81 (MQ=255)
                                                          cAGAGTGATTAACTAAACATAGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGcgccg                  >  1:165053/1‑81 (MQ=255)
                                                            gagTGATTAACTAAACATAGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCt                >  1:165066/1‑81 (MQ=255)
                                                                    aaCTAAACATAGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCt         <  1:165021/80‑1 (MQ=255)
                                                                         aaCATAGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCTTCATCg   <  1:165022/81‑1 (MQ=255)
                                                                           cATAGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGGCGCTGGCACCTTCATCGa  >  2:165070/1‑80 (MQ=255)
                                                                           cATAGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCTTCATCGa  >  2:165072/1‑80 (MQ=255)
                                                                              |                                                                            
AGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATCGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCTTCATCGA  >  NC_000913/807367‑807521

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