Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,155,146:1 +A 23.4% coding (783/870 nt) prkB ← predicted phosphoribulokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,155,1461.A23.4% 105.9 / 45.5 47coding (783/870 nt)prkBpredicted phosphoribulokinase
Reads supporting (aligned to +/- strand):  ref base . (9/27);  new base A (5/6);  total (14/33)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

AATCGTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCCAGCAGCCAGGGGAAAT  >  W3110S.gb/4155084‑4155207
                                                               |                                                             
aaTCGTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACg                                                        <  1:2200207/70‑1 (MQ=255)
aaTCGTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACg                                                        <  1:1781522/70‑1 (MQ=255)
aaTCGTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACg                                                        <  1:1305306/70‑1 (MQ=255)
    gTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTAtt                                                    <  1:2611828/71‑1 (MQ=255)
    gTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTAtt                                                    <  1:1147243/71‑1 (MQ=255)
    gTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTAtt                                                    <  1:3115468/71‑1 (MQ=255)
    gTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTAtt                                                    <  1:3149094/71‑1 (MQ=255)
    gTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTAtt                                                    <  1:1070145/70‑1 (MQ=255)
    gTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTAtt                                                    <  1:280912/70‑1 (MQ=255)
    gTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTAtt                                                    <  1:3050156/70‑1 (MQ=255)
    gTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTAtt                                                    <  1:1237777/70‑1 (MQ=255)
    gTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTAtt                                                    <  1:483325/70‑1 (MQ=255)
    gTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTAtt                                                    <  1:1177152/70‑1 (MQ=255)
    gTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTAtt                                                    <  1:800172/70‑1 (MQ=255)
    gTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTAtt                                                    <  1:957039/70‑1 (MQ=255)
     ttGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTAtt                                                    <  1:45124/70‑1 (MQ=255)
     ttGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTAtt                                                    <  1:218085/70‑1 (MQ=255)
     ttGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTAtt                                                    <  1:1946945/69‑1 (MQ=255)
                gcagcaTAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTTGGAaatt                                       <  1:2437199/71‑2 (MQ=255)
                gcagcaTAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGTGAaatg                                       <  1:1916034/71‑1 (MQ=255)
                gcagcaTAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAaatg                                       <  1:3263087/71‑1 (MQ=255)
                gcagcaTAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAaatg                                       <  1:3026114/71‑1 (MQ=255)
                gcagcaTAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAaatg                                       <  1:512471/71‑1 (MQ=255)
                gcagcaTAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAaatg                                       <  1:2453003/71‑1 (MQ=255)
                gcagcaTAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAaatg                                       <  1:251856/71‑1 (MQ=255)
                gcagcaTAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAaatg                                       <  1:1296804/71‑1 (MQ=255)
                gcagcaTAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAaatg                                       <  1:1388452/71‑1 (MQ=255)
                 cagcaTAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGt                                                       <  1:3030038/54‑1 (MQ=255)
                 cagcaTAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGt                                                       <  1:1150544/54‑1 (MQ=255)
                 cagcaTAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGt                                                       <  1:2756926/54‑1 (MQ=255)
                 cagcaTAATTAATTCCATTGCCAGACCCATTTTGCCGCCAGGTACCACTAACGt                                                       <  1:3305102/54‑1 (MQ=255)
                     aTAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAaatgaat                                    >  1:374258/1‑69 (MQ=255)
                           aattCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGCCTTg                            >  1:2653473/1‑71 (MQ=255)
                           aattCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGACTTg                            >  1:1262876/1‑71 (MQ=255)
                            attCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGCCTTGc                           <  1:199251/71‑1 (MQ=255)
                             ttacATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTg                            >  1:11614/4‑70 (MQ=255)
                             ttCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAGCGTATTGATGTGGGAAATGAATGAGCCTTg                            >  1:3020148/1‑70 (MQ=255)
                             ttCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTg                            >  1:1205760/1‑70 (MQ=255)
                             ttCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTg                            >  1:1385478/1‑70 (MQ=255)
                                cATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACa                       <  1:2188058/71‑1 (MQ=255)
                                    gCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGc                   >  1:623210/1‑71 (MQ=255)
                                    gCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGc                   >  1:1073862/1‑71 (MQ=255)
                                    gCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGc                   >  1:3147952/1‑71 (MQ=255)
                                    gCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGc                   >  1:2004712/1‑71 (MQ=255)
                                    gCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGc                   >  1:2917868/1‑71 (MQ=255)
                                    gCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGc                   >  1:2305346/1‑71 (MQ=255)
                                                      cccGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCCAGCAGCCAGGGGAAAt  >  1:3007594/1‑71 (MQ=255)
                                                               |                                                             
AATCGTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCCAGCAGCCAGGGGAAAT  >  W3110S.gb/4155084‑4155207

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