Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,739,327 C→T 34.1% V195I (GTT→ATT)  yieH ← predicted hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,739,3270CT34.1% 96.7 / 78.4 82V195I (GTT→ATT) yieHpredicted hydrolase
Reads supporting (aligned to +/- strand):  ref base C (25/29);  new base T (11/17);  total (36/46)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

TCCCAACCACGCGCTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATACCAGATTGTG  >  W3110S.gb/3739257‑3739395
                                                                      |                                                                    
tCCCAACCACGCGCTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAAc                                                                      <  1:4954509/71‑1 (MQ=255)
    aaCCACGCGCTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAATGATc                                                                  <  1:3795774/71‑1 (MQ=255)
    aaCCACGCGCTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAATGATc                                                                  <  1:2026402/71‑1 (MQ=255)
         cgcgCTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCtt                                                             <  1:421774/71‑1 (MQ=255)
         cgcgCTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCt                                                              >  1:4347094/1‑70 (MQ=255)
          gcgcTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTa                                                            <  1:4018240/71‑1 (MQ=255)
           cgcTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAATGATCgg                                                                <  1:3165684/66‑1 (MQ=255)
           cgcTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAATGATCgg                                                                <  1:2061885/66‑1 (MQ=255)
            gcTTTCCACAGTTCATGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTAt                                                           <  1:3290003/70‑1 (MQ=255)
             cTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGAGTTTCGGGTGAATGATCGGCTTATtg                                                         >  1:4006821/1‑71 (MQ=255)
             cTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAATGATCGGCt                                                              >  1:4454205/1‑66 (MQ=255)
             cTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAATGATCGGCt                                                              >  1:555762/1‑66 (MQ=255)
             cTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAATGATCGGCTTATtg                                                         >  1:1674404/1‑71 (MQ=255)
               ttCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACg                                                                     <  1:2840624/57‑1 (MQ=255)
               ttCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACg                                                                     <  1:3359261/57‑1 (MQ=255)
                  cacaGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAATGATCGGCTTATTGTGCgg                                                    >  1:119845/1‑71 (MQ=255)
                    caGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTa                                                            >  1:611647/1‑61 (MQ=255)
                     aGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTa                                                            >  1:5001524/1‑60 (MQ=255)
                        tCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGt                                                  >  1:4522134/1‑67 (MQ=255)
                        tCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGt                                                  >  1:3502510/1‑67 (MQ=255)
                        tCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGgc                                              >  1:2377084/1‑71 (MQ=255)
                            gTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAATGATCGGCTTATTGTGCGGGTCGgcg                                             >  1:4404530/1‑68 (MQ=255)
                            gTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAATGATCGGCTTATTGTGCGGGTCGGCGCAg                                          <  1:2018421/71‑1 (MQ=255)
                            gTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGa                                                                    <  1:3806659/45‑1 (MQ=255)
                            gTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGa                                                                    <  1:806017/45‑1 (MQ=255)
                            gTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTAt                                                           <  1:4829394/54‑1 (MQ=255)
                              aaCTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAg                                          >  1:3355950/1‑69 (MQ=255)
                              aaCTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAg                                          >  1:4573915/1‑69 (MQ=255)
                              aaCTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAg                                          >  1:3229113/1‑69 (MQ=255)
                                cTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTAt                                                           <  1:2847052/50‑1 (MQ=255)
                                 tGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAATGATCGGCTTa                                                            <  1:980293/48‑1 (MQ=255)
                                  gCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGCACGATCGGCTTATTGTGCGGGTCGGCGCAGAAg                                       >  1:3175587/1‑68 (MQ=255)
                                  gCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAATGATCGGCt                                                              <  1:858317/45‑1 (MQ=255)
                                  gCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAATGATCGGCt                                                              <  1:4775077/45‑1 (MQ=255)
                                  gCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAATGATCGGCTTATTGTGCGGGTCGGCGCAGAAGt                                      >  1:3502148/1‑69 (MQ=255)
                                  gCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAATGATCGGCTTATTGTGCGGGTCGGCGCAGAAGt                                      >  1:4549430/1‑69 (MQ=255)
                                  gCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGc                                                      >  1:2051147/1‑53 (MQ=255)
                                  gCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGa                                         >  1:942303/1‑66 (MQ=255)
                                  gCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAg                                       >  1:236602/1‑68 (MQ=255)
                                  gCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAg                                       >  1:353443/1‑68 (MQ=255)
                                   cGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGa                                   <  1:1591870/71‑1 (MQ=255)
                                       aGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAacac                               <  1:2542534/71‑1 (MQ=255)
                                       aGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAacac                               <  1:2536651/71‑1 (MQ=255)
                                        gATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTAt                                                           <  1:3154281/42‑1 (MQ=255)
                                            ggTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTa                                                            <  1:2128616/37‑1 (MQ=255)
                                            ggTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCa                          <  1:3707515/71‑1 (MQ=255)
                                             gTAAAGGTGGTGACTTTCGGGTGAATGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCAt                         <  1:1135385/71‑1 (MQ=255)
                                             gTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCAt                         <  1:4660084/71‑1 (MQ=255)
                                              tAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATa                        <  1:3872767/71‑1 (MQ=255)
                                              tAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATa                        <  1:1054181/71‑1 (MQ=255)
                                               aaaGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGaa                                        >  1:3432435/1‑54 (MQ=255)
                                               aaaGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGaa                                        >  1:1615159/1‑54 (MQ=255)
                                               aaaGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATAc                       >  1:2903318/1‑71 (MQ=255)
                                                aaGGTGGTGACTTTCGGGTGAATGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCa                          <  1:1419969/67‑1 (MQ=255)
                                                aaGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATAcc                      <  1:964211/71‑1 (MQ=255)
                                                   gtggtgACTTTCGGGTGAATGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGc                   >  1:777503/1‑71 (MQ=255)
                                                    tggtgACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCg                  <  1:2897630/71‑1 (MQ=255)
                                                      gtgACTTTCGGGTGAATGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTc                <  1:2235753/71‑1 (MQ=255)
                                                      gtgACTTTCGGGTGAATGATCGGCTTATTGTGCGGGGCGGCGCAGAAGTAGAACACTTCCATACCTGCGTc                <  1:4695324/71‑1 (MQ=255)
                                                      gtgACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAg                                    <  1:2319152/51‑1 (MQ=255)
                                                       tgACTTTCGGGTGAATGATCGGCTTATTGTGCGGGTCGGCGCAg                                          <  1:1766542/44‑1 (MQ=255)
                                                       tgACTTTCGGGTGAATGATCGGCTTATTGTGCGGGTCGGCGCAg                                          <  1:4148757/44‑1 (MQ=255)
                                                       tgACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTc                <  1:1452236/70‑1 (MQ=255)
                                                       tgACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCg               <  1:1728924/71‑1 (MQ=255)
                                                         aCTTTCGGGTGAATGATCGGCTTATTGTGCGGGTCGGCGCAg                                          <  1:1778922/42‑1 (MQ=255)
                                                         aCTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCa                          >  1:1673123/1‑58 (MQ=255)
                                                            ttCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTc                            <  1:1281238/53‑1 (MQ=255)
                                                            ttCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTc                >  1:3553433/1‑65 (MQ=255)
                                                            ttCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTc                >  1:1332273/1‑65 (MQ=255)
                                                             tCGGGTGAATGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGa              >  1:3521379/1‑66 (MQ=255)
                                                             tCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATAc           >  1:4428038/1‑69 (MQ=255)
                                                             tCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATAc           <  1:3571571/69‑1 (MQ=255)
                                                             tCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATAc           >  1:877864/1‑69 (MQ=255)
                                                              cGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATACCa         <  1:978137/70‑1 (MQ=255)
                                                               gggTGAATGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATACCAg        <  1:1517006/70‑1 (MQ=255)
                                                               gggTGAATGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATACCAg        <  1:373491/70‑1 (MQ=255)
                                                                  tGAATGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATACCAGATtg    >  1:3309894/1‑71 (MQ=255)
                                                                  tGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACt                              <  1:1173335/45‑1 (MQ=255)
                                                                  tGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTc                            <  1:805866/47‑1 (MQ=255)
                                                                  tGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTc                            <  1:4537432/47‑1 (MQ=255)
                                                                   gAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCa                          >  1:4959605/1‑48 (MQ=255)
                                                                    aaCGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATACCAGATtgtg  >  1:4243228/1‑71 (MQ=255)
                                                                     aCGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATACCAGATtgtg  >  1:4358160/1‑70 (MQ=255)
                                                                      |                                                                    
TCCCAACCACGCGCTTTCCACAGTTCAGGTAACTGCGAAAGATGGGTAAAGGTGGTGACTTTCGGGTGAACGATCGGCTTATTGTGCGGGTCGGCGCAGAAGTAGAACACTTCCATACCTGCGTCGATACCAGATTGTG  >  W3110S.gb/3739257‑3739395

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: