Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	387,209	T→C	L80P (CTG→CCG)	tauC →	taurine ABC transporter membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	387,209	0	T	C	100.0%	141.0 / NA	41	L80P (CTG→CCG)	tauC	taurine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/23);  total (18/23)

GCCGCAACAGGTACTGGAAAAACTACTCACCATTGCCGGACCGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCTGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTATTTGCAGCGGTGTTGTTCGGTATTCCGGTCGGGATC  >  NC_000913/387132‑387279                                                                              |                                                                       gccgcAACAGGTACTGGAAAAACTACTCACCATTGCCGGACCGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGc                                                                     <  2:87817/81‑1 (MQ=255)    gcAACAGGTACTGGAAAAACTACTCACCATTGCCGGACCGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGc                                                                  >  1:87867/1‑81 (MQ=255)      aaCAGGTACTGGAAAAACTACTCACCATTGCCGGACCGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCa                                                                <  1:87818/81‑1 (MQ=255)        cAGGTACTGGAAAAACTACTCACCATTGCCGGACCGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGt                                                              >  2:87841/1‑81 (MQ=255)          ggTACTGGAAAAACTACTCACCATTGCCGGACCGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCt                                                            <  1:87819/81‑1 (MQ=255)              cTGGAAAAACTACTCACCATTGCCGGACCGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGAcg                                                        <  2:87821/81‑1 (MQ=255)              cTGGAAAAACTACTCACCATTGCCGGACCGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGAcg                                                        >  1:87862/1‑81 (MQ=255)              cTGGAAAAACTACTCACCATTGCCGGACCGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGAc                                                         <  1:87820/80‑1 (MQ=255)                    aaaCTACTCACCATTGCCGGACCGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCAtt                                                  <  2:87822/81‑1 (MQ=255)                    aaaCTACTCACCATTGCCGGACCGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCAtt                                                  <  1:87823/81‑1 (MQ=255)                      actactCACCATTGCCGGACCGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTAt                                                <  1:87824/81‑1 (MQ=255)                      actactCACCATTGCCGGACCGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTAt                                                <  1:87825/81‑1 (MQ=255)                       ctactCACCATTGCCGGACCGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATg                                               >  1:87865/1‑81 (MQ=255)                            cACCATTGCCGGACCGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGgc                                          >  2:87872/1‑81 (MQ=255)                               cATTGCCGGACCGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCt                                       <  1:87826/81‑1 (MQ=255)                               cATTGCCGGACCGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCt                                       <  1:87827/81‑1 (MQ=255)                                aTTGCCGGACCGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTa                                      >  1:386941/1‑81 (MQ=255)                                   gCCGGACCGCAAGTCTTTTTGGACGCCACGCTGTCGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCCCTAttt                                   >  2:87871/1‑81 (MQ=255)                                   gCCGGACCGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTAttt                                   >  2:87875/1‑81 (MQ=255)                                    ccGGACCGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTAttt                                   <  2:87828/80‑1 (MQ=255)                                     cGGACCGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTATTTGc                                 <  2:87829/81‑1 (MQ=255)                                         ccGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTATTTGCAGCg                             >  2:87858/1‑81 (MQ=255)                                           gCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTATTTGCAGCGgt                           <  1:87830/81‑1 (MQ=255)                                           gCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTATTTGCAGCGgt                           >  2:87870/1‑81 (MQ=255)                                              aGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTATTTGCAGCGGtgtt                        <  2:87831/81‑1 (MQ=255)                                              aGGCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTATTTGCAGCGGtgtt                        <  2:87832/81‑1 (MQ=255)                                                gCTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCTCATTATGCTGGCGCTATTTGCAGCGGtgttgt                      >  1:87861/1‑81 (MQ=255)                                                 cTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTATTTGCAGCGGTgttgtt                     <  1:87833/81‑1 (MQ=255)                                                 cTTTATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTATTTGCAGCGGTgttgtt                     <  1:87834/81‑1 (MQ=255)                                                   ttATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTATTTGCAGCGGTGTTGTTCg                   >  1:87868/1‑81 (MQ=255)                                                    tATGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTATTTGCAGCGGTGTTGTTc                    <  2:87835/79‑1 (MQ=255)                                                      tGGACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTATTTGCAGCGGTGTTGTTCGGTa                >  2:87901/1‑81 (MQ=255)                                                        gACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTATTTGCAGCGGTGTTGTTCGGTAtt              >  1:87876/1‑81 (MQ=255)                                                        gACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTATTTGCAGCGGTGTTGTTCGGTAtt              >  2:87880/1‑81 (MQ=255)                                                        gACGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTATTTGCAGCGGTGTTGTTCGGTAt               <  2:87836/80‑1 (MQ=255)                                                         aCGCCACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTATTTGCAGCGGTGTTGTTCGGTATTc             <  1:87837/81‑1 (MQ=255)                                                            ccACGCTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTATTTGCAGCGGTGTTGTTCGGTATTCCgg          <  2:87838/81‑1 (MQ=255)                                                                 cTGTGGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTATTTGCAGCGGTGTTGTTCGGTATTCCGGTCgg      >  2:87877/1‑80 (MQ=255)                                                                   gtgGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTATTTGCAGCGGTGTTGTTCGGTATTCCGGTCGGGAt   >  1:87887/1‑81 (MQ=255)                                                                   gtgGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTATTTGCAGCGGTGTTGTTCGGTATTCCGGTCGGGAt   >  1:87878/1‑81 (MQ=255)                                                                    tgGCAGCATCCGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTATTTGCAGCGGTGTTGTTCGGTATTCCGGTCGGGATc  <  1:87839/81‑1 (MQ=255)                                                                              |                                                                       GCCGCAACAGGTACTGGAAAAACTACTCACCATTGCCGGACCGCAAGGCTTTATGGACGCCACGCTGTGGCAGCATCTGGCAGCCAGTCTGACGCGCATTATGCTGGCGCTATTTGCAGCGGTGTTGTTCGGTATTCCGGTCGGGATC  >  NC_000913/387132‑387279

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