Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	155,463	A→T	*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.gb	155,463	0	A	T	100.0%	299.5 / NA	99	*247K (TAA→AAA)	ecpD	predicted periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (41/58);  total (41/58)

TGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTT  >  W3110S.gb/155395‑155522                                                                     |                                                            tgatgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                           <  1:3524807/71‑1 (MQ=255) tgatgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                           <  1:3605350/71‑1 (MQ=255)   atgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAg                                                          <  1:3436092/70‑1 (MQ=255)   atgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                         <  1:4187694/71‑1 (MQ=255)      atAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTg                                                      >  1:4402246/1‑71 (MQ=255)       tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                     <  1:2145031/71‑1 (MQ=255)       tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                     <  1:1859383/71‑1 (MQ=255)       tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                     <  1:2034955/71‑1 (MQ=255)       tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                     <  1:3011025/71‑1 (MQ=255)       tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                     <  1:1096910/71‑1 (MQ=255)       tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                     <  1:4549281/71‑1 (MQ=255)         aTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGc                                                    <  1:2790962/70‑1 (MQ=255)          tttttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCATGATTTCAGCCTgg                                                     <  1:472500/68‑1 (MQ=255)            tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                 <  1:2055712/70‑1 (MQ=255)            tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                 <  1:4001246/70‑1 (MQ=255)            tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                 <  1:21502/70‑1 (MQ=255)            tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                 <  1:2876597/70‑1 (MQ=255)            tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                 <  1:1540270/70‑1 (MQ=255)            tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                 <  1:1200488/70‑1 (MQ=255)            tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                 <  1:3854131/70‑1 (MQ=255)              tAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAc                                               <  1:4370229/70‑1 (MQ=255)                 tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                           <  1:4441725/71‑1 (MQ=255)                 tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                           <  1:2441412/71‑1 (MQ=255)                 tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                           <  1:294657/71‑1 (MQ=255)                 tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                           <  1:812123/71‑1 (MQ=255)                  atatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                           <  1:1473828/70‑1 (MQ=255)                   tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                           <  1:2064464/53‑1 (MQ=255)                   tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                           <  1:3463186/53‑1 (MQ=255)                   tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                           <  1:3901474/69‑1 (MQ=255)                   tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                         >  1:2243204/1‑71 (MQ=255)                   tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCa                                          <  1:2130584/70‑1 (MQ=255)                   tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCa                                          <  1:2599738/70‑1 (MQ=255)                   tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCa                                          <  1:2889824/70‑1 (MQ=255)                   tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCa                                          <  1:3548042/70‑1 (MQ=255)                    atTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                 <  1:889105/62‑1 (MQ=255)                     ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                 >  1:2071574/1‑61 (MQ=255)                     ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                       >  1:2054793/1‑71 (MQ=255)                     ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                       >  1:4060136/1‑71 (MQ=255)                     ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                       >  1:102090/1‑71 (MQ=255)                     ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                       >  1:2762452/1‑71 (MQ=255)                     ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                       >  1:4147468/1‑71 (MQ=255)                     ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                       >  1:3819844/1‑71 (MQ=255)                     ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                       >  1:2880111/1‑71 (MQ=255)                     ttCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt                                        <  1:2284497/70‑1 (MQ=255)                      tCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTg                                      >  1:3141880/1‑71 (MQ=255)                      tCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTg                                      >  1:1316507/1‑71 (MQ=255)                        tataGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                 <  1:3390022/58‑1 (MQ=255)                          taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                           <  1:1042865/62‑1 (MQ=255)                          taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                   >  1:2071019/1‑70 (MQ=255)                          taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACATTCAATTgcgc                                   >  1:1016089/1‑70 (MQ=255)                            gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccac                                >  1:4298053/1‑71 (MQ=255)                            gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccac                                >  1:4102919/1‑71 (MQ=255)                             tCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                           <  1:3453934/43‑1 (MQ=255)                             tCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCa                                                           <  1:4659266/43‑1 (MQ=255)                               aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                              >  1:4618145/1‑70 (MQ=255)                               aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacca                              >  1:4628867/1‑70 (MQ=255)                               aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGcaacca                              >  1:57240/1‑70 (MQ=255)                                  cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                          <  1:2000995/71‑1 (MQ=255)                                    aTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                        >  1:3950096/1‑71 (MQ=255)                                     tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAt                                                  <  1:142959/44‑1 (MQ=255)                                     tGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAt                                                  <  1:1120435/44‑1 (MQ=255)                                     tGCTTCCTCCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAt                                                  <  1:282521/44‑1 (MQ=38)                                       cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAc                                               >  1:1717088/1‑45 (MQ=255)                                       cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAc                                               >  1:2729018/1‑45 (MQ=255)                                       cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGccacc                               >  1:3920820/1‑61 (MQ=255)                                       cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                      <  1:714340/70‑1 (MQ=255)                                       cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                      <  1:1625528/70‑1 (MQ=255)                                       cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                      <  1:183880/70‑1 (MQ=255)                                       cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                      <  1:454411/70‑1 (MQ=255)                                       cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                      >  1:4071623/1‑70 (MQ=255)                                       cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                      >  1:3508421/1‑70 (MQ=255)                                       cTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAt                       >  1:4541251/1‑69 (MQ=255)                                         tCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAAt                   >  1:3199410/1‑71 (MQ=255)                                         tCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAAt                   >  1:2169446/1‑71 (MQ=255)                                           cTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAcc                                              <  1:4241280/42‑1 (MQ=255)                                             gCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt                                        >  1:4393666/1‑46 (MQ=255)                                              cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAcc                                              >  1:2026200/1‑39 (MQ=255)                                              cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAcc                                              >  1:4432844/1‑39 (MQ=255)                                                 aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                      <  1:4160307/60‑1 (MQ=255)                                                 aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                      <  1:436306/60‑1 (MQ=255)                                                 aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa           <  1:1017987/71‑1 (MQ=255)                                                 aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa           <  1:3875100/71‑1 (MQ=255)                                                 aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa           <  1:1438087/71‑1 (MQ=255)                                                 aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa           <  1:2275137/71‑1 (MQ=255)                                                 aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa           <  1:2942476/71‑1 (MQ=255)                                                  aTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa           >  1:35420/1‑70 (MQ=255)                                                   ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                          <  1:2455122/54‑1 (MQ=255)                                                   ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                          <  1:4378504/54‑1 (MQ=255)                                                     atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                            >  1:2026976/1‑50 (MQ=255)                                                     atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                            >  1:3683414/1‑50 (MQ=255)                                                     atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa                        >  1:495987/1‑54 (MQ=255)                                                     atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg       >  1:502406/1‑71 (MQ=255)                                                     atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg       >  1:3745169/1‑71 (MQ=255)                                                      tatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg       <  1:3584097/70‑1 (MQ=255)                                                       atCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCa     >  1:3579001/1‑71 (MQ=255)                                                       atCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCa     >  1:4675748/1‑71 (MQ=255)                                                        tCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCAc    <  1:2724399/71‑1 (MQ=255)                                                          cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtt  >  1:2551542/1‑71 (MQ=255)                                                          cTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtt  >  1:1405377/1‑71 (MQ=255)                                                                     |                                                            TGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTT  >  W3110S.gb/155395‑155522

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: ‑