Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,155,146:1 +A 60.6% 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.A60.6% 45.5 / 47.4 33coding (783/870 nt)prkBpredicted phosphoribulokinase
Reads supporting (aligned to +/- strand):  ref base . (2/11);  new base A (15/5);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 1.27e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

CTTCCATCAATCGTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCCAGCAGCCAGGG  >  W3110S.gb/4155076‑4155202
                                                                       |                                                        
cttccATCAATCGTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTAcc                                                           <  1:448732/71‑1 (MQ=255)
       cAATCGTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAAc                                                    >  1:715258/1‑71 (MQ=255)
        aaTCGTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACg                                                   >  1:1185112/1‑71 (MQ=255)
        aaTCGTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACg                                                   >  1:660099/1‑71 (MQ=255)
        aaTCGTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACg                                                   >  1:1445463/1‑71 (MQ=255)
        aaTCGTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACg                                                   >  1:416366/1‑71 (MQ=255)
                        gcagcaTAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAaatg                                  <  1:1458070/71‑1 (MQ=255)
                        gcagcaTAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAaatg                                  <  1:3597383/71‑1 (MQ=255)
                        gcagcaTAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAaatg                                  <  1:3247065/71‑1 (MQ=255)
                                attaattCCATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGc                           <  1:860926/71‑1 (MQ=255)
                                        cATTGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCaa                    <  1:2780283/70‑1 (MQ=255)
                                         aTTGCCAGACCCATTTTGCCGCCCGGTACCACTAACg                                                   <  1:2225539/37‑1 (MQ=255)
                                         aTTGCCAGACCCATTTTGCCGCCCGGTACCACTAACg                                                   <  1:2530494/37‑1 (MQ=255)
                                         aTTGCCAGACCCATTTTGCCGCCCGGTACCACTAACg                                                   <  1:2776301/37‑1 (MQ=255)
                                         aTTGCCAGACCCATTTTGCCGCCCGGTACCACTAACg                                                   <  1:3665501/37‑1 (MQ=255)
                                         aTTGCCAGACCCATTTTGCCGCCCGGTACCACTAACg                                                   <  1:1743395/37‑1 (MQ=255)
                                          ttGCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTg                                              <  1:723203/42‑1 (MQ=38)
                                          ttGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATc                <  1:2599589/71‑1 (MQ=255)
                                          ttGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATc                <  1:1205200/71‑1 (MQ=255)
                                            gCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCg               >  1:3070956/1‑71 (MQ=255)
                                            gCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCg               >  1:1105664/1‑71 (MQ=255)
                                            gCCAGACCCATTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCg               >  1:2799682/1‑71 (MQ=255)
                                                  cccATTTTGCCGCCCGGTACCTCTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATc                <  1:294007/63‑1 (MQ=255)
                                                     aTTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGcc                          >  1:1216364/1‑51 (MQ=255)
                                                     aTTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGcc                          >  1:11573/1‑51 (MQ=255)
                                                     aTTTTGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGcc                          >  1:2139322/1‑51 (MQ=255)
                                                     aTTTTGCCGCCCCGTACCAACTAACGTAATGATGTGGGAAATGAATGAGcc                          >  1:2968821/1‑51 (MQ=255)
                                                      ttttGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCCAGCAGCCa     <  1:3598475/71‑1 (MQ=255)
                                                      ttttGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCCAGCAGCCa     <  1:3605418/71‑1 (MQ=255)
                                                        ttGCCGCCCGGTGCCACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCCAGCAGCCAggg  >  1:1916094/1‑71 (MQ=255)
                                                        ttGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCCAGCAGCCAgg   >  1:2955862/1‑71 (MQ=255)
                                                        ttGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCCAGCAGCCAgg   >  1:1675318/1‑71 (MQ=255)
                                                        ttGCCGCCCGGTACCAACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCCAGCAGCCAgg   >  1:1013627/1‑71 (MQ=255)
                                                        ttGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCCAGCAGCCAggg  >  1:1578459/1‑71 (MQ=255)
                                                                       |                                                        
CTTCCATCAATCGTTGCACCAGCGGCAGCATAATTAATTCCATTGCCAGACCCATTTTGCCGCCCGGTACCACTAACGTATTGATGTGGGAAATGAATGAGCCTTGCAACATCGCCAGCAGCCAGGG  >  W3110S.gb/4155076‑4155202

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: