Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	385,737	A→G	D169G (GAT→GGT)	tauA →	taurine ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	385,737	0	A	G	100.0%	78.3 / NA	23	D169G (GAT→GGT)	tauA	taurine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/9);  total (14/9)

ATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGC  >  NC_000913/385619‑385859                                                                                                                       |                                                                                                                           aTCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCt                                                                                                        >  1:498378/1‑139 (MQ=255)       accaccCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCAc                                                                                                  >  2:36100/1‑139 (MQ=255)        ccacccacTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCAcc                                                                                                 >  1:652269/1‑139 (MQ=255)                         ctgGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCt                                                                                                               <  1:260297/108‑1 (MQ=255)                         ctgGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCt                                                                                                               >  2:260297/1‑108 (MQ=255)                              ggcACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGaaaaa                                                                           >  2:525612/1‑139 (MQ=255)                                   tGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACgg                                                                      >  1:317860/1‑139 (MQ=255)                                               gCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                          <  1:525612/139‑1 (MQ=255)                                                        ccGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGa                                                 <  2:743384/139‑1 (MQ=255)                                                          gggCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATg                                                                                                                    >  2:81422/1‑70 (MQ=255)                                                          gggCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATg                                                                                                                    <  1:81422/70‑1 (MQ=255)                                                                        tgtgAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAg                                    >  1:318539/1‑136 (MQ=255)                                                                        tgtgAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAg                                    <  2:318539/136‑1 (MQ=255)                                                                        tgtgAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTgg                                 >  1:460158/1‑139 (MQ=255)                                                                          tgAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTgggg                               >  1:711004/1‑139 (MQ=255)                                                                               cTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGgcgcgc                          >  1:271660/1‑139 (MQ=255)                                                                                tGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAg                                                                  >  2:636792/1‑98 (MQ=255)                                                                                tGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAg                                                                  <  1:636792/98‑1 (MQ=255)                                                                                  cAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCaa                       <  2:652269/139‑1 (MQ=255)                                                                                      cgcccgcGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCt                   >  2:798715/1‑139 (MQ=255)                                                                                       gcccgcGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTg                  <  1:798715/139‑1 (MQ=255)                                                                                         ccgcGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGa                <  1:36100/139‑1 (MQ=255)                                                                                                       gCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTgc  >  2:462267/1‑139 (MQ=255)                                                                                                                       |                                                                                                                           ATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGC  >  NC_000913/385619‑385859

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