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

ATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAA  >  W3110S.gb/155397‑155530
                                                                  |                                                                   
atgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                                 <  1:1497365/71‑1 (MQ=255)
     aaTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                             >  1:1294676/1‑70 (MQ=255)
         tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                         <  1:160001/70‑1 (MQ=255)
                tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                                 >  1:269405/1‑71 (MQ=255)
                tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                                 >  1:1249903/1‑71 (MQ=255)
                tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                                 >  1:1741610/1‑71 (MQ=255)
                  ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                               >  1:565351/1‑71 (MQ=255)
                  ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                               >  1:23913/1‑71 (MQ=255)
                  ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                               >  1:662803/1‑71 (MQ=255)
                  ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                               >  1:884250/1‑71 (MQ=255)
                    cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgc                                             >  1:105584/1‑71 (MQ=255)
                     tataGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcg                                            >  1:907893/1‑71 (MQ=255)
                       taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                               >  1:1172970/1‑50 (MQ=255)
                       taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                               >  1:81874/1‑50 (MQ=255)
                       taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                               >  1:1045328/1‑50 (MQ=255)
                             cGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                                    <  1:1251147/71‑1 (MQ=255)
                                 aTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                                >  1:1655225/1‑71 (MQ=255)
                                 aTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                                >  1:434082/1‑71 (MQ=255)
                                  tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                   <  1:1032193/51‑1 (MQ=255)
                                  tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                   <  1:1382333/51‑1 (MQ=255)
                                    cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                                 >  1:1026909/1‑35 (MQ=255)
                                    cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                                 >  1:1132632/1‑35 (MQ=255)
                                    cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                                 >  1:1202683/1‑35 (MQ=255)
                                    cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                                 >  1:260693/1‑35 (MQ=255)
                                    cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                              <  1:1220942/70‑1 (MQ=255)
                                    cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                              <  1:642107/70‑1 (MQ=255)
                                    cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                              <  1:1837689/70‑1 (MQ=255)
                                    cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                              <  1:1732807/70‑1 (MQ=255)
                                       ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                          >  1:1598281/1‑71 (MQ=255)
                                        cTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATgg                         <  1:597848/71‑1 (MQ=255)
                                        cTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATgg                         <  1:782091/71‑1 (MQ=255)
                                              aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                   <  1:1643337/71‑1 (MQ=255)
                                              aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                   <  1:347064/71‑1 (MQ=255)
                                              aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                   <  1:278372/71‑1 (MQ=255)
                                              aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                   <  1:689147/71‑1 (MQ=255)
                                              aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                   <  1:1768036/71‑1 (MQ=255)
                                              aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                   <  1:1057904/71‑1 (MQ=255)
                                              aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                   <  1:1683553/71‑1 (MQ=255)
                                               aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                   <  1:1094254/70‑1 (MQ=255)
                                                ttATATCCTTGCCTGATTTCAGCCTGGCCTTACCTTc                                                   <  1:1672376/37‑1 (MQ=37)
                                                ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa                 >  1:1701481/1‑71 (MQ=255)
                                                ttATATCCTTCCTGATTTCAGCCTGGCATTACCTTc                                                   <  1:418825/36‑1 (MQ=37)
                                                    atCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCa             >  1:7909/1‑71 (MQ=255)
                                                    atCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCa             >  1:826460/1‑71 (MQ=255)
                                                               aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGaa  >  1:568352/1‑71 (MQ=255)
                                                                  |                                                                   
ATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAA  >  W3110S.gb/155397‑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: