Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 4,643,405 T→C H182H (CAT→CAC creD → inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,643,4050TC100.0% 150.7 / NA 41H182H (CAT→CACcreDinner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/28);  total (13/28)

CGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATAT  >  W3110S.gb/4643350‑4643455
                                                       |                                                  
cggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATccc                                               >  1:3716134/1‑61 (MQ=255)
cggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATccc                                               >  1:1601688/1‑61 (MQ=255)
cggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCt                                              <  1:1017897/62‑1 (MQ=255)
    gcTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTa                                           >  1:2066490/1‑61 (MQ=255)
    gcTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTAc                                          >  1:3103525/1‑62 (MQ=255)
    gcTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTAc                                          >  1:4077855/1‑62 (MQ=255)
    gcTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTAc                                          >  1:2171620/1‑62 (MQ=255)
    gcTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTAc                                          >  1:4532405/1‑62 (MQ=255)
    gcTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTAc                                          >  1:1842160/1‑62 (MQ=255)
    gcTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTAc                                          >  1:1423580/1‑62 (MQ=255)
    gcTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTAc                                          >  1:1088311/1‑62 (MQ=255)
     cTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTAcc                                         <  1:7010/62‑1 (MQ=255)
     cTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTAcc                                         <  1:508923/62‑1 (MQ=255)
     cTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTAcc                                         <  1:742033/62‑1 (MQ=255)
     cTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTAcc                                         <  1:774777/62‑1 (MQ=255)
     cTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTAcc                                         <  1:4332646/62‑1 (MQ=255)
     cTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTAcc                                         <  1:3763232/62‑1 (MQ=255)
     cTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTAcc                                         <  1:3121468/62‑1 (MQ=255)
     cTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTAcc                                         <  1:1047005/62‑1 (MQ=255)
     cTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTAcc                                         <  1:2096225/62‑1 (MQ=255)
     cTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTAcc                                         <  1:2189533/62‑1 (MQ=255)
     cTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTAcc                                         <  1:2261844/62‑1 (MQ=255)
     cTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTAc                                          >  1:425145/1‑61 (MQ=255)
      tGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTACCt                                        <  1:3145462/62‑1 (MQ=255)
                                  gCAAGGCGGGCAGGGCGGGCACATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGa            <  1:568551/62‑1 (MQ=255)
                                   cAAGGCGGGCAGGGCGTGCACATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGa            <  1:840181/61‑1 (MQ=255)
                                     aGGCGGGCAGGGCGTGCACATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGa            <  1:2643141/59‑1 (MQ=255)
                                     aGGCGGGCAGGGCGTGCACATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGa            <  1:408424/59‑1 (MQ=255)
                                     aGGCGGGCAGGGCGTGCACATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGa            <  1:2986371/59‑1 (MQ=255)
                                     aGGCGGGCAGGGCGTGCACATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGa            <  1:3714920/59‑1 (MQ=255)
                                     aGGCGGGCAGGGCGTGCACATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGa            <  1:4521928/59‑1 (MQ=255)
                                     aGGCGGGCAGGGCGTGCACATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGa            <  1:2011543/59‑1 (MQ=255)
                                     aGGCGGGCAGGGCGTGCACATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGa            <  1:510526/59‑1 (MQ=255)
                                     aGGCGGGCAGGGCGTGCACATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGa            <  1:1669148/59‑1 (MQ=255)
                                     aGGCGGGCAGGGCGTGCACATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGa            <  1:1663679/59‑1 (MQ=255)
                                      ggCGGGCAGGGCGTGCACATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGa            <  1:3090253/58‑1 (MQ=255)
                                         gggCAGGGCGTGCACATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGa            <  1:1940182/55‑1 (MQ=255)
                                            cAGGGCGTGCACATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGaa     >  1:3347964/1‑59 (MQ=255)
                                            cAGGGCGTGCACATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGaa     >  1:499271/1‑59 (MQ=255)
                                            cAGGGCGTGCACATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAatat  <  1:4480986/62‑1 (MQ=255)
                                            cAGGGCGTGCACATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAatat  <  1:3015344/62‑1 (MQ=255)
                                                       |                                                  
CGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATAT  >  W3110S.gb/4643350‑4643455

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