Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,207,066 A→G 13.7% E440E (GAA→GAG ymfN → predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,207,0660AG13.7% 129.4 / 13.3 51E440E (GAA→GAGymfNpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (21/23);  new base G (3/4);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATT  >  W3110S.gb/1206998‑1207125
                                                                    |                                                           
gAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGagaaagaa                                                             <  1:1362391/69‑1 (MQ=255)
gAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGagaaagaa                                                             <  1:1596962/69‑1 (MQ=255)
  aaaCAAAATCTATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAagagga                                                           <  1:890615/69‑1 (MQ=255)
  aaaCAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAagagga                                                           <  1:2935727/69‑1 (MQ=255)
  aaaCAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAagagga                                                           >  1:2567979/1‑69 (MQ=255)
  aaaCAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAagagga                                                           >  1:255458/1‑69 (MQ=255)
  aaaCAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAagagga                                                           <  1:1508886/69‑1 (MQ=255)
  aaaCAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAagagga                                                           <  1:2009819/69‑1 (MQ=255)
  aaaCAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAagaaga                                                           >  1:854344/1‑69 (MQ=255)
  aaaCAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAagaaga                                                           <  1:1768168/69‑1 (MQ=255)
   aaCAAAATCGCTGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGac                                                          >  1:1811839/1‑69 (MQ=255)
   aaCAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGac                                                          >  1:921168/1‑69 (MQ=255)
   aaCAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGac                                                          >  1:2527164/1‑69 (MQ=255)
   aaCAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGac                                                          >  1:1268613/1‑69 (MQ=255)
     cAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGacac                                                        <  1:969468/69‑1 (MQ=255)
     cAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGacac                                                        <  1:2302733/69‑1 (MQ=255)
         aTCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTg                                                    <  1:1073331/69‑1 (MQ=255)
          tCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGaca                                                         >  1:3416288/1‑63 (MQ=255)
               ggTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGacac                                                        >  1:1705844/1‑59 (MQ=255)
               ggTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTg                                                    >  1:873829/1‑63 (MQ=255)
               ggTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGAt                                              >  1:2257140/1‑69 (MQ=255)
               ggTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGAt                                              <  1:2300323/69‑1 (MQ=255)
                  gCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATcac                                           <  1:1090329/69‑1 (MQ=255)
                  gCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATcac                                           <  1:1725077/69‑1 (MQ=255)
                  gCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATcac                                           <  1:3622588/69‑1 (MQ=255)
                  gCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATca                                            <  1:930801/68‑1 (MQ=255)
                  gCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATca                                            <  1:854947/68‑1 (MQ=255)
                      ttGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTg                                       >  1:2955615/1‑69 (MQ=255)
                        gcgcTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAGGACACGCTGTCTGATCACATTGAg                                     <  1:2064427/69‑1 (MQ=255)
                        gcgcTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAg                                     >  1:3377311/1‑69 (MQ=255)
                        gcgcTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAg                                     >  1:2322285/1‑69 (MQ=255)
                          gcTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGt                                    <  1:1661574/68‑1 (MQ=255)
                          gcTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGt                                    <  1:1198254/68‑1 (MQ=255)
                          gcTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGt                                    <  1:1202550/68‑1 (MQ=255)
                          gcTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGt                                    <  1:2518936/68‑1 (MQ=255)
                          gcTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGt                                    <  1:13224/68‑1 (MQ=255)
                           cTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTcc                                  <  1:1250534/69‑1 (MQ=255)
                           cTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTcc                                  <  1:2628795/69‑1 (MQ=255)
                           cTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTcc                                  <  1:295130/69‑1 (MQ=255)
                             gATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTa                                <  1:3061684/69‑1 (MQ=255)
                              aTTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAagaaga                                                           >  1:2743107/1‑41 (MQ=255)
                              aTTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAagaaga                                                           <  1:240032/41‑1 (MQ=255)
                              aTTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAagaaga                                                           <  1:2137621/41‑1 (MQ=255)
                              aTTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAagaaga                                                           <  1:1181153/41‑1 (MQ=255)
                              aTTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTaa                               >  1:283312/1‑68 (MQ=255)
                                      ggTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACg                              <  1:369372/62‑1 (MQ=255)
                                      ggTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACg                              <  1:1636501/62‑1 (MQ=255)
                                        ttGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATca                                            <  1:1997261/46‑1 (MQ=255)
                                              gagCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCttt               <  1:412626/69‑1 (MQ=255)
                                               agCCATGCTGTACGAGAAAGAGGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTa              <  1:1115460/69‑1 (MQ=255)
                                                    tGCTGTACGAGAAAGAGGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGa         <  1:2991795/69‑1 (MQ=255)
                                                    tGCTGTACGAGAAAGAGGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGa         <  1:960623/69‑1 (MQ=255)
                                                      cTGTACGAGAAAGAGGACACTCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAg        >  1:240350/1‑68 (MQ=255)
                                                      cTGTACGAGAAAGAGGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAg        >  1:2909583/1‑68 (MQ=255)
                                                      cTGTACGAGAAAGAGGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAg        >  1:1590065/1‑68 (MQ=255)
                                                      cTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAg        >  1:3698546/1‑68 (MQ=255)
                                                       tGTACGCGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCttt               >  1:478353/1‑60 (MQ=255)
                                                       tGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCttt               >  1:2852286/1‑60 (MQ=255)
                                                       tGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCttt               >  1:1139987/1‑60 (MQ=255)
                                                        gTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGt                                    <  1:1641663/38‑1 (MQ=255)
                                                           cGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAAtt  >  1:70470/1‑69 (MQ=255)
                                                           cGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAAtt  >  1:395315/1‑69 (MQ=255)
                                                           cGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAAtt  >  1:928259/1‑69 (MQ=255)
                                                           cGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAAtt  >  1:1905394/1‑69 (MQ=255)
                                                           cGAGAAAGAAGACACGCTGACTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAAt   >  1:2583461/1‑68 (MQ=255)
                                                                 agaagaCACGCTGTCTGATCACATTGAGTCCTACGGGATCcgctcgc                  <  1:281295/47‑1 (MQ=255)
                                                                 agaagaCACGCTGTCTGATCACATTGAGTCCTACGGGATCcgctcgc                  <  1:2860077/47‑1 (MQ=255)
                                                                    |                                                           
GAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATT  >  W3110S.gb/1206998‑1207125

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