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

GATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAG  >  W3110S.gb/155399‑155528
                                                                |                                                                 
gatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAg                                                                >  1:3706315/1‑68 (MQ=255)
gatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAg                                                                <  1:1208120/68‑1 (MQ=255)
gatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAg                                                                <  1:3716062/68‑1 (MQ=255)
gatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                               >  1:3317053/1‑69 (MQ=255)
gatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                               >  1:1978346/1‑69 (MQ=255)
gatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                               >  1:1848063/1‑69 (MQ=255)
gatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                               >  1:469244/1‑69 (MQ=255)
gatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                               >  1:1689594/1‑69 (MQ=255)
   aaTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                             <  1:2760884/68‑1 (MQ=255)
   aaTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                             <  1:2307917/68‑1 (MQ=255)
   aaTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                             <  1:2765530/68‑1 (MQ=255)
      ttttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCa                                                         <  1:2224552/69‑1 (MQ=255)
       tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAt                                                        <  1:1413265/69‑1 (MQ=255)
              tatTCTATAGTCACGCTATGCTTCCTGCTGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                 >  1:3293211/1‑69 (MQ=255)
              tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                 >  1:2891163/1‑69 (MQ=255)
              tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                 >  1:2604716/1‑69 (MQ=255)
              tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                 <  1:1577351/69‑1 (MQ=255)
              tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                 >  1:28963/1‑69 (MQ=255)
              tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                 >  1:60654/1‑69 (MQ=255)
              tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                 <  1:625699/69‑1 (MQ=255)
              tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                 >  1:2109270/1‑69 (MQ=255)
                ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                 <  1:900304/51‑1 (MQ=255)
                ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                               >  1:247301/1‑69 (MQ=255)
                ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                               >  1:1624374/1‑69 (MQ=255)
                ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                               >  1:2781655/1‑69 (MQ=255)
                ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                               >  1:2952042/1‑69 (MQ=255)
                ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCa                                                <  1:2117484/68‑1 (MQ=255)
                ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCa                                                <  1:2312807/68‑1 (MQ=255)
                ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCATTGCCTGATTTCAGCCTGGCATTACCTTCaa                                               >  1:2824959/1‑69 (MQ=255)
                  cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                             >  1:1792785/1‑69 (MQ=255)
                  cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                             >  1:2375002/1‑69 (MQ=255)
                                 gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgc                                           <  1:1914385/56‑1 (MQ=255)
                                 gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                                  <  1:2166465/65‑1 (MQ=255)
                                 gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                                  <  1:2658620/65‑1 (MQ=255)
                                  cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                       >  1:358547/1‑43 (MQ=255)
                                  cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                       >  1:496820/1‑43 (MQ=255)
                                  cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                       >  1:1171716/1‑43 (MQ=255)
                                  cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                 >  1:1576618/1‑49 (MQ=255)
                                  cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                              >  1:3357021/1‑68 (MQ=255)
                                   ttCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                                  <  1:3377694/63‑1 (MQ=255)
                                    tCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            <  1:2404893/68‑1 (MQ=255)
                                    tCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            <  1:1429025/68‑1 (MQ=255)
                                    tCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            <  1:306959/68‑1 (MQ=255)
                                    tCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            <  1:3181849/68‑1 (MQ=255)
                                    tCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                           <  1:942887/69‑1 (MQ=255)
                                    tCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                           <  1:1140868/69‑1 (MQ=255)
                                    tCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                           <  1:277250/69‑1 (MQ=255)
                                        gCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            >  1:3371196/1‑64 (MQ=255)
                                                atatCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                 >  1:1511330/1‑35 (MQ=255)
                                                atatCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                 >  1:2744080/1‑35 (MQ=255)
                                                atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa               >  1:644287/1‑69 (MQ=255)
                                                atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa               >  1:1608222/1‑69 (MQ=255)
                                                atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa               >  1:3390662/1‑69 (MQ=255)
                                                atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa               >  1:364716/1‑69 (MQ=255)
                                                    ccTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            <  1:1052750/52‑1 (MQ=255)
                                                    ccTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            <  1:3301550/52‑1 (MQ=255)
                                                    ccTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            <  1:2911078/52‑1 (MQ=255)
                                                    ccTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaa                <  1:3476886/64‑1 (MQ=255)
                                                     cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            <  1:2875400/51‑1 (MQ=255)
                                                      ttGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaa                <  1:1812124/62‑1 (MQ=255)
                                                        gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACttt       >  1:2804159/1‑69 (MQ=255)
                                                             aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            >  1:2114711/1‑43 (MQ=255)
                                                             aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            >  1:386487/1‑43 (MQ=255)
                                                             aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg  >  1:3449844/1‑69 (MQ=255)
                                                             aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg  >  1:2044584/1‑69 (MQ=255)
                                                             aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg  >  1:147482/1‑69 (MQ=255)
                                                             aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg  >  1:3115071/1‑69 (MQ=255)
                                                                |                                                                 
GATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAG  >  W3110S.gb/155399‑155528

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