Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	869,180	A→G	N554S (AAC→AGC)	gsiA →	glutathione ABC transporter ATP binding subunit GsiA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	869,180	0	A	G	100.0%	105.4 / NA	34	N554S (AAC→AGC)	gsiA	glutathione ABC transporter ATP binding subunit GsiA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/22);  total (12/22)

CAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTG‑TGGCGGTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGC  >  NC_000913/869047‑869312                                                                                                                                         |                                                                                                                                     cAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTG‑TGGCGGTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCgca                                                                                                                                <  1:1262232/140‑1 (MQ=255) cAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTG‑TGGCGGTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCgca                                                                                                                                <  1:804570/140‑1 (MQ=255) cAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTG‑TGGCGGTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCgca                                                                                                                                <  1:537187/140‑1 (MQ=255) cAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGGGGGGTTGTCATCGTG‑TGGCGGGGGTGGTTTTTTGGGGAAATTTTTTAAATTGGTCC‑CCCCGCGCGGGCCTTGGAAAAccccc                                                                                                                                 >  2:1199967/1‑138 (MQ=255)  gacaggATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTG‑TGGCGGTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCgcag                                                                                                                               <  1:1417541/135‑1 (MQ=255)  aGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTG‑TGGCGGTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCgcag                                                                                                                               <  1:524353/140‑1 (MQ=255)  aGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTG‑TGGCGGTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCgcag                                                                                                                               <  1:429390/140‑1 (MQ=255)       gATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTAGTCATCGTGTTGGGGGGGGG‑GTTTTTCGGGGAAATTTTTTAATTTGTTCCCCGGGGCCC‑GGGCTTGTAAAACCCcccgcacc                                                                                                                           >  2:1003067/1‑138 (MQ=255)                              tATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTG‑TGGCGGTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTggcga                                                                                                   <  1:1003067/140‑2 (MQ=255)                                   cccACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTG‑TGGCGGTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGtt                                                                                              <  1:1617934/140‑1 (MQ=255)                                     cACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTG‑TGGCGGTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTcc                                                                                            <  1:1199967/140‑1 (MQ=255)                                                ggtggtAGAGCGGATTAGTCATCGTG‑TGGCGGTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAAc                                                                                 <  1:895596/140‑1 (MQ=255)                                                        agCGGATTAGTCATCGTG‑TGGCGGTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGa                                                                         >  2:697859/1‑140 (MQ=255)                                                         gCGGATTAGTCATCGTG‑TGGCGGTGAT‑GTATCTCGGGGAAAATTTTTAAATTTGGCCCCGGCGCGC‑GGGCTTCGAAAACCCCCCGCATCCTTATACGCGTAAAATACTTGCGGGCGTTTCCGGCGCCGGAACGTcccccc                                                                        >  2:1557843/1‑137 (MQ=255)                                                            gATTAGTCATCGTG‑TGGCGGTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAAc                                                                     >  1:1379458/1‑140 (MQ=255)                                                            gATTAGTCATCGTG‑TGGCGGTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGGCGCTGAACCGTCCCGACAAc                                                                     >  1:187547/1‑140 (MQ=255)                                                            gATTAGTCATCGTG‑TGGCGGTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGGCGCTGAACCCTCCCGACaaa                                                                     >  2:682858/1‑139 (MQ=255)                                                            gATTAGTCATCGTG‑TGGCGGCGCT‑GTATCTCGGGCCAATTGTTGAAATTTGGCCACCGCGCCC‑GGTCTTCGGAAGCCCGCAGCCTCCTTATACGCGTAAATTACTTGCGGCAGTTCCCGTGGCTGAAACCTcccctaaaa                                                                     >  2:697417/1‑134 (MQ=255)                                                               tAGTCATCGTG‑TGGCGGTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGAc                                                                  >  1:1041514/1‑140 (MQ=255)                                                                 gTCATCGTG‑TGGCGGTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCg                                                                <  1:1557843/140‑1 (MQ=255)                                                                 gTCATCGTG‑TGGCGGTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCg                                                                <  1:934451/140‑1 (MQ=255)                                                                 gTCATCGTG‑TGGCGGTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCg                                                                <  1:697417/140‑1 (MQ=255)                                                                  tCATCGTG‑TGGCGGTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGc                                                               >  2:748315/1‑140 (MQ=255)                                                                           tgGCGGTGAT‑GTATCTCGGGCAAATCGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTa                                                       >  1:1357514/1‑140 (MQ=255)                                                                             ccGGTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTAct                                                     >  2:1450818/2‑140 (MQ=255)                                                                                gTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCGCGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTActgct                                                  >  2:1309593/1‑140 (MQ=255)                                                                                  gAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGt                                                <  1:697859/140‑1 (MQ=255)                                                                                      gTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCgg                                             <  2:646918/140‑1 (MQ=255)                                                                                      gTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCgg                                             <  2:1105666/140‑1 (MQ=255)                                                                                                 aaaTTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTccc                                  >  1:865861/1‑140 (MQ=255)                                                                                                   aTTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAg                                <  1:101697/140‑1 (MQ=255)                                                                                                   aTTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAg                                <  1:748315/140‑1 (MQ=255)                                                                                                   aTTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAg                                <  1:699963/140‑1 (MQ=255)                                                                                                   aTTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAg                                <  1:1450818/140‑1 (MQ=255)                                                                                                     tgttgAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCa                              <  2:1041514/140‑1 (MQ=255)                                                                                                             ttGGTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCAt                      <  1:682858/140‑1 (MQ=255)                                                                                                               ggTCCACGGCGCGC‑GGTCTTCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCt                    >  2:945874/1‑140 (MQ=255)                                                                                                                                  ttCGAAAGCCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGc  >  1:970930/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                     CAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTG‑TGGCGGTGAT‑GTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGC‑GGTCTTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGC  >  NC_000913/869047‑869312

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