Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,264,229 C→T 28.3% L172L (CTG→CTA fruA ← fused fructose‑specific PTS enzyme IIB'BC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,264,2290CT28.3% 85.1 / 36.7 53L172L (CTG→CTAfruAfused fructose‑specific PTS enzyme IIB'BC components
Reads supporting (aligned to +/- strand):  ref base C (20/18);  new base T (9/6);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

TGCGCGGTTTTCTTCAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGT  >  W3110S.gb/2264164‑2264291
                                                                 |                                                              
tGCGCGGTTTTCTTCAGCGCCAGACCGGTAGAGGTACTATACATCGGTTTACCAGCAAATTTCGCTAg                                                              <  1:504509/68‑1 (MQ=255)
    cgGTTTTCTTCAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCTAGATc                                                           <  1:3313342/67‑1 (MQ=255)
    cgGTTTTCTTCAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCTAGATc                                                           <  1:1253396/67‑1 (MQ=255)
       ttttCTTCAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCTAGATCCACt                                                       <  1:1230023/68‑1 (MQ=255)
        tttCTTCAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAGATCCACTTc                                                     >  1:3541376/1‑69 (MQ=255)
         ttcttcAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAg                                                              >  1:118838/1‑59 (MQ=255)
         ttcttcAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAg                                                              >  1:3597383/1‑59 (MQ=255)
         ttcttcAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAg                                                              >  1:2128426/1‑59 (MQ=255)
         ttcttcAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAg                                                              >  1:3722890/1‑59 (MQ=255)
            ttcAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCTAGATCCACTTCGat                                                  <  1:740843/68‑1 (MQ=255)
                    cAGACCGGTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAGATCCACTTGGATATCTgccg                                          <  1:964632/68‑1 (MQ=255)
                          ggTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGata                                                 >  1:406797/1‑55 (MQ=255)
                          ggTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGata                                                 <  1:2374017/55‑1 (MQ=255)
                          ggTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGata                                                 >  1:1048026/1‑55 (MQ=255)
                            tAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATc                                  >  1:320646/1‑68 (MQ=255)
                            tAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATc                                  >  1:288984/1‑68 (MQ=255)
                                ggTACGATACATCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATc                                  >  1:1243971/1‑64 (MQ=255)
                                ggTACGATACATCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATc                                  >  1:3430980/1‑64 (MQ=255)
                                ggTACGATACATCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATc                                  >  1:3212394/1‑64 (MQ=255)
                                   aCGATACATCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGATATc                                               >  1:2291587/1‑48 (MQ=255)
                                   aCGATACATCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATc                          >  1:1112475/1‑69 (MQ=255)
                                   aCGATACATCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATc                          >  1:3126288/1‑69 (MQ=255)
                                    cGATACATCGGTTTACCAGCAAATTTCGCTAGATCCACTTCGATATCTGCCGCCACAATCACCAGATcc                         >  1:120285/1‑69 (MQ=255)
                                    cGATACATCGGTTTACCAGCAAATTTCGCTAGATCCACTTCGATATCTGCCGCCACAATCACCAGATcc                         >  1:809541/1‑69 (MQ=255)
                                    cGATACATCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATcc                         <  1:1040790/69‑1 (MQ=255)
                                      aTACATCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCgc                       >  1:496929/1‑69 (MQ=255)
                                          aTCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCcac                                      <  1:53037/50‑1 (MQ=255)
                                          aTCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCcac                                      <  1:432736/50‑1 (MQ=255)
                                          aTCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCcac                                      <  1:726426/50‑1 (MQ=255)
                                          aTCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCcac                                      <  1:2465943/50‑1 (MQ=255)
                                          aTCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCcac                                      <  1:249172/50‑1 (MQ=255)
                                           tCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCgctgct                   >  1:3365248/1‑68 (MQ=255)
                                            cGGTTTACCAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCtgctgc                 <  1:512391/69‑1 (MQ=255)
                                             ggTTTACCAGCAAATTTCGCTAGATCCACTTTGata                                                 >  1:2269579/1‑36 (MQ=37)
                                             ggTTTACCAGCAAATTTCGCTAGATCCACTTGGata                                                 >  1:633051/1‑36 (MQ=37)
                                             ggTTTACCAGCAAATTTCGCTAGATCCACTTCGata                                                 >  1:2737634/1‑36 (MQ=255)
                                             ggTTTACCAGCAAATTTCGCTAGATCCACTTCGata                                                 >  1:3644739/1‑36 (MQ=255)
                                             ggTTTACCAGCAAATTTCGCTAGATCCACTTCGata                                                 >  1:1084387/1‑36 (MQ=255)
                                             ggTTTACCAGCAAATTTCGCTAGATCCACTTCGATATCTGCCGCCACAATCACCAg                             >  1:2597912/1‑56 (MQ=255)
                                             ggTTTACCAGCAAATTTCGCTAGATCCACTTCGATATCTGCCGCCACAATCACCAg                             >  1:385826/1‑56 (MQ=255)
                                             ggTTTACCAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCtgctgc                 <  1:1029983/68‑1 (MQ=255)
                                                ttACCAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGAct             <  1:2065716/69‑1 (MQ=255)
                                                ttACCAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGAct             <  1:1249165/69‑1 (MQ=255)
                                                     aGCAAATTTCGCTAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGActtctt         <  1:3050423/68‑1 (MQ=255)
                                                           tttCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGt  <  1:1256536/69‑1 (MQ=255)
                                                           tttCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGt  <  1:2901592/69‑1 (MQ=255)
                                                           tttCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGt  <  1:3072822/69‑1 (MQ=255)
                                                            ttCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGAct             <  1:578421/57‑1 (MQ=255)
                                                            ttCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGAct             <  1:1432226/57‑1 (MQ=255)
                                                            ttCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGAct             >  1:287002/1‑57 (MQ=255)
                                                            ttCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGAct             <  1:81341/57‑1 (MQ=255)
                                                            ttCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGAct             >  1:894621/1‑57 (MQ=255)
                                                            ttCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGAct             >  1:3000252/1‑57 (MQ=255)
                                                                 |                                                              
TGCGCGGTTTTCTTCAGCGCCAGACCGGTAGAGGTACGATACATCGGTTTACCAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGT  >  W3110S.gb/2264164‑2264291

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