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% 275.6 / NA 93*247K (TAA→AAA) ecpDpredicted periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (41/52);  total (41/52)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAG  >  W3110S.gb/155396‑155528
                                                                   |                                                                 
gatgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAg                                                                >  1:1926018/1‑71 (MQ=255)
 atgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                               >  1:2691026/1‑71 (MQ=255)
  tgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGcc                                                              <  1:3498766/71‑1 (MQ=255)
  tgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGcc                                                              <  1:3942974/71‑1 (MQ=255)
  tgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                               <  1:1339529/70‑1 (MQ=255)
    atAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                             <  1:500813/70‑1 (MQ=255)
     tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                           >  1:3938113/1‑71 (MQ=255)
      aaTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                           >  1:2343506/1‑70 (MQ=255)
      aaTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTg                                                            >  1:2552369/1‑69 (MQ=255)
       aTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                             <  1:1998896/67‑1 (MQ=255)
       aTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGc                                                          <  1:3606994/70‑1 (MQ=255)
        tttttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAt                                                        <  1:3279533/71‑1 (MQ=255)
          tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                       <  1:2789861/70‑1 (MQ=255)
          tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                       <  1:3339487/70‑1 (MQ=255)
          tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                       <  1:4204251/70‑1 (MQ=255)
          tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                       <  1:1933117/70‑1 (MQ=255)
             aGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                             >  1:1491982/1‑61 (MQ=255)
              gTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCtt                                                  <  1:585745/71‑1 (MQ=255)
              gTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCt                                                   <  1:3658011/70‑1 (MQ=255)
              gTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCt                                                   <  1:495260/70‑1 (MQ=255)
               tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                 <  1:974482/71‑1 (MQ=255)
                 tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                 >  1:2008346/1‑53 (MQ=255)
                 tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                               >  1:2959542/1‑71 (MQ=255)
                  atTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAcc                                                    <  1:720695/65‑1 (MQ=255)
                  atTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAcc                                                    <  1:4235578/65‑1 (MQ=255)
                   ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                             <  1:1230917/71‑1 (MQ=255)
                   ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                             <  1:554502/71‑1 (MQ=255)
                   ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                             <  1:3207486/71‑1 (MQ=255)
                    tCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCa                                                         >  1:1593883/1‑58 (MQ=255)
                    tCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCa                                                         >  1:4150586/1‑58 (MQ=255)
                    tCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCa                                                         >  1:515488/1‑58 (MQ=255)
                    tCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTg                                            >  1:3602366/1‑71 (MQ=255)
                     cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                       <  1:304383/59‑1 (MQ=255)
                     cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTg                                            <  1:419268/70‑1 (MQ=255)
                     cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTg                                            <  1:624065/70‑1 (MQ=255)
                        taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                         >  1:445007/1‑70 (MQ=255)
                        taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgc                                           <  1:1484913/68‑1 (MQ=255)
                             aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                    >  1:3269153/1‑70 (MQ=255)
                              cGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                    <  1:671357/69‑1 (MQ=255)
                              cGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                    >  1:22050/1‑69 (MQ=255)
                               gCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAg                                 <  1:2595580/71‑1 (MQ=255)
                               gCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAg                                 <  1:3252481/71‑1 (MQ=255)
                               gCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAg                                 <  1:3247473/71‑1 (MQ=255)
                               gCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAg                                 >  1:2512103/1‑71 (MQ=255)
                               gCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCACCCTGGCAt                                                        >  1:2433942/1‑48 (MQ=38)
                                cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                    <  1:2076433/67‑1 (MQ=255)
                                  aTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                 <  1:2373569/52‑1 (MQ=255)
                                  aTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                              >  1:2430755/1‑71 (MQ=255)
                                   tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAt                             >  1:3197199/1‑71 (MQ=255)
                                   tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAt                             >  1:1692516/1‑71 (MQ=255)
                                    gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCa                                                         >  1:3080722/1‑42 (MQ=255)
                                    gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAACTTCAATTGCGCCACCAAAGTCAt                             >  1:3473748/1‑70 (MQ=255)
                                     cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                               >  1:87051/1‑35 (MQ=255)
                                     cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGcca                                       >  1:176932/1‑59 (MQ=255)
                                     cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            <  1:1527209/70‑1 (MQ=255)
                                     cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            <  1:316538/70‑1 (MQ=255)
                                     cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                           <  1:966103/71‑1 (MQ=255)
                                     cTTCCTGCGGAATTATATCCTTGCCTGATTTAAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                           <  1:1166189/71‑1 (MQ=255)
                                        ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAAt                         <  1:1309137/70‑1 (MQ=255)
                                         cTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                        <  1:2032957/70‑1 (MQ=255)
                                           gCGGAATTTATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCg                     >  1:2565075/1‑70 (MQ=255)
                                            cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCg                     >  1:3756945/1‑70 (MQ=255)
                                              gAATTATATCCTTGCCTGATTTCAGCCTGGGATTACCTTCAATTGCGCCACCAAAGTCATTa                           <  1:2601861/62‑1 (MQ=255)
                                                aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaa                <  1:2149242/71‑1 (MQ=255)
                                                aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaa                <  1:3944760/71‑1 (MQ=255)
                                                aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                 >  1:1675499/1‑70 (MQ=255)
                                                 ttATATCCTTTCCTGATTTCAGCCTGGCATTACCTTc                                                 <  1:1404393/37‑1 (MQ=38)
                                                 ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacc                                     <  1:3664381/49‑1 (MQ=255)
                                                   atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                             <  1:2570109/39‑1 (MQ=255)
                                                   atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAAt                         >  1:2640621/1‑59 (MQ=255)
                                                    tatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCGACCAAAGTCATTAATGGCGTAAAAATg             <  1:3494942/70‑1 (MQ=255)
                                                    tatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg             >  1:198436/1‑70 (MQ=255)
                                                       ccTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                             <  1:1584536/35‑1 (MQ=255)
                                                       ccTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                             <  1:2034999/35‑1 (MQ=255)
                                                       ccTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                              >  1:995937/1‑50 (MQ=255)
                                                       ccTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACt         >  1:3072296/1‑71 (MQ=255)
                                                       ccTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACt         >  1:2003358/1‑71 (MQ=255)
                                                       ccTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACt         >  1:3067456/1‑71 (MQ=255)
                                                        cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtt        <  1:4204971/71‑1 (MQ=255)
                                                        cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtt        <  1:3119816/71‑1 (MQ=255)
                                                         ttGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            <  1:3569052/50‑1 (MQ=255)
                                                         ttGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            <  1:1355925/50‑1 (MQ=255)
                                                          tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                          >  1:1817744/1‑51 (MQ=255)
                                                          tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                          >  1:2716292/1‑51 (MQ=255)
                                                          tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                          >  1:3563484/1‑51 (MQ=255)
                                                          tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtt        >  1:2156488/1‑69 (MQ=255)
                                                           gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATgg                       >  1:2396614/1‑53 (MQ=255)
                                                           gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg     >  1:1162593/1‑71 (MQ=255)
                                                              tGATTTCAGCCTGGCATTACCTTCAATTGCGccacc                                     <  1:4017127/36‑1 (MQ=255)
                                                              tGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                                <  1:4175720/41‑1 (MQ=255)
                                                              tGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                                <  1:1590620/41‑1 (MQ=255)
                                                              tGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTa                   >  1:3685414/1‑54 (MQ=255)
                                                               gATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg  >  1:260000/1‑70 (MQ=255)
                                                                   |                                                                 
GATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAG  >  W3110S.gb/155396‑155528

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