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

TTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAA  >  W3110S.gb/155406‑155495                                                          |                                 tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                              >  1:2285296/1‑62 (MQ=255) tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                              >  1:3923515/1‑62 (MQ=255) tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                              >  1:3287852/1‑62 (MQ=255) tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                              >  1:1915878/1‑62 (MQ=255)   tAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGcc                             <  1:1883403/61‑1 (MQ=255)   tAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGcc                             <  1:2201151/61‑1 (MQ=255)   tAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGcc                             <  1:2735784/61‑1 (MQ=255)   tAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                            >  1:4079360/1‑62 (MQ=255)   tAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                            >  1:2846214/1‑62 (MQ=255)   tAGTATATTCTATAGTAACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                            >  1:4674235/1‑62 (MQ=255)                 gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa              >  1:2498263/1‑62 (MQ=255)                 gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa              >  1:2602569/1‑62 (MQ=255)                 gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa              >  1:493015/1‑62 (MQ=255)                  tCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt             >  1:4595983/1‑62 (MQ=255)                  tCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt             >  1:1582952/1‑62 (MQ=255)                  tCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt             >  1:3932790/1‑62 (MQ=255)                  tCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt             >  1:53521/1‑62 (MQ=255)                  tCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt             >  1:2711347/1‑62 (MQ=255)                  tCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt             >  1:2245323/1‑62 (MQ=255)                  tCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAt             >  1:1890402/1‑62 (MQ=255)                           gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgc          <  1:3429719/56‑1 (MQ=255)                           gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgc          <  1:3425007/56‑1 (MQ=255)                           gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgc          <  1:230611/56‑1 (MQ=255)                           gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgc          <  1:4351745/56‑1 (MQ=255)                           gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgc          <  1:481270/56‑1 (MQ=255)                           gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgc          <  1:1843406/56‑1 (MQ=255)                           gCTTCCTGCGGAATTATATCCTTGCATGATTTCAGCCTGGCATTACCTTCAATTgc          <  1:810103/56‑1 (MQ=255)                              tCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaa  <  1:3573322/61‑1 (MQ=255)                              tCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaa  <  1:1498605/61‑1 (MQ=255)                                                          |                                 TTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAA  >  W3110S.gb/155406‑155495

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