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

TGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAG  >  W3110S.gb/155395‑155528
                                                                    |                                                                 
tgatgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                                 <  1:2918390/71‑1 (MQ=255)
 gatgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAg                                                                >  1:1274116/1‑71 (MQ=255)
  atgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                               >  1:2577082/1‑71 (MQ=255)
   tgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGcc                                                              <  1:1565091/71‑1 (MQ=255)
     atAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                             <  1:1966310/70‑1 (MQ=255)
     atAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                             <  1:2624123/70‑1 (MQ=255)
      tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                           <  1:530569/71‑1 (MQ=255)
       aaTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                           >  1:1050359/1‑70 (MQ=255)
        aTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGc                                                          <  1:1564499/70‑1 (MQ=255)
         tttttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAt                                                        <  1:3025655/71‑1 (MQ=255)
              aGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAcc                                                    <  1:3119881/70‑1 (MQ=255)
               gTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCt                                                   <  1:1883124/70‑1 (MQ=255)
                  tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                               >  1:138091/1‑71 (MQ=255)
                  tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                               >  1:2588547/1‑71 (MQ=255)
                  tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCa                                                <  1:1257826/70‑1 (MQ=255)
                  tatTCTATAGTCACGCTATGCTTCATGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                               >  1:1198873/1‑71 (MQ=255)
                   atTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCa                                                         <  1:2106277/60‑1 (MQ=255)
                    ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                 >  1:546954/1‑67 (MQ=255)
                    ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                                 >  1:543049/1‑67 (MQ=255)
                     tCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTa                                                      >  1:2000382/1‑61 (MQ=255)
                     tCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAc                                                     >  1:2681208/1‑62 (MQ=255)
                        ataGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCa                                                <  1:1412099/64‑1 (MQ=255)
                        ataGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                         >  1:2632306/1‑71 (MQ=255)
                           gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCa                                                         >  1:853065/1‑52 (MQ=255)
                            tCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacc                                     <  1:1457223/71‑1 (MQ=255)
                            tCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccac                                      >  1:2663123/1‑70 (MQ=255)
                              aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                               >  1:324385/1‑43 (MQ=255)
                              aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                    <  1:3143706/70‑1 (MQ=255)
                              aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                    <  1:2976221/70‑1 (MQ=255)
                               cGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                    >  1:3142720/1‑69 (MQ=255)
                               cGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                    >  1:237692/1‑69 (MQ=255)
                                gCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAg                                 >  1:1738670/1‑71 (MQ=255)
                                    tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTc                               <  1:1725560/69‑1 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            >  1:2451661/1‑70 (MQ=255)
                                      cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            <  1:438031/70‑1 (MQ=255)
                                       ttCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                              <  1:1827688/67‑1 (MQ=255)
                                        tCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATAAAt                         >  1:1844608/1‑71 (MQ=255)
                                         ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                          >  1:958666/1‑69 (MQ=255)
                                         ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTaa                          <  1:2959712/69‑1 (MQ=255)
                                         ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                        >  1:1289464/1‑71 (MQ=255)
                                         ccTGCGGAATTATATCCTTGCCAGATTTCAGCCTGGCATTACCTTCAATTgc                                           >  1:2736834/1‑52 (MQ=39)
                                          cTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAc                                                     >  1:1695785/1‑41 (MQ=255)
                                          cTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccac                                      >  1:81341/1‑56 (MQ=255)
                                          cTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccac                                      >  1:454718/1‑56 (MQ=255)
                                             cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCg                     >  1:1406397/1‑70 (MQ=255)
                                              ggAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGt                    >  1:1245374/1‑70 (MQ=255)
                                               gAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaa                  <  1:646445/71‑1 (MQ=255)
                                                aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgc                                           <  1:1608729/45‑1 (MQ=255)
                                                 aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                                >  1:2224722/1‑55 (MQ=255)
                                                 aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                                >  1:2186662/1‑55 (MQ=255)
                                                 aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaa                <  1:2600872/71‑1 (MQ=255)
                                                  ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                            <  1:1681597/58‑1 (MQ=255)
                                                  ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaa                  <  1:1282274/68‑1 (MQ=255)
                                                   tataTCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAAt              >  1:337201/1‑71 (MQ=255)
                                                   tataTCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAAt              >  1:1164166/1‑71 (MQ=255)
                                                       tCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                    <  1:2040157/45‑1 (MQ=255)
                                                            gCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                         <  1:128634/35‑1 (MQ=255)
                                                            gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                                    <  1:2220613/40‑1 (MQ=255)
                                                            gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtttt      <  1:1061172/70‑1 (MQ=255)
                                                            gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtttt      <  1:3149981/70‑1 (MQ=255)
                                                            gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtttt      <  1:1141992/70‑1 (MQ=255)
                                                               tGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg     >  1:2290450/1‑68 (MQ=255)
                                                                gATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg  >  1:2721250/1‑70 (MQ=255)
                                                                 aTTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg             >  1:493528/1‑58 (MQ=255)
                                                                    |                                                                 
TGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAG  >  W3110S.gb/155395‑155528

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: