Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,207,066 A→G 20.6% 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,0660AG20.6% 67.2 / 14.4 34E440E (GAA→GAGymfNpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (15/12);  new base G (3/4);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCT  >  W3110S.gb/1206996‑1207136
                                                                      |                                                                      
cGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGagaaagaa                                                                        <  1:1296738/71‑1 (MQ=255)
cGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGagaaagaa                                                                        <  1:3453319/71‑1 (MQ=255)
cGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGagaaagaa                                                                        <  1:1593971/71‑1 (MQ=255)
   aaaaCAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGagaaagaa                                                                        <  1:1032600/68‑1 (MQ=255)
   aaaaCAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGagaaagaa                                                                        <  1:1172347/68‑1 (MQ=255)
   aaaaCAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGagaaagaa                                                                        <  1:1291908/68‑1 (MQ=255)
   aaaaCAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGagaaagaa                                                                        <  1:3627874/68‑1 (MQ=255)
   aaaaCAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGagaaagaa                                                                        <  1:2403904/68‑1 (MQ=255)
        aaaaTCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCGTGCTGTACGAGAAAGAAGACACGCt                                                                <  1:782536/71‑1 (MQ=255)
              gATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAGGACACGCTGTCTGa                                                          >  1:168316/1‑71 (MQ=255)
                 ggTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTATGAGAAAGAGGACACGCTGTCTGATca                                                       >  1:452218/1‑71 (MQ=255)
                 ggTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAGGACACGCTGTCTGATca                                                       >  1:45726/1‑71 (MQ=255)
                  gTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATcac                                                      <  1:2473335/71‑1 (MQ=255)
                  gTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATca                                                       <  1:470031/70‑1 (MQ=255)
                   tGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATca                                                       <  1:2615614/69‑1 (MQ=255)
                     cAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAGGACACGCTGTCTGATCACAtt                                                   <  1:1980648/71‑1 (MQ=255)
                     cAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAGGACACGCTGTCTGATCACAtt                                                   <  1:453294/71‑1 (MQ=255)
                        ttGCGCTGATTATGGCGGTTGGCAGAGTCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAg                                                <  1:1365078/71‑1 (MQ=255)
                        ttGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAg                                                <  1:2723560/71‑1 (MQ=255)
                         tGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAg                                                >  1:3934709/1‑70 (MQ=255)
                         tGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGt                                               >  1:3124349/1‑71 (MQ=255)
                         tGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGt                                               >  1:3816543/1‑71 (MQ=255)
                         tGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGt                                               >  1:1400452/1‑71 (MQ=255)
                         tGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGt                                               >  1:921271/1‑71 (MQ=255)
                          gcgcTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAagaaga                                                                      >  1:419201/1‑47 (MQ=255)
                          gcgcTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTc                                              >  1:16451/1‑71 (MQ=255)
                          gcgcTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTc                                              >  1:1537249/1‑71 (MQ=255)
                             cTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTg                                                           <  1:3938200/55‑1 (MQ=255)
                               gATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACg                                         >  1:3845953/1‑71 (MQ=255)
                                  tATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGac                                                                     >  1:2503139/1‑40 (MQ=255)
                                  tATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGac                                                                     >  1:915979/1‑40 (MQ=255)
                                  tATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGac                                                                     >  1:3901359/1‑40 (MQ=255)
                                       cGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACGGGATCcgc                                 <  1:534503/71‑1 (MQ=255)
                                          ttGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACggg                                       <  1:738304/62‑1 (MQ=255)
                                          ttGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACggg                                       <  1:2455167/62‑1 (MQ=255)
                                           tGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACGGGATCcgctcgc                             >  1:3472316/1‑71 (MQ=255)
                                              cAGAGCCATGCTGTACGAGAAAGAGGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCttt                          <  1:1038574/71‑1 (MQ=255)
                                              cAGAGCCATGCTGTACGAGAAAGAGGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCttt                          <  1:2480356/71‑1 (MQ=255)
                                              cAGAGCCATGATGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCttt                          <  1:1904194/71‑1 (MQ=255)
                                                   ccATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTg                     <  1:3416879/71‑1 (MQ=255)
                                                        cTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACGGGATCcgctc                               >  1:261896/1‑56 (MQ=255)
                                                        cTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAgg                  >  1:739341/1‑69 (MQ=255)
                                                        cTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGt                 >  1:3262065/1‑70 (MQ=255)
                                                                gaaagagGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAAtt             <  1:3816095/66‑1 (MQ=255)
                                                                gaaagagGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAAtt             <  1:919758/66‑1 (MQ=255)
                                                                   agaggaCACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCAt     >  1:3254284/1‑71 (MQ=255)
                                                                      agaCACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCt  <  1:2688312/71‑1 (MQ=255)
                                                                      agaCACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCt  <  1:2537400/71‑1 (MQ=255)
                                                                      agaCACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCt  <  1:895347/71‑1 (MQ=255)
                                                                      |                                                                      
CGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGCTGTACGAGAAAGAAGACACGCTGTCTGATCACATTGAGTCCTACGGGATCCGCTCGCTTTAACTGAGGTAATTATGATCATGCT  >  W3110S.gb/1206996‑1207136

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