Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	871,477	C→T	Q171* (CAG→TAG)	gsiC →	glutathione ABC transporter membrane subunit GsiC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	871,477	0	C	T	100.0%	107.4 / NA	30	Q171* (CAG→TAG)	gsiC	glutathione ABC transporter membrane subunit GsiC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/7);  total (23/7)

TTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTT  >  NC_000913/871400‑871553                                                                              |                                                                             ttGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGgtag                                                                            >  1:177344/1‑80 (MQ=255) ttGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGgtag                                                                            >  1:177359/1‑80 (MQ=255) ttGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGgtag                                                                            >  2:177368/1‑80 (MQ=255)   gCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGTAGCAc                                                                         >  1:177366/1‑81 (MQ=255)      ctGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGTAGCACTa                                                                       >  2:177354/1‑80 (MQ=255)          ggATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGTAGCACTACAttt                                                                  >  1:177358/1‑81 (MQ=255)          ggATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGTAGCACTACAttt                                                                  <  2:177321/81‑1 (MQ=255)            aTGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGTAGCACTACATTTTa                                                                >  2:177353/1‑81 (MQ=255)                     aTTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGTAGCACTACATTTTAccctccct                                                        >  1:177365/1‑80 (MQ=255)                       tCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGTAGCACTACATTTTACCCTCCCTGAc                                                     >  2:177367/1‑81 (MQ=255)                         aGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGTAGCACTACATTTTACCCTCCCTGAc                                                     <  2:177322/79‑1 (MQ=255)                           gTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGTAGCACTACATTTTACCCTCCCTGACGCTc                                                 <  1:177323/81‑1 (MQ=255)                           gTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGTAGCACTACATTTTACCCTCCCTGACGCTc                                                 >  2:177364/1‑81 (MQ=255)                                  ccGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGTAGCACTACATTTTACCCTCCCTGACGCTCGGCGCgg                                          >  1:177361/1‑81 (MQ=255)                                   cGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGTAGCACTACATTTTACCCTCCCTGACGCTCGGCGCgg                                          >  2:177352/1‑80 (MQ=255)                                       gAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGTAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTg                                     >  1:177375/1‑81 (MQ=255)                                         aCTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGTAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGgc                                   >  1:177360/1‑81 (MQ=255)                                          cTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGTAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGgcc                                  >  2:177374/1‑81 (MQ=255)                                             ggctggctGCCTACCGTGGGAGCAGACAGTTGGTAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGgccgcc                               >  1:177377/1‑81 (MQ=255)                                              gctggctgCCTACCGTGGGAGCAGACAGTTGGTAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGccgccg                              >  1:177372/1‑81 (MQ=255)                                                tggctgCCTACCGTGGGAGCAGACAGTTGGTAGCACTACATTTTACCCTTCCTGACGCTCGGCGCGGCAGTGGCCGCCGt                             >  2:177378/1‑80 (MQ=255)                                                     gCCTACCGTGGGAGCAGACAGTTGGTAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGgc                       >  1:177379/1‑81 (MQ=255)                                                         aCCGTGGGAGCAGACAGTTGGTAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGgcgcgc                   >  1:177371/1‑81 (MQ=255)                                                            gTGGGAGCAGACAGTTGGTAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCttt                <  1:177324/81‑1 (MQ=255)                                                              gggAGCAGACAGTTGGTAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTAc              >  2:177386/1‑81 (MQ=255)                                                               ggAGCAGACAGTTGGTAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTAcc             >  1:177363/1‑81 (MQ=255)                                                               ggAGCAGACAGTTGGTAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTAcc             >  2:177370/1‑81 (MQ=255)                                                                gAGCAGACAGTTGGTAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTAccc            <  1:177325/81‑1 (MQ=255)                                                                 aGCAGACAGTTGGTAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCcg           >  1:177362/1‑81 (MQ=255)                                                                  gCAGACAGTTGGTAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCcgc          >  1:177393/1‑81 (MQ=255)                                                                    aGACAGTTGGTAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCcgcgc        <  2:177326/81‑1 (MQ=255)                                                                    aGACAGTTGGTAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCcgcgc        >  2:177389/1‑81 (MQ=255)                                                                          ttGGTAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGtt  <  2:177327/81‑1 (MQ=255)                                                                              |                                                                             TTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCGTT  >  NC_000913/871400‑871553

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