Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	386,360	C→T	P52S (CCT→TCT)	tauB →	taurine ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	386,360	0	C	T	100.0%	152.6 / NA	44	P52S (CCT→TCT)	tauB	taurine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/23);  total (21/23)

TTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGC  >  NC_000913/386226‑386488                                                                                                                                       |                                                                                                                                 tttACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTAt                                                                                                                             <  1:436791/140‑1 (MQ=255)   tACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCa                                                                                                                           >  2:529444/1‑140 (MQ=255)         gATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATgg                                                                                                                     <  1:225807/140‑1 (MQ=255)          aTTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGgc                                                                                                                    >  1:258715/1‑140 (MQ=255)          aTTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACCGGTGGTTCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGgc                                                                                                                    >  2:753382/1‑140 (MQ=255)                    cAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAAc                                                                                                          >  2:307500/1‑140 (MQ=255)                    cAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAAc                                                                                                          >  1:627649/1‑140 (MQ=255)                           gCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCggg                                                                                                   >  1:100436/1‑140 (MQ=255)                               tGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAg                                                                                               <  1:529444/140‑1 (MQ=255)                                             cTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGatc                                                                                 <  1:190100/140‑3 (MQ=255)                                              tGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCg                                                                                <  2:627649/140‑1 (MQ=255)                                               gACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCgg                                                                               >  2:444687/1‑140 (MQ=255)                                                  gCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAg                                                                            <  1:73474/140‑1 (MQ=255)                                                    tGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCa                                                                          <  1:546543/140‑1 (MQ=255)                                                      gAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCaga                                                                        >  2:671490/1‑140 (MQ=255)                                                                   tACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGtt                                                           >  1:745804/1‑140 (MQ=255)                                                                           gtgCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGaa                                                   >  2:762514/1‑140 (MQ=255)                                                                              cTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGa                                                <  1:16747/140‑1 (MQ=255)                                                                              cTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGa                                                <  2:258715/140‑1 (MQ=255)                                                                                ggggCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAg                                              <  2:663003/140‑1 (MQ=255)                                                                                ggggCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAg                                              <  2:783619/140‑1 (MQ=255)                                                                                  ggCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAggg                                            >  1:257/1‑140 (MQ=255)                                                                                  ggCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAggg                                            >  2:473641/1‑140 (MQ=255)                                                                                   gCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGc                                           <  1:307500/140‑1 (MQ=255)                                                                                      gTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGctac                                        >  1:560583/1‑140 (MQ=255)                                                                                      gTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAAGGctac                                        >  1:582358/1‑140 (MQ=255)                                                                                       tCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGa                                                >  2:510208/1‑131 (MQ=255)                                                                                       tCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGa                                                <  1:510208/131‑1 (MQ=255)                                                                                       tCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGa                                                <  2:177835/131‑1 (MQ=255)                                                                                       tCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGa                                                >  1:177835/1‑131 (MQ=255)                                                                                       tCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGctact                                       >  1:629583/1‑140 (MQ=255)                                                                                             tGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTg                                 <  1:202300/140‑1 (MQ=255)                                                                                                  tAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCa                            <  2:100436/140‑1 (MQ=255)                                                                                                  tAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCa                            <  2:372060/140‑1 (MQ=255)                                                                                                    aaaCCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAAt                          <  1:444687/140‑1 (MQ=255)                                                                                                     aaCCACCCTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATg                         <  1:762514/140‑1 (MQ=255)                                                                                                          cccTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCa                            >  1:814906/1‑132 (MQ=255)                                                                                                          cccTGCTGAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCa                            <  2:814906/132‑1 (MQ=255)                                                                                                             tgctgAATCTGATTGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGAc                 <  1:1733/140‑1 (MQ=255)                                                                                                                         ttGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTc     >  1:647556/1‑140 (MQ=255)                                                                                                                         ttGCCGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGGGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTc     >  1:416491/1‑140 (MQ=255)                                                                                                                            ccGGTTTTGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGc  <  1:473641/140‑1 (MQ=255)                                                                                                                                  ttGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGt       >  2:121171/1‑129 (MQ=255)                                                                                                                                  ttGTGTCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGt       <  1:121171/129‑1 (MQ=255)                                                                                                                                       |                                                                                                                                 TTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGC  >  NC_000913/386226‑386488

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: ‑