Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	868,422	(C)8→7	coding (903/1872 nt)	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	868,415	0	C	.	100.0%	98.7 / NA	24	coding (896/1872 nt)	gsiA	glutathione ABC transporter ATP binding subunit GsiA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (10/14);  total (10/14)

TTTTCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGT  >  NC_000913/868287‑868547                                                                                                                                 |                                                                                                                                     ttttCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGa                                                                                                                          <  1:63160/140‑1 (MQ=255)         cACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGaaaa                                                                                                                  <  2:567700/140‑1 (MQ=255)              cAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACggtg                                                                                                             >  1:636160/1‑140 (MQ=255)              cAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACggtg                                                                                                             >  1:610055/1‑140 (MQ=255)                        acacCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCg                                                                                                   >  1:152335/1‑140 (MQ=255)                            ccGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAAcc                                                                                               <  2:636160/140‑1 (MQ=255)                            ccGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAAcc                                                                                               <  2:610055/140‑1 (MQ=255)                              gTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGgtggtg                                                                                                          <  1:585769/127‑1 (MQ=255)                              gTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGgtggtg                                                                                                          >  2:585769/1‑127 (MQ=255)                               tGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGt                                                                                            >  2:349927/1‑140 (MQ=255)                               tGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCAAGCAGTAATCGGTGGCGGATGGCGAACCTGt                                                                                            >  1:318071/1‑140 (MQ=255)                                    tGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTAc                                                                                       >  1:610128/1‑140 (MQ=255)                                     gTTAGCTGCTGTTTCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACg                                                                                      <  2:390615/140‑1 (MQ=255)                                      ttAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGa                                                                                     <  1:573548/140‑1 (MQ=255)                                      ttAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGa                                                                                     <  1:145604/140‑1 (MQ=255)                                              tgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTaa                                                                             <  2:526184/140‑1 (MQ=255)                                                     cAACTTGGCGCGATGATAGGGTTAGATTATTCCCGTCGTTTCCCGTTGTTATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCAAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTc                                                                      <  2:318071/140‑1 (MQ=255)                                                               cGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTcc                                                            <  1:349927/140‑1 (MQ=255)                                                                           tAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCg                                                <  2:152335/140‑1 (MQ=255)                                                                           tAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCg                                                <  1:483590/140‑1 (MQ=255)                                                                                  tCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGttgttg                                         >  1:611207/1‑136 (MQ=255)                                                                                        aCGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATcg                                   >  1:619667/1‑140 (MQ=255)                                                                                                             tGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGGGGAGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGGTTGTTGAATCGGGTATCGCGGGAAGTGCATGc              >  2:601310/1‑140 (MQ=255)                                                                                                                         gAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGt  <  1:516638/140‑1 (MQ=255)                                                                                                                                 |                                                                                                                                     TTTTCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGT  >  NC_000913/868287‑868547

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