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

ATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATT  >  NC_000913/385619‑385868                                                                                                                       |                                                                                                                                    aTCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCt                                                                                                                 >  1:61363/1‑139 (MQ=255)        ccacccacTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCAcc                                                                                                          >  1:456108/1‑139 (MQ=255)         cacccacTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACcg                                                                                                         >  1:596047/1‑139 (MQ=255)              acTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGt                                                                                                    <  2:456108/139‑1 (MQ=255)                tACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTa                                                                                                  <  1:333049/139‑1 (MQ=255)                     ccTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCa                                                                                                            <  1:1100260/124‑1 (MQ=255)                     ccTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCa                                                                                                            >  2:1100260/1‑124 (MQ=255)                     ccTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGcc                                                                                             <  2:239712/139‑1 (MQ=255)                        gctgGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGa                                                                                                                               <  2:905523/102‑1 (MQ=255)                        gctgGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGa                                                                                                                               >  1:905523/1‑102 (MQ=255)                         ctgGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTgg                                                                                         >  1:258525/1‑139 (MQ=255)                          tgGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGAt                                                                                                                              >  2:448344/1‑101 (MQ=255)                          tgGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGAt                                                                                                                              <  1:448344/101‑1 (MQ=255)                                 aCTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGAc                                                                                 >  2:39842/1‑139 (MQ=255)                                          cTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGgtg                                                                        <  1:39842/139‑1 (MQ=255)                                               gCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                                   <  2:596047/139‑1 (MQ=255)                                               gCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                                   <  1:10140/139‑1 (MQ=255)                                               gCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                                   <  2:1060747/139‑1 (MQ=255)                                               gCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                                   <  2:61363/139‑1 (MQ=255)                                               gCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                                   >  1:238281/1‑139 (MQ=255)                                                         cGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                                   >  2:24246/1‑129 (MQ=255)                                                         cGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                                   <  1:24246/129‑1 (MQ=255)                                                            gCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                                   <  1:237860/126‑1 (MQ=255)                                                            gCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                                   >  2:237860/1‑126 (MQ=255)                                                                        tgtgAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTgg                                          >  2:342937/1‑139 (MQ=255)                                                                         gtgAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTggg                                         >  1:273641/1‑139 (MQ=255)                                                                          tgAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATg                                                                                                                             <  1:42445/54‑1 (MQ=255)                                                                          tgAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATg                                                                                                                             >  2:42445/1‑54 (MQ=255)                                                                          tgAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTgggg                                        <  2:273641/139‑1 (MQ=255)                                                                                tGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGa                                                                                        <  1:208078/85‑1 (MQ=255)                                                                                tGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGa                                                                                        >  2:208078/1‑85 (MQ=255)                                                                                         ccgcGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAgg                                                     >  1:435971/1‑111 (MQ=255)                                                                                         ccgcGATTATCGCTGCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAgg                                                     <  2:435971/111‑1 (MQ=255)                                                                                                       gCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCaa                                >  2:999808/1‑118 (MQ=255)                                                                                                       gCCTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCaa                                <  1:999808/118‑1 (MQ=255)                                                                                                        ccTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGc                                                                                              <  1:77617/55‑1 (MQ=255)                                                                                                        ccTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGc                                                                                              >  2:77617/1‑55 (MQ=255)                                                                                                        ccTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                                   >  2:170942/1‑82 (MQ=255)                                                                                                        ccTGGCAGCGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                                   <  1:170942/82‑1 (MQ=255)                                                                                                                cGGGGAGGTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGAtt  >  1:566915/1‑139 (MQ=255)                                                                                                                       |                                                                                                                                    ATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATT  >  NC_000913/385619‑385868

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