Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,272,642	C→T	P46S (CCG→TCG)	chaC →	glutathione‑specific gamma‑glutamylcyclotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,272,642	0	C	T	100.0%	113.3 / NA	37	P46S (CCG→TCG)	chaC	glutathione‑specific gamma‑glutamylcyclotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/17);  total (20/17)

TAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGCCTGCGCCTGACCGCCGGGCGCG  >  NC_000913/1272511‑1272765                                                                                                                                    |                                                                                                                            tAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAA‑‑‑‑CGGCTTCGCTGGCGGCGACGCTGGCTTATCGACCTGATGAGCGATCGGTGTGGATcc                                                                                                                 >  1:339551/1‑139 (MQ=255)   aCGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTAGCTGGCGGCGACACTGTCTTGTCGAACTGACGAGGGATCGGTGTGGAt                                                                                                                   >  2:290347/1‑140 (MQ=255)      cgTGATTTCTTGATGAATGCCGCTTGTAAAACGGCGTTTTGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCtt                                                                                                                <  2:340520/140‑1 (MQ=255)      cgTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGATGCTAGATTGTCGACCTTATGAGTGATCGGTGTGGATcct                                                                                                                >  2:389227/1‑138 (MQ=255)      cgTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCtt                                                                                                                <  1:715641/140‑1 (MQ=255)       gTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTc                                                                                                               >  1:71114/1‑140 (MQ=255)        tGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCtt                                                                                                                <  1:45931/138‑1 (MQ=255)        tGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCtt                                                                                                                >  2:45931/1‑138 (MQ=255)        tGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCg                                                                                                              >  1:514228/1‑140 (MQ=255)           tttCTTGATGAATGCGGATTGTGAATCGGGATTTGGTTCCATTGAAGAATCAGTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCt                                                                                                           <  1:315808/140‑1 (MQ=255)                 gatgaATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGAt                                                                                                     >  2:138266/1‑140 (MQ=255)                   tgaATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGGCCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCg                                                                                                   <  1:138266/140‑1 (MQ=255)                            aTTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTgg                                                                                          >  2:606463/1‑140 (MQ=255)                            aTTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTgg                                                                                          >  2:37512/1‑140 (MQ=255)                                     cGGCATTTGGTGTCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCa                                                                                 <  1:818633/140‑1 (MQ=255)                                          tttGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGa                                                                            <  1:37512/140‑1 (MQ=255)                                                    ttGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGaa                                                                  >  1:800822/1‑140 (MQ=255)                                                         gaaTCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTg                                                             <  2:71114/140‑1 (MQ=255)                                                                     tGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACt                                                 <  1:48965/140‑1 (MQ=255)                                                                            cagAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGGCCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGa                                          >  2:818966/1‑140 (MQ=255)                                                                            cagAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTgg                                               <  2:472912/135‑1 (MQ=255)                                                                            cagAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTgg                                               >  1:472912/1‑135 (MQ=255)                                                                            cagAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGAAATCCGGCACTGGAGTGTACCGAATGCTGCACCGGTACACTGGTTGGa                                          >  2:240152/1‑140 (MQ=255)                                                                               aacaacGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATgg                                       >  2:25379/1‑140 (MQ=255)                                                                                    cGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATcc                                  >  2:738590/1‑139 (MQ=255)                                                                                    cGGGCGGCTTCGCTGCCGGCGATGCGGGCATGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATcg                                  <  1:290347/140‑1 (MQ=255)                                                                                      ggcggcTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTCCCGCATCTTGCACCGGAACACAGGTTGGATGGCATcgcg                                >  2:299774/1‑140 (MQ=255)                                                                                      ggcggcTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACGGAATCGTCCACGGGGACACTGTTTGGATGGCATcgcg                                >  1:255523/1‑140 (MQ=255)                                                                                      ggcggcTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATTGCATcgcg                                >  2:769331/1‑140 (MQ=255)                                                                                      ggcggcTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATcgcg                                >  1:635672/1‑140 (MQ=255)                                                                                         ggcTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCAt                             >  1:677528/1‑140 (MQ=255)                                                                                         ggcTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCAt                             >  2:763314/1‑140 (MQ=255)                                                                                                 tGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGTTCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGCCTg                     <  1:606463/140‑1 (MQ=255)                                                                                                        gACGCTGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGCCTGCGCCTGa              <  1:818966/140‑1 (MQ=255)                                                                                                             tGGCTTGTCGACCTGGTGAGGGATCGGTGGGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGCCTGCGCCTGAccgcc         <  1:299774/140‑1 (MQ=255)                                                                                                             tGGCTTGTCGACCTGATGAGGGGTCGGTGTGGATCGTCGACTATGGATGGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGCCTGCGCCTGAccgcc         <  2:255523/140‑1 (MQ=255)                                                                                                             tGGCTTGTCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGCCTGCGCCTGAccgcc         <  1:769331/140‑1 (MQ=255)                                                                                                                    tCGACCTGATGAGGGATCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGCCTGCGCCTGACCGCCGGGcgcg  <  2:635672/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                            TAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGCCTGCGCCTGACCGCCGGGCGCG  >  NC_000913/1272511‑1272765

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