Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,938,519 A→G 70.4% 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,5190AG70.4% 80.4 / 41.4 54I166V (ATC→GTC) fucRDNA‑binding transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (13/3);  new base G (18/20);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 3.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GCACACTTGAGCGCAAATATGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCACTGTGGGACTCC  >  W3110S.gb/2938451‑2938578
                                                                    |                                                           
gcACACTTGAGCGCAAATATGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTc                                                           <  1:3187052/71‑1 (MQ=255)
 cacaCTTGAGCGCAAATATGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCg                                                          <  1:695226/71‑1 (MQ=255)
   cacTTGAGCGCAAATATGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGAt                                                        >  1:1284888/1‑71 (MQ=255)
          gcgcAAATATGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGttt                                                  <  1:1608387/70‑1 (MQ=255)
            gcAAATATGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTAtt                                               <  1:1144367/71‑1 (MQ=255)
            gcAAATATGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTAtt                                               <  1:1439698/71‑1 (MQ=255)
            gcAAATATGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTAtt                                               <  1:1315228/71‑1 (MQ=255)
                 tatGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTc                                          >  1:1674791/1‑71 (MQ=255)
                  atGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCt                                         >  1:2769876/1‑71 (MQ=255)
                      cTGTTACGGCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTtgtg                                     <  1:842903/71‑1 (MQ=255)
                         ttACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTtgtg                                     <  1:766160/68‑1 (MQ=255)
                         ttACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTtgtg                                     <  1:2395520/68‑1 (MQ=255)
                         ttACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGa                                    <  1:122633/69‑1 (MQ=255)
                          tACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATc                                                       >  1:3513302/1‑49 (MQ=255)
                          tACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATc                                                       >  1:1133335/1‑49 (MQ=255)
                          tACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATc                                                       >  1:2141472/1‑49 (MQ=255)
                          tACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATc                                                       >  1:1842498/1‑49 (MQ=255)
                           aCGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGAt                                                        >  1:2315874/1‑47 (MQ=255)
                           aCGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGAt                                                        >  1:2711303/1‑47 (MQ=255)
                              tCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATc                             <  1:1239033/71‑1 (MQ=255)
                              tCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATc                             <  1:3040324/71‑1 (MQ=255)
                              tCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATc                             <  1:2685311/71‑1 (MQ=255)
                                aaTCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGaa                                   >  1:2003488/1‑63 (MQ=255)
                                aaTCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGaa                                   >  1:1612104/1‑63 (MQ=255)
                                aaTCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTAGTGaa                                   >  1:3666674/1‑63 (MQ=255)
                                    ccTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGaa                                   >  1:1902696/1‑59 (MQ=255)
                                    ccTCGCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGaa                                   >  1:2808671/1‑59 (MQ=255)
                                    ccTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATag                        >  1:3567514/1‑70 (MQ=255)
                                    ccTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATag                        >  1:282391/1‑70 (MQ=255)
                                     cTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATa                         >  1:698071/1‑68 (MQ=255)
                                     cTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATa                         >  1:3417023/1‑68 (MQ=255)
                                     cTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATa                         >  1:2345708/1‑68 (MQ=255)
                                     cTCGCTGATTTCCCAACTTAAATCGCTCGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATa                         >  1:2738465/1‑68 (MQ=255)
                                        gCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGAt                              >  1:770969/1‑60 (MQ=255)
                                        gCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGAt                              >  1:2264802/1‑60 (MQ=255)
                                        gCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGAt                              >  1:1346331/1‑60 (MQ=255)
                                        gCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGAt                              >  1:1784434/1‑60 (MQ=255)
                                        gCTGATTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGCGAt                              >  1:1115801/1‑60 (MQ=255)
                                            aTTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAgcggcg                <  1:3575806/70‑1 (MQ=255)
                                            aTTTCCCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAgcggcg                <  1:57621/70‑1 (MQ=255)
                                            aTTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGgcgc               <  1:1280946/71‑1 (MQ=255)
                                              ttCGCAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCa              <  1:3306333/70‑1 (MQ=255)
                                               tCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCAc             <  1:3350229/70‑1 (MQ=255)
                                               tCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCAc             <  1:21829/70‑1 (MQ=255)
                                                cccAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCACtg           >  1:714267/1‑71 (MQ=255)
                                                 ccAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCACtgt          <  1:833323/71‑1 (MQ=255)
                                                 ccAACTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCACtgt          <  1:3051751/71‑1 (MQ=255)
                                                   aaCTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCACTGTgg        >  1:1785328/1‑71 (MQ=255)
                                                   aaCTTAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCACTGTgg        >  1:1714568/1‑71 (MQ=255)
                                                       tAAATCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATa                         >  1:1524719/1‑50 (MQ=255)
                                                       tAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGaa                                   >  1:1190856/1‑40 (MQ=255)
                                                       tAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGaa                                   >  1:859906/1‑40 (MQ=255)
                                                        aaaTCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCACTGTGGGACt    <  1:723540/70‑1 (MQ=255)
                                                         aaTCGCTGGAAGTCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCACTGTGGGACTcc  <  1:1864026/71‑1 (MQ=255)
                                                                    |                                                           
GCACACTTGAGCGCAAATATGGCTGTTACGTCAATCCCTCGCTGATTTCCCAACTTAAATCGCTGGAAATCGATCTGTTTATTTTTTCTTGTGAAGGGATCGATAGCAGCGGCGCACTGTGGGACTCC  >  W3110S.gb/2938451‑2938578

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: