Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,938,519 A→G 56.3% I166V (ATC→GTC)  fucR → DNA‑binding transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,938,5190AG56.3% 30.6 / 90.6 71I166V (ATC→GTC) fucRDNA‑binding transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (9/22);  new base G (27/13);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.83e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

CTTGAGCGCAAATATGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCACTGTGGGACTCCA  >  W3110S.gb/2938456‑2938579
                                                               |                                                            
cTTGAGCGCAAATATGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGAt                                                         >  1:3512731/1‑69 (MQ=255)
 ttGAGCGCAAATATGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATc                                                        >  1:444901/1‑69 (MQ=255)
  tGAGCGCAAATATGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATc                                                        >  1:673518/1‑68 (MQ=255)
  tGAGCGCAAATATGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATc                                                        >  1:3162249/1‑68 (MQ=255)
  tGAGCGCAAATATGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATc                                                        >  1:2329568/1‑68 (MQ=255)
  tGAGCGCAAATATGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATc                                                        >  1:1598580/1‑68 (MQ=255)
   gAGCGCAAATATGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGa                                                          <  1:1674337/65‑1 (MQ=255)
           atatGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTAtttt                                              <  1:1005505/69‑1 (MQ=255)
           atatGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTAtttt                                              <  1:2759651/69‑1 (MQ=255)
           atatGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTAtttt                                              <  1:1001216/69‑1 (MQ=255)
                gCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGa                                                          <  1:244582/52‑1 (MQ=255)
                  tGTTAAGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTATTTTTTCTtgtg                                      <  1:1239265/69‑1 (MQ=255)
                   gTTACGTTAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTtgtg                                      <  1:2523290/69‑1 (MQ=255)
                   gTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTtgtg                                      <  1:1391275/69‑1 (MQ=255)
                   gTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTtgtg                                      <  1:3354327/69‑1 (MQ=255)
                      aCGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGaa                                    >  1:2489011/1‑68 (MQ=255)
                      aCGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGaa                                    >  1:2981021/1‑68 (MQ=255)
                      aCGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGaa                                    >  1:3263068/1‑68 (MQ=255)
                      aCGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGa                                     >  1:802670/1‑67 (MQ=255)
                      aCGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCg                                                           <  1:3699424/45‑1 (MQ=255)
                      aCGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCg                                                           <  1:2336285/45‑1 (MQ=255)
                      aCGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCt                                          <  1:2282208/62‑1 (MQ=255)
                      aCGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCt                                          <  1:2348262/62‑1 (MQ=255)
                       cGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAgg                                  >  1:1151290/1‑69 (MQ=255)
                         tCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGa                                >  1:3049054/1‑69 (MQ=255)
                           aaTCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATc                              <  1:2917173/69‑1 (MQ=255)
                           aaTCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATc                              <  1:2442537/69‑1 (MQ=255)
                           aaTCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAg                                   >  1:2381401/1‑64 (MQ=255)
                           aaTCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATc                              <  1:2567662/69‑1 (MQ=255)
                           aaTCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATc                              >  1:1216987/1‑69 (MQ=255)
                           aaTCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATc                              <  1:1644996/69‑1 (MQ=255)
                             tCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCg                             <  1:2655486/68‑1 (MQ=255)
                              cccTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGAt                           <  1:1333951/69‑1 (MQ=255)
                              cccTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGAt                           <  1:469355/69‑1 (MQ=255)
                               ccTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGAt                           >  1:447727/1‑68 (MQ=255)
                               ccTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGAt                           >  1:2462818/1‑68 (MQ=255)
                               ccTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTc                                           >  1:503839/1‑52 (MQ=255)
                               ccTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATc                              <  1:458503/65‑1 (MQ=255)
                               ccTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATc                              >  1:1534792/1‑65 (MQ=255)
                                 tCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATagc                        <  1:1589239/69‑1 (MQ=255)
                                   gCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTtgtg                                      <  1:162278/53‑1 (MQ=255)
                                   gCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTtgtg                                      >  1:673295/1‑53 (MQ=255)
                                   gCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATagcag                      >  1:1652269/1‑69 (MQ=255)
                                   gCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATagcag                      >  1:2340991/1‑69 (MQ=255)
                                         ttCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATagc                        <  1:2857073/61‑1 (MQ=255)
                                         ttCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATagc                        <  1:172774/61‑1 (MQ=255)
                                          tCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATc                              >  1:3046229/1‑54 (MQ=255)
                                          tCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATc                              >  1:2002363/1‑54 (MQ=255)
                                          tCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATc                              >  1:2596209/1‑54 (MQ=255)
                                          tCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCa               <  1:128093/69‑1 (MQ=255)
                                          tCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCa               <  1:3118582/69‑1 (MQ=255)
                                              aaCTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCACtgt           >  1:1380851/1‑69 (MQ=255)
                                              aaCTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCACtg            >  1:3153993/1‑68 (MQ=255)
                                                cTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGa                            >  1:2012330/1‑50 (MQ=255)
                                                cTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGa                            <  1:877940/50‑1 (MQ=255)
                                                cTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGa                            <  1:2049395/50‑1 (MQ=255)
                                                cTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGa                            <  1:2225359/50‑1 (MQ=255)
                                                cTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGa                            >  1:273148/1‑50 (MQ=255)
                                                 ttAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCACTGTggg        >  1:1834263/1‑69 (MQ=255)
                                                  tAAATCGCTGGAAGTCGATCTGTTTATTTTTTATTGTGAAGGGATCGATAGCAGCGGCGCACTGTGGGa       >  1:3133074/1‑69 (MQ=255)
                                                  tAAATCGCTGCAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCACTGTGGGa       >  1:483309/1‑69 (MQ=255)
                                                  tAAATCGCAGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATa                          >  1:2113033/1‑50 (MQ=39)
                                                   aaaTCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATc                              <  1:3507319/45‑1 (MQ=255)
                                                   aaaTCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGAt                           >  1:1915826/1‑48 (MQ=255)
                                                   aaaTCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGAt                           >  1:1701791/1‑48 (MQ=255)
                                                   aaaTCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGAt                           >  1:2913760/1‑48 (MQ=255)
                                                      tCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCACTGTGGGACTcc   <  1:908778/69‑1 (MQ=255)
                                                       cGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATGGCAGCGGCGCACTGTGGGACTCCa  >  1:684269/1‑69 (MQ=255)
                                                            gAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCACTGTgg         <  1:126567/57‑1 (MQ=255)
                                                            gAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCACTGTgg         <  1:1698956/57‑1 (MQ=255)
                                                            gAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCACTGTgg         <  1:1279055/57‑1 (MQ=255)
                                                               |                                                            
CTTGAGCGCAAATATGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCACTGTGGGACTCCA  >  W3110S.gb/2938456‑2938579

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: