Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,577,704	C→T	V547I (GTC→ATC)	yddA ←	ABC transporter family protein YddA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,577,704	0	C	T	100.0%	67.5 / NA	21	V547I (GTC→ATC)	yddA	ABC transporter family protein YddA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/8);  total (13/8)

AGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATA  >  NC_000913/1577578‑1577820                                                                                                                               |                                                                                                                     aGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGATACCggg                                                                                                                <  1:307233/133‑1 (MQ=255)   cgcCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGATACCGGGtt                                                                                                              <  2:332605/133‑1 (MQ=255)         aTAAAATGCCGCCAGTAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGATACCGGGTTGAtgtg                                                                                                        <  1:187018/133‑1 (MQ=255)                  cGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGATACCGGGTTGATGTGTAACCATaa                                                                                               >  2:31641/1‑133 (MQ=255)                                         atttatttATATCTCCGTACAGCATGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGATACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGtttt                                                                        >  2:347289/1‑133 (MQ=255)                                             atttatATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGATACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCAc                                                                    <  1:347289/133‑1 (MQ=255)                                                tataTCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGATCCAGATACCGGGTTGATTTGTAACCACAATGACGTCGATTGTGGGTAGGTTTTCACGCa                                                                 >  2:202701/1‑133 (MQ=255)                                                   atCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGATACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTa                                                              >  2:324692/1‑133 (MQ=255)                                                   atCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGATACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCAATa                                                              >  2:299230/1‑133 (MQ=255)                                                        cGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGATACCGGGTTGATGTGTAACCTTAATGACGACGCATGTTGGTAGTTTATCACTCACTAAACGc                                                         >  2:179342/1‑133 (MQ=255)                                                         gTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGATACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCa                                                        >  1:9876/1‑133 (MQ=255)                                                               cagGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGATACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGc                                                  >  2:27004/1‑133 (MQ=255)                                                                    cTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGATACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATa                                             <  1:324692/133‑1 (MQ=255)                                                                    cTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGATACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGATTTTCACACACTAATCGCAGTAAGCGGATa                                             >  1:178632/1‑133 (MQ=255)                                                                              ttttATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGATACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGtt                                   >  2:128459/1‑133 (MQ=255)                                                                                  ataAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGATACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTc                               <  1:128459/133‑1 (MQ=255)                                                                                       aCCGCGCTTATGTCACAACTATAATCGGCCAGGTTCCAGATACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGAt                          <  1:179342/133‑1 (MQ=255)                                                                                       aCCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGATACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTc                                  <  1:54996/125‑1 (MQ=255)                                                                                       aCCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGATACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTc                                  >  2:54996/1‑125 (MQ=255)                                                                                                       aaaTATCATCCGCCACGTTCCAGAAACCAGGTAGATGTGTAACCATAATGACGCCGCCTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTc          <  2:178632/133‑1 (MQ=255)                                                                                                        aaTATCATCGGCCAGGTTCCAGATACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCCTGTTCCTAAAGTTGAGAGGTAGTTTCGTCa         >  1:170505/1‑133 (MQ=255)                                                                                                            tcatcGGCCAGGTTCCAGATACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGaa     >  2:458669/1‑133 (MQ=255)                                                                                                               tcGGCCAGGTTCCAGATACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAata  >  2:60753/1‑133 (MQ=255)                                                                                                                               |                                                                                                                     AGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATA  >  NC_000913/1577578‑1577820

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