Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,695,708 C→T 33.8% L372F (CTT→TTT)  yfhD → predicted transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,695,7080CT33.8% 86.6 / 65.0 71L372F (CTT→TTT) yfhDpredicted transglycosylase
Reads supporting (aligned to +/- strand):  ref base C (21/26);  new base T (11/13);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

CCCGCGCCCTGACGGCAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACATGAAGCCTACGCT  >  W3110S.gb/2695637‑2695773
                                                                       |                                                                 
cccGCGCCCTGACGGCAAAAACCAAAGGGAATCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTAc                                                                   >  1:4597873/1‑71 (MQ=255)
   gcgcCCTGACGGCAAAAACCAACGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTAttt                                                                 <  1:1011928/71‑1 (MQ=255)
    cgcCCTGACGGCAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTa                                                                <  1:2755317/71‑1 (MQ=255)
     gcCCTGACGGCAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTa                                                                <  1:5126408/70‑1 (MQ=255)
       ccTGACGGCAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGc                                                              >  1:1911463/1‑70 (MQ=255)
        cTGACGGCAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACtt                                                                 >  1:3432992/1‑66 (MQ=255)
        cTGACGGCAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGc                                                              >  1:3331874/1‑69 (MQ=255)
          gACGGCAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTTGCCAg                                                           >  1:2547891/1‑70 (MQ=255)
          gACGGCAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAg                                                           >  1:161212/1‑70 (MQ=255)
           aCGGCAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGa                                                          <  1:4341567/70‑1 (MQ=255)
              gCAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAAcc                                                      <  1:21572/71‑1 (MQ=255)
               cAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTATTTAGCCAGAAAccc                                                     >  1:874980/1‑71 (MQ=255)
                 aaaaCCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTATTTAGCCAGAAACCCTa                                                   >  1:148547/1‑71 (MQ=255)
                 aaaaCCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTa                                                   >  1:1476423/1‑71 (MQ=255)
                 aaaaCCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTa                                                   <  1:2103305/71‑1 (MQ=255)
                  aaaCCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTATTTAGc                                                              <  1:1235602/59‑1 (MQ=255)
                  aaaCCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTAt                                                  <  1:4901684/71‑1 (MQ=255)
                  aaaCCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTAt                                                  <  1:4791759/71‑1 (MQ=255)
                   aaCCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTATTTAGCCAGAAACCCTAt                                                  >  1:3697597/1‑70 (MQ=255)
                     ccAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAAccc                                                     <  1:4852661/65‑1 (MQ=255)
                     ccAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAAccc                                                     <  1:3180293/65‑1 (MQ=255)
                     ccAAAGCGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAAccc                                                     <  1:2433966/65‑1 (MQ=255)
                      cAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGaa                                                         >  1:4189318/1‑60 (MQ=255)
                         aGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCa                                           <  1:394680/71‑1 (MQ=255)
                         aGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCa                                           <  1:3431341/71‑1 (MQ=255)
                            gAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTATTTAGCCAGAAACCCTATTACAGCAAGc                                        <  1:4127183/71‑1 (MQ=255)
                                 cTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACt                                   >  1:3953291/1‑71 (MQ=255)
                                 cTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACt                                   >  1:104117/1‑71 (MQ=255)
                                            gCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCt                                                    <  1:476740/43‑1 (MQ=255)
                                            gCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCt                                                    <  1:2879714/43‑1 (MQ=255)
                                            gCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCt                                                    <  1:1784421/43‑1 (MQ=255)
                                             cTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCt                       <  1:1974823/71‑1 (MQ=255)
                                                aCGTAAAACAGCGTCTGCCTTTATTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCg                     <  1:5123547/70‑1 (MQ=255)
                                                aCGTAAAACAGCGTCTGCCTTTATTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGt                    >  1:439158/1‑71 (MQ=255)
                                                aCGTAAAACAGCGTCTGCCTTTATTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGt                    >  1:112804/1‑71 (MQ=255)
                                                aCGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCt                                       <  1:3768346/52‑1 (MQ=255)
                                                 cGTAAAACAGCGTCTGCCTTTATTTAGCCAGAAACCCTATTa                                                >  1:4702826/1‑42 (MQ=255)
                                                 cGTAAAACAGCGTCTGCCTTTATTTAGCCAGAAACCCTATTa                                                >  1:5042091/1‑42 (MQ=255)
                                                 cGTAAAACAGCGTCTGCCTTTATTTAGCCAGAAACCCTATTa                                                >  1:998783/1‑42 (MQ=255)
                                                 cGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTg                   <  1:3205072/71‑1 (MQ=255)
                                                  gTAAAACAGCGTCTGCCTTTATTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCg                     >  1:3193131/1‑68 (MQ=255)
                                                  gTAAAACAGCGTCTGCCTTTATTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCg                     >  1:2412447/1‑68 (MQ=255)
                                                  gTAAAACAGCGTCTGCCTTTATTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTgg                  <  1:173888/71‑1 (MQ=255)
                                                  gTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTa                                                   >  1:3332929/1‑38 (MQ=255)
                                                  gTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGa                                     >  1:1308051/1‑52 (MQ=255)
                                                  gTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGa                                     >  1:3923019/1‑52 (MQ=255)
                                                   tAAAACAGCGTCTGCCTTTATTTAGCCAGAAACCCTATTACAGCAAGCTg                                      <  1:4578695/50‑1 (MQ=255)
                                                   tAAAACAGCGTCTGCCTTTATTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGa                 <  1:2929498/71‑1 (MQ=255)
                                                   tAAAACAGCGTCTGCCTTTATTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGa                 <  1:2142563/71‑1 (MQ=255)
                                                   tAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCaa                                          <  1:5072481/46‑1 (MQ=255)
                                                   tAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGa                 <  1:4628694/71‑1 (MQ=255)
                                                   tAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGa                 <  1:632859/71‑1 (MQ=255)
                                                   tAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGa                 <  1:3155116/71‑1 (MQ=255)
                                                     aaaCAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACa               >  1:1714094/1‑71 (MQ=255)
                                                      aaCAGCGTCTGCCTTTATTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCt                       <  1:4797070/62‑1 (MQ=255)
                                                       aCAGCGTCTGCCTTTATTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACAt              <  1:2360951/70‑1 (MQ=255)
                                                       aCAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACATg             <  1:2638913/71‑1 (MQ=255)
                                                       aCAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACATg             <  1:62010/71‑1 (MQ=255)
                                                           cGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACATGAAGc         >  1:4400998/1‑71 (MQ=255)
                                                            gTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACATGAAGcc        >  1:4495582/1‑71 (MQ=255)
                                                             tCTGCCTTTATTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACa               >  1:4026104/1‑63 (MQ=255)
                                                              cTGCCTTTATTTAGCCAGAAACCCTATTACAGCAAGCTGACt                                   <  1:4392436/42‑1 (MQ=255)
                                                                gCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACATGAAGCCTACg    >  1:2499719/1‑71 (MQ=255)
                                                                  cTTTATTTAGCCAGAAGCCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACATGAAGCCTACGCt  <  1:2166947/71‑1 (MQ=255)
                                                                  cTTTATTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACATGAAGCCTACGCt  <  1:2936279/71‑1 (MQ=255)
                                                                  cTTTATTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACATGAAGCCTACGCt  <  1:586390/71‑1 (MQ=255)
                                                                  cTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACATGAAGCCTACGCt  <  1:1311093/71‑1 (MQ=255)
                                                                    ttACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACATGAAGCCTACGCt  <  1:2746567/69‑1 (MQ=255)
                                                                    ttACTAAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTc                      >  1:3295468/1‑49 (MQ=255)
                                                                     tACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTa                           >  1:2823186/1‑43 (MQ=255)
                                                                     tACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTa                           >  1:4910324/1‑43 (MQ=255)
                                                                     tACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTa                           >  1:2422944/1‑43 (MQ=255)
                                                                     tACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTa                           >  1:2645765/1‑43 (MQ=255)
                                                                       |                                                                 
CCCGCGCCCTGACGGCAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACATGAAGCCTACGCT  >  W3110S.gb/2695637‑2695773

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: