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

GATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTT  >  W3110S.gb/155396‑155524
                                                                   |                                                             
gatgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAg                                                            >  1:2521516/1‑71 (MQ=255)
 atgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                           <  1:674020/71‑1 (MQ=255)
 atgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                           <  1:693444/71‑1 (MQ=255)
  tgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGcc                                                          >  1:4435239/1‑71 (MQ=255)
     tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                       <  1:1498045/71‑1 (MQ=255)
     tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                       >  1:4155126/1‑71 (MQ=255)
       aTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCa                                                     <  1:4091220/71‑1 (MQ=255)
       aTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCa                                                     >  1:2900144/1‑71 (MQ=255)
        tttttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                           >  1:2643306/1‑64 (MQ=255)
        tttttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                           >  1:2614721/1‑64 (MQ=255)
        tttttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAt                                                    >  1:1135638/1‑71 (MQ=255)
          tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTg                                                        >  1:457341/1‑65 (MQ=255)
           ttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTa                                                  <  1:1729412/70‑1 (MQ=255)
              gTATATTCTATAGTCACGCTATGCTTCCTGCCGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCt                                               <  1:819599/70‑1 (MQ=255)
                 tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                             <  1:2336675/69‑1 (MQ=255)
                 tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                           >  1:1598097/1‑71 (MQ=255)
                  atTCTATAGTCACGCTATGCTTCCCGCGGAATTATATCCTTGCCTGATTTCAGc                                                           <  1:2773178/54‑1 (MQ=39)
                   ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCa                                                     >  1:1763712/1‑59 (MQ=255)
                      tataGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgc                                       >  1:80926/1‑70 (MQ=255)
                       ataGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                     >  1:2297089/1‑71 (MQ=255)
                       ataGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                     >  1:3000383/1‑71 (MQ=255)
                        taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                     <  1:1537417/70‑1 (MQ=255)
                          gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGAGccac                                  >  1:534104/1‑71 (MQ=255)
                           tCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTg                                                        <  1:4602688/48‑1 (MQ=255)
                           tCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTg                                                        <  1:5143838/48‑1 (MQ=255)
                             aCGCTATGCTTCCTGCGGAATTATATCCTTGCGTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                >  1:5136529/1‑70 (MQ=255)
                             aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                >  1:5020874/1‑70 (MQ=255)
                              cGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                >  1:2370162/1‑69 (MQ=255)
                              cGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTAAGCCTGGCATTACCTTCAATTGCGccacca                                <  1:3821806/69‑1 (MQ=255)
                                cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAc                                                 >  1:1767534/1‑50 (MQ=255)
                                cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCa                                            <  1:4032550/55‑1 (MQ=255)
                                cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                            <  1:4452590/71‑1 (MQ=255)
                                cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                            >  1:2880745/1‑71 (MQ=255)
                                cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                            <  1:705697/71‑1 (MQ=255)
                                  aTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                   <  1:4598986/46‑1 (MQ=255)
                                  aTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCa                                            <  1:1256995/53‑1 (MQ=255)
                                  aTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                          >  1:4801029/1‑71 (MQ=255)
                                   tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                          <  1:717853/70‑1 (MQ=255)
                                    gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAt                                                    >  1:4568452/1‑43 (MQ=255)
                                    gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAt                         >  1:1156278/1‑70 (MQ=255)
                                     cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                        <  1:4654939/70‑1 (MQ=255)
                                     cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAt                         >  1:1561823/1‑69 (MQ=255)
                                       tCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAAt                     >  1:2467740/1‑71 (MQ=255)
                                        ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                      <  1:862938/69‑1 (MQ=255)
                                        ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAAt                     <  1:3582557/70‑1 (MQ=255)
                                        ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAAt                     <  1:4684531/70‑1 (MQ=255)
                                        ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                    >  1:645133/1‑71 (MQ=255)
                                        ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                    >  1:818280/1‑71 (MQ=255)
                                          tGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAt                         >  1:1145141/1‑64 (MQ=255)
                                           gCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                         >  1:3719795/1‑47 (MQ=255)
                                              gAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaa                               <  1:452846/54‑1 (MQ=255)
                                               aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGt                >  1:2128903/1‑68 (MQ=255)
                                                aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaa            <  1:512657/71‑1 (MQ=255)
                                                 ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa           >  1:926218/1‑71 (MQ=255)
                                                 ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa           >  1:2840583/1‑71 (MQ=255)
                                                   atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg         >  1:4128136/1‑71 (MQ=255)
                                                    tatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                <  1:292710/47‑1 (MQ=255)
                                                    tatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                        <  1:1943830/55‑1 (MQ=255)
                                                      tCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCAc      >  1:3586355/1‑71 (MQ=255)
                                                       ccTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACt     >  1:504924/1‑71 (MQ=255)
                                                       ccTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACt     >  1:1540939/1‑71 (MQ=255)
                                                        cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtt    >  1:4823852/1‑71 (MQ=255)
                                                        cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtt    <  1:4659914/71‑1 (MQ=255)
                                                          tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtttt  >  1:2745331/1‑71 (MQ=255)
                                                           gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                       >  1:2233230/1‑49 (MQ=255)
                                                             cTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg         <  1:1891678/61‑1 (MQ=255)
                                                                   |                                                             
GATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTT  >  W3110S.gb/155396‑155524

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: