Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,211,170	(T)8→7	coding (1246/2187 nt)	espL4 ←	type III secretion system effector EspL4

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,211,163	0	T	.	100.0%	125.2 / NA	30	coding (1253/2187 nt)	espL4	type III secretion system effector EspL4
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (14/16);  total (14/16)

TACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTAT  >  NZ_CP009273/4211078‑4211248                                                                                      |                                                                                      tACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGAttttt                                                                                   >  2:234195/1‑90 (MQ=255) tACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGAttttt                                                                                   >  1:694596/1‑90 (MQ=255) tACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGAttttt                                                                                   >  1:239293/1‑90 (MQ=255)        tttttATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGA‑TTTTTTTCgaga                                                                           <  1:121943/90‑1 (MQ=255)            tATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTaa                                                                       <  1:15087/90‑1 (MQ=255)             aTCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAAt                                                                      <  1:347971/90‑1 (MQ=255)              tCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAATc                                                                     <  1:402698/90‑1 (MQ=255)                  tACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAATCCTTc                                                                 >  2:681989/1‑90 (MQ=255)                  tACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAATCCTTc                                                                 >  2:26703/1‑90 (MQ=255)                  tACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAATCCTTc                                                                 >  1:95193/1‑90 (MQ=255)                    cGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAATCCTTCat                                                               <  2:611069/90‑1 (MQ=255)                             gCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAATCCTTCATATCCTGttt                                                      >  2:519720/1‑90 (MQ=255)                             gCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAATCCTTCATATCCTGttt                                                      >  2:306559/1‑90 (MQ=255)                               cagaaaTAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCt                                                    >  1:650399/4‑90 (MQ=255)                                  aaaTAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGaa                                                 >  1:20/1‑89 (MQ=255)                                      aatCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATg                                             <  2:699899/90‑1 (MQ=255)                                          cAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAAtt                                         <  1:250211/90‑1 (MQ=255)                                          cAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAAtt                                         <  1:320731/90‑1 (MQ=255)                                               aaaCCATTTTTATCTTTTGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGaaaa                                    <  1:238662/90‑1 (MQ=255)                                                 aCCATTTTTATCTTTTGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATa                                  <  2:348890/90‑1 (MQ=255)                                                        ttATCTTTTGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCaa                           <  2:522249/90‑1 (MQ=255)                                                             ttttGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGata                      <  2:126856/90‑1 (MQ=255)                                                             ttttGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGata                      <  2:287813/90‑1 (MQ=255)                                                             ttttGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGata                      >  2:570519/1‑90 (MQ=255)                                                              tttGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGatat                     >  1:85243/1‑90 (MQ=255)                                                              tttGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGatat                     >  2:308082/1‑90 (MQ=255)                                                               ttGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATc                    >  1:282560/1‑90 (MQ=255)                                                                tGCCTCCAGAATATGCATGAGA‑TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATcc                   <  1:681989/90‑1 (MQ=255)                                                                           tatGCATGAGA‑TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATa        <  2:63638/90‑1 (MQ=255)                                                                            atGCATGAGA‑TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATaa       >  1:432368/1‑90 (MQ=255)                                                                               cATGAGA‑TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAAttt    >  2:499104/1‑90 (MQ=255)                                                                                 tGAGA‑TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTtat  <  2:92960/90‑1 (MQ=255)                                                                                 tGAGA‑TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTtat  >  1:99859/1‑90 (MQ=255)                                                                                      |                                                                                      TACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTAT  >  NZ_CP009273/4211078‑4211248

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