Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,998,457 T→C 6.5% E171E (GAA→GAG sdiA ← DNA‑binding transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,998,4570TC 6.5% 171.3 / 3.7 61E171E (GAA→GAGsdiADNA‑binding transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (27/30);  new base C (2/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

CCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTT  >  W3110S.gb/1998393‑1998524
                                                                |                                                                   
ccgccgTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTtcatc                                                                <  1:293552/70‑1 (MQ=255)
    cgTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTa                                                           <  1:1088096/71‑1 (MQ=255)
    cgTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTa                                                           <  1:2520002/71‑1 (MQ=255)
        cACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACg                                                       >  1:2998948/1‑71 (MQ=255)
        cACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACg                                                       >  1:3106725/1‑71 (MQ=255)
          ccTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTa                                                           <  1:84951/65‑1 (MQ=255)
             cAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCAt                                                    >  1:1037191/1‑69 (MQ=255)
             cAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCAt                                                    >  1:1764422/1‑69 (MQ=255)
                aaTTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTa                                                           <  1:320204/59‑1 (MQ=255)
                aaTTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTa                                                           <  1:1302146/59‑1 (MQ=255)
                aaTTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTa                                                           <  1:565812/59‑1 (MQ=255)
                aaTTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTa                                                           <  1:1089327/59‑1 (MQ=255)
                aaTTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCagag                                               <  1:2299657/71‑1 (MQ=255)
                  tttCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGcc                                             <  1:2149823/71‑1 (MQ=255)
                  tttCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGc                                              >  1:3350406/1‑70 (MQ=255)
                    tCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCa                                                  >  1:3163839/1‑64 (MQ=255)
                         ttCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGAc                                      >  1:1100355/1‑71 (MQ=255)
                         ttCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGAc                                      >  1:2572961/1‑71 (MQ=255)
                           ggcgcTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACt                                     <  1:3317878/69‑1 (MQ=255)
                           ggcgcTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACt                                     <  1:3313340/69‑1 (MQ=255)
                           cgcgCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACtt                                    >  1:3401194/1‑71 (MQ=255)
                           cgcgCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACtt                                    >  1:91137/1‑71 (MQ=255)
                           cgcgCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACtt                                    >  1:2557784/1‑71 (MQ=255)
                           cgcgCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACtt                                    >  1:2251742/1‑71 (MQ=255)
                           cgcgCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACt                                     <  1:943621/70‑1 (MQ=255)
                           cgcgCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACt                                     <  1:2992885/70‑1 (MQ=255)
                           cgcgCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACt                                     <  1:1300958/70‑1 (MQ=255)
                           cgcgCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCACTTAAACGCATCAGAGCCATCAGACt                                     <  1:2866675/70‑1 (MQ=255)
                                  gCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTcgcg                               <  1:1363061/69‑1 (MQ=255)
                                               tctGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCAt                >  1:3580911/1‑71 (MQ=255)
                                                   gCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCAttttt            <  1:1145921/71‑1 (MQ=255)
                                                   gCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCAttttt            <  1:805550/71‑1 (MQ=255)
                                                   gCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCAttttt            <  1:1553249/71‑1 (MQ=255)
                                                      tcatcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGt                        <  1:214704/56‑1 (MQ=255)
                                                      tcatcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGt                        <  1:3336046/56‑1 (MQ=255)
                                                      tcatcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGt                        <  1:2805331/56‑1 (MQ=255)
                                                      tcatcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCAttttt            <  1:2869550/68‑1 (MQ=255)
                                                      tcatcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCAttttt            <  1:1132905/68‑1 (MQ=255)
                                                      tcatcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCAttttt            <  1:942681/68‑1 (MQ=255)
                                                      tcatcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCAttttt            <  1:61293/68‑1 (MQ=255)
                                                       catcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGt                        <  1:182932/55‑1 (MQ=255)
                                                        atcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGAc                                      >  1:3303163/1‑40 (MQ=255)
                                                        atcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTaa                      >  1:3141385/1‑56 (MQ=255)
                                                        atcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTaa                      >  1:2213177/1‑56 (MQ=255)
                                                        atcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTaa                      >  1:1936355/1‑56 (MQ=255)
                                                        atcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTaa                      >  1:57559/1‑56 (MQ=255)
                                                        atcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTaa                      >  1:1406174/1‑56 (MQ=255)
                                                        atcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTaa                      >  1:1191808/1‑56 (MQ=255)
                                                        atcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTg       >  1:1938933/1‑71 (MQ=255)
                                                         tcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGc      <  1:385932/71‑1 (MQ=255)
                                                         tcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGc      <  1:2543297/71‑1 (MQ=255)
                                                         tcaCTATCTCATCATTTAAACGCATCAGAGCCATCAg                                        <  1:1259545/37‑1 (MQ=255)
                                                         tcaCTATCTCATCATTTAAACGCATCAGAGCCATCAg                                        <  1:2549889/37‑1 (MQ=255)
                                                           aCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAg    >  1:2108631/1‑71 (MQ=255)
                                                           aCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAg    >  1:2104386/1‑71 (MQ=255)
                                                           aCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAg    >  1:3141961/1‑71 (MQ=255)
                                                           aCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAg    >  1:197865/1‑71 (MQ=255)
                                                           aCTATCTCATCATTTAAACGCATCAGAGCCATCAGAc                                      >  1:32270/1‑37 (MQ=255)
                                                           aCTATCTCATCATTTAAACGCATCAGAGCCATCAGAc                                      >  1:1699545/1‑37 (MQ=255)
                                                              aTTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGtt  >  1:2115919/1‑70 (MQ=255)
                                                              aTTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGtt  >  1:1166180/1‑70 (MQ=255)
                                                               tttCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGtt  <  1:211382/69‑1 (MQ=255)
                                                               tttCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGtt  <  1:550519/69‑1 (MQ=255)
                                                                |                                                                   
CCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTT  >  W3110S.gb/1998393‑1998524

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