Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,249,660 A→C 100% Y344* (TAT→TAG)  dhaH ← fused predicted dihydroxyacetone‑specific PTS enzyme HPr component and EI component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,249,6600AC85.8% 79.2 / 1.6 28Y344* (TAT→TAG) dhaHfused predicted dihydroxyacetone‑specific PTS enzyme HPr component and EI component
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base C (8/16);  total (11/17)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTC  >  W3110S.gb/1249601‑1249716
                                                           |                                                        
tatTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCat                                                         >  1:635837/1‑61 (MQ=255)
tatTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCat                                                         >  1:496980/1‑61 (MQ=255)
 atTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCCTATTCTGCCGTg                                              >  1:571733/1‑71 (MQ=255)
    catcatCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCCTATTCTGCCGTGCa                                            <  1:2804445/70‑1 (MQ=255)
    catcatCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCaa                                           >  1:1876121/1‑71 (MQ=255)
        atcCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCCTATTCTGCCGTGCAATGtt                                       <  1:1722187/71‑1 (MQ=255)
        atcCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCCTATTCTGCCGTGCAATGtt                                       <  1:2513828/71‑1 (MQ=255)
        atcCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCCTATTCTGCCGTGCAATGtt                                       <  1:2419996/71‑1 (MQ=255)
        atcCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCCTATTCTGCCGTGCAATGtt                                       <  1:2160833/71‑1 (MQ=255)
         tcCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCCTATTCTGCCGTGCAATGtt                                       <  1:1377135/70‑1 (MQ=255)
           cAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCAt                                    >  1:1741148/1‑71 (MQ=255)
                gCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCCTATTCTGCCGTGCAATGTTCATg                                   >  1:319566/1‑67 (MQ=255)
                 cTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGa                               <  1:1061129/70‑1 (MQ=255)
                  tGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCCTATTCTGCCGTGCAATGTTCa                                     <  1:878707/63‑1 (MQ=255)
                    gTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCCTATTCTGCCGTGCAATGTTCATGCTGAAGGa                           <  1:1874956/71‑1 (MQ=255)
                     tATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCCTATTCTGCCGTGCAATGtt                                       <  1:2423519/58‑1 (MQ=255)
                       ttGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCCTATTCTGCCGTGCAATGTTCAt                                    <  1:1860300/59‑1 (MQ=255)
                         gCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCCTATTCTGCCGTGCAATGTTCa                                     >  1:221134/1‑56 (MQ=255)
                         gCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCCTATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTc                       <  1:2601467/70‑1 (MQ=255)
                         gCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCa                                     >  1:2573853/1‑56 (MQ=255)
                                 aGTTCTTTAAGAACTTGCTGCCAGGCCTATTCTGCCGTGCAATGTTCa                                     >  1:1590063/1‑48 (MQ=255)
                                 aGTTCTTTAAGAACTTGCTGCCAGGCCTATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTgcc                <  1:724262/69‑1 (MQ=255)
                                   ttCTTTAAGAACTTGCTGCCAGGCCTATTCTGCCGTGCa                                            <  1:2026028/39‑1 (MQ=255)
                                   ttCTTTAAGAACTTGCTGCCAGGCCTATTCTGCCGTGCAATGTTCAt                                    <  1:870642/47‑1 (MQ=255)
                                   ttCTTTAAGAACTTGCTGCCAGGCCTATTCTGCCGTGCAATCTTCATGCTGAAg                             <  1:1855439/54‑1 (MQ=255)
                                     cTTTAAGAACTTGCTGCCAGGCCTATTCTGCCTTGCAATGTTCATGCTGAAGGa                           >  1:1542199/1‑54 (MQ=255)
                                     cTTTAAGAACTTGCTGCCAGGCCTATTCTGCCGTGCAATGTTCATGCTGAAGGa                           >  1:1770315/1‑54 (MQ=255)
                                     cTTTAAGAACTTGCTGCCAGGCCTATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTgccgccgc           <  1:321046/70‑1 (MQ=255)
                                                tgctgcCAGGCCTATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTc  >  1:2446250/1‑68 (MQ=255)
                                                tgctgcCAGGCCTATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTc  >  1:1179052/1‑68 (MQ=255)
                                                           |                                                        
TATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTC  >  W3110S.gb/1249601‑1249716

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