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%	213.5 / NA	65	*304* (TAG→TAA)	gsiD	glutathione ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/31);  total (34/31)

TGAGGCTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTG  >  NC_000913/872666‑872934                                                                                                                                        |                                                                                                                                      tGAGGCTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTa                                                                                                                                    <  1:634209/139‑1 (MQ=255) tGAGGCTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTa                                                                                                                                    <  1:1777027/139‑1 (MQ=255)      cTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGttt                                                                                                                               >  1:869316/1‑139 (MQ=255)               atatGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAata                                                                                                                          >  2:1303626/1‑135 (MQ=255)               atatGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGCACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAata                                                                                                                          <  1:1303626/135‑1 (MQ=255)                     ttATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGaaa                                                                                                                >  1:1134885/1‑139 (MQ=255)                     ttATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGaaa                                                                                                                <  2:20480/139‑1 (MQ=255)                       aTCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAgg                                                                                                              <  1:65546/139‑1 (MQ=255)                         cgcgCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGt                                                                                                            >  2:1409749/1‑139 (MQ=255)                         cgcgCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGt                                                                                                            >  2:186057/1‑139 (MQ=255)                         cgcgCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGt                                                                                                            >  1:1645789/1‑139 (MQ=255)                             ccgcATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATc                                                                                                        >  2:408573/1‑139 (MQ=255)                                          tttttCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGaaa                                                                                           >  2:322974/1‑139 (MQ=255)                                               ccGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATAtt                                                                                                                        >  1:668004/1‑105 (MQ=255)                                               ccGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATAtt                                                                                                                        <  2:668004/105‑1 (MQ=255)                                                    ccTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAAt                                                                                                         >  2:413421/1‑115 (MQ=255)                                                    ccTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAAt                                                                                                         <  1:413421/115‑1 (MQ=255)                                                      tGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTaa                                                                               >  1:897231/1‑139 (MQ=255)                                                         cTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAGAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCa                                                                            <  2:869316/139‑1 (MQ=255)                                                           aTTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAg                                                                          <  2:1134885/139‑1 (MQ=255)                                                             ttttCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCa                                                                            >  2:537381/1‑135 (MQ=255)                                                             ttttCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCa                                                                            <  1:537381/135‑1 (MQ=255)                                                                      cGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTa                                                                                                                                    <  2:755055/70‑1 (MQ=255)                                                                      cGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTa                                                                                                                                    >  1:755055/1‑70 (MQ=255)                                                                      cGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGaa                                                                                                                 <  1:717787/89‑1 (MQ=255)                                                                      cGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGaa                                                                                                                 >  1:1964731/1‑89 (MQ=255)                                                                      cGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGaa                                                                                                                 <  2:1964731/89‑1 (MQ=255)                                                                      cGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGaa                                                                                                                 >  2:717787/1‑89 (MQ=255)                                                                      cGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAAt                                                                                                         <  2:1058020/97‑1 (MQ=255)                                                                      cGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAAt                                                                                                         >  1:1058020/1‑97 (MQ=255)                                                                      cGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTca                                                                                                  >  1:46242/1‑104 (MQ=255)                                                                      cGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTca                                                                                                  <  2:46242/104‑1 (MQ=255)                                                                              cGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGa                                                       <  2:1645789/139‑1 (MQ=255)                                                                              cGTTCAATTAGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGa                                                       <  1:1409749/139‑1 (MQ=255)                                                                               gTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGcc                                                                                   >  1:838129/1‑110 (MQ=255)                                                                               gTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGcc                                                                                   <  2:838129/110‑1 (MQ=255)                                                                                ttCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTctgc                                                          >  2:986120/1‑134 (MQ=255)                                                                                ttCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTctgc                                                          <  1:986120/134‑1 (MQ=255)                                                                                 tCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACG‑TTGAATATTGCTTGAAAGGGTAAt                                                                                                         >  1:1484974/1‑85 (MQ=255)                                                                                 tCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACG‑TTGAATATTGCTTGAAAGGGTAAt                                                                                                         <  2:1484974/85‑1 (MQ=255)                                                                                  cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGaa                                                                                                                            <  2:1695445/66‑1 (MQ=255)                                                                                  cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGaa                                                                                                                            >  1:1695445/1‑66 (MQ=255)                                                                                  cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAAt                                                                                                         >  1:1565794/1‑85 (MQ=255)                                                                                  cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAAt                                                                                                         <  2:1565794/85‑1 (MQ=255)                                                                                  cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGAtttt                                                   >  2:116955/1‑139 (MQ=255)                                                                                  cAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGAtttt                                                   >  2:1052166/1‑139 (MQ=255)                                                                                       tgttgGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGACCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGa                                              >  2:763645/1‑139 (MQ=255)                                                                                        gttgGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGaa                                             <  1:1052166/139‑1 (MQ=255)                                                                                             gCGATGGTTTACGCAATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGt                                                                                                            <  1:932390/71‑1 (MQ=255)                                                                                             gCGATGGTTTACGCAATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGt                                                                                                            >  2:932390/1‑71 (MQ=255)                                                                                                aTGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGaaa                                                                                                                >  2:1155780/1‑64 (MQ=255)                                                                                                aTGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGaaa                                                                                                                <  1:1155780/64‑1 (MQ=255)                                                                                                 tGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCa                                                                            >  1:442604/1‑99 (MQ=255)                                                                                                 tGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCa                                                                            <  2:442604/99‑1 (MQ=255)                                                                                                           gATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAatta                                                                                        <  2:1095874/77‑1 (MQ=255)                                                                                                           gATGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAatta                                                                                        >  1:1095874/1‑77 (MQ=255)                                                                                                             tGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCt                                                                                                    <  1:1354230/63‑1 (MQ=255)                                                                                                             tGCGCTGGATCCGAAAATTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAg                        >  1:1273957/1‑139 (MQ=255)                                                                                                             tGCGCTGGATCCGAAAATAAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCt                                                                                                    >  2:1354230/1‑63 (MQ=255)                                                                                                                       ccGAAAATTAAAGGATAATTACGTTTGAATATTGCTTg                                                                                                                   >  2:1824966/1‑38 (MQ=255)                                                                                                                       ccGAAAATTAAAGGATAATTACGTTTGAATATTGCTTg                                                                                                                   <  1:1824966/38‑1 (MQ=255)                                                                                                                         gAAAATTTAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAACGCCGCAGCAAggg            >  1:1349253/1‑139 (MQ=255)                                                                                                                           aaaTTAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGAc          <  1:186057/139‑1 (MQ=255)                                                                                                                              ttAAAGGATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAAt       >  1:1497639/1‑139 (MQ=255)                                                                                                                                   ggATAATTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTg  <  2:1497639/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                      TGAGGCTCGAGCGGATATGGTTATCGCGCCGCATGTCGCTGTTTTTCCGGCCCTGGCTATTTTTCTGACCGTACTGGCGTTCAATTTGTTGGGCGATGGTTTACGCGATGCGCTGGATCCGAAAATTAAAGGATAGTTACGTTTGAATATTGCTTGAAAGGGTAATCACCTCACAGGAAATTATTGCCCTAAGCAAGTGTTGTAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTG  >  NC_000913/872666‑872934

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