Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	869,338	C→T	P607S (CCG→TCG)	gsiA →	glutathione ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	869,338	0	C	T	100.0%	95.9 / NA	31	P607S (CCG→TCG)	gsiA	glutathione ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/19);  total (12/19)

TAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGC  >  NC_000913/869202‑869472                                                                                                                                         |                                                                                                                                       tAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCg                                                                                                                                      >  2:362256/1‑139 (MQ=255)      tACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCa                                                                                                                                 <  2:268294/139‑1 (MQ=255)                            tGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCa                                                                                                           >  2:284558/1‑139 (MQ=255)                                   tCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACg                                                                                                    <  2:309301/139‑1 (MQ=255)                                           aCGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATg                                                                                            <  2:933664/139‑1 (MQ=255)                                            cGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGc                                                                                           >  1:460153/1‑139 (MQ=255)                                              aCCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGt                                                                                         <  2:132077/139‑1 (MQ=255)                                                               tgTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCgg                                                                        >  2:486214/1‑139 (MQ=255)                                                                 tCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACAACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGgag                                                                      >  2:744783/1‑139 (MQ=255)                                                                            tCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGc                                                           <  1:284558/139‑1 (MQ=255)                                                                            tCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGc                                                           <  2:275696/139‑1 (MQ=255)                                                                            tCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGc                                                           <  1:530848/139‑1 (MQ=255)                                                                            tCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGc                                                           <  1:486214/139‑1 (MQ=255)                                                                            tCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGc                                                           <  1:362256/139‑1 (MQ=255)                                                                                 gCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAagag                                                      <  2:830105/139‑1 (MQ=255)                                                                                  cAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCg                                                                         <  2:341956/119‑1 (MQ=255)                                                                                  cAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCg                                                                         >  1:341956/1‑119 (MQ=255)                                                                                   aaTATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCt                                                    <  2:397404/139‑1 (MQ=255)                                                                                          aTCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACAcc                                             <  2:460153/139‑1 (MQ=255)                                                                                             tGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTa                                          <  1:744783/139‑1 (MQ=255)                                                                                              gCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAg                                         <  1:169779/139‑1 (MQ=255)                                                                                               cGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGt                                        >  2:332340/1‑139 (MQ=255)                                                                                                 tGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTgg                                      >  2:127609/1‑139 (MQ=255)                                                                                                       agagGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTaca                                               >  1:379994/1‑124 (MQ=255)                                                                                                       agagGTGGCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTaca                                               <  2:379994/124‑1 (MQ=255)                                                                                                              gCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAAc                                                                                                                   >  2:983963/1‑49 (MQ=255)                                                                                                              gCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAAc                                                                                                                   >  2:718081/1‑49 (MQ=255)                                                                                                              gCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAAc                                                                                                                   <  1:718081/49‑1 (MQ=255)                                                                                                              gCAGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAAc                                                                                                                   <  1:983963/49‑1 (MQ=255)                                                                                                                aGCCGTCTCGTTGCAATGCGTCGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTggg                       <  1:332340/139‑1 (MQ=255)                                                                                                                                   gtcgggtcggGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATg    <  2:523073/139‑1 (MQ=255)                                                                                                                                    tcgggtcgggGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATgg   >  1:319257/1‑139 (MQ=255)                                                                                                                                     cGGGTCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGc  >  2:52263/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                       TAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGC  >  NC_000913/869202‑869472

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