Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,527,342 C→T 26.4% E479K (GAG→AAG)  yjhG ← predicted dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,527,3420CT26.4% 160.7 / 62.8 91E479K (GAG→AAG) yjhGpredicted dehydratase
Reads supporting (aligned to +/- strand):  ref base C (31/36);  new base T (14/10);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

CGCCAATAATGACCAGAATATCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGTC  >  W3110S.gb/4527273‑4527408
                                                                     |                                                                  
cGCCAATAATGACCAGAATATCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCg                                                                   <  1:531364/71‑1 (MQ=255)
cGCCAATAATGACCAGAATATCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCg                                                                   <  1:1818770/71‑1 (MQ=255)
cGCCAATAATGACCAGAATATCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTATCTCg                                                                   <  1:2776481/71‑1 (MQ=255)
   cAATAATGACCAGAATATCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGa                                                                 <  1:752990/70‑1 (MQ=255)
       aatGACCAGAATATCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGa                                                             >  1:1074704/1‑70 (MQ=255)
        atGACCAGAATATCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTtctc                                                                    >  1:473742/1‑62 (MQ=255)
        atGACCAGAATATCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGAt                                                            >  1:2066699/1‑70 (MQ=255)
         tGACCAGAATATCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTTGGACAGATaa                                                          <  1:512816/71‑1 (MQ=255)
         tGACCAGAATATCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTTGGACAGATaa                                                          <  1:298322/71‑1 (MQ=255)
         tGACCAGAATATCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATaa                                                          <  1:4032461/71‑1 (MQ=255)
         tGACCAGAATATCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATaa                                                          <  1:2479057/71‑1 (MQ=255)
         tGACCAGAATATCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATaa                                                          <  1:200019/71‑1 (MQ=255)
          gACCAGAATATCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATaaa                                                         >  1:884550/1‑71 (MQ=255)
                aaTATCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCt                                                      <  1:1747864/68‑1 (MQ=255)
                 atatCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTTGGACAGATAAACCTTCGc                                                  >  1:353731/1‑71 (MQ=255)
                 atatCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGa                                                                 >  1:1286569/1‑56 (MQ=255)
                 atatCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGa                                                                 >  1:2463568/1‑56 (MQ=255)
                  tatCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTTGGACAGATAAACCTTCGc                                                  <  1:3234690/70‑1 (MQ=255)
                  tatCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGc                                                  >  1:2177611/1‑70 (MQ=255)
                    tCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCcac                                               <  1:3542901/71‑1 (MQ=255)
                    tCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCca                                                >  1:357309/1‑70 (MQ=255)
                    tCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCca                                                >  1:2943053/1‑70 (MQ=255)
                     cgcccgccTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCcaca                                              >  1:4065738/1‑71 (MQ=255)
                      gcccgccTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTTGGACAGATAAACCTTCGCcaca                                              <  1:3432740/70‑1 (MQ=255)
                      gcccgccTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTTGGACAGATAAACCTTCGCcaca                                              <  1:247566/70‑1 (MQ=255)
                      gcccgccTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTTGGACAGATAAACCTTCGCCacac                                             <  1:1969288/71‑1 (MQ=255)
                      gcccgccTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATaa                                                          >  1:320447/1‑58 (MQ=255)
                         cgccTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTTGGACAGATAAACCTTCGCCACACCtt                                          <  1:707499/71‑1 (MQ=255)
                         cgccTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCtt                                          >  1:3527200/1‑71 (MQ=255)
                         cgccTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCtt                                          <  1:1517576/71‑1 (MQ=255)
                         cgccTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCtt                                          <  1:2637492/71‑1 (MQ=255)
                          gccTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTTGGACAGATAAACCTTCGCCACACCttt                                         <  1:788910/71‑1 (MQ=255)
                            cTTGATCTTGTCATGTTTGATATCTTAAATCGCACTTTTCTCGGACAGATaa                                                          <  1:3345677/52‑1 (MQ=255)
                            cTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATaa                                                          <  1:44246/52‑1 (MQ=255)
                               gATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATgg                                     >  1:424324/1‑70 (MQ=255)
                                aTCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTTGGACAGATAAACCTTCGCCACACCTTTATGGt                                    <  1:2886267/70‑1 (MQ=255)
                                aTCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTTGGACAGATAAACCTTCGCCACACCTTTATGGt                                    <  1:1564171/70‑1 (MQ=255)
                                 tCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTTGGACAGATAAACCTTCGCca                                                >  1:1380066/1‑57 (MQ=255)
                                     gTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATa                                <  1:3434874/69‑1 (MQ=255)
                                     gTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATa                                <  1:1343493/69‑1 (MQ=255)
                                     gTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATa                                <  1:3544911/69‑1 (MQ=255)
                                     gTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATa                                <  1:241253/69‑1 (MQ=255)
                                     gTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATa                                <  1:431495/69‑1 (MQ=255)
                                     gTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATa                                <  1:1664461/69‑1 (MQ=255)
                                      tCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGa                              <  1:2280931/70‑1 (MQ=255)
                                      tCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGa                              <  1:2276285/70‑1 (MQ=255)
                                        aTGTTTGATATCGTAAATCGCACTTTTCTTGGACAGATAAACCTTCGCCACACCTTTa                                        >  1:601360/1‑58 (MQ=255)
                                          gTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATgg                                     <  1:2777484/59‑1 (MQ=255)
                                          gTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGt                                    <  1:55386/60‑1 (MQ=255)
                                          gTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCt                         >  1:1107289/1‑71 (MQ=255)
                                           tttGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCtt                        <  1:882638/71‑1 (MQ=255)
                                            ttGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTa                                   >  1:214974/1‑59 (MQ=255)
                                            ttGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTa                                   >  1:3732772/1‑59 (MQ=255)
                                            ttGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTa                                   >  1:238387/1‑59 (MQ=255)
                                               atatCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCt                     <  1:1979866/70‑1 (MQ=255)
                                                tatCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTa                                        <  1:2224157/50‑1 (MQ=255)
                                                tatCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCtca                   >  1:1277963/1‑71 (MQ=255)
                                                tatCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCtca                   <  1:670180/71‑1 (MQ=255)
                                                    gTAAATCGCACTTTTCTTGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCatcaa               >  1:1738320/1‑71 (MQ=255)
                                                    gTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTa                                        >  1:575817/1‑46 (MQ=255)
                                                    gTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTa                                        >  1:87813/1‑46 (MQ=255)
                                                    gTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTa                                        >  1:3101790/1‑46 (MQ=255)
                                                    gTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATa                            >  1:3930310/1‑58 (MQ=255)
                                                    gTAAATCGCACTTTTCTCGGACAGATAAACATTCGCCACAACTTTATGGTAATAGATACCTTGCTCatcaa               >  1:2024610/1‑71 (MQ=255)
                                                       aaTCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCatcaatc             <  1:1954579/70‑1 (MQ=255)
                                                         tCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAatcatc          <  1:1379958/71‑1 (MQ=255)
                                                            cACTTTTCTTGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATAcc                          >  1:426076/1‑52 (MQ=255)
                                                            cACTTTTCTTGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAg       >  1:186238/1‑71 (MQ=255)
                                                            cACTTTTCTTGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAg       >  1:2025888/1‑71 (MQ=255)
                                                              cTTTTCTCGGACAGATAAACCTTCGCCACACGTTTATg                                      >  1:1427196/1‑38 (MQ=255)
                                                              cTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGt                                    <  1:1253929/40‑1 (MQ=255)
                                                              cTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGt                                    <  1:2570108/40‑1 (MQ=255)
                                                              cTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAggg     >  1:1541001/1‑71 (MQ=255)
                                                              cTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAggg     >  1:3111489/1‑71 (MQ=255)
                                                              cTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAggg     <  1:2966709/71‑1 (MQ=255)
                                                              cTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTACTCATCAATCATCGAggg     >  1:1921820/1‑71 (MQ=255)
                                                              cTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTACTCATCAATCATCGAggg     >  1:2042122/1‑71 (MQ=255)
                                                               ttttCTCGGACAGATAAACCTTCGCCACACCTTTATGGt                                    <  1:1026825/39‑1 (MQ=255)
                                                               ttttCTCGGACAGATAAACCTTCGCCACACCTTTATGGt                                    <  1:2727121/39‑1 (MQ=255)
                                                               ttttCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGAt                             <  1:2912128/46‑1 (MQ=255)
                                                               ttttCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCtt                        <  1:2388750/51‑1 (MQ=255)
                                                                tttCTTGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGt   >  1:1737226/1‑71 (MQ=255)
                                                                tttCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAg       >  1:2126539/1‑67 (MQ=255)
                                                                tttCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAg       >  1:1586221/1‑67 (MQ=255)
                                                                 ttCTTGGACAGATAAACCTTCGCCACACCTTTATGGTaa                                  >  1:644109/1‑39 (MQ=255)
                                                                 ttCTTGGACAGATAAACCTTCGCCACACCTTTATGGTaa                                  >  1:2244656/1‑39 (MQ=255)
                                                                 ttCTTGGACAGATAAACCTTCGCCACACCTTTATGGTaa                                  >  1:1680932/1‑39 (MQ=255)
                                                                 ttCTTGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGt   >  1:664009/1‑70 (MQ=255)
                                                                 ttCTTGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGt   >  1:1482047/1‑70 (MQ=255)
                                                                 ttCTTGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGt   >  1:2044663/1‑70 (MQ=255)
                                                                 ttCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGt   >  1:3562819/1‑70 (MQ=255)
                                                                 ttCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGt   >  1:707263/1‑70 (MQ=255)
                                                                  tctcGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGTc  <  1:3192698/70‑1 (MQ=255)
                                                                     |                                                                  
CGCCAATAATGACCAGAATATCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGTC  >  W3110S.gb/4527273‑4527408

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