Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 287 C→T intergenic (+32/‑50) thrL → / → thrA thr operon leader peptide/fused aspartokinase I and homoserine dehydrogenase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2870CT100.0% 131.8 / NA 40intergenic (+32/‑50)thrL/thrAthr operon leader peptide/fused aspartokinase I and homoserine dehydrogenase I
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/13);  total (27/13)

TTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAA  >  W3110S.gb/230‑331
                                                         |                                            
ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAg                                       >  1:439690/1‑65 (MQ=255)
ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAg                                       >  1:2546911/1‑65 (MQ=255)
ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAg                                       >  1:1508514/1‑65 (MQ=255)
ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAg                                       >  1:1430600/1‑65 (MQ=255)
      caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCggg                                 >  1:1183347/1‑65 (MQ=255)
      caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCggg                                 >  1:953626/1‑65 (MQ=255)
      caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCggg                                 >  1:802810/1‑65 (MQ=255)
      caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCggg                                 >  1:589232/1‑65 (MQ=255)
      caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCggg                                 >  1:2678861/1‑65 (MQ=255)
      caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCggg                                 >  1:2469718/1‑65 (MQ=255)
      caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCggg                                 >  1:239794/1‑65 (MQ=255)
      caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCggg                                 >  1:1295962/1‑65 (MQ=255)
      caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCggg                                 >  1:1308522/1‑65 (MQ=255)
      caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCggg                                 >  1:1417743/1‑65 (MQ=255)
      caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCgg                                  >  1:430212/1‑64 (MQ=255)
      caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCg                                   >  1:1125132/1‑63 (MQ=255)
      caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCg                                   >  1:736400/1‑63 (MQ=255)
               gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCtttttttt                        >  1:2142439/1‑65 (MQ=255)
               gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCtttttttt                        >  1:2253127/1‑65 (MQ=255)
               gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCtttttttt                        >  1:832672/1‑65 (MQ=255)
               gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCtttttttt                        >  1:1410595/1‑65 (MQ=255)
               gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCtttttttt                        >  1:1710278/1‑65 (MQ=255)
               gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCtttttttt                        >  1:2717056/1‑65 (MQ=255)
               gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCtttttttt                        >  1:851805/1‑65 (MQ=255)
               gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCtttttttt                        >  1:453220/1‑65 (MQ=255)
               gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCttcttttt                        >  1:2516791/1‑65 (MQ=255)
                     gCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGAc                   >  1:491470/1‑64 (MQ=255)
                                 aGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGa       <  1:1367739/64‑1 (MQ=255)
                                 aGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGa       <  1:1944291/64‑1 (MQ=255)
                                  ggAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCACAGGTAACGAgg     <  1:2428908/65‑1 (MQ=39)
                                  ggAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCAAATGTAACGAgg     <  1:2902680/65‑1 (MQ=39)
                                  ggAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAgg     <  1:2204621/65‑1 (MQ=255)
                                  ggAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAgg     <  1:1286013/65‑1 (MQ=255)
                                  ggAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAgg     <  1:1599546/65‑1 (MQ=255)
                                  ggAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAgg     <  1:2804994/65‑1 (MQ=255)
                                  ggAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTCGACCAAAGGTAACGAgg     <  1:1326534/64‑1 (MQ=38)
                                     aaCACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaa  <  1:1265998/65‑1 (MQ=255)
                                     aaCACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaa  <  1:2302765/65‑1 (MQ=255)
                                     aaCACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaa  <  1:1645211/65‑1 (MQ=255)
                                     aaCACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaa  <  1:1510272/65‑1 (MQ=255)
                                                         |                                            
TTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAA  >  W3110S.gb/230‑331

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: