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%	148.8 / NA	42	D169G (GAT→GGT)	tauA	taurine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/21);  total (21/21)

GCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTT  >  NC_000913/385607‑385870                                                                                                                                   |                                                                                                                                      gCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATg                                                                                                                               <  2:1561233/139‑1 (MQ=255) gCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATg                                                                                                                               <  2:502171/139‑1 (MQ=255)          tttATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATg                                                                                                                      >  2:390807/1‑139 (MQ=255)           ttATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGt                                                                                                                     <  2:79763/139‑1 (MQ=255)             aTCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCt                                                                                                                   >  2:471047/1‑139 (MQ=255)                  caccacCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCa                                                                                                              >  2:533509/1‑139 (MQ=255)                             aCAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTaa                                                                                                   <  2:1297903/139‑1 (MQ=255)                                gCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGc                                                                                                <  1:763042/139‑1 (MQ=255)                                   tgctgGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTa                                                                                                                        >  2:203343/1‑112 (MQ=255)                                   tgctgGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTa                                                                                                                        <  1:203343/112‑1 (MQ=255)                                      tgGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGa                                                                                          >  1:1404788/1‑139 (MQ=255)                                                      cTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGaaaaa                                                                                      <  1:1640153/127‑1 (MQ=255)                                                      cTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGaaaaa                                                                                      >  2:1640153/1‑127 (MQ=255)                                                         gggCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATg                                                                                                                               <  1:345184/83‑1 (MQ=255)                                                         gggCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATg                                                                                                                               >  2:345184/1‑83 (MQ=255)                                                           gCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                                     <  1:533509/139‑1 (MQ=255)                                                                       ggCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAAc                                                                                                  >  2:961814/1‑98 (MQ=255)                                                                       ggCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAAc                                                                                                  <  1:961814/98‑1 (MQ=255)                                                                           aGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTc                                                     >  2:1476561/1‑139 (MQ=255)                                                                           aGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTc                                                     <  1:471047/139‑1 (MQ=255)                                                                                  aTTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGt                                              >  2:429111/1‑139 (MQ=255)                                                                                  aTTGAGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGt                                              >  1:971584/1‑139 (MQ=255)                                                                                      tgAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGaaaaa                                                                                      <  1:754736/95‑1 (MQ=255)                                                                                      tgAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGaaaaa                                                                                      >  2:754736/1‑95 (MQ=255)                                                                                       gAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCt                                                                                                                          <  2:1071441/58‑1 (MQ=255)                                                                                       gAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCt                                                                                                                          >  1:1071441/1‑58 (MQ=255)                                                                                       gAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGccc                                                                                              <  2:451333/86‑1 (MQ=255)                                                                                       gAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGccc                                                                                              >  1:451333/1‑86 (MQ=255)                                                                                         aCCTGCAGCCGCCCGCGATTGTCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                                     >  2:1272032/1‑109 (MQ=255)                                                                                         aCCTGCAGCCGCCCGCGATTGTCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                                     <  1:1272032/109‑1 (MQ=255)                                                                                         aCCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGgcgc                                       >  2:226689/1‑139 (MQ=255)                                                                                          ccTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGgcgcg                                      >  2:161162/1‑139 (MQ=255)                                                                                             gCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGaa                                                                                         >  2:1215365/1‑85 (MQ=255)                                                                                             gCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGaa                                                                                         <  1:1215365/85‑1 (MQ=255)                                                                                                  cgcccgcGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGaaaa                                                                                       <  2:1202328/82‑1 (MQ=255)                                                                                                  cgcccgcGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGaaaa                                                                                       >  1:1202328/1‑82 (MQ=255)                                                                                                     ccgcGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGa                           <  1:429111/139‑1 (MQ=255)                                                                                                     ccgcGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGa                           <  1:226689/139‑1 (MQ=255)                                                                                                                         cAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                                     >  1:834291/1‑77 (MQ=255)                                                                                                                         cAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                                     <  2:834291/77‑1 (MQ=255)                                                                                                                            cGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGAtt    <  1:1476561/139‑1 (MQ=255)                                                                                                                              gggAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGAtttt  >  2:928030/1‑139 (MQ=255)                                                                                                                                   |                                                                                                                                      GCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTT  >  NC_000913/385607‑385870

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