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

GATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAA  >  W3110S.gb/155399‑155530
                                                                |                                                                   
gatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                               >  1:550574/1‑71 (MQ=255)
    aTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGc                                                            <  1:904011/70‑1 (MQ=255)
    aTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCa                                                           >  1:851457/1‑71 (MQ=255)
       tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTa                                                        >  1:10532/1‑71 (MQ=255)
                ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                               >  1:1532775/1‑71 (MQ=255)
                  cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCGGCCTGGCATTACCTTCAATTgc                                             >  1:1295090/1‑71 (MQ=255)
                           cGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                           <  1:1239675/64‑1 (MQ=255)
                            gCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAg                                   <  1:865763/71‑1 (MQ=255)
                             cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                                  >  1:320637/1‑71 (MQ=255)
                                tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                                <  1:1755169/70‑1 (MQ=255)
                                 gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGTCATTACCTTCaa                                                 >  1:1451522/1‑52 (MQ=39)
                                  cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                                 >  1:564344/1‑51 (MQ=255)
                                  cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                                 >  1:1507939/1‑51 (MQ=255)
                                  cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt                                                <  1:256822/52‑1 (MQ=255)
                                  cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt                                                <  1:61760/52‑1 (MQ=255)
                                  cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                              >  1:935442/1‑70 (MQ=255)
                                      cTGCGGAATTATATCTTTGCCTGATTTCAGCCTGGCAt                                                          >  1:460499/1‑38 (MQ=37)
                                      cTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                      <  1:1284982/58‑1 (MQ=255)
                                      cTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATgg                         >  1:282967/1‑71 (MQ=255)
                                                atatCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                   >  1:141381/1‑35 (MQ=255)
                                                 tatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg               >  1:1385006/1‑70 (MQ=255)
                                                    ccTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACt           >  1:286753/1‑71 (MQ=255)
                                                         ccTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtt          >  1:1281575/1‑67 (MQ=255)
                                                             aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg    <  1:443083/69‑1 (MQ=255)
                                                             aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg    <  1:675749/69‑1 (MQ=255)
                                                             aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGaa  >  1:1654303/1‑71 (MQ=255)
                                                             aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGaa  >  1:1254219/1‑71 (MQ=255)
                                                                |                                                                   
GATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAA  >  W3110S.gb/155399‑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: