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

GATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTG  >  W3110S.gb/155396‑155525
                                                                   |                                                              
gatgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAg                                                             >  1:379840/1‑71 (MQ=255)
gatgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAg                                                             >  1:1515420/1‑71 (MQ=255)
     tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                        <  1:2436563/71‑1 (MQ=255)
      aaTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                        >  1:1419465/1‑70 (MQ=255)
       aTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGc                                                       <  1:1908180/70‑1 (MQ=255)
        tttttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAt                                                     >  1:1040153/1‑71 (MQ=255)
          tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                    <  1:2563787/70‑1 (MQ=255)
          tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                    <  1:2514592/70‑1 (MQ=255)
           ttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                    <  1:3759671/69‑1 (MQ=255)
           ttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTa                                                   >  1:3658464/1‑70 (MQ=255)
              gTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                    <  1:1658217/66‑1 (MQ=255)
               tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                              <  1:960271/71‑1 (MQ=255)
               tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                              <  1:581242/71‑1 (MQ=255)
               tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                              <  1:723191/71‑1 (MQ=255)
                 tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                            >  1:2279024/1‑71 (MQ=255)
                 tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                            >  1:2083190/1‑71 (MQ=255)
                 tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATACTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                            >  1:2021402/1‑71 (MQ=255)
                 tatTCTATAGTCACGCTATGCTTCCTCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt                                           <  1:3203706/71‑1 (MQ=255)
                   ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                            >  1:857208/1‑53 (MQ=255)
                   ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                          >  1:1875488/1‑71 (MQ=255)
                   ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                          >  1:3732519/1‑71 (MQ=255)
                   ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                          >  1:824709/1‑71 (MQ=255)
                   ttCTATAGTCACGCTATGATTCCTGCGGAATTATATCCTTGCCTGATTTCAGACTGGAAt                                                     >  1:1468639/1‑60 (MQ=255)
                              cGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                              <  1:3148805/56‑1 (MQ=255)
                                cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                             <  1:2831345/71‑1 (MQ=255)
                                cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                             <  1:1254874/71‑1 (MQ=255)
                                     cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                            >  1:4039706/1‑35 (MQ=255)
                                     cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt                                           >  1:3713280/1‑52 (MQ=255)
                                     cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt                                           >  1:3117059/1‑52 (MQ=255)
                                     cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                         <  1:123235/70‑1 (MQ=255)
                                     cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                         <  1:3279009/70‑1 (MQ=255)
                                     cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                         <  1:590205/70‑1 (MQ=255)
                                     cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                        <  1:4014700/71‑1 (MQ=255)
                                        ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGcca                                    >  1:1720210/1‑56 (MQ=255)
                                        ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                     >  1:1336833/1‑71 (MQ=255)
                                          tGCGGAATTATATCCTTGCCTGATTTCAGCCTGGc                                                       <  1:3576188/35‑1 (MQ=255)
                                          tGCGGAATTATATCCTTGCCTGATTTCAGCCTGGc                                                       <  1:4063322/35‑1 (MQ=255)
                                            cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                            >  1:2842344/1‑44 (MQ=255)
                                            cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                            >  1:207281/1‑44 (MQ=255)
                                               aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                         <  1:728036/60‑1 (MQ=255)
                                               aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa              <  1:2951894/71‑1 (MQ=255)
                                               aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa              <  1:2899882/71‑1 (MQ=255)
                                               aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa              <  1:1391509/71‑1 (MQ=255)
                                                 ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                       >  1:1256159/1‑60 (MQ=255)
                                                 ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                       >  1:704474/1‑60 (MQ=255)
                                                  tataTCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                              <  1:388866/36‑1 (MQ=255)
                                                  tataTCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                              <  1:539921/36‑1 (MQ=255)
                                                  tataTCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                              <  1:3577732/36‑1 (MQ=255)
                                                  tataTCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                              <  1:1427559/36‑1 (MQ=255)
                                                   atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg          >  1:2918837/1‑71 (MQ=255)
                                                   atatCCTAGCCTGATTTCAGCCTGCCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg          >  1:727370/1‑71 (MQ=255)
                                                      tCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                       >  1:2557593/1‑55 (MQ=255)
                                                      tCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                       >  1:662140/1‑55 (MQ=255)
                                                      tCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCAc       >  1:3296875/1‑71 (MQ=255)
                                                        cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtt     >  1:3073793/1‑71 (MQ=255)
                                                        cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtt     <  1:1412655/71‑1 (MQ=255)
                                                        cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtt     <  1:1997999/71‑1 (MQ=255)
                                                         ttGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                         <  1:11131/50‑1 (MQ=255)
                                                          tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtttt   >  1:1997487/1‑71 (MQ=255)
                                                          tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtttt   >  1:1122031/1‑71 (MQ=255)
                                                          tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtttt   >  1:3288777/1‑71 (MQ=255)
                                                            ccTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg  <  1:2549630/70‑1 (MQ=255)
                                                             cTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                        <  1:2383443/47‑1 (MQ=255)
                                                             cTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                        <  1:2378774/47‑1 (MQ=255)
                                                             cTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                        <  1:2828851/47‑1 (MQ=255)
                                                             cTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                        <  1:2927942/47‑1 (MQ=255)
                                                             cTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                        <  1:3031154/47‑1 (MQ=255)
                                                             cTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                        <  1:880115/47‑1 (MQ=255)
                                                                   |                                                              
GATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTG  >  W3110S.gb/155396‑155525

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: