Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,181,866	(T)5→4	coding (964/1047 nt)	potD ←	spermidine preferential ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,181,862	0	T	.	100.0%	94.6 / NA	23	coding (968/1047 nt)	potD	spermidine preferential ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (5/18);  total (5/18)

ACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAA  >  NC_000913/1181786‑1181937                                                                             |                                                                            aCGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCAtttt                                                                          >  1:213126/1‑80 (MQ=255)         cTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCA‑TTTTAATGGTTTc                                                                <  1:213103/81‑1 (MQ=255)               gCTTCTGATAATACTCTTCATAAATGCTGCT‑GCTGCGCCAACGTCATTCTGCCATTCGCCA‑TTTTAATGGTTTCAGCATcc                                                         >  1:213140/1‑81 (MQ=255)                    tGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCA‑TTTTAATGGTTTCAGCATCCggg                                                      >  2:213132/1‑80 (MQ=255)                      ataataCTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCA‑TTTTAATGGTTTCAGCATCCGGGTag                                                   <  2:213104/81‑1 (MQ=255)                         ataCTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCA‑TTTTAATGGTTTCAGCATCCGGGTAGAgt                                                <  1:213106/81‑1 (MQ=255)                         ataCTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCA‑TTTTAATGGTTTCAGCATCCGGGTAGAgt                                                <  2:213105/81‑1 (MQ=255)                            ctctTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCA‑TTTTAATGGTTTCAGCATCCGGGTAGAGTGtt                                             <  1:213107/81‑1 (MQ=255)                                 cATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCA‑TTTTAATGGTTTCAGCATCCGGGTAGAGTGtttt                                           <  2:213108/78‑1 (MQ=255)                                 cATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCA‑TTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTAt                                         <  1:213109/80‑1 (MQ=255)                                     aaTGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCA‑TTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTc                                    <  2:213110/81‑1 (MQ=255)                                            ctggctgCGCCAACGTCATTCTGCCATTCGCCA‑TTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCAc                               <  2:213111/79‑1 (MQ=255)                                            ctggctgCGCCAACGTCATTCTGCCATTCGCCA‑TTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACtt                             >  1:213133/1‑81 (MQ=255)                                            ctggctgCGCCAACGTCATTCTGCCATTCGCCA‑TTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACtt                             <  1:213113/81‑1 (MQ=255)                                            ctggctgCGCCAACGTCATTCTGCCATTCGCCA‑TTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACtt                             <  2:213112/81‑1 (MQ=255)                                                 tgCGCCAACGTCATTCTGCCATTCGCCA‑TTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGa                        >  1:213134/1‑81 (MQ=255)                                                    gcCAACGTCATTCTGCCATTCGCCA‑TTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACtt                     <  2:213114/81‑1 (MQ=255)                                                      cAACGTCATTCTGCCATTCGCCA‑TTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTaa                   <  1:213115/81‑1 (MQ=255)                                                         cGTCATTCTGCCATTCGCCA‑TTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAg                <  2:213116/81‑1 (MQ=255)                                                                cTGCCATTCGCCA‑TTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTAcgc         >  2:213135/1‑81 (MQ=255)                                                                   ccATTCGCCA‑TTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGcgcc      <  1:213117/81‑1 (MQ=255)                                                                   ccATTCGCCA‑TTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGcgcc      <  2:213119/81‑1 (MQ=255)                                                                   ccATTCGCCA‑TTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGcgcc      <  1:213118/81‑1 (MQ=255)                                                                        cGCCA‑TTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCaa  <  1:213120/80‑1 (MQ=255)                                                                             |                                                                            ACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAA  >  NC_000913/1181786‑1181937

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