Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,689,625 C→T A213T (GCC→ACC)  bcsZ ← endo‑1,4‑D‑glucanase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,689,6250CT100.0% 92.4 / NA 30A213T (GCC→ACC) bcsZendo‑1,4‑D‑glucanase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/14);  total (16/14)

GTCAGGCATCATGCCTACCCACATGTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCA  >  NC_000913/3689494‑3689736
                                                                                                                                   |                                                                                                               
gTCAGGCATCATGCCTACCCACATGTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTcc                                                                                                          <  1:736250/139‑1 (MQ=255)
           tGCCTACCCACATGTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACg                                                                                               >  1:309740/1‑139 (MQ=255)
                    acaTGTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGt                                                                                      >  2:1049347/1‑139 (MQ=255)
                            aCGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGAt                                                                                          <  1:279992/127‑1 (MQ=255)
                            aCGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGAt                                                                                          >  2:279992/1‑127 (MQ=255)
                                           tAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGgcg                                                               >  1:147811/1‑139 (MQ=255)
                                                gctgATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCa                                                          >  1:194167/1‑139 (MQ=255)
                                                gctgATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCa                                                          >  2:581246/1‑139 (MQ=255)
                                                  tgATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCaaa                                                        <  2:194167/139‑1 (MQ=255)
                                                    aTCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGc                                                      <  1:1049347/139‑1 (MQ=255)
                                                              tttCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGAt                                                                                          <  2:605024/93‑1 (MQ=255)
                                                              tttCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGAt                                                                                          >  1:605024/1‑93 (MQ=255)
                                                                          gCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGt                                <  2:147811/139‑1 (MQ=255)
                                                                              ccAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACg                                                                                               >  2:928782/1‑72 (MQ=255)
                                                                              ccAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACg                                                                                               <  1:928782/72‑1 (MQ=255)
                                                                                    tttGTCTTTCTCATAGCGTACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAgg                      <  1:1135621/139‑1 (MQ=255)
                                                                                             ctcATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTc                                                                                                           <  1:736827/45‑1 (MQ=255)
                                                                                             ctcATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTc                                                                                                           >  2:736827/1‑45 (MQ=255)
                                                                                                  aGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCaa                                                         >  2:430293/1‑90 (MQ=255)
                                                                                                  aGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCaa                                                         <  1:430293/90‑1 (MQ=255)
                                                                                                   gcgcACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCgg                                                    <  1:889016/94‑1 (MQ=255)
                                                                                                   gcgcACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCgg                                                    >  2:889016/1‑94 (MQ=255)
                                                                                                    cgcACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAAc      >  2:1139440/1‑139 (MQ=255)
                                                                                                      cACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATa                                                                                                   >  1:498248/1‑44 (MQ=255)
                                                                                                      cACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATa                                                                                                   <  2:498248/44‑1 (MQ=255)
                                                                                                      cACCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGc    >  2:910026/1‑139 (MQ=255)
                                                                                                       aCCCAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGcc   <  1:581246/139‑1 (MQ=255)
                                                                                                        cccAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCa  >  2:211723/1‑139 (MQ=255)
                                                                                                         ccAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTcgcg                                                                                >  2:1313771/1‑60 (MQ=255)
                                                                                                         ccAGTCTGGCGAAAAACCTTTCGGGGTGGTTTCCAGCAATAAACGTTGATTGGTTTcgcg                                                                                <  1:1313771/60‑1 (MQ=255)
                                                                                                                                   |                                                                                                               
GTCAGGCATCATGCCTACCCACATGTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCA  >  NC_000913/3689494‑3689736

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: