Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 386,844 T→A 11.4% L213Q (CTG→CAG)  tauB → taurine ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913386,8440TA11.4% 82.6 / 2.7 37L213Q (CTG→CAG) tauBtaurine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (16/15);  new base A (2/2);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01

AGATGCAAACCCTGCTGCTGAAACTCTGGCAGGAGACGGGCAAGCAGGTGCTGTTGATTACCCACGATATAGAAGAAGCGGTGTTTATGGCGACTGAACTGGTTCTGCTTTCATCCGGCCCTGGCCGTGTGCTGGAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGGAGAGTCGAGCCGCAGCATCAAGTCCGATCCACAATTCATCGCCATGCGCGAATATGTTTTAAGCCGCGTATTTGAGCAACGGGAGGCGTTCTCATGAGTGTG  >  NC_000913/386703‑386979
                                                                                                                                             |                                                                                                                                       
aGATCCAAGCCCTGCTGCTGAAACTCTGGCAGGAGACGGGCAAGCAGGTGCTGTTGATTACCCACGATATAGAAGAAGCGGTGTTTATGGCGACTGAACTGGTTCTGCTTTCATCCGGCCCTGGCCGTGTGCTGGAGCGGCTg                                                                                                                                        <  1:735328/143‑1 (MQ=255)
 gATGCAAACCCTGCTGCTGAAACTCTGGCAGGAGACGGGCAAGCAGGTGCTGTTGATTACCCACGATATAGAAGAAGCGGTGTTTATGGCGACTGAACTGGTTCTGCTTTCATCCGGCCCTGGCCGTGTGCTGGAGCGGCTGc                                                                                                                                       <  1:608461/143‑1 (MQ=255)
      aaaCCCTGCTGCTGAAACTCTGGCAGGAGACGGGCAAGCAGGTGCTGTTGATTACCCACGATATAGAGGAGGCGGCGTTTATGGCGCCTGAACCTGCTCTGCGGTCGCCCGGCGCTGGCCGTGTGGTCGGGCGGCCCCCGCTc                                                                                                                                  >  1:76666/1‑143 (MQ=255)
       aaCCCTGCTGCTGAAGCTCTGGCAGGAGACGGGCAAGCAGGTGCTGTTGACTACCCACGATATAGAAGAAGCGGTGTTTATGGCGACTGAACTGGTTCTGCTTTCATCCGGCCCTGGCCGTGTGCTGGAGCGGCAGCCGCTcc                                                                                                                                 <  2:774962/143‑2 (MQ=255)
       aaCCCTGCTGCTGAAACTCTGGCAGGAGCCGGGCAAGCCGCCGCTGTTGATTACCCCCGAGATAGGAGAACCGGTGCTTCTGGCGACTCGACTGGCACCGGGGTCATCCGGCCCTTGCCGTGTGCTTGTGCACCTGCCGCTcc                                                                                                                                 >  2:548292/1‑142 (MQ=255)
         tccTGCTGCTGAGCGGCTGTGAGGTGTCGGGCACGCGGGTGCTGTTGATTACCTACGATATAGAAGAAGCGGTGTTTATGGCGACTGACCTGGTTCTGCTTTCATCCGGCCCTGGCCGGGTGCTGGAGCGGCAGCCGCTccag                                                                                                                               <  2:687936/142‑4 (MQ=255)
                gcggAAACTCTGGCAGGAGACGGGCAAGCAGGTGCTGTTGATTACCCACGATATAGAAGAAGCGGTGTTTATGGCGACTGAACTGGTTCTGCTTTCATCCGGCCCTGGCCGTGTGCTGGAGCGGCTGCCGCTCAACTTTGCTc                                                                                                                        <  1:247333/140‑1 (MQ=255)
                  tgAAACTCTGGCAGGAGACGGGCAAGCAGGTGCTGTTGATTACCCACGATATAGAAGACGCGGTGTTTATGGCGACTGAACTCGTTCCGCCCCCATCCGGCCCGGGCCGCGTGCTGGCGCGGCAGCCGCTCACCCtcgctcgc                                                                                                                      >  1:579151/1‑143 (MQ=255)
                       cacgggCAGGACACGGGCAAGCAGGTGCTGTTGATTACCCGCGATATAGAAGATGCGGTGTTTATGGCGACTGAACTGGTTCTGCTTTCATCCGGCCCTGGCCGTGTGCTGGAGCGGCTGCCGCTCAACTTTGCTCGCCGCtt                                                                                                                 <  2:490422/139‑1 (MQ=255)
                               ggggaCGGGCACGCAGGGTCTGTTGATGACCCACGATATAGAAGAAGCGGTGTTTATGGCGACTGAACTGGTTCTGCTTTCATCCGGCCCTGGCCGTGTGCTGGAGAGGCTGCCGCTCAACTGTGCTCGACGCTTTGTTGCgg                                                                                                         <  2:112337/140‑1 (MQ=255)
                                gagaCGGGCCAGCAGGTGCTGTTGATTACCCACGATATAGAACCAGCGGTGTTTCTGGCGACTGAACTCGCCCTGCTTTCACCCGGCCCTGGCCGCCTGCTGGCGCGGCTCCCCGACAACTATGCTCGCCGCTTCGTTGCgag                                                                                                        >  2:288550/1‑141 (MQ=255)
                                 ggaCGGGCGAGCAGGTGCTGTTGATTACCCACGATATAGAAGAAGCGGTGTTTATGGCGACTGAACTGGTTCTGCTTTCATCCGGCCCTGGCCGTGTGCTGGAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGga                                                                                                       <  2:177325/142‑1 (MQ=255)
                                     gggCAAGCAGGTGCTGTTGATTACCCACGATATAGAAGAATCGGTGTTTATGGCGACTGAACTGGTTATGCTTTCGTCCGGCCCTGGCGGTGTGCTGCAGCGGCTGCCGCTCACCGTTGCTCGCCGATTGGTTGCGGgggact                                                                                                   >  1:477176/1‑141 (MQ=255)
                                             cggcgCTCTTGATTACCCACGACATAGAAGAAGCGGTGTTTATGGCGACGGAACTGGTTCTGCTTTCATCCGGCCCTGGCCGTGTGCTGGAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGGAGAGTCGAGCCgc                                                                                           <  1:410607/139‑1 (MQ=255)
                                                  cTGTTGATTACCCACGATATAGAAGAAGCGGTGTTTATGGCGACTGAACTGGTTCTGCTTTCATCCGGCCCTGGCCGTGTGCTGGAGCGGCTGCCGCTCCACTTTGCTCGCCGCTTTGTGGGGGGAGAGGCGAGCCGCAGCAt                                                                                      >  1:342893/1‑143 (MQ=255)
                                                  cTGTTGATTACCCACGATATAGAAGAAGCGGTGTTTATGGCGACTGAACTGGTTCTGCTTTCATCCGGCCCTGGCCGTGTGCTGGAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGGAGCGGCGCGCCGCGGCAt                                                                                      >  2:526984/1‑143 (MQ=255)
                                                  cTGTGGATTACCCACGATATAGAAGAAGCGGTGTTTATGGCGACTGAACTGGTTCTGCTTTCATCCGGCCCTGGCCGTGTGCTGGAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGGAGAGTCGAGCCGCAGCAt                                                                                      <  1:543312/143‑1 (MQ=255)
                                                     ttGATTACCCACGATATAGAAGAAGCGGTGTTTATGGCGACTGAACTGGTTCTGCTTTCATCCGGCCCTGGCCGTGTGCTGGAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGGAGAGTCGAGCCgcagcggcca                                                                                   >  2:267103/1‑138 (MQ=255)
                                                       gactACCCACGATATCGAAGAAGCGGTGTTTATGGCGCCTGACCTGGTTCTGCTTTCATCCGGCCCTCGCCGTGCACTGCAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGGAGAGTCGAGCCGCAGCAGCaaat                                                                                 <  1:152740/140‑3 (MQ=255)
                                                              cACGATATAGAAGAAGCGGTGTTTATGGCGACTGAACTGGTTCTGCTTTCATCCGGCCCTGGCCGTGTGCTGGAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGGCGAGTCGAGCCGCAGCATCCAGTCCGCTcc                                                                          >  1:472902/1‑143 (MQ=255)
                                                                cGATATAGAAGAAGCGGTGTTTATGGCGACTGAACTGGTTCTGCTTTCATCCGGCCCTGGCCGTGTGCTGGAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGGAGAGTCGAGCCGCAGCATCAAGTCCGATCcac                                                                        <  1:267103/143‑1 (MQ=255)
                                                                  atatAGAAGAAGCGGGGTTTAGGGCGACTGAACTGGGTCGGCGTTCAGCCGGCCCTGGCCGTGCGCTGGAGCGGCGGCCGCTCAACTTTGCGCGCCGCTTTGTTGCGGGAGAGTCGAGCCGCAGCATCAAGTCCGATCCACaa                                                                      <  1:82425/143‑1 (MQ=255)
                                                                       gaagaaGCGGTGTTTATGGCGACTGAACTGGTTCTGCTTTCATCCGGCCCTGGCCGTGTGCTGGAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGGAGAGTCGAGCCGCAGCCTCAAGTACGATCCACGATTGAt                                                                 >  1:212636/1‑143 (MQ=255)
                                                                                gtgtgtATGGCGGCTGAACTGGTTCTGCTTTCATCCGGCCCTGGCCGTGTGCTTGAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGGAGAGTCGAGCCGCAGCATCAAGTCCGATCCACAATTCATCGCCATgcg                                                        <  2:72585/143‑1 (MQ=255)
                                                                                    ttATGGCGACTGAACTGGTTCTGCTTTCATCCGGCCCTGGCCGTGTGCTGGAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGGAGAGTCGAGCCGCAGCATCAAGTCCGATCCACAATTCATCGCCATGCGCGaa                                                    >  2:13313/1‑143 (MQ=255)
                                                                                          cGACTGAACGGGTTCTGCTTTCATCCGGCCCTGGCCGTGTGCTGGAGCGGCTGCCGCTCAACCTTGCTCGCCGCTTTGTTGCGGGAGAGTCGAGCCGCAGCATCAAGTCCGATCCACAATTCATCGCCATGCGCGAATATGtt                                              <  2:579151/143‑1 (MQ=255)
                                                                                           gccTGCACTGGTTCTGCTTTCTTCCGGCCGTGGCCGTGTGCTGGAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGGAGAGTCGAGCCGCAGCATCAAGTCCGATCCACAATTCATCGCCATGCGCGAATATGttt                                             <  2:325294/141‑1 (MQ=255)
                                                                                                aaCTGGTTCTGCTTTCATCCGGCCCTGGCCGTGTGCTGGAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGGAGAGTCGAGCCGCAGCATCAAGTACGATCCACAATTCCTCGCCATGCGCGAATATGTTTTAAGc                                        >  1:771548/1‑143 (MQ=255)
                                                                                                   gtgTTCTGCGTACATCCGGCCCTGGCCGTGTGCTGGAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGGAGAGTCGAGCCGCAGCATCAAGTCCGATCCACAATTCATCGCCATGCGCGAATATGTTTTAAGCcgc                                     <  1:646650/141‑1 (MQ=255)
                                                                                                     gTTCTGCTTTCATCCGGCCCTGGCCGTGTGCTGGAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGGCGAGTCGAGCCGCAGCATCAAGTCCGATCCACAATTCATCGCCATGCGCGAACCGGTTTTAAGCcgcgg                                   >  1:516856/1‑142 (MQ=255)
                                                                                                        cTGCTTTCATCCGGCCCTGGCCGTGTGCTGGAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGGAGAGTCGAGCCGCAGCATCAAGTCCGATCCACAATTTATCGCCATGCGCGAATATGGTTTAAGCGGCGTGtt                                >  1:119335/1‑143 (MQ=255)
                                                                                                            tttCATCCGGCCCTGGCCGTGTGCTGGAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGGAGAGTCGAGCCGCAGCATCAAGTCCGATCCACAATTCATCGCCATGCGCGAATATGTTTTAAGCCGCGGATTTGAg                            >  1:48315/1‑143 (MQ=255)
                                                                                                                  ccGGCCCTGGCCGTGTGCTGGAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGGAGAGTCGAGCCGCAGCATCAAGTCCGATCCACAATTCATCGCCATGCGCGAATATGTTTTAAGCCGCGCATCTGGGCAACgg                      >  2:83500/1‑143 (MQ=255)
                                                                                                                    ggCCCTGGCCGTGTGCTGGAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGGAGAGTAGAGCCGCAGCATCAAGTCCGATCCACAATTCATCGCCATGCGCGAATATGTTTTAAGCCGCGTATTTGAGCAACGGGa                    >  1:351802/1‑143 (MQ=255)
                                                                                                                             cGTGTGCGGGAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTCCGGGAGAGTCGAGCCGCAGCATCAAGTCCGATCCACAATTCATCGCCATGCGCGAATATGTTTTAAGCCGCGTATTTGAGCAACGGGAGGCGTtctc           <  1:354590/143‑1 (MQ=255)
                                                                                                                                    ttgAGCGGCAGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGGAGAGTCGAGCCGCAGCATCAAGTCCGATCCACAATTCATCGCCATGCGCGAATATGTTTTAAGCCGCGTATTTGAGCAACGGGAGGCGTTCTCATGCgtg    >  1:155787/3‑143 (MQ=255)
                                                                                                                                      gAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGGAGAGTCGAGCCGCAGCATCAAGAACGAGGGACAATTCATCGCCATGCGCGACTATGTTTTAAGCCGCGTATTTGAGCTACGCGAGGCGTTCTCATGAGtgtg  >  1:110758/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                       
AGATGCAAACCCTGCTGCTGAAACTCTGGCAGGAGACGGGCAAGCAGGTGCTGTTGATTACCCACGATATAGAAGAAGCGGTGTTTATGGCGACTGAACTGGTTCTGCTTTCATCCGGCCCTGGCCGTGTGCTGGAGCGGCTGCCGCTCAACTTTGCTCGCCGCTTTGTTGCGGGAGAGTCGAGCCGCAGCATCAAGTCCGATCCACAATTCATCGCCATGCGCGAATATGTTTTAAGCCGCGTATTTGAGCAACGGGAGGCGTTCTCATGAGTGTG  >  NC_000913/386703‑386979

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: