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

Only 100 of 123 total aligned reads displayed.

TGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAG  >  W3110S.gb/155395‑155528
                                                                    |                                                                 
tgatgatAATTTTTAGTATATTCTATAGTCACTCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                 <  1:491132/71‑1 (MQ=255)
tgatgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                 <  1:3957087/71‑1 (MQ=255)
 gatgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAg                                                                >  1:3375160/1‑71 (MQ=255)
     atAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                 <  1:4557090/66‑1 (MQ=255)
     atAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                 <  1:491434/66‑1 (MQ=255)
     atAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                 <  1:581186/66‑1 (MQ=255)
     atAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                 <  1:1143234/66‑1 (MQ=255)
      tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                 <  1:1891690/65‑1 (MQ=255)
      tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                           >  1:4690675/1‑71 (MQ=255)
      tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                           >  1:4309994/1‑71 (MQ=255)
      tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                           >  1:1561039/1‑71 (MQ=255)
      tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                           >  1:4853907/1‑71 (MQ=255)
      tAATTTTTAGTATATTCTATAGTCAAGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                           >  1:3195654/1‑71 (MQ=255)
        aTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGc                                                          <  1:1281182/70‑1 (MQ=255)
        aTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCa                                                         >  1:4579893/1‑71 (MQ=255)
                tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                 <  1:3203189/71‑1 (MQ=255)
                tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                 <  1:3923280/71‑1 (MQ=255)
                    ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACTTTCAAt                                              <  1:722101/70‑1 (MQ=255)
                    ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt                                              <  1:3117371/70‑1 (MQ=255)
                    ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt                                              <  1:1166439/70‑1 (MQ=255)
                     tCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                 >  1:1590879/1‑50 (MQ=255)
                         taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCa                                                <  1:434055/63‑1 (MQ=255)
                         taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCa                                                <  1:268520/63‑1 (MQ=255)
                         taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCa                                                <  1:252733/63‑1 (MQ=255)
                         taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCa                                                <  1:292414/63‑1 (MQ=255)
                         taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCa                                                <  1:5242948/63‑1 (MQ=255)
                           gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                               >  1:4883242/1‑46 (MQ=255)
                           gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                               >  1:1396495/1‑46 (MQ=255)
                           gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                               >  1:1331432/1‑46 (MQ=255)
                           gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt                                              >  1:3766431/1‑63 (MQ=255)
                           gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt                                              >  1:3759905/1‑63 (MQ=255)
                           gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt                                              >  1:3585143/1‑63 (MQ=255)
                           gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt                                              >  1:5445602/1‑63 (MQ=255)
                           gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccac                                      <  1:2379712/71‑1 (MQ=255)
                           gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccac                                      <  1:2283755/71‑1 (MQ=255)
                           gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccac                                      <  1:5170070/71‑1 (MQ=255)
                           gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccac                                      <  1:1246122/71‑1 (MQ=255)
                           gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccac                                      <  1:661836/71‑1 (MQ=255)
                                gCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAg                                 >  1:259932/1‑71 (MQ=255)
                                gCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAg                                 >  1:2916227/1‑71 (MQ=255)
                                gCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAg                                 >  1:1703012/1‑71 (MQ=255)
                                gCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAg                                 >  1:5440112/1‑71 (MQ=255)
                                  tATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTg                                                            >  1:5317237/1‑42 (MQ=255)
                                  tATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                                >  1:638653/1‑70 (MQ=255)
                                  tATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                                >  1:481086/1‑70 (MQ=255)
                                  tATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                                >  1:4919983/1‑70 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            <  1:3154455/70‑1 (MQ=255)
                                       ttCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                                <  1:846659/65‑1 (MQ=255)
                                       ttCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                                <  1:1294750/65‑1 (MQ=255)
                                       ttCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                                <  1:288402/65‑1 (MQ=255)
                                          cTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                           <  1:2623052/35‑1 (MQ=255)
                                          cTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                        <  1:5132803/70‑1 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                          >  1:3219456/1‑65 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                          >  1:3164492/1‑65 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                          >  1:1422354/1‑65 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                          >  1:2776512/1‑65 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                        >  1:4702062/1‑67 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                        >  1:5278675/1‑67 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                        >  1:5165576/1‑67 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                        >  1:2431151/1‑67 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                        >  1:3413013/1‑67 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                        >  1:4549073/1‑67 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                        >  1:4544844/1‑67 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                        >  1:4499564/1‑67 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                        >  1:4307884/1‑67 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                        >  1:3216816/1‑67 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                        >  1:650764/1‑67 (MQ=255)
                                                  ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTTCGccacca                                    <  1:2835048/50‑1 (MQ=255)
                                                  ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                    <  1:2781493/50‑1 (MQ=255)
                                                  ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                    <  1:5441896/50‑1 (MQ=255)
                                                  ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                    <  1:5397307/50‑1 (MQ=255)
                                                  ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                    <  1:732611/50‑1 (MQ=255)
                                                  ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                    <  1:2111967/50‑1 (MQ=255)
                                                  ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                    <  1:4874259/50‑1 (MQ=255)
                                                  ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                    <  1:4626054/50‑1 (MQ=255)
                                                  ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                    <  1:2370554/50‑1 (MQ=255)
                                                  ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                    <  1:26063/50‑1 (MQ=255)
                                                    atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                                  >  1:1824126/1‑50 (MQ=255)
                                                    atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                                  >  1:1415510/1‑50 (MQ=255)
                                                    atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                                  >  1:3823556/1‑50 (MQ=255)
                                                    atatCCTTGCCTGATTTCAGCCTGGCATTACCGTCAATTGCGCCACCaaa                                  >  1:5518190/1‑50 (MQ=255)
                                                    atatCCTTGCCTGATTTAAGCCTGGCATTACCTTCAATTGCGCCACCaaa                                  >  1:5203063/1‑50 (MQ=255)
                                                       tCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                    <  1:4964109/45‑1 (MQ=255)
                                                         cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCTCCACCAAAGTCATTAATGGCGTAAAAATGCACt         <  1:3000691/70‑1 (MQ=255)
                                                         cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAt                             >  1:2260040/1‑50 (MQ=255)
                                                         cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAt                             >  1:2729758/1‑50 (MQ=255)
                                                         cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACt         <  1:3865139/70‑1 (MQ=255)
                                                         cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACt         <  1:1048670/70‑1 (MQ=255)
                                                         cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACt         <  1:2070846/70‑1 (MQ=255)
                                                               tGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            <  1:797044/45‑1 (MQ=255)
                                                                gATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCa   >  1:3598042/1‑69 (MQ=255)
                                                                gATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg  >  1:5344147/1‑70 (MQ=255)
                                                                gATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg  >  1:4947000/1‑70 (MQ=255)
                                                                gATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg  >  1:3854326/1‑70 (MQ=255)
                                                                gATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg  >  1:1052947/1‑70 (MQ=255)
                                                                gATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGAACTTTTGCAg  >  1:2523629/1‑70 (MQ=255)
                                                                 aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg  <  1:3981676/69‑1 (MQ=255)
                                                                 aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg  <  1:3591907/69‑1 (MQ=255)
                                                                 aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg  <  1:1729508/69‑1 (MQ=255)
                                                                 aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg  <  1:1091532/69‑1 (MQ=255)
                                                                    |                                                                 
TGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAG  >  W3110S.gb/155395‑155528

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: