Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 4,155,147 +A coding (782/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.A85.7% 124.9 / 19.7 49coding (783/870 nt)prkBpredicted phosphoribulokinase
Reads supporting (aligned to +/- strand):  ref base . (6/1);  new base A (23/19);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCCAGCAGCCAGGGGAAATCG  >  W3110S.gb/4155091‑4155209
                                                        |                                                               
gCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGa                                                    <  1:2409174/70‑1 (MQ=255)
                 attaattCCATTGCCAGACCCATTTTGGCGCCCGGTACCAACTAACGTATTGATGTGGGa                                             <  1:1413494/60‑1 (MQ=255)
                 attaattCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGa                                             <  1:705734/60‑1 (MQ=255)
                 attaattCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGa                                             <  1:50144/60‑1 (MQ=255)
                 attaattCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGcc                                 >  1:3575909/1‑71 (MQ=255)
                 attaattCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGcc                                 >  1:3202360/1‑71 (MQ=255)
                 attaattCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGcc                                 >  1:2339539/1‑71 (MQ=255)
                   taattCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGaaa                                           >  1:1119027/1‑60 (MQ=255)
                   taattCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGaaa                                           >  1:3619951/1‑60 (MQ=255)
                   taattCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGaaa                                           >  1:2050843/1‑60 (MQ=255)
                    aattCCATTGCCAGACCCATTTTGCCTCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCt                                <  1:582501/70‑1 (MQ=255)
                    aattCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCt                                <  1:500502/70‑1 (MQ=255)
                    aattCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCt                                <  1:4007617/70‑1 (MQ=255)
                    aattCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCt                                <  1:1334307/70‑1 (MQ=255)
                    aattCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCt                                <  1:2283685/70‑1 (MQ=255)
                    aattCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCt                                <  1:3403899/70‑1 (MQ=255)
                    aattCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTAATGATGTGGGAAATGAATGAGCCt                                <  1:1685739/70‑1 (MQ=255)
                     attCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCt                                <  1:2029384/69‑1 (MQ=255)
                      ttCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGc                             <  1:3507575/71‑1 (MQ=255)
                      ttCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGc                             <  1:1941337/71‑1 (MQ=255)
                           ttGCCAGACCCATTTTGCCGCCCGGTACCAACTAAc                                                           >  1:1389844/1‑36 (MQ=37)
                           ttGCCAGACCCATTTTGCCGCCCGGTACCAACTAAGGTATTGATGTGGGAAATGAATGAGCCTTGCAACAt                        >  1:564716/1‑71 (MQ=255)
                           ttGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAAc                          >  1:3351476/1‑69 (MQ=255)
                           ttGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACa                         >  1:419458/1‑70 (MQ=255)
                           ttGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACAt                        >  1:1051918/1‑71 (MQ=255)
                           ttGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACAt                        >  1:1743599/1‑71 (MQ=255)
                           ttGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACAt                        >  1:1139396/1‑71 (MQ=255)
                           ttGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACAt                        >  1:451458/1‑71 (MQ=255)
                           ttGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACAt                        >  1:326547/1‑71 (MQ=255)
                               cAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGcc                    <  1:560091/71‑1 (MQ=255)
                                  aCCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGATTGAGCCTTGCAACATCGCcag                  >  1:3972043/1‑70 (MQ=255)
                                  aCCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCcag                  >  1:339028/1‑70 (MQ=255)
                                  aCCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCcag                  >  1:1472342/1‑70 (MQ=255)
                                  aCCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCcag                  >  1:2883878/1‑70 (MQ=255)
                                  aCCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCcag                  >  1:1366085/1‑70 (MQ=255)
                                  aCCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCcag                  >  1:1824025/1‑70 (MQ=255)
                                  aCCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCcag                  >  1:4016784/1‑70 (MQ=255)
                                  aCCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCcag                  >  1:1291055/1‑70 (MQ=255)
                                  aCCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCcag                  >  1:3583410/1‑70 (MQ=255)
                                  aCCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCcag                  >  1:1657884/1‑70 (MQ=255)
                                  aCCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCcag                  >  1:1684801/1‑70 (MQ=255)
                                       ttttGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGc                                  <  1:668319/49‑1 (MQ=39)
                                         ttGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGATTGAGCCTTGCat                           >  1:3896121/1‑52 (MQ=255)
                                         ttGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCaa                           >  1:2286318/1‑53 (MQ=255)
                                         ttGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCaa                           >  1:306052/1‑53 (MQ=255)
                                            ccgccCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCCAGCAGCCAGGGGa       <  1:3995763/71‑1 (MQ=255)
                                            ccgccCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCCAGCAGCCAGGGGa       <  1:3311553/71‑1 (MQ=255)
                                            ccgccCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCCAGCAGCCAGGGGa       <  1:3565396/71‑1 (MQ=255)
                                                ccGGTACCACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCCAGCAGCCAGGGGAAATCg  <  1:3613623/71‑1 (MQ=255)
                                                        |                                                               
GCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCCAGCAGCCAGGGGAAATCG  >  W3110S.gb/4155091‑4155209

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: