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→GAAthiBthiamine/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

GCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTGAGAAGCGAAACGCCATCTTCCAGCGCCACCAGTTTCAGTT  >  NC_000913/74875‑75329
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
gCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGTCGGGGTTTTACGTTTTTTTTTTTCTTTAACGaaatgcaaggag                                                                                                                                                                                                                              >  1:751377/1‑228 (MQ=255)
gCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTTGTTATAAAAAAAgggaaagaag                                                                                                                                                                                                                              >  2:244461/1‑231 (MQ=255)
gCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTTTTCTTTTTCTTAACAAAAGCGCAGTAg                                                                                                                                                                                                                              >  1:68849/1‑238 (MQ=255)
        cTTCCCCCCAGTTTTTGTTGCCCGTGCCCGTTTTTTCCGCCATTTTCGCCCGGGTTGGGGGGGCGTCTTCGCCATAGATTTTTTGCACCAATAGCAACAGCCCCAGCCCGGGTGTACTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCa                                                                                                                                                                                                                     <  2:689352/239‑1 (MQ=255)
         ttCGCTCCAGCTTTGGGTTACTGGCACCGTTTTCGTCGCCAGTTGCTGCCGGGCTGGGGGGGCGTCATGACCAAAGACTTTTTGCATCCATAGCAACAGCCCCATCCCCGGTTTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCa                                                                                                                                                                                                                     <  1:62529/238‑1 (MQ=255)
          tCGCTCCCGCCTTGGGGCACGTTGCCCTTTTTCTTCGCCAGTTGCTGCAGGGCTTGGGGGGCGTCATCGCCTTAGACTTTTTGCACCCATAGCAACAGCCCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCaaa                                                                                                                                                                                                                   <  1:377130/239‑1 (MQ=255)
          tCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGCCGGGGTTTTCCTTTTTTTATTTTTATAAACAAAAGCGAAGTGGtcatcatcaca                                                                                                                                                                                                                   >  1:321225/1‑237 (MQ=255)
                        ggTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGCGTTTTTAGTTTTTTTTTTTCATATACAAAAAGGCAGTAGTCAATATTACACTGTACGaaacggg                                                                                                                                                                                                     >  1:20953/1‑235 (MQ=255)
                             ccgtgaccgtTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTCAGTTTGTTTTTTTCTTAAAACAAGGCGGAGTAGCCCTAATCAACAGGGAAGAAAAAGTCatta                                                                                                                                                                                                >  1:823010/1‑239 (MQ=255)
                                                                     ggCGTCAGCGCCATAGCCTTTTGGCATCCATAGCAACAGCCCCAGCCCCGGTGTACGGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCcaca                                                                                                                                                        <  1:244461/238‑1 (MQ=255)
                                                                                                   tAGCAACAGCCCCGGCCCGGGTGTACTGGTGCGCGGATCCTGATAAACCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTgg                                                                                                                       <  1:743231/239‑1 (MQ=255)
                                                                                                       acacgACCCAGCCCGGGTGTATGGGTGCCCGGATCCTGATAAACCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTggcggc                                                                                                                   <  1:264522/235‑1 (MQ=255)
                                                                                                            aCCCAGCCCCGGTGTACTGGTGCTCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGATTTTTCAGTTTGTTCTTTTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCCGGAAAGTTAACGGCATCCGCTGCCACACCGCTTTTGGGCAACAGTCCGGTTTTACTTGAGGCGACTa                                                                                                              >  1:697796/1‑239 (MQ=255)
                                                                                                                 gCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTCAGTTTTTTTTTGTCTTAAACGAAGGCGAAGTAGCCCTAATCAAACGGTACGAAAATGTCATTATTCCCGCCGCCGGGGACGCTAACGGCATCCCCCGCCACACCCCTTTTTGCCAAAAACCCGGCTTTATTGGCGGCGGCTTACAgg                                                                                                         >  1:571791/1‑239 (MQ=255)
                                                                                                                     gggTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAg                                                                                                                                      <  1:657677/206‑1 (MQ=255)
                                                                                                                               ggtgggggATTCTTATAAATCCCCCGCAAGTTTTGATCGCTCCCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGccc                                                                                           <  1:233420/233‑1 (MQ=255)
                                                                                                                                  cgcgGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTTTTTTTTTCCTAAAACAAAGGGGAGGAGCCCATAATCAAACGTGAGGCAAGGTCATTTTTTCAGACCCCCGGGAAGGTAAACGCCACCGCCTCCACCCCCCCTTTGGTCAAAAAGCCGGGTTTTATTGTGCCGCGCTAAAAGATTGTTATACAcccaa                                                                                          >  2:231752/1‑236 (MQ=255)
                                                                                                                                          cttaTAACTCCCCCTCCATTTTTGACCGCTCTCAACCAGTTTTTTCAGGCTTGGTGGCGGTTTTTTCAGTTTGTTTTTGTCATAAACGAAGGCGAAGTAGCCAAAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACAt                                                                                <  2:20953/236‑1 (MQ=255)
                                                                                                                                                                tGATCGCTCTCACCCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCAAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTc                                                          <  1:439629/239‑1 (MQ=255)
                                                                                                                                                                                              tttgGGGGTTTTTTCTGTTTTTTTTTTTCATAACCAAAGGGGAAGTAGCAATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTGAGAAGCGAAACGcc                            <  2:823010/237‑1 (MQ=255)
                                                                                                                                                                                                        tttttCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAATTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTGAGAAGCGAAACGCCATCTTCCa                    <  1:574853/238‑1 (MQ=255)
                                                                                                                                                                                                           ttCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACCGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCCTCCGTCGACGGTTGATAAGCGCAAAGCCATCTTCCAgcg                 >  2:753735/1‑238 (MQ=255)
                                                                                                                                                                                                            tCAGTTTGTTCTGTTAAAAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTTTTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTGAGAAGCGAAACGCCATCTTCCAGCGcc               <  1:231752/239‑1 (MQ=255)
                                                                                                                                                                                                                    tttttGTCATAAAGTAAGGCAAGTTAGGTATAATCAAACGGTACGAAATTTGCATTATCCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTGAGAAGCGAAACGCCATCTTCCAGCGCCACCAGttt       <  2:279749/236‑1 (MQ=255)
                                                                                                                                                                                                                         gTCATGAACGAAGGCGAAGTGGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTGAGAAGCGAAACGCCATCTTCCAGCGCCACCAGTTTCAGtt  <  1:594890/239‑1 (MQ=255)
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
GCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTGAGAAGCGAAACGCCATCTTCCAGCGCCACCAGTTTCAGTT  >  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: