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,9670CA96.6% 199.5 / ‑0.8 59A171S (GCG→TCG) xapRtranscriptional activator of xapAB
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (32/25);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

GGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTT  >  NC_000913/2521870‑2522065
                                                                                                 |                                                                                                  
ggcggcATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAacgacg                                                                                                  >  1:518084/1‑100 (MQ=255)
 gcggcATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAAcgacga                                                                                                 >  1:6765/1‑100 (MQ=255)
  cggcATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGAt                                                                                                <  2:92214/100‑1 (MQ=255)
    gcATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTc                                                                                              >  2:796561/1‑100 (MQ=255)
    gcATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTc                                                                                              <  1:1319952/100‑1 (MQ=255)
           aaCAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAAc                                                                                       >  2:1430252/1‑100 (MQ=255)
             cAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCg                                                                                     >  1:1119809/1‑100 (MQ=255)
              aaaaTACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGt                                                                                    >  2:728035/1‑100 (MQ=255)
               aaaTACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTa                                                                                   <  2:1014340/100‑1 (MQ=255)
                 aTACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAg                                                                                 <  1:636208/100‑1 (MQ=255)
                   aCTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCt                                                                               <  1:1428770/100‑1 (MQ=255)
                    cTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTg                                                                              >  1:1781526/1‑100 (MQ=255)
                     tCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTgg                                                                             <  1:1281959/100‑1 (MQ=255)
                      cGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGt                                                                            >  1:41425/1‑100 (MQ=255)
                       gTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTa                                                                           >  2:354440/1‑100 (MQ=255)
                       gTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTa                                                                           >  2:353019/1‑100 (MQ=255)
                           cGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAAc                                                                       <  1:401852/100‑1 (MQ=255)
                            gTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAAcc                                                                      <  2:477730/100‑1 (MQ=255)
                            gTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTGTAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAAcc                                                                      <  1:495954/100‑1 (MQ=255)
                               gcgcTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGt                                                                   >  1:903033/1‑100 (MQ=255)
                                      cggCGGGACGGTGAAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCggtggt                                                            >  2:1320260/3‑100 (MQ=255)
                                      cAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCggtggt                                                            <  1:433921/100‑1 (MQ=255)
                                      cAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCggtggt                                                            <  1:441901/100‑1 (MQ=255)
                                        gCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTc                                                          <  1:1384687/100‑1 (MQ=255)
                                            gACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTc                                                      >  1:1072234/1‑100 (MQ=255)
                                               ggTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGcc                                                   <  2:55010/100‑1 (MQ=255)
                                                gTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCa                                                  >  2:1939517/1‑100 (MQ=255)
                                                 tGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCAt                                                 >  1:399264/1‑100 (MQ=255)
                                                    aaaaTGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCg                                              >  1:1447165/1‑100 (MQ=255)
                                                      aaTGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGcc                                            >  2:1749163/1‑100 (MQ=255)
                                                      aaTGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGcc                                            <  2:872538/100‑1 (MQ=255)
                                                       atgatgAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCa                                           <  2:6765/100‑1 (MQ=255)
                                                         gatgaGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGa                                         >  2:906673/1‑100 (MQ=255)
                                                          atgaGAGATGATGCTCTTCAGGCATCACCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGAt                                        >  2:1187492/1‑100 (MQ=255)
                                                          atgaGAGATGATGCTCCTCAGGAATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAAGCAGTCTGTGGTTCTGTCGCCAGGCTCCAGa                                         >  1:646213/1‑99 (MQ=255)
                                                           tgaGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATc                                       <  2:1119810/100‑1 (MQ=255)
                                                               agaTGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCgg                                   <  2:300858/100‑1 (MQ=255)
                                                                gatgatGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGc                                  >  1:984378/1‑100 (MQ=255)
                                                                gatgatGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGc                                  <  2:795535/100‑1 (MQ=255)
                                                                gatgatGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGc                                  <  1:1589171/100‑1 (MQ=255)
                                                                gatgatGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCAATTCGCCAGATCCCg                                    >  1:1495462/1‑98 (MQ=255)
                                                                   gttgCTCTTCATGCATCTCCACCAGAAACGACGATTCATTCAACCGTAAGCTGGTAAAACCAGTCGGTCGTTCTGTCGCCCTTCCCCAGACCCCGGCa                                 <  1:362074/96‑1 (MQ=255)
                                                                    atgCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCa                              >  1:1606269/1‑100 (MQ=255)
                                                                     tgCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCaa                             <  2:554301/100‑1 (MQ=255)
                                                                      gCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAg                            >  1:514109/1‑100 (MQ=255)
                                                                       ctctTCAGGCATCGCCACCAGAAACGACAATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAg                            <  1:165042/99‑1 (MQ=255)
                                                                         ctTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTc                         >  1:1552271/1‑100 (MQ=255)
                                                                         ctTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTc                         >  2:1282930/1‑100 (MQ=255)
                                                                             aGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCgcgg                     <  2:41425/100‑1 (MQ=255)
                                                                              ggCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCgcggc                    >  2:686331/1‑100 (MQ=255)
                                                                                   cGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTcc               >  2:709837/1‑100 (MQ=255)
                                                                                    gCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCa              >  1:144922/1‑100 (MQ=255)
                                                                                     ccaccaGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAg             >  2:374482/1‑100 (MQ=255)
                                                                                       accaGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCa           >  1:228459/1‑100 (MQ=255)
                                                                                        ccaGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCaa          >  2:1780918/1‑100 (MQ=255)
                                                                                         caGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAg         >  1:443137/1‑100 (MQ=255)
                                                                                          aGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAgg        <  2:325756/100‑1 (MQ=255)
                                                                                            aaaCGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTCCGCGGCGTTCCAGCAAGGcc      <  1:1187077/100‑1 (MQ=255)
                                                                                             aaCGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCa     >  2:1528164/1‑100 (MQ=255)
                                                                                             aaCGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCa     >  1:786268/1‑100 (MQ=255)
                                                                                              acgacgATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCAt    >  2:1326219/1‑100 (MQ=255)
                                                                                                gCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCAttt  <  1:1640205/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
GGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTT  >  NC_000913/2521870‑2522065

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: