Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 155,463 A→T 100% *247K (TAA→AAA)  ecpD ← predicted periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb155,4630AT100.0% 234.5 / NA 79*247K (TAA→AAA) ecpDpredicted periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (33/46);  total (33/46)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAA  >  W3110S.gb/155395‑155530
                                                                    |                                                                   
tgatgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                   <  1:3795324/71‑1 (MQ=255)
tgatgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                   <  1:1749199/71‑1 (MQ=255)
 gatgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGATTTATATCCTTGCCTGATTTCAg                                                                  >  1:529519/1‑71 (MQ=255)
   tgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGcc                                                                <  1:1405532/71‑1 (MQ=255)
     atAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                               <  1:4994271/70‑1 (MQ=255)
      tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                             <  1:2118229/71‑1 (MQ=255)
       aaTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                             >  1:1860094/1‑70 (MQ=255)
        aTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTg                                                              >  1:1120936/1‑68 (MQ=255)
        aTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGc                                                            <  1:3878336/70‑1 (MQ=255)
        aTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGc                                                            <  1:1512880/70‑1 (MQ=255)
          ttttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                   >  1:886809/1‑61 (MQ=255)
           tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                         <  1:4235123/70‑1 (MQ=255)
           tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                         <  1:2651461/70‑1 (MQ=255)
           tttAGTATATTATATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                         >  1:1273824/1‑70 (MQ=255)
              aGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCt                                                     >  1:2555737/1‑71 (MQ=255)
              aGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCt                                                     >  1:1439550/1‑71 (MQ=255)
               gTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                   <  1:3565966/56‑1 (MQ=255)
               gTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCt                                                     <  1:391336/70‑1 (MQ=255)
                tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                   <  1:2449578/71‑1 (MQ=255)
                  tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAt                                                          >  1:4866837/1‑62 (MQ=255)
                    ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                               >  1:300028/1‑71 (MQ=255)
                    ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt                                                <  1:2441448/70‑1 (MQ=255)
                        ataGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                           >  1:1899197/1‑71 (MQ=255)
                           gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccac                                        >  1:3923181/1‑71 (MQ=255)
                              aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                      >  1:5093934/1‑70 (MQ=255)
                              aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                      <  1:3457743/70‑1 (MQ=255)
                              aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                      >  1:3932178/1‑70 (MQ=255)
                                gCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAt                                                          >  1:3451324/1‑48 (MQ=255)
                                gCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAg                                   >  1:2544487/1‑71 (MQ=255)
                                 cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                                  <  1:2092860/71‑1 (MQ=255)
                                  tATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGcc                                                                <  1:4547126/40‑1 (MQ=255)
                                  tATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                      <  1:1451731/66‑1 (MQ=255)
                                  tATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                                  >  1:1379764/1‑70 (MQ=255)
                                   aTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                   <  1:3888676/36‑1 (MQ=255)
                                   aTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                   <  1:1664576/36‑1 (MQ=255)
                                    tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                   <  1:366973/35‑1 (MQ=255)
                                    tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                                <  1:501794/70‑1 (MQ=255)
                                     gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                               >  1:3246726/1‑38 (MQ=255)
                                     gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAt                               >  1:404227/1‑70 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                              <  1:2582180/70‑1 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                              >  1:1719722/1‑70 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                              <  1:5100117/70‑1 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAt                               >  1:3853735/1‑69 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                             <  1:1486098/71‑1 (MQ=255)
                                       ttCCTGCGGTATTATATCCTTGCCTGATTTCAGCCt                                                               <  1:3971630/36‑1 (MQ=37)
                                       ttCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                            >  1:914611/1‑71 (MQ=255)
                                        tCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAAt                           >  1:3705699/1‑71 (MQ=255)
                                         ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTat                            <  1:3883925/69‑2 (MQ=255)
                                         ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                            <  1:3340115/69‑1 (MQ=255)
                                         ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                            <  1:3334363/69‑1 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTc                                 <  1:566715/60‑1 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTc                                 <  1:3284154/60‑1 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCg                       >  1:373851/1‑70 (MQ=255)
                                               gAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTa                     >  1:3629296/1‑70 (MQ=255)
                                                aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                               <  1:2906494/43‑1 (MQ=255)
                                                aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                               <  1:3415031/43‑1 (MQ=255)
                                                aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCGTTAATGGCGTaaa                   <  1:2609502/71‑1 (MQ=255)
                                                 aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATTGCGTaaaa                  <  1:1214276/71‑1 (MQ=255)
                                                 aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaa                  <  1:2250392/71‑1 (MQ=255)
                                                 aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaa                  <  1:2445593/71‑1 (MQ=255)
                                                  ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                                 >  1:3194288/1‑39 (MQ=255)
                                                  ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                                 >  1:3174595/1‑39 (MQ=255)
                                                  ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgc                                             <  1:2005621/43‑1 (MQ=255)
                                                    atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt                                                >  1:2524913/1‑38 (MQ=255)
                                                    atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt                                                >  1:871753/1‑38 (MQ=255)
                                                    atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                           >  1:2155611/1‑43 (MQ=255)
                                                     tatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg               >  1:4819312/1‑70 (MQ=255)
                                                       tCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacc                                       <  1:1883289/44‑1 (MQ=255)
                                                       tCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaa                  <  1:932968/65‑1 (MQ=255)
                                                         cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtt          >  1:4741448/1‑71 (MQ=255)
                                                         cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACt           <  1:1276265/70‑1 (MQ=255)
                                                            gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACttt         <  1:4944172/69‑1 (MQ=255)
                                                             ccTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg       <  1:2469896/70‑1 (MQ=255)
                                                              cTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGc              <  1:2961288/62‑1 (MQ=255)
                                                              cTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtttt        >  1:939998/1‑68 (MQ=255)
                                                              cTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCa     <  1:4167159/71‑1 (MQ=255)
                                                              cTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCa     <  1:2091880/71‑1 (MQ=255)
                                                               tGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                                >  1:3333398/1‑43 (MQ=255)
                                                                 aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGaa  <  1:4174351/71‑1 (MQ=255)
                                                                    |                                                                   
TGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAA  >  W3110S.gb/155395‑155530

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: