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% 249.4 / NA 81*247K (TAA→AAA) ecpDpredicted periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (43/38);  total (43/38)
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:2903838/1‑71 (MQ=255)
gatgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAg                                                                  >  1:3894850/1‑71 (MQ=255)
   gatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                               >  1:3746812/1‑71 (MQ=255)
    atAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                               <  1:1726362/70‑1 (MQ=255)
    atAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTg                                                              >  1:4436296/1‑71 (MQ=255)
     tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                             <  1:4063191/71‑1 (MQ=255)
     tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                             >  1:2553194/1‑71 (MQ=255)
     tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                             <  1:883501/71‑1 (MQ=255)
      aaTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                   >  1:3088639/1‑64 (MQ=255)
      aaTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                   >  1:4656298/1‑64 (MQ=255)
       aTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                                 >  1:4022287/1‑65 (MQ=255)
        tttttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAt                                                          <  1:3458213/71‑1 (MQ=255)
          tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                         <  1:4159906/70‑1 (MQ=255)
               tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                   <  1:53672/71‑1 (MQ=255)
               tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                   <  1:4511473/71‑1 (MQ=255)
               tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                   <  1:1743787/71‑1 (MQ=255)
               tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                   <  1:2561604/71‑1 (MQ=255)
               tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                   <  1:4145933/71‑1 (MQ=255)
                 tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                                 >  1:3252829/1‑71 (MQ=255)
                 tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                                 <  1:3094850/71‑1 (MQ=255)
                 tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                                 >  1:530037/1‑71 (MQ=255)
                 tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCAGATTTCAGCCTGGCATTACCTTCaa                                                 >  1:172111/1‑71 (MQ=255)
                   ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCa                                                           >  1:833204/1‑59 (MQ=255)
                   ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                               <  1:701570/71‑1 (MQ=255)
                   ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                               >  1:1318710/1‑71 (MQ=255)
                   ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt                                                <  1:780869/70‑1 (MQ=255)
                     cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgc                                             >  1:2680007/1‑71 (MQ=255)
                     cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTg                                              <  1:904142/70‑1 (MQ=255)
                        taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                   <  1:1725619/46‑1 (MQ=255)
                        taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                           <  1:2704052/70‑1 (MQ=255)
                          gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGcca                                         <  1:2651813/70‑1 (MQ=255)
                          gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGcca                                         <  1:538281/70‑1 (MQ=255)
                           tCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                             <  1:2945045/49‑1 (MQ=255)
                           tCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAc                                                       >  1:3263411/1‑55 (MQ=255)
                           tCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAc                                                       >  1:4645071/1‑55 (MQ=255)
                            cACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacc                                       >  1:2799179/1‑70 (MQ=255)
                            cACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacc                                       >  1:1456970/1‑70 (MQ=255)
                             aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCa                                                           >  1:1257268/1‑49 (MQ=255)
                             aCGCTATGCTTCATGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                      >  1:3380494/1‑70 (MQ=255)
                                  aTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                                >  1:3603302/1‑71 (MQ=255)
                                   tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTg                                              <  1:3849528/56‑1 (MQ=255)
                                   tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                                <  1:1252181/70‑1 (MQ=255)
                                   tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                                <  1:1674317/70‑1 (MQ=255)
                                    gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                   >  1:3537042/1‑50 (MQ=255)
                                     cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                                 >  1:2163055/1‑35 (MQ=255)
                                     cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAATTCAtt                              <  1:431014/70‑1 (MQ=255)
                                      ttCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                             >  1:2462438/1‑70 (MQ=255)
                                        ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAAt                           <  1:1833758/70‑1 (MQ=255)
                                        ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAAt                           <  1:1828250/70‑1 (MQ=255)
                                           gCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCg                       >  1:2468364/1‑71 (MQ=255)
                                            cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                           >  1:2418085/1‑50 (MQ=255)
                                            cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCg                       >  1:2847484/1‑70 (MQ=255)
                                            cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAGTa                             >  1:70870/1‑64 (MQ=255)
                                              gAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                                  <  1:259228/57‑1 (MQ=255)
                                                aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaa                  <  1:4337485/71‑1 (MQ=255)
                                                aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaa                  <  1:3399263/71‑1 (MQ=255)
                                                aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaa                  <  1:3960124/71‑1 (MQ=255)
                                                 ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCa                                                  <  1:2844469/38‑1 (MQ=255)
                                                 ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                            >  1:2854635/1‑60 (MQ=255)
                                                 ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaa                  >  1:3903215/1‑70 (MQ=255)
                                                   atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGcca                                         <  1:3088932/45‑1 (MQ=255)
                                                   atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTc                                 >  1:373779/1‑53 (MQ=255)
                                                   atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAAt                >  1:4585692/1‑70 (MQ=255)
                                                   atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg               >  1:3424655/1‑71 (MQ=255)
                                                   atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCCCCACCaaa                                    >  1:538915/1‑50 (MQ=255)
                                                    tatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg               >  1:661814/1‑70 (MQ=255)
                                                     atCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt                                                >  1:1504643/1‑36 (MQ=255)
                                                     atCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                           >  1:263243/1‑41 (MQ=255)
                                                     atCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                           >  1:832622/1‑41 (MQ=255)
                                                        cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACt           <  1:429345/70‑1 (MQ=255)
                                                          tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtttt        >  1:4698981/1‑71 (MQ=255)
                                                             cTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtttt        >  1:1761655/1‑68 (MQ=255)
                                                              tGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg               <  1:2627030/60‑1 (MQ=255)
                                                                aTTTCAGCCTGGCATTACCTTCTATTGCGCCACCAAAGTCa                                >  1:605498/1‑41 (MQ=38)
                                                                aTTTCAGCCTGGCATTACCTTCAATTGCGccacca                                      <  1:4212700/35‑1 (MQ=255)
                                                                aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                                >  1:3449427/1‑41 (MQ=255)
                                                                aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                   <  1:719463/54‑1 (MQ=255)
                                                                aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg    <  1:3059820/69‑1 (MQ=255)
                                                                aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGaa  <  1:615990/71‑1 (MQ=255)
                                                                aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGaa  <  1:4439803/71‑1 (MQ=255)
                                                                aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGaa  >  1:950440/1‑71 (MQ=255)
                                                                   |                                                                   
GATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAA  >  W3110S.gb/155396‑155530

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