Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	826,223	G→T	G207V (GGT→GTT)	ychP →	predicted invasin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	826,223	0	G	T	100.0%	98.8 / NA	31	G207V (GGT→GTT)	ychP	predicted invasin
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (15/16);  total (15/16)

ACAGCTCGGATGCGCATGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGG  >  minE/826156‑826285                                                                    |                                                               aCAGCTCGGATGCGCATGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTttgttg                                                              <  1:2837238/70‑1 (MQ=255) aCAGCTCGGATGCGCATGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGTTGa                                                             >  1:247133/1‑71 (MQ=255) aCAGCTCGGATGCGCATGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGTTGa                                                             >  1:250905/1‑71 (MQ=255)  cAGCTCGGATGCGCATGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGTTGAt                                                            <  1:2607736/71‑1 (MQ=255)   aGCTCGGATGCGCATGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGTTGATc                                                           >  1:1536721/1‑71 (MQ=255)        ggATGCGCATGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGTTGATCGTGtt                                                      <  1:1122266/71‑1 (MQ=255)        ggATGCGCATGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGTTGATCGTGtt                                                      <  1:3125471/71‑1 (MQ=255)          aTGCGCATGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGTTGATCGTGTTGa                                                    >  1:2412881/1‑71 (MQ=255)              gcATGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTttgttg                                                              <  1:1141585/57‑1 (MQ=255)              gcATGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTttgttg                                                              <  1:343728/57‑1 (MQ=255)              gcATGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTttgttg                                                              <  1:1285753/57‑1 (MQ=255)              gcATGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTttgttg                                                              <  1:931392/57‑1 (MQ=255)              gcATGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTttgttg                                                              <  1:117196/57‑1 (MQ=255)                aTGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGTTGATCGTGTTGATTTGtt                                              >  1:1540205/1‑71 (MQ=255)                 tGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTAATTTGTTGATCGTGTTGATTTGtt                                              >  1:1381468/1‑70 (MQ=255)                   ccGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGTTGATCGTGTTGATTTGTTTaa                                           <  1:2391089/71‑1 (MQ=255)                        cTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGTTGATCGTGTTGAtt                                                  >  1:607339/1‑59 (MQ=255)                        cTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGTTGATCGTGTTGATTTGTTTAACTCTg                                      >  1:648684/1‑71 (MQ=255)                             aaCACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGTTGATCGTGTTGATTTGTTTAACTCTGGTAc                                  <  1:2444701/70‑1 (MQ=255)                             aaCACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGTTGATCGTGTTGATTTGTTTAACTCTGGTAc                                  <  1:1299559/70‑1 (MQ=255)                             aaCACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGTTGATCGTGTTGATTTGTTTAACTCTGGTACg                                 >  1:2012164/1‑71 (MQ=255)                                     aaTACCAGTGTCAGCCTAGAACAGTATATTGTTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATc                          >  1:2839069/1‑70 (MQ=255)                                          cAGTGTCAGCCTAGAACAGTATTTTGTTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTAt                           <  1:2240252/64‑1 (MQ=255)                                          cAGTGTCAGCCTAGAACAGTATTTTGTTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAAt                     >  1:1575717/1‑70 (MQ=255)                                          cAGTGTCAGCCTAGAACAGTATTTTGTTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATc                    >  1:999003/1‑71 (MQ=255)                                          cAGTGTCAGCCTAGAACAGTATTTTGTTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATc                    <  1:1951119/71‑1 (MQ=255)                                           aGTGTCAGCCTAGAACAGTATTTTGTTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATcc                   <  1:752278/71‑1 (MQ=255)                                             tgtCAGCCTAGAACAGTATTTTGTTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCg                 <  1:584779/71‑1 (MQ=255)                                             tgtCAGCCTAGAACAGTATTTTGTTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCg                 <  1:77902/71‑1 (MQ=255)                                              gtCAGCCTAGAACAGTATTTTGTTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGt                >  1:1181149/1‑71 (MQ=255)                                               tCAGCCTAGAACAGTATTTTGTTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGt                >  1:2241652/1‑70 (MQ=255)                                                  gCCTAGAACAGTATTTTGTTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTAt                           <  1:2020660/56‑1 (MQ=255)                                                  gCCTAGAACAGTATTTTGTTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTcgcg            <  1:3109248/71‑1 (MQ=255)                                                    cTAGAACAGTATTTTGTTGATCGTGTTGATTTGttt                                             <  1:908656/36‑1 (MQ=255)                                                          cAGTATTTTGTTGATCGTGTTGATTTGTTTAACTCTGGTAc                                  <  1:325960/41‑1 (MQ=255)                                                           aGTATTTTGTTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTg   <  1:2500726/71‑1 (MQ=255)                                                            gTATTTTGTTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTgg  >  1:1383295/1‑71 (MQ=255)                                                                    |                                                               ACAGCTCGGATGCGCATGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGG  >  minE/826156‑826285

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