Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,521,967 C→A A171S (GCG→TCG)  xapR ← transcriptional activator of xapAB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,521,9670CA87.5% 179.4 / 21.0 65A171S (GCG→TCG) xapRtranscriptional activator of xapAB
Reads supporting (aligned to +/- strand):  ref base C (4/4);  new base A (31/25);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.00e-01

GGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCA  >  NC_000913/2521870‑2522062
                                                                                                 |                                                                                               
ggcggcATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAacgacg                                                                                               <  1:1454786/100‑1 (MQ=255)
    gcATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTc                                                                                           >  1:497419/1‑100 (MQ=255)
    gcATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTc                                                                                           <  2:623139/100‑1 (MQ=255)
     cATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCa                                                                                          >  2:1013917/1‑100 (MQ=255)
       ttGTAACAAAATTCTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAAGCATCGCCACCAGAAACGACGATTCATg                                                                                        <  1:1963993/100‑1 (MQ=255)
       ttGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATg                                                                                        <  1:1756678/100‑1 (MQ=255)
        tGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGc                                                                                       <  2:78342/100‑1 (MQ=255)
          tAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCaa                                                                                     >  2:1571858/1‑100 (MQ=255)
           aaCAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGTCGATTCATGCagc                                                                                    <  1:1384955/100‑3 (MQ=255)
            aCAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAAcc                                                                                   >  2:821875/1‑100 (MQ=255)
            aCAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAAcc                                                                                   >  1:415636/1‑100 (MQ=255)
                       gTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTa                                                                        >  2:123425/1‑100 (MQ=255)
                        tCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTaa                                                                       >  2:1397914/1‑100 (MQ=255)
                        tCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTaa                                                                       >  1:1778470/1‑100 (MQ=255)
                         cACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTaaa                                                                      >  1:847082/1‑100 (MQ=255)
                          aCGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTaaaa                                                                     >  1:1543918/1‑100 (MQ=255)
                           cGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAAc                                                                    <  2:2041502/100‑1 (MQ=255)
                           cGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAAc                                                                    >  1:586051/1‑100 (MQ=255)
                           cGTAGCGCGTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAAc                                                                    <  1:1408961/100‑1 (MQ=255)
                            gTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAAcc                                                                   <  2:1778455/100‑1 (MQ=255)
                            gTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAAcc                                                                   <  1:1271189/100‑1 (MQ=255)
                              aGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAg                                                                 <  2:497418/100‑1 (MQ=255)
                              aGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAg                                                                 >  1:831553/1‑100 (MQ=255)
                                cgcTTCAAGCGGGATGGAGGATAATGATGAGAGATGCTGATCTTCAGGCAGCTCCACCAGAAACGACGATGCATTCAACCGGAAGCTGGTAAAACCAGTc                                                               >  2:885149/1‑100 (MQ=255)
                                 gcTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCg                                                              >  2:426191/1‑100 (MQ=255)
                                  cTTCCAGCGGGACGGTGGAAAATGATGAGAGATGGTGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCgg                                                             <  1:1546331/100‑1 (MQ=255)
                                  cTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCgg                                                             <  1:1957102/100‑1 (MQ=255)
                                   ttCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCggt                                                            <  1:1076853/100‑1 (MQ=255)
                                    tCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCggtg                                                           <  2:299701/100‑1 (MQ=255)
                                     ccAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCggtgg                                                          <  2:146300/100‑1 (MQ=255)
                                      cAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCggtggt                                                         <  1:1530447/100‑1 (MQ=255)
                                              cGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGc                                                 >  1:1769384/1‑100 (MQ=255)
                                               ggTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGcc                                                <  1:1771706/100‑1 (MQ=255)
                                               ggTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGcc                                                <  1:1571868/100‑1 (MQ=255)
                                                  ggAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCAtt                                             >  2:287646/1‑100 (MQ=255)
                                                   gAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTc                                            >  1:1978251/1‑100 (MQ=255)
                                                   gAAAACGATGAGAGCTGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCAAGCAACCGTAAGCTGGTAAAACCAGCCGCTGGTTCTTTCGCCATTc                                            >  2:144912/1‑100 (MQ=255)
                                                    aaaaTGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCg                                           >  1:1366244/1‑100 (MQ=255)
                                                       atgatgAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCa                                        <  1:1833442/100‑1 (MQ=255)
                                                           tgaGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATc                                    >  2:1767357/1‑100 (MQ=255)
                                                            gagagaTGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATcc                                   <  1:1013922/100‑1 (MQ=255)
                                                            gagagaTGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATcc                                   <  2:1363811/100‑1 (MQ=255)
                                                             agagaTGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATccc                                  >  1:1380529/1‑100 (MQ=255)
                                                             agagaTGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATccc                                  >  2:1924662/1‑100 (MQ=255)
                                                               agaTGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCgg                                <  1:632384/100‑1 (MQ=255)
                                                                gatgatGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGc                               >  2:528904/1‑100 (MQ=255)
                                                                  tgatgCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCAt                             >  1:1918181/1‑100 (MQ=255)
                                                                   gatgCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATc                            <  1:990931/100‑1 (MQ=255)
                                                                     tgCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCaa                          <  2:415635/100‑1 (MQ=255)
                                                                      gCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAg                         <  2:586049/100‑1 (MQ=255)
                                                                         ctTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTc                      >  2:174353/1‑100 (MQ=255)
                                                                           tCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTcgc                    >  2:741717/1‑100 (MQ=255)
                                                                           tCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTcgc                    >  2:1833194/1‑100 (MQ=255)
                                                                            cAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTcgcg                   >  2:1292396/1‑100 (MQ=255)
                                                                              ggCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCgcggc                 <  2:2025604/100‑1 (MQ=255)
                                                                              ggCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCgcggc                 >  1:2025622/1‑100 (MQ=255)
                                                                                cATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCgcggc                 <  1:885152/98‑1 (MQ=255)
                                                                                    gCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTc             <  2:692609/98‑1 (MQ=255)
                                                                                     ccaccaGACACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCCGATCCCGGCATCCAGTTCGCGGCGTTCCAg          >  2:1749691/1‑100 (MQ=255)
                                                                                      caccaGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGc         >  2:1717166/1‑100 (MQ=255)
                                                                                       accaGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCa        <  1:1767371/100‑1 (MQ=255)
                                                                                       accaGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCa        >  2:186664/1‑100 (MQ=255)
                                                                                        ccaGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCaa       >  1:1768817/1‑100 (MQ=255)
                                                                                         caGAAACGACGATTCATGCAACCGTAAGCTTGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAg      <  2:239933/100‑1 (MQ=255)
                                                                                           gAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGc    >  2:1619844/1‑100 (MQ=255)
                                                                                             aaCGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCa  >  1:1501096/1‑100 (MQ=255)
                                                                                                 |                                                                                               
GGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCA  >  NC_000913/2521870‑2522062

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