Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,066,112 C→T 5.1% V209M (GTG→ATG)  yggE ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,066,1120CT 5.1% 257.7 / 2.5 78V209M (GTG→ATG) yggEconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (46/28);  new base T (2/2);  total (48/30)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

AACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCAC  >  W3110S.gb/3066042‑3066181
                                                                      |                                                                     
aaCTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGacac                                                                       >  1:5949065/1‑71 (MQ=255)
 aCTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACAcc                                                                      >  1:5447762/1‑71 (MQ=255)
   tGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCgg                                                                    >  1:2167098/1‑71 (MQ=255)
     gAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGgcg                                                                  >  1:4483357/1‑71 (MQ=255)
          ccACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCa                                                             <  1:138852/71‑1 (MQ=255)
             caTCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGAGACCGGCGCGGCATCg                                                          >  1:2684668/1‑71 (MQ=255)
             caTCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCg                                                          >  1:7354735/1‑71 (MQ=255)
             caTCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCg                                                          >  1:3708566/1‑71 (MQ=255)
              aTCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCgg                                                         >  1:1729876/1‑71 (MQ=255)
               tCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATc                                                           >  1:7167843/1‑68 (MQ=255)
               tCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGc                                                        <  1:536930/71‑1 (MQ=255)
                cGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCt                                                       >  1:1874731/1‑71 (MQ=255)
                cGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCt                                                       >  1:452080/1‑71 (MQ=255)
                     tGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCttt                                                     <  1:3033223/68‑1 (MQ=255)
                     tGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTtc                                                    <  1:7320798/69‑1 (MQ=255)
                      gATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCTGACACCGGCGCGGCATCGGCTTtc                                                    <  1:6039943/68‑1 (MQ=255)
                       atcatcAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCttt                                                     >  1:5397582/1‑66 (MQ=255)
                       atcatcAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTtcatc                                                 <  1:4046131/70‑1 (MQ=255)
                        tcatcaAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTcatcat                                               >  1:6514024/1‑71 (MQ=255)
                         catcaAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTtcatca                                                <  1:6848660/69‑1 (MQ=255)
                         catcaAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCatcatc                                              >  1:3576331/1‑71 (MQ=255)
                          atcaAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACATCGGCGCGGCATCGGCt                                                       <  1:2450602/61‑1 (MQ=255)
                          atcaAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACATCGGCGCGGCATCGGCt                                                       <  1:1895903/61‑1 (MQ=255)
                          atcaAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATcc                                             <  1:310311/71‑1 (MQ=255)
                           tcaAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGc                                                        <  1:4595001/59‑1 (MQ=255)
                            caAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGacac                                                                       >  1:6894521/1‑43 (MQ=255)
                            caAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACATCGGCGCGGCATCGGCttt                                                     >  1:5889973/1‑61 (MQ=255)
                            caAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACATCGGCGCGGCATCGGCttt                                                     >  1:4082449/1‑61 (MQ=255)
                              aaCTGAATAGCGGCCTGCTCGTATGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCAc                                         <  1:4757592/71‑1 (MQ=255)
                              aaCTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCAt                                                            <  1:6918812/52‑1 (MQ=255)
                              aaCTGAATAGCGGCCTGCTCGTAAGTTTCATGGGCGGACACCGGCGCGGCAt                                                            <  1:1804362/52‑1 (MQ=255)
                                 tGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCAt                                      >  1:2355072/1‑71 (MQ=255)
                                  gAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCAtt                                     >  1:3288532/1‑71 (MQ=255)
                                  gAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCAtt                                     >  1:4328047/1‑71 (MQ=255)
                                    aTAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATtgg                                   >  1:7683322/1‑71 (MQ=255)
                                     tAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATtggg                                  <  1:5767750/71‑1 (MQ=255)
                                      aGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATtgggc                                 >  1:242075/1‑71 (MQ=255)
                                      aGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATtgggc                                 >  1:4367993/1‑71 (MQ=255)
                                       gCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCttt                                                     >  1:5875236/1‑50 (MQ=255)
                                       gCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCttt                                                     >  1:3231404/1‑50 (MQ=255)
                                       gCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATtgggct                                >  1:1427028/1‑71 (MQ=255)
                                        cGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATtgggctg                               <  1:37517/71‑1 (MQ=255)
                                          gCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATtgggctgg                              <  1:2692694/70‑1 (MQ=255)
                                           ccTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATtgggctggg                             >  1:889246/1‑70 (MQ=255)
                                              gCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCgg                                                         <  1:3628511/39‑1 (MQ=255)
                                               cTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATc                                                           >  1:5295292/1‑36 (MQ=255)
                                               cTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCgg                                                         <  1:4953917/38‑1 (MQ=255)
                                               cTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGa                         >  1:6636344/1‑70 (MQ=255)
                                                  gTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATtg                                    >  1:2637943/1‑56 (MQ=255)
                                                  gTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGt                     >  1:3331620/1‑71 (MQ=255)
                                                   tAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATcc                                             >  1:3212297/1‑46 (MQ=255)
                                                   tAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCAt                                      >  1:5023387/1‑53 (MQ=255)
                                                   tAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGt                     >  1:3488566/1‑70 (MQ=255)
                                                   tAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGt                     <  1:531436/70‑1 (MQ=255)
                                                    aaGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTg                   >  1:7644070/1‑71 (MQ=255)
                                                    aaGTTTCCTGGGCGGACACCGCAGCGGCATCGGCTTTCATCATCCGc                                           >  1:3453240/1‑47 (MQ=39)
                                                    aaGTTTCCTGG‑‑GGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGtt                    >  1:329318/1‑68 (MQ=255)
                                                     aGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTtca                                                   >  1:5768664/1‑38 (MQ=255)
                                                     aGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTtc                                                    >  1:3931305/1‑37 (MQ=255)
                                                     aGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTg                   >  1:1434307/1‑70 (MQ=255)
                                                     aGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTg                   >  1:67108/1‑70 (MQ=255)
                                                      gTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGa                 <  1:6377449/71‑1 (MQ=255)
                                                       tttCCTGGGCGGACACCGGTGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATa                       <  1:6445514/64‑1 (MQ=255)
                                                       tttCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATa                       <  1:3823751/64‑1 (MQ=255)
                                                        ttCCTGGGCGGACACCGGCGCGGCATCGGCTTTcatcat                                               >  1:7397183/1‑39 (MQ=255)
                                                        ttCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGaaa               >  1:7463104/1‑71 (MQ=255)
                                                        ttCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGa                 <  1:484046/69‑1 (MQ=255)
                                                         tCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGaaa               >  1:948689/1‑70 (MQ=255)
                                                         tCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGaaa               >  1:2404069/1‑70 (MQ=255)
                                                         tCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGaaa               >  1:2021206/1‑70 (MQ=255)
                                                          ccTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAAc              <  1:7762477/70‑1 (MQ=255)
                                                           cTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATcc                                             >  1:806671/1‑38 (MQ=255)
                                                           cTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATcc                                             >  1:108098/1‑38 (MQ=255)
                                                           cTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAAc              <  1:294461/69‑1 (MQ=255)
                                                             gggCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCAATGGGCTGGGCTGATAGTTgg                  >  1:7148665/1‑63 (MQ=255)
                                                              ggCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATgg          <  1:3619193/70‑1 (MQ=255)
                                                               gCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGt                     <  1:6117803/58‑1 (MQ=255)
                                                               gCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGt                     <  1:4573813/58‑1 (MQ=255)
                                                               gCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTgg                  >  1:3573719/1‑61 (MQ=255)
                                                                 ggACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATg           <  1:2083924/66‑1 (MQ=255)
                                                                  gACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAgcg     >  1:1223150/1‑71 (MQ=255)
                                                                     acCGGCGCGGCATCGGCTTTCATCATCCGCACCATtgggctgggc                            <  1:948562/45‑1 (MQ=255)
                                                                     acCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATa                       <  1:268568/50‑1 (MQ=255)
                                                                     acCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCAc  >  1:3791740/1‑71 (MQ=255)
                                                                      |                                                                     
AACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCAC  >  W3110S.gb/3066042‑3066181

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: