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%	211.6 / NA	50	coding (896/1872 nt)	gsiA	glutathione ABC transporter ATP binding subunit GsiA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (21/29);  total (21/29)

TTTCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGT  >  NC_000913/868288‑868547                                                                                                                                |                                                                                                                                     tttCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAg                                                                                                                         <  1:453201/140‑1 (MQ=255) tttCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAg                                                                                                                         <  2:1074444/140‑1 (MQ=255)        cACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGaaaa                                                                                                                  <  2:885907/140‑1 (MQ=255)           cGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACgg                                                                                                               >  2:401095/1‑140 (MQ=255)             cAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGaa                                                                                                                    <  2:614222/133‑1 (MQ=255)             cAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGaa                                                                                                                    >  1:976666/1‑133 (MQ=255)             cAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGaa                                                                                                                    <  2:976666/133‑1 (MQ=255)             cAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGaa                                                                                                                    >  1:614222/1‑133 (MQ=255)                     ttACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATgg                                                                                                     <  2:511517/140‑1 (MQ=255)                         acCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCTGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGaa                                                                                                 <  2:61329/140‑1 (MQ=255)                         acCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGaa                                                                                                 <  2:250317/140‑1 (MQ=255)                         acCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGaa                                                                                                 <  1:997126/140‑1 (MQ=255)                              tGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGa                                                                                                        >  2:974525/1‑128 (MQ=255)                              tGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGa                                                                                                        <  1:974525/128‑1 (MQ=255)                              tGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGt                                                                                            >  2:993213/1‑140 (MQ=255)                              tGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGt                                                                                            >  1:531461/1‑140 (MQ=255)                              tGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGt                                                                                            >  1:530100/1‑140 (MQ=255)                              tGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGt                                                                                            >  2:601662/1‑140 (MQ=255)                                   tGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTAc                                                                                       <  2:657641/140‑1 (MQ=255)                                    gTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTCTTACg                                                                                      <  2:318219/140‑1 (MQ=255)                                          tgctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCg                                                                                >  1:995250/1‑140 (MQ=255)                                           gctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGt                                                                               <  2:30215/140‑1 (MQ=255)                                            ctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTa                                                                              >  1:964996/1‑140 (MQ=255)                                            ctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTa                                                                              <  1:37251/140‑1 (MQ=255)                                            ctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTa                                                                              <  2:537522/140‑1 (MQ=255)                                                    cAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTc                                                                      <  1:564884/140‑1 (MQ=255)                                                    cAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTc                                                                      <  1:401095/140‑1 (MQ=255)                                                           gTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGtt                                                               <  2:530100/140‑1 (MQ=255)                                                                aTGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCt                                                          <  2:826825/140‑1 (MQ=255)                                                                     aGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTgcg                                                     <  2:988531/140‑1 (MQ=255)                                                                            gATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGt                                              <  1:649860/140‑1 (MQ=255)                                                                                 tCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTttgtt                                         >  1:1103965/1‑140 (MQ=255)                                                                                 tCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTttgtt                                         >  2:460006/1‑140 (MQ=255)                                                                                   cccGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGa                                       >  1:129791/1‑140 (MQ=255)                                                                                      gACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCg                                                <  2:459645/128‑1 (MQ=255)                                                                                      gACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCg                                                >  1:459645/1‑128 (MQ=255)                                                                                       aCGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATcg                                   >  2:819850/1‑140 (MQ=255)                                                                                          tttCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGt                                <  1:928787/140‑1 (MQ=255)                                                                                                 ttGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTCGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCgg                         <  2:995250/140‑1 (MQ=255)                                                                                                  tGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCggg                        <  2:63508/140‑1 (MQ=255)                                                                                                  tGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCggg                        <  2:531461/140‑1 (MQ=255)                                                                                                   gATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGAGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGa                       >  1:172004/1‑140 (MQ=255)                                                                                                      atCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGt                    >  1:904720/1‑140 (MQ=255)                                                                                                          cTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCAt                >  1:768359/1‑140 (MQ=255)                                                                                                            tGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGc              >  2:664324/1‑140 (MQ=255)                                                                                                              aaCATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCg            >  2:253227/1‑140 (MQ=255)                                                                                                                        gAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGt  <  2:1103965/140‑1 (MQ=255)                                                                                                                        gAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGt  <  1:664324/140‑1 (MQ=255)                                                                                                                           aCAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGt  >  2:835979/1‑137 (MQ=255)                                                                                                                           aCAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGt  <  1:835979/137‑1 (MQ=255)                                                                                                                                |                                                                                                                                     TTTCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGT  >  NC_000913/868288‑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: ‑