Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,007,539 G→A 14.3% pseudogene (367/393 nt) yhiK → ECK3474:JW5674:b3489; hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,007,5390GA14.3% 71.6 / 4.7 28pseudogene (367/393 nt)yhiKECK3474:JW5674:b3489; hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (11/13);  new base A (2/2);  total (13/15)
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.98e-01

GTATTATGGGGCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAGGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACGGTGGTATGTCTTTTTAATCCTGC  >  W3110S.gb/4007474‑4007603
                                                                 |                                                                
gTATTATGGGGCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAagg                                                                >  1:3492676/1‑68 (MQ=255)
     atGGGGCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAGGGCTTAgt                                                          >  1:3067472/1‑69 (MQ=255)
     atGGGGCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAGGGCTTAgt                                                          >  1:3700502/1‑69 (MQ=255)
        gggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAGGGCt                                                              >  1:2327769/1‑62 (MQ=255)
        gggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAGGGCTTAgttgt                                                       >  1:2979706/1‑69 (MQ=255)
        gggCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAGGGCTTAgttg                                                        <  1:760277/68‑1 (MQ=255)
                 tCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAGTTGTTTTCAACTa                                              <  1:3373342/69‑1 (MQ=255)
                 tCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAGTTGTTTTCAACTa                                              <  1:3376719/69‑1 (MQ=255)
                        aggaagCACGTGATAAAAATAATAAAATAATTTATTCCTCAAGGCTTAGTTGTTTTCAACTACATTaaa                                       >  1:3107564/1‑69 (MQ=255)
                            agCACGTGATAAAAATAATAAAATAATTTATTCCTCAGGGCTTAGTTGTTTTCAACTACATTAAAACTc                                   <  1:1672196/69‑1 (MQ=255)
                            agCACGTGATAAAAATAATAAAATAATTTATTCCTCAGGGCTTAGTTGTTTTCAACTACATTAAAACTc                                   <  1:3101942/69‑1 (MQ=255)
                            agCACGTGATAAAAATAATAAAATAATTTATTCCTCAGGGCTTAGTTGTTTTAACTACATTAAAACTcc                                  <  1:3147617/69‑1 (MQ=255)
                             gCACGTGATAAAAATAATAAAATAATTTATTCCTCAGGGCTTAGTTGTTTTCAACTACATTAAAACTcc                                  <  1:1605658/69‑1 (MQ=255)
                                cGTGATAAAAATAATAAAATAATTTATTCCTCAGGGCTTAGTTGTTTTCAACTACATTAAAACTCCATg                               <  1:877650/69‑1 (MQ=255)
                                cGTGATAAAAATAATAAAATAATTTATTCCTCAGGGCTTAGTTGTTTTCAACTACATTAAAACTCCATg                               <  1:1006456/69‑1 (MQ=255)
                                    aTAAAAATAATAAAATAATTTATTCCTCAGGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATa                           >  1:1797345/1‑69 (MQ=255)
                                    aTAAAAATAATAAAATAATTTATTCCTCAAGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATa                           >  1:3178248/1‑69 (MQ=255)
                                          ataataaaataatTTATTCCTCAGGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACggtgg                     >  1:366390/1‑69 (MQ=255)
                                          ataataaaataatTTATTCCTCAGGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACggtgg                     <  1:3356414/69‑1 (MQ=255)
                                             ataaaataatTTATTCCTCAGGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACGGTGGTAt                  >  1:999338/1‑69 (MQ=255)
                                             ataaaataatTTATTCCTCAGGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACGGTGGTAt                  >  1:3020047/1‑69 (MQ=255)
                                                aaataatTTATTCCTCAGGGCTTAGTTGTTTTCAACTACATTAAAACt                                    >  1:3128579/1‑48 (MQ=255)
                                                  ataatTTATTCCTCAGGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACggtggt                    <  1:3423219/62‑1 (MQ=255)
                                                  ataatTTATTCCTCAGGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACggtggt                    <  1:231332/62‑1 (MQ=255)
                                                  ataatTTATTCCTCAGGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACggtggt                    <  1:1209672/62‑1 (MQ=255)
                                                              tCAGGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACGGTGGTATGTCTTTTTAATCCTGc  <  1:652575/68‑1 (MQ=255)
                                                              tCAGGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACGGTGGTATGTCTTTTTAATCCTGc  <  1:138982/68‑1 (MQ=255)
                                                               cAGGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACGGTGGTATGTCTTTTTAAt       >  1:2650972/1‑62 (MQ=255)
                                                               cAGGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACGGTGGTATGTCTTTTTAAt       >  1:1051378/1‑62 (MQ=255)
                                                                 |                                                                
GTATTATGGGGCTAAATTCCCTGAAGGAAGCACGTGATAAAAATAATAAAATAATTTATTCCTCAGGGCTTAGTTGTTTTCAACTACATTAAAACTCCATGAATACGGTGGTATGTCTTTTTAATCCTGC  >  W3110S.gb/4007474‑4007603

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: