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

TGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAA  >  W3110S.gb/155398‑155530
                                                                 |                                                                   
tgatAATTTTTAGTATATTCTACAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGcc                                                                >  1:2413171/1‑71 (MQ=255)
   tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTg                                                              >  1:1235777/1‑70 (MQ=255)
    aaTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAg                                                                  >  1:2489735/1‑65 (MQ=255)
    aaTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                             >  1:2397559/1‑70 (MQ=255)
    aaTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                             >  1:1380622/1‑70 (MQ=255)
     aTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAg                                                                  >  1:4631114/1‑64 (MQ=255)
     aTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAg                                                                  >  1:800503/1‑64 (MQ=255)
        tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                         <  1:2908947/70‑1 (MQ=255)
         ttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                         <  1:1783575/69‑1 (MQ=255)
          tAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                         <  1:1867638/68‑1 (MQ=255)
           aGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAcc                                                      <  1:263936/70‑1 (MQ=255)
             tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                             <  1:2488969/61‑1 (MQ=255)
               tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCa                                                  <  1:4572253/70‑1 (MQ=255)
                atTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt                                                >  1:3833913/1‑71 (MQ=255)
                 ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                               <  1:2264903/71‑1 (MQ=255)
                  tCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTg                                              >  1:4374226/1‑71 (MQ=255)
                   cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                             >  1:3245575/1‑55 (MQ=255)
                   cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAc                                                       >  1:855663/1‑61 (MQ=255)
                   cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAc                                                       >  1:1083325/1‑61 (MQ=255)
                   cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                                 >  1:2656473/1‑67 (MQ=255)
                   cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgc                                             >  1:101798/1‑71 (MQ=255)
                   cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgc                                             >  1:930222/1‑71 (MQ=255)
                     ataGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                             <  1:3667962/53‑1 (MQ=255)
                      taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                           <  1:3313278/70‑1 (MQ=255)
                       aGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                           <  1:2277692/69‑1 (MQ=255)
                        gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                         <  1:1688259/54‑1 (MQ=255)
                        gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                         <  1:2213932/54‑1 (MQ=255)
                        gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt                                                <  1:2448016/63‑1 (MQ=255)
                         tCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacc                                       <  1:2122587/71‑1 (MQ=255)
                         tCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacc                                       <  1:3670010/71‑1 (MQ=255)
                          cACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTAAATTGCGccac                                        >  1:828141/1‑69 (MQ=255)
                           aCGCTATGCTTCCTGCGTAATTATATCCTTGCCTGATTTCAGCCTGGCAt                                                          >  1:2331375/1‑50 (MQ=39)
                           aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCa                                                           >  1:1724001/1‑49 (MQ=255)
                           aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                      <  1:3912939/70‑1 (MQ=255)
                           aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                      >  1:3367704/1‑70 (MQ=255)
                           aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                      >  1:2614887/1‑70 (MQ=255)
                           aCGCTATGCTTCCTGCGGAATTATATCATTGCCTGATTTCAGCCTGGCAt                                                          >  1:526191/1‑50 (MQ=255)
                             gCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAg                                   >  1:805902/1‑71 (MQ=255)
                              cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                                  >  1:1437087/1‑71 (MQ=255)
                              cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                                  <  1:3353167/71‑1 (MQ=255)
                                aTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTg                                                              <  1:4421764/41‑1 (MQ=255)
                                 tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAttg                                                        <  1:278719/46‑2 (MQ=255)
                                  gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                                 >  1:3818935/1‑36 (MQ=255)
                                  gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAt                               >  1:3571513/1‑70 (MQ=255)
                                   cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                                 >  1:4332547/1‑35 (MQ=255)
                                   cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                                 >  1:4337984/1‑35 (MQ=255)
                                   cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                                 >  1:2764480/1‑35 (MQ=255)
                                   cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                              <  1:1606311/70‑1 (MQ=255)
                                   cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                              <  1:1903243/70‑1 (MQ=255)
                                   cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCCGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                              <  1:3255844/70‑1 (MQ=255)
                                    ttCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                              <  1:3167338/69‑1 (MQ=255)
                                     tCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                             <  1:3031076/37‑1 (MQ=255)
                                     tCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAAt                           >  1:4152128/1‑71 (MQ=255)
                                     tCCTGCGGAATTATATCCTTCCTGATTTCAGCCTgg                                                             <  1:4223481/36‑1 (MQ=37)
                                       cTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                          <  1:334598/70‑1 (MQ=255)
                                         gCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                                 >  1:3016281/1‑45 (MQ=255)
                                         gCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCg                       >  1:168670/1‑71 (MQ=255)
                                                 atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                                 >  1:3198616/1‑37 (MQ=255)
                                                 atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                                >  1:943627/1‑54 (MQ=255)
                                                 atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg               >  1:4199102/1‑71 (MQ=255)
                                                 atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg               >  1:712525/1‑71 (MQ=255)
                                                  tatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg               <  1:2162736/70‑1 (MQ=255)
                                                   atCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTc                                 <  1:1174575/51‑1 (MQ=255)
                                                   atCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                             <  1:755352/55‑1 (MQ=255)
                                                   atCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                             <  1:750106/55‑1 (MQ=255)
                                                   atCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCa             >  1:979731/1‑71 (MQ=255)
                                                    tCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGTGCCACCAAAGTCAt                               <  1:2579843/52‑1 (MQ=255)
                                                    tCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCAc            >  1:73553/1‑71 (MQ=255)
                                                     ccTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGc                        >  1:4356164/1‑58 (MQ=255)
                                                             gATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                                  <  1:1447798/40‑1 (MQ=255)
                                                             gATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg    >  1:3732125/1‑70 (MQ=255)
                                                              aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                            >  1:2476288/1‑45 (MQ=255)
                                                              aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtt          >  1:2132699/1‑63 (MQ=255)
                                                              aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg    <  1:214006/69‑1 (MQ=255)
                                                              aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg    <  1:1098264/69‑1 (MQ=255)
                                                              aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGaa  <  1:1763667/71‑1 (MQ=255)
                                                                 |                                                                   
TGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAA  >  W3110S.gb/155398‑155530

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: