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%	200.1 / NA	64	*247K (TAA→AAA)	ecpD	predicted periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (27/37);  total (27/37)

ATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTG  >  W3110S.gb/155397‑155525                                                                   |                                                               atgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                            >  1:885824/1‑71 (MQ=255) atgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                            >  1:1275866/1‑71 (MQ=255)     tAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                        >  1:28041/1‑71 (MQ=255)          tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                    <  1:1935625/70‑1 (MQ=255)          tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                    <  1:1990981/70‑1 (MQ=255)           ttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                    <  1:3400844/69‑1 (MQ=255)           ttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTa                                                   <  1:228969/70‑1 (MQ=255)           ttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTa                                                   <  1:2347185/70‑1 (MQ=255)           ttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTa                                                   <  1:281602/70‑1 (MQ=255)           ttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTa                                                   <  1:3230935/70‑1 (MQ=255)           ttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTa                                                   <  1:1448404/70‑1 (MQ=255)           ttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTa                                                   <  1:56042/70‑1 (MQ=255)            tAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTa                                                   <  1:2340136/69‑1 (MQ=255)            tAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTa                                                   <  1:1964475/69‑1 (MQ=255)               tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGcc                                                           <  1:565087/58‑1 (MQ=255)               tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGcc                                                           <  1:2803741/58‑1 (MQ=255)               tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                              <  1:1801597/71‑1 (MQ=255)               tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                              <  1:2216256/71‑1 (MQ=255)               tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                              <  1:208590/71‑1 (MQ=255)               tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                              <  1:1922575/71‑1 (MQ=255)               tataTTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                              <  1:1861665/71‑1 (MQ=255)                atatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTc                                              <  1:3307978/70‑1 (MQ=255)                 tatTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                          >  1:2704037/1‑57 (MQ=255)                        taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                      >  1:1128460/1‑70 (MQ=255)                        taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                      >  1:2997064/1‑70 (MQ=255)                        taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                      >  1:59578/1‑70 (MQ=255)                        taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                      >  1:948381/1‑70 (MQ=255)                        taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                      >  1:1318661/1‑70 (MQ=255)                        taGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCAGGCATTACCTTCAATTgcgc                                      >  1:1884709/1‑70 (MQ=255)                             aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                      >  1:1246484/1‑65 (MQ=255)                              cGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                               <  1:445265/71‑1 (MQ=255)                              cGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                               <  1:1852567/71‑1 (MQ=255)                              cGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                               <  1:2880256/71‑1 (MQ=255)                              cGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                               <  1:1578621/71‑1 (MQ=255)                                cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                             <  1:2115124/71‑1 (MQ=255)                                cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                             <  1:676634/71‑1 (MQ=255)                                cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                             <  1:1376807/71‑1 (MQ=255)                                cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                             <  1:483559/71‑1 (MQ=255)                                cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                             <  1:4013919/71‑1 (MQ=255)                                cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                             <  1:345667/71‑1 (MQ=255)                                cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                             <  1:1536302/71‑1 (MQ=255)                                cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                             <  1:1550901/71‑1 (MQ=255)                                cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGt                             <  1:265743/71‑1 (MQ=255)                                cTATGCTTCCTGCGGAATTATATCCTTGCATGATTTCAGCATGGCATTACCTTCAATTGCGCCACCAAAGt                             <  1:1853894/71‑1 (MQ=255)                                    gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAt                          >  1:4111616/1‑70 (MQ=255)                                           gCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                               <  1:3581107/58‑1 (MQ=255)                                            cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTa                                                   <  1:558257/37‑1 (MQ=255)                                                  tataTCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCATAGt                             <  1:1114868/53‑1 (MQ=255)                                                   atatCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaaa                               <  1:2306087/50‑1 (MQ=255)                                                      tCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCAc       >  1:43318/1‑71 (MQ=255)                                                      tCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCAc       >  1:2330315/1‑71 (MQ=255)                                                      tCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCAc       >  1:446962/1‑71 (MQ=255)                                                      tCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCAc       >  1:1746258/1‑71 (MQ=255)                                                           gacTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg  >  1:3560240/3‑71 (MQ=255)                                                           gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACttt    >  1:1416338/1‑69 (MQ=255)                                                           gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg  >  1:3026730/1‑71 (MQ=255)                                                           gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg  >  1:3889507/1‑71 (MQ=255)                                                           gCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg  >  1:1299955/1‑71 (MQ=255)                                                            ccTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg  >  1:1908457/1‑70 (MQ=255)                                                            ccTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg  >  1:1986189/1‑70 (MQ=255)                                                            ccTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg  >  1:3236016/1‑70 (MQ=255)                                                            ccTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg  >  1:2017969/1‑70 (MQ=255)                                                            ccTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg  >  1:2022754/1‑70 (MQ=255)                                                            ccTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg  >  1:211360/1‑70 (MQ=255)                                                                   |                                                               ATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTG  >  W3110S.gb/155397‑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: ‑