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% 231.8 / NA 75*247K (TAA→AAA) ecpDpredicted periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (42/33);  total (42/33)
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:1745173/71‑1 (MQ=255)
    gatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                               >  1:3905416/1‑71 (MQ=255)
    gatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                               >  1:3884298/1‑71 (MQ=255)
      tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                             <  1:1972248/71‑1 (MQ=255)
       aaTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                             >  1:4334024/1‑70 (MQ=255)
       aaTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                             >  1:1741719/1‑70 (MQ=255)
           tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                         <  1:1278662/70‑1 (MQ=255)
           tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                         <  1:1457848/70‑1 (MQ=255)
           tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTa                                                        >  1:2540814/1‑71 (MQ=255)
            ttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAc                                                       >  1:3475535/1‑71 (MQ=255)
                  tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                         <  1:4198179/63‑1 (MQ=255)
                  tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                         <  1:2685215/63‑1 (MQ=255)
                  tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                         <  1:3230620/63‑1 (MQ=255)
                  tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                         <  1:4493971/63‑1 (MQ=255)
                  tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                         <  1:3677346/63‑1 (MQ=255)
                  tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCa                                                  <  1:1010511/70‑1 (MQ=255)
                    ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                                 >  1:2592828/1‑53 (MQ=255)
                    ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                               >  1:1890067/1‑71 (MQ=255)
                    ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                               >  1:73918/1‑71 (MQ=255)
                    ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                               >  1:1408646/1‑71 (MQ=255)
                    ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                               >  1:3778880/1‑71 (MQ=255)
                     tCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                                 >  1:4846355/1‑52 (MQ=255)
                     tCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                                 >  1:702901/1‑52 (MQ=255)
                     tCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                         <  1:2812942/60‑1 (MQ=255)
                     tCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                               >  1:1697683/1‑70 (MQ=255)
                      cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                               >  1:1695708/1‑53 (MQ=255)
                      cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                               >  1:3732344/1‑53 (MQ=255)
                      cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                               >  1:1523180/1‑53 (MQ=255)
                      cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgc                                             >  1:488872/1‑71 (MQ=255)
                         taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcg                                            <  1:2489581/69‑1 (MQ=255)
                         taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcg                                            <  1:1307909/69‑1 (MQ=255)
                           gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGc                                                            >  1:4433781/1‑51 (MQ=255)
                             cACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                   <  1:2058608/42‑1 (MQ=255)
                              aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                   <  1:4049443/57‑1 (MQ=255)
                              aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                      <  1:4801423/70‑1 (MQ=255)
                                gCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAg                                   >  1:4357885/1‑71 (MQ=255)
                                gCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAg                                   >  1:4375340/1‑71 (MQ=255)
                                gCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAg                                   >  1:3950448/1‑71 (MQ=255)
                                  tATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                   <  1:679667/53‑1 (MQ=255)
                                   aTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                                >  1:4885478/1‑71 (MQ=255)
                                   aTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                                >  1:3796880/1‑71 (MQ=255)
                                    tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                                    <  1:4412866/66‑1 (MQ=255)
                                     gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAt                               >  1:5067066/1‑70 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                              <  1:885303/70‑1 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                              >  1:3414494/1‑70 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                              <  1:1390697/70‑1 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                             <  1:1277287/71‑1 (MQ=255)
                                         ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                            <  1:2777737/69‑1 (MQ=255)
                                         ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                          >  1:740349/1‑71 (MQ=255)
                                          cTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGcca                                         <  1:100766/55‑1 (MQ=255)
                                          cTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                          <  1:1346708/70‑1 (MQ=255)
                                            gCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacc                                       >  1:2042841/1‑55 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGt                      <  1:2373605/71‑1 (MQ=255)
                                              ggAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTa                     <  1:3829310/71‑1 (MQ=255)
                                                aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccac                                        <  1:3631280/50‑1 (MQ=255)
                                                 aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                   >  1:1007841/1‑70 (MQ=255)
                                                  ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTc                                 <  1:3712420/55‑1 (MQ=255)
                                                  ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa                 <  1:5116695/71‑1 (MQ=255)
                                                    atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                                 >  1:1122016/1‑37 (MQ=255)
                                                     tatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg               <  1:619343/70‑1 (MQ=255)
                                                     tatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg               >  1:2194849/1‑70 (MQ=255)
                                                       tCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTg                                              >  1:2272993/1‑37 (MQ=255)
                                                           tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtttt        >  1:3109086/1‑71 (MQ=255)
                                                           tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtttt        >  1:5030008/1‑71 (MQ=255)
                                                           tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtttt        >  1:2543299/1‑71 (MQ=255)
                                                            gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACttt         <  1:2036185/69‑1 (MQ=255)
                                                            gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg       >  1:5100526/1‑71 (MQ=255)
                                                            gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg       >  1:649651/1‑71 (MQ=255)
                                                                gATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                                >  1:4820095/1‑42 (MQ=255)
                                                                 aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATgg                         >  1:3893255/1‑48 (MQ=255)
                                                                 aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCa             >  1:2584329/1‑60 (MQ=255)
                                                                 aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg    <  1:801987/69‑1 (MQ=255)
                                                                 aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGaa  <  1:1887016/71‑1 (MQ=255)
                                                                 aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGaa  >  1:1258608/1‑71 (MQ=255)
                                                                 aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGaa  >  1:4592472/1‑71 (MQ=255)
                                                                    |                                                                   
TGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAA  >  W3110S.gb/155395‑155530

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