Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,902,804 C→A 64.3% intergenic (+774/+599) ygcE → / ← ygcF predicted kinase/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,902,8040CA64.3% 2.7 / 34.1 28intergenic (+774/+599)ygcE/ygcFpredicted kinase/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (6/4);  new base A (11/7);  total (17/11)
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

GCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGC  >  W3110S.gb/2902742‑2902869
                                                              |                                                                 
gCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCgg                                                           >  1:2548044/1‑71 (MQ=255)
    gcgcGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGATGGCGCgggg                                                         <  1:4515203/69‑1 (MQ=255)
    gcgcGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGATGGCGCgggg                                                         <  1:565064/69‑1 (MQ=255)
    gcgcGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGATGGCGCgggg                                                         <  1:571088/69‑1 (MQ=255)
       cgGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAacac                                                    <  1:2292442/71‑1 (MQ=255)
            aaCACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACgg                                                  >  1:1846478/1‑68 (MQ=255)
            aaCACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGgcg                                               <  1:4715414/71‑1 (MQ=255)
             acacCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGATGGCGCGGGGAACACGGgcgc                                              <  1:1583386/71‑1 (MQ=255)
                ccTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACgg                                                  >  1:3317915/1‑64 (MQ=255)
                    ttttCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGaa                                        >  1:2053455/1‑70 (MQ=255)
                    ttttCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGATGGCGCGGGGaa                                                       >  1:3852853/1‑55 (MQ=255)
                     tttCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGATGGCGCGGGGAACACggg                                                 >  1:4898218/1‑60 (MQ=255)
                      ttCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGATGGCGCGGGGAACACgg                                                  >  1:2506830/1‑58 (MQ=255)
                      ttCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGATGGCGCGGGGAACACGGGCGCACGGAATa                                      >  1:988065/1‑70 (MQ=255)
                        cGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGATGGCGCggg                                                          <  1:1422832/48‑1 (MQ=255)
                            aaTCTATGGACTATTGCTATTCGGTTTATCCCCGATGGCGCGGGGAACACGGGCGCACGGAATACAAAGcc                               >  1:3528975/1‑71 (MQ=255)
                                    gACTATTGCTATTCGGTTTATCCCCGATGGCGCGGGGAaca                                                     >  1:4951673/1‑41 (MQ=255)
                                      cTATTGCTATTCGGTTTATCCCCGATGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGc                     >  1:1628719/1‑71 (MQ=255)
                                       tATTGCTATTCGGTTTATCCCCGATGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCt                    >  1:859214/1‑71 (MQ=255)
                                        aTTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGgcg                                               <  1:4885843/43‑1 (MQ=255)
                                         ttGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGc                                <  1:4117152/57‑1 (MQ=255)
                                          tGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCgg                 >  1:1059053/1‑71 (MQ=255)
                                               ttCGGTTTATCCCCGATGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCgt                             <  1:4946004/54‑1 (MQ=255)
                                               ttCGGTTTATCCCCGATGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCgt                             <  1:1026899/54‑1 (MQ=255)
                                                 cGGTTTATCCCCGATGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTAt                         >  1:3542086/1‑56 (MQ=255)
                                                 cGGTTTATCCCCGATGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATcc          >  1:3127517/1‑71 (MQ=255)
                                                 cGGTTTATCCCCGATGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATcc          >  1:2314450/1‑71 (MQ=255)
                                                         ccccGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTgga  >  1:4282570/1‑70 (MQ=255)
                                                              |                                                                 
GCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGC  >  W3110S.gb/2902742‑2902869

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: