Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,536,106	T→C	E261G (GAA→GGA)	cysP ←	thiosulfate/sulfate ABC transporter substrate‑binding protein CysP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,536,106	0	T	C	100.0%	202.8 / NA	62	E261G (GAA→GGA)	cysP	thiosulfate/sulfate ABC transporter substrate‑binding protein CysP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (43/19);  total (43/19)

GTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAA  >  NZ_CP009273/2536021‑2536190                                                                                      |                                                                                     gTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGt                                                                                  <  1:1108267/90‑1 (MQ=255) gTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGt                                                                                  <  2:284678/90‑1 (MQ=255)   tATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTAc                                                                                >  1:360737/1‑90 (MQ=255)    aTTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTAcc                                                                               >  2:817188/1‑90 (MQ=255)    aTTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTAcc                                                                               >  1:311430/1‑90 (MQ=255)    aTTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTAcc                                                                               >  1:95075/1‑90 (MQ=255)                   aTAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCTGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGtt                                                                <  2:89463/90‑1 (MQ=255)                     aGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGtttt                                                              >  1:349952/1‑90 (MQ=255)                     aGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGtttt                                                              >  2:156637/1‑90 (MQ=255)                     aGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGtttt                                                              >  2:668124/1‑90 (MQ=255)                       tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGGTTTGGCGGGTTTTCCCGTACCGTTGGCCTGCACGTTTTTa                                                            >  1:773505/1‑90 (MQ=255)                       tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa                                                            >  2:1216396/1‑90 (MQ=255)                       tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa                                                            >  2:1004895/1‑90 (MQ=255)                       tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa                                                            >  2:238483/1‑90 (MQ=255)                       tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa                                                            >  1:773235/1‑90 (MQ=255)                       tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa                                                            >  2:942184/1‑90 (MQ=255)                       tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa                                                            >  1:1128730/1‑90 (MQ=255)                       tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa                                                            >  1:288854/1‑90 (MQ=255)                         ggTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGtt                                                                <  2:210649/84‑1 (MQ=255)                         ggTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGtt                                                                >  1:210649/1‑84 (MQ=255)                         ggTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATc                                                          <  2:1159300/90‑1 (MQ=255)                          gTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCa                                                         >  1:1320281/1‑90 (MQ=255)                          gTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCa                                                         >  1:361068/1‑90 (MQ=255)                          gTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCa                                                         >  2:1249546/1‑90 (MQ=255)                          gTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCa                                                         >  2:1325740/1‑90 (MQ=255)                               gatgGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGt                                                                                  >  2:349547/1‑60 (MQ=255)                               gatgGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGt                                                                                  <  1:349547/60‑1 (MQ=255)                                    tttGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCa                                               >  2:727989/1‑90 (MQ=255)                                    tttGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCa                                               >  2:275029/1‑90 (MQ=255)                                              cGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTc                                     <  2:549952/90‑1 (MQ=255)                                                         gCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCGCGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGt                          >  2:609189/1‑90 (MQ=255)                                                         gCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGt                          >  2:621906/1‑90 (MQ=255)                                                         gCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGt                          >  1:334135/1‑90 (MQ=255)                                                         gCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGt                          >  1:182038/1‑90 (MQ=255)                                                            aGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAGTGTTgg                       >  1:1416639/1‑90 (MQ=255)                                                            aGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTgg                       >  1:314718/1‑90 (MQ=255)                                                            aGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTgg                       >  2:1011031/1‑90 (MQ=255)                                                            aGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTgg                       >  2:572928/1‑90 (MQ=255)                                                            aGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTgg                       >  2:348791/1‑90 (MQ=255)                                                            aGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTgg                       >  1:934819/1‑90 (MQ=255)                                                               tCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGttt                    <  1:1131899/90‑1 (MQ=255)                                                                cAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGtttt                   <  2:885275/90‑1 (MQ=255)                                                                 aGATAGGCTTTGGCGGC‑TTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTg                        >  2:482414/1‑83 (MQ=255)                                                                 aGATAGGCTTTGGCGGC‑TTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTg                        <  1:482414/83‑1 (MQ=255)                                                                   aTAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCgg                >  2:453126/1‑90 (MQ=255)                                                                       gCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATc            >  2:361255/1‑90 (MQ=255)                                                                       gCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATc            >  1:1404802/1‑90 (MQ=255)                                                                          ttGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATcacc         >  2:405154/1‑90 (MQ=255)                                                                          ttGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATcacc         <  2:216100/90‑1 (MQ=255)                                                                          ttGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATcacc         <  1:668124/90‑1 (MQ=255)                                                                          ttGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATcacc         <  2:349952/90‑1 (MQ=255)                                                                          ttGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATcacc         <  2:597614/90‑1 (MQ=255)                                                                          ttGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATCCCGCCAGAATGTTGGTTTTCGGAATcacc         <  2:1480931/90‑1 (MQ=255)                                                                           tGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATcacca        >  2:15465/1‑90 (MQ=255)                                                                                gcTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGtt                         <  2:1517047/68‑1 (MQ=255)                                                                                gcTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGtt                         >  1:1517047/1‑68 (MQ=255)                                                                                gcTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTc                  <  2:583087/75‑1 (MQ=255)                                                                                gcTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTc                  >  1:583087/1‑75 (MQ=255)                                                                                gcTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCa   >  1:908230/1‑90 (MQ=255)                                                                                gcTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCa   >  2:686437/1‑90 (MQ=255)                                                                                 cTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCaa  <  1:1299046/90‑1 (MQ=255)                                                                                 cTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCaa  <  1:1257826/90‑1 (MQ=255)                                                                                      |                                                                                     GTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAA  >  NZ_CP009273/2536021‑2536190

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