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

GATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAA  >  W3110S.gb/155396‑155530
                                                                   |                                                                   
gatgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAg                                                                  >  1:2001070/1‑71 (MQ=255)
gatgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAg                                                                  >  1:1932264/1‑71 (MQ=255)
  tgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGcc                                                                <  1:1442288/71‑1 (MQ=255)
       aTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCa                                                           >  1:938110/1‑71 (MQ=255)
                atatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGc                                                            <  1:231630/61‑1 (MQ=255)
                 tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                   <  1:701021/53‑1 (MQ=255)
                 tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                                 >  1:1997191/1‑71 (MQ=255)
                 tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                                 <  1:731367/71‑1 (MQ=255)
                   ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                               <  1:1309162/71‑1 (MQ=255)
                           tCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                             <  1:1732424/49‑1 (MQ=255)
                             aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaa                                     <  1:1907929/71‑1 (MQ=255)
                                   tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAt                                                          <  1:295570/44‑1 (MQ=255)
                                    gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGc                                                            <  1:708551/41‑1 (MQ=255)
                                    gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAt                                                          <  1:1571419/43‑1 (MQ=255)
                                     cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTa                                                        >  1:2070880/1‑44 (MQ=255)
                                       tCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAAt                           >  1:1240499/1‑71 (MQ=255)
                                        ccTGCGGAATTATATCCTTGCCTTATTTCAGCCTGGCATTAc                                                       >  1:842595/1‑42 (MQ=255)
                                        ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAAt                           <  1:443647/70‑1 (MQ=255)
                                            cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCg                       >  1:936642/1‑70 (MQ=255)
                                             ggAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                                >  1:1444646/1‑60 (MQ=255)
                                                aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaa                  <  1:308144/71‑1 (MQ=255)
                                                   atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg               >  1:128575/1‑71 (MQ=255)
                                                    tatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg               <  1:1997514/70‑1 (MQ=255)
                                                       ccTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACAAAAGTCATTAATgg                         >  1:591681/1‑57 (MQ=255)
                                                           gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACttt         <  1:1596287/69‑1 (MQ=255)
                                                           gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg       >  1:1784957/1‑71 (MQ=255)
                                                            ccTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAt                               <  1:1418638/46‑1 (MQ=255)
                                                             cTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg               <  1:373359/61‑1 (MQ=255)
                                                               gATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg    >  1:1408443/1‑70 (MQ=255)
                                                               gATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg    >  1:370654/1‑70 (MQ=255)
                                                                aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                                <  1:1967245/41‑1 (MQ=255)
                                                                aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                            >  1:802111/1‑45 (MQ=255)
                                                                aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCg                       >  1:1361176/1‑50 (MQ=255)
                                                                aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg    <  1:386/69‑1 (MQ=255)
                                                                aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGaa  >  1:2012262/1‑71 (MQ=255)
                                                                aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGaa  <  1:1088920/71‑1 (MQ=255)
                                                                   |                                                                   
GATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAA  >  W3110S.gb/155396‑155530

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