Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 75,091 G→T 37.4% D130E (GAC→GAA)  thiB ← thiamine/thiamine pyrophosphate/thiamine monophosphate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091375,0910GT37.4% 25.8 / 3.4 24D130E (GAC→GAA) thiBthiamine/thiamine pyrophosphate/thiamine monophosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (3/12);  new base T (7/2);  total (10/14)
Fisher's exact test for biased strand distribution p-value = 1.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAG‑CCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGG‑CAAACAGTCCGGTTTTACTG‑GCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTGAGAAGCGAAACGCCATCTTCCAGCGCCACCAGTTTCAGTT  >  NC_000913/74875‑75329
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
gCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGTC‑GGGGTTTTACGTTTTTTTTTTTCTT‑TAACGaaatgcaaggag                                                                                                                                                                                                                              >  1:751377/1‑228 (MQ=255)
gCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTTGTTAT‑AAAAAAAgggaaagaag                                                                                                                                                                                                                              >  2:244461/1‑231 (MQ=255)
gCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGTTTTTTCTTTTTCT‑TAACAAAAGCGCAGTAg                                                                                                                                                                                                                              >  1:68849/1‑238 (MQ=255)
        cTTCCCCCCAGTTTTTGTTGCCCGTGCCCGTTTTTTCCGCCATTTTCGCCCGGGTTGGGGGGGCGTCTTCGCCATAGATTTTTTGCACCAATAGCAACAGCCCCAGCCCGGGTGTACTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCa                                                                                                                                                                                                                     <  2:689352/239‑1 (MQ=255)
         ttCGCTCCAGCTTTGGGTTACTGGCACCGTTTTCGTCGCCAGTTGCTGCCGGGCTGGGGGGGCGTCATGACCAAAGACTTTTTGCATCCATAGCAACAGCCCCATCCCCGGTTTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCa                                                                                                                                                                                                                     <  1:62529/238‑1 (MQ=255)
          tCGCTCCCGCCTTGGGGCACGTTGCCCTTTTTCTTCGCCAGTTGCTGCAGGGCTTGGGGGGCGTCATCGCCTTAGACTTTTTGCACCCATAGCAACAGCCCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCaaa                                                                                                                                                                                                                   <  1:377130/239‑1 (MQ=255)
          tCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGCC‑GGGGTTTTCCTTTTTTTATTTTTAT‑AAACAAAAGCGAAGTGGtcatcatcaca                                                                                                                                                                                                                   >  1:321225/1‑237 (MQ=255)
                        ggTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGCGTTTTTAGTTTTTTTTTTTCAT‑ATACAAAAAGGCAGTAGTCAATATTACACTGTACGaaacggg                                                                                                                                                                                                     >  1:20953/1‑235 (MQ=255)
                             ccgtgaccgtTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTCAGTTTGTTTTTTTCTT‑AAAACAAGGCGGAGTAGCCCTAATCAACAGGGAAGAAAAAGTCatta                                                                                                                                                                                                >  1:823010/1‑239 (MQ=255)
                                                                     ggCGTCAGCGCCATAGCCTTTTGGCATCCATAGCAACAGCCCCAGCCCCGGTGTACGGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAG‑CCGCCGGGAACGTTAACGGCATCCGCTGCcaca                                                                                                                                                        <  1:244461/238‑1 (MQ=255)
                                                                                                   tAGCAACAGCCCCGGCCCGGGTGTACTGGTGCGCGGATCCTGATAAACCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAG‑CCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGG‑CAAACAGTCCGGTTTTACTg‑g                                                                                                                       <  1:743231/239‑1 (MQ=255)
                                                                                                       acacgACCCAGCCCGGGTGTATGGGTGCCCGGATCCTGATAAACCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAG‑CCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGG‑CAAACAGTCCGGTTTTACTg‑gcggc                                                                                                                   <  1:264522/235‑1 (MQ=255)
                                                                                                            aCCCAGCCCCGGTGTACTGGTGCTCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGATTTTTCAGTTTGTTCTTTTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAG‑CCGCCCGGAAAGTTAACGGCATCCGCTGCCACACCGCTTTTGG‑GCAACAGTCCGGTTTTACTT‑GAGGCGACTa                                                                                                              >  1:697796/1‑239 (MQ=255)
                                                                                                                 gCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTCAGTTTTTTTTTGTCTT‑AAACGAAGGCGAAGTAGCCCTAATCAAACGGTACGAAAATGTCATTATTCCCG‑CCGCCGGGGACGCTAACGGCATCCCCCGCCACACCCCTTTTTG‑CCAAAAACCCGGCTTTATTG‑GCGGCGGCTTACAgg                                                                                                         >  1:571791/1‑239 (MQ=255)
                                                                                                                     gggTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAG‑CCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGG‑CAAACAg                                                                                                                                      <  1:657677/206‑1 (MQ=255)
                                                                                                                               ggtgggggATTCTTATAAATCCCCCGCAAGTTTTGATCGCTCCCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAG‑CCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGG‑CAAACAGTCCGGTTTTACTG‑GCGGCGTCTAACAGGTTGTTATCCAGccc                                                                                           <  1:233420/233‑1 (MQ=255)
                                                                                                                                  cgcgGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGTTTTTTTTTTTCCTAAAACAAAGGGGAGGAGCCCATAATCAAACGTGAGGCAAG‑GTCATTTTTTCAGACCCCCGGGAAGGTAAACGCCA‑CCGCCTCCACCCCCCCTTTGGTCAAAAAGCCGGGTTTTATTGTGCCGCG‑CTAAAAGATTGTTATACAcccaa                                                                                          >  2:231752/1‑236 (MQ=255)
                                                                                                                                          cttaTAACTCCCCCTCCATTTTTGACCGCTCTCAACCAGTTTTTTCAGGCTTGGTGGC‑GGTTTTTTCAGTTTGTTTTTGTCAT‑AAACGAAGGCGAAGTAGCCAAAATCAAACGGTACGAAAGTGTCATTATTCCAG‑CCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGG‑CAAACAGTCCGGTTTTACTG‑GCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACAt                                                                                <  2:20953/236‑1 (MQ=255)
                                                                                                                                                                tGATCGCTCTCACCCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCAAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAG‑CCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGG‑CAAACAGTCCGGTTTTACTG‑GCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTc                                                          <  1:439629/239‑1 (MQ=255)
                                                                                                                                                                                              tttgGG‑GGTTTTTTCTGTTTTTTTTTTTCAT‑AACCAAAGGGGAAGTAGCAATAATCAAACGGTACGAAAGTGTCATTATTCCAG‑CCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGG‑CAAACAGTCCGGTTTTACTG‑GCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTGAGAAGCGAAACGcc                            <  2:823010/237‑1 (MQ=255)
                                                                                                                                                                                                        tttttCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAATTGTCATTATTCCAG‑CCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGG‑CAAACAGTCCGGTTTTACTG‑GCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTGAGAAGCGAAACGCCATCTTCCa                    <  1:574853/238‑1 (MQ=255)
                                                                                                                                                                                                           ttCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAG‑CCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGG‑CAAACAGTCCGGTTTTACTG‑GCGGCGTCTAACCGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCCTCCGTCGACGGTTGATAAGCGCAAAGCCATCTTCCAgcg                 >  2:753735/1‑238 (MQ=255)
                                                                                                                                                                                                            tCAGTTTGTTCTGTTAAA‑AAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTTTTCCAG‑CCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGG‑CAAACAGTCCGGTTTTACTG‑GCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTGAGAAGCGAAACGCCATCTTCCAGCGcc               <  1:231752/239‑1 (MQ=255)
                                                                                                                                                                                                                    tttttGTCAT‑AAAGTAAGGCAAGTTAGGTATAATCAAACGGTACGAAATTTGCATTATCCCAG‑CCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGG‑CAAACAGTCCGGTTTTACTG‑GCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTGAGAAGCGAAACGCCATCTTCCAGCGCCACCAGttt       <  2:279749/236‑1 (MQ=255)
                                                                                                                                                                                                                         gTCAT‑GAACGAAGGCGAAGTGGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAG‑CCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGG‑CAAACAGTCCGGTTTTACTG‑GCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTGAGAAGCGAAACGCCATCTTCCAGCGCCACCAGTTTCAGtt  <  1:594890/239‑1 (MQ=255)
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
GCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAT‑AAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAG‑CCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGG‑CAAACAGTCCGGTTTTACTG‑GCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTGAGAAGCGAAACGCCATCTTCCAGCGCCACCAGTTTCAGTT  >  NC_000913/74875‑75329

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