Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 810,589 T→G 84.1% M275R (ATG→AGG)  bioB → biotin synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb810,5890TG84.1% 69.0 / 2.4 25M275R (ATG→AGG) bioBbiotin synthase
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base G (8/13);  total (10/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 = 9.41e-01

AACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGA  >  W3110S.gb/810530‑810648
                                                           |                                                           
aaCCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTa                                                     <  1:127090/68‑1 (MQ=255)
aaCCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTa                                                     <  1:2146451/68‑1 (MQ=255)
aaCCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTa                                                     <  1:2604344/68‑1 (MQ=255)
aaCCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTAt                                                    <  1:3007515/69‑1 (MQ=255)
        aCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGgc                                             >  1:2506815/1‑68 (MQ=255)
            gcgcCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCtt                                                       >  1:1882143/1‑54 (MQ=255)
            gcgcCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCtt                                                       >  1:3664210/1‑54 (MQ=255)
                cTTTCTGCCGGAGGCTATCAGATGAACGAACAGACTGAGGCGATGTGCtt                                                       <  1:3302699/50‑1 (MQ=37)
                   tCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTAt                                                    >  1:1009114/1‑50 (MQ=255)
                   tCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGAt                                  >  1:365348/1‑68 (MQ=255)
                    cTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGAttt                                >  1:256456/1‑69 (MQ=255)
                         ggACGCGAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTc                              <  1:2795513/66‑1 (MQ=255)
                         ggACGCGAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTa                            <  1:825212/68‑1 (MQ=255)
                            cgcgAGCAGATGAACGAACAGACTCAGGCGATGTGCtt                                                       <  1:707133/38‑1 (MQ=255)
                                 gCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCaa                   <  1:1052613/69‑1 (MQ=255)
                                 gCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCaa                   <  1:371915/69‑1 (MQ=255)
                                    gATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAAct                >  1:2369794/1‑69 (MQ=255)
                                    gATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAAct                >  1:820340/1‑69 (MQ=255)
                                       gaacgaacAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTc                                     >  1:474221/1‑45 (MQ=255)
                                       gaacgaacAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTc                                     >  1:831657/1‑45 (MQ=255)
                                       gaacgaacAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCa                    <  1:1630526/62‑1 (MQ=255)
                                       gaacgaacAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCa                    <  1:1485639/62‑1 (MQ=255)
                                              cAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCAc       <  1:2103405/68‑1 (MQ=255)
                                              cAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCAc       <  1:1789677/68‑1 (MQ=255)
                                                  cTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGa  <  1:2384925/69‑1 (MQ=255)
                                                           |                                                           
AACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGA  >  W3110S.gb/810530‑810648

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: