Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	522,457	T→G	L681R (CTG→CGG)	ybbP →	putative ABC transporter membrane subunit YbbP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	522,457	0	T	G	100.0%	69.0 / NA	20	L681R (CTG→CGG)	ybbP	putative ABC transporter membrane subunit YbbP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (12/8);  total (12/8)

TCCGCTGGGAGAATGGCAACGGCATGTTGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGTGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCT  >  NC_000913/522331‑522586                                                                                                                               |                                                                                                                                  tCCGCTGGGAGAATGGCAACGGCATGTTGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCGGGAAGt                                                                                                                             <  2:443001/133‑1 (MQ=255)      tGGGAGAATGACAACGGCATGTTGTCACTACCCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCGGGAAGTGAtgg                                                                                                                        <  2:219410/133‑1 (MQ=255)                    cGGCATGTTGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCGGGAAGTGATGGTGGTACTGGTCAcc                                                                                                          <  1:329740/133‑1 (MQ=255)                         tttgaacacacCTAAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCGGGAAGTGGTGGTGGTACTGGTCACCGCCTg                                                                                                     <  2:157421/128‑1 (MQ=255)                          ttcgACACAACCCAAGCGCCAGTTCCCGACAATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCGGGAAGTGATGGTGGTACTGGTCACCGCCTGc                                                                                                    <  1:236512/130‑1 (MQ=255)                          gTTGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCGGGAAGTGATGGTGGTACTGGTCACCGCCTGc                                                                                                    >  2:45608/1‑133 (MQ=255)                           ttGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCGGGAAGTGATGGTGGTACTGGTCACCGCCTGCg                                                                                                   <  1:516745/133‑1 (MQ=255)                            tGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCGGGAAGTGATGGTGGTACTGGTCACCGCCTGCgg                                                                                                  <  1:45608/133‑1 (MQ=255)                                            ccAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCGGGAAGTGATGGTGGTACTGGTCACCGCCTGCGGTATGttgctgttgctg                                                                                  >  2:323161/1‑133 (MQ=255)                                                          ttAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCGGGAAGTGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCAcaggcgcaggt                                                                    >  2:36026/1‑133 (MQ=255)                                                                  ttAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCGGGAAGTGATGGTGGTACTGGTCACCGCCTCCGGTATGTAGCTGTTGGTGGCACAGGTGCATGTGGGAATgc                                                            >  2:240501/1‑133 (MQ=255)                                                                      atatTGGCGCCATTTTAAAACAGGTCGGTCAGGTGCTGGAGAAGGTAAGTCGGGCGCGGGAAGTGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATgcgtca                                                        <  2:188358/133‑1 (MQ=255)                                                                            gcgcGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCGGGAAGTGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATgcgtcagcgtca                                                  >  2:426362/1‑133 (MQ=255)                                                                            gcgcGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCGGGAAGTGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATgcgtcagcgtca                                                  >  1:20351/1‑133 (MQ=255)                                                                            gcgcGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCGGGAAGTGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATgcgtcagcgtca                                                  >  2:88372/1‑133 (MQ=255)                                                                                            ggtcggtcAGGTGCTGGAGCAGGTAAGTCGGGCGCGGGAAGTGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGgtggtg                                  >  1:586108/1‑133 (MQ=255)                                                                                            ggtcggtcAGGTGCTGGAGCAGGTAAGTCGGGCGCGGGAAGTGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGACATCAGGAGCTGgtggtg                                  >  1:257036/1‑133 (MQ=255)                                                                                             gtcggtcAGGTGCTGGAGCAGGTAAGTCGGGCGCGGGAAGTGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGgtgt                                 >  2:582688/1‑133 (MQ=255)                                                                                                 gtcAGGTGCTGGAGCAGGTAAGTCAGGCGCGGGAAGTGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGgcg                             >  2:102896/1‑133 (MQ=255)                                                                                                      gTGCTGGAGCAGGTAAGTCGGGCGCGGGAAGTGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCacac                        >  2:143127/1‑133 (MQ=255)                                                                                                                            cgcgGGAAGTGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAActgct  <  1:582688/133‑1 (MQ=255)                                                                                                                               |                                                                                                                                  TCCGCTGGGAGAATGGCAACGGCATGTTGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGTGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCT  >  NC_000913/522331‑522586

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