Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	872,430	Δ1 bp	coding (541/912 nt)	gsiD →	glutathione ABC transporter membrane subunit GsiD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	872,428	0	C	.	100.0%	100.3 / NA	24	coding (539/912 nt)	gsiD	glutathione ABC transporter membrane subunit GsiD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (5/19);  total (5/19)

GATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGC  >  NC_000913/872348‑872504                                                                                 |                                                                             gATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGc                                                                              <  1:165539/81‑1 (MQ=255)  aTCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGcc                                                                             <  2:165540/81‑1 (MQ=255)     gCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGccgcc                                                                          >  2:165585/1‑78 (MQ=255)      cTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGccgcct                                                                         <  1:165541/81‑5 (MQ=255)         ttgttgCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTG‑CCGCCTgg                                                                      <  2:165542/80‑1 (MQ=255)                 gtgtTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTG‑CCGCCTGGTGCGCGGCa                                                             <  1:165543/81‑1 (MQ=255)                    ttGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTG‑CCGCCTGGTGCGCGGCAac                                                           <  1:165544/80‑1 (MQ=255)                     tGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTG‑CCGCCTGGTGCGCGGCAacac                                                         >  2:165593/1‑81 (MQ=255)                        gAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTG‑CCGCCTGGTGCGCGGCAACACGCt                                                      <  2:165545/81‑1 (MQ=255)                           gCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTG‑CCGCCTGGTGCGCGGCAACACGCTGgt                                                   <  2:165546/81‑1 (MQ=255)                               cATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTG‑CCGCCTGGTGCGCGGCAACACGCTGGTGTTg                                               <  2:165547/81‑1 (MQ=255)                                 ttGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTG‑CCGCCTGGTGCGCGGCAACACGCTGGTGTTGaa                                             <  2:165548/81‑1 (MQ=255)                                  tGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTG‑CCGCCTGGTGCGCGGCAACACGCTGGTGTTGaaa                                            >  2:165583/1‑81 (MQ=255)                                   gCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTG‑CCGCCTGGTGCGCGGCAACACGCTGGTGTTGaaa                                            >  1:165582/1‑80 (MQ=255)                                       aCGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTG‑CCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCa                                       <  2:165549/81‑1 (MQ=255)                                                ttGCAGTCGCCATTTTTTCCATCCCCGCGTTTG‑CCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTAt                              <  2:165550/81‑1 (MQ=255)                                                      tCGCCATTTTTTCCATCCCCGCGTTTG‑CCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTc                        <  2:165551/81‑1 (MQ=255)                                                          cATTTTTTCCATCCCCGCGTTTG‑CCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCa                    <  2:165552/81‑1 (MQ=255)                                                            ttttttCCATCCCCGCGTTTG‑CCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACg                  <  2:165553/81‑1 (MQ=255)                                                               tttCCATCCCCGCGTTTG‑CCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAg               >  1:165591/1‑81 (MQ=255)                                                                ttCCATCCCCGCGTTTG‑CCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAg               <  2:165554/80‑1 (MQ=255)                                                                 tCCATCCCCGCGTTTG‑CCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGt              <  2:165555/80‑1 (MQ=255)                                                                   cATCCCCGCGTTTG‑CCGCCTGGTGCGCGGCAACACGCTGGTATTGAAGCAGCAAACCTTTATTGAGTCGGCACGCAGTAtt           <  2:165556/81‑1 (MQ=255)                                                                       cccGCGTTTG‑CCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTg       >  1:165589/1‑81 (MQ=255)                                                                         cgcgTTTG‑CCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGcc     <  2:165557/81‑1 (MQ=255)                                                                         cgcgTTTG‑CCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGcc     <  1:165558/81‑1 (MQ=255)                                                                            gTTTG‑CCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGc  <  2:165560/81‑1 (MQ=255)                                                                            gTTTG‑CCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGc  <  1:165559/81‑1 (MQ=255)                                                                                 |                                                                             GATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGC  >  NC_000913/872348‑872504

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: ‑