Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	872,801	G→A	*304* (TAG→TAA)	gsiD →	glutathione ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	872,801	0	G	A	100.0%	165.2 / NA	51	*304* (TAG→TAA)	gsiD	glutathione ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/28);  total (23/28)

TGAGGCTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGG  >  NC_000913/872666‑872924                                                                                                                                        |                                                                                                                            tGAGGCTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTa                                                                                                                          <  2:384666/139‑1 (MQ=255)    ggCTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGt                                                                                                                       >  2:802876/1‑139 (MQ=255)    ggCTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGt                                                                                                                       <  2:1090437/139‑1 (MQ=255)      cTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGATATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGttt                                                                                                                     >  1:906577/1‑139 (MQ=255)         gAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGaa                                                                                                                  <  2:1333082/139‑1 (MQ=255)           gCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAata                                                                                                                <  1:525933/139‑1 (MQ=255)                    gTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGaa                                                                                                       <  2:507254/139‑1 (MQ=255)                     ttATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGaaa                                                                                                      <  2:1387956/139‑1 (MQ=255)                         cgcgCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACGGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGt                                                                                                  >  2:823603/1‑139 (MQ=255)                            gccgcATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAAt                                                                                               <  1:1469390/139‑1 (MQ=255)                             ccgcATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATc                                                                                              >  2:852948/1‑139 (MQ=255)                              cgcATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCa                                                                                             >  2:1012910/1‑139 (MQ=255)                                   gTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTca                                                                                        <  1:1012910/139‑1 (MQ=255)                                   gTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTca                                                                                        <  1:832943/139‑1 (MQ=255)                                                 ggCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGc                                                                          >  2:378627/1‑139 (MQ=255)                                                  gCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGcc                                                                         <  1:1250624/139‑1 (MQ=255)                                                  gCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGcc                                                                         <  1:161249/139‑1 (MQ=255)                                                        gCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGc                                                                   >  2:1292051/1‑139 (MQ=255)                                                             ttttCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGa                                                                                   <  1:1354997/118‑1 (MQ=255)                                                             ttttCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGa                                                                                   >  2:1354997/1‑118 (MQ=255)                                                                    aCCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTca                                                                                        >  1:1455895/1‑106 (MQ=255)                                                                    aCCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTca                                                                                        <  2:1455895/106‑1 (MQ=255)                                                                      cGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGaaa                                                                                                      >  1:40929/1‑90 (MQ=255)                                                                      cGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGaaa                                                                                                      <  2:40929/90‑1 (MQ=255)                                                                      cGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAAt                                                                                               >  1:583773/1‑97 (MQ=255)                                                                      cGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAAt                                                                                               <  2:583773/97‑1 (MQ=255)                                                                       gTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACttt                                                    >  1:1273847/1‑139 (MQ=255)                                                                           tGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTctgc                                                <  1:802876/139‑1 (MQ=255)                                                                                ttCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGAtt                                           >  1:49565/1‑139 (MQ=255)                                                                                  cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTa                                                                                                                          >  2:1386359/1‑58 (MQ=255)                                                                                  cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTa                                                                                                                          <  1:1386359/58‑1 (MQ=255)                                                                                  cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGaa                                                                                                       >  2:25414/1‑77 (MQ=255)                                                                                  cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGaa                                                                                                       <  1:25414/77‑1 (MQ=255)                                                                                  cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGAtttt                                         >  2:1476388/1‑139 (MQ=255)                                                                                         ttgGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAAt                                  >  1:1441820/1‑139 (MQ=255)                                                                                           gggCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAg                                                                >  2:88714/1‑107 (MQ=255)                                                                                           gggCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAg                                                                <  1:88714/107‑1 (MQ=255)                                                                                                 tGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAAt                          <  1:1476388/139‑1 (MQ=255)                                                                                                   gTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAAttt                        <  1:1292051/139‑1 (MQ=255)                                                                                                      tACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCa                                                                  <  2:1476248/94‑1 (MQ=255)                                                                                                      tACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCa                                                                  >  1:1476248/1‑94 (MQ=255)                                                                                                             tGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCt                                                                                          >  2:751872/1‑63 (MQ=255)                                                                                                             tGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCt                                                                                          <  1:751872/63‑1 (MQ=255)                                                                                                             tGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGAtttt                                         >  2:343481/1‑112 (MQ=255)                                                                                                             tGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGAtttt                                         <  1:343481/112‑1 (MQ=255)                                                                                                              gcgcTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGc             <  2:1361893/139‑1 (MQ=255)                                                                                                              gcgcTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGc             <  2:49565/139‑1 (MQ=255)                                                                                                              gcgcTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGc             >  1:848833/1‑139 (MQ=255)                                                                                                                         gAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGTCGCAGCAAcgg  <  1:852948/139‑1 (MQ=255)                                                                                                                              ttAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGa                                             <  1:85473/91‑1 (MQ=255)                                                                                                                              ttAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGa                                             >  2:85473/1‑91 (MQ=255)                                                                                                                                        |                                                                                                                            TGAGGCTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGG  >  NC_000913/872666‑872924

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