Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,138,489 T→C 6.1% L44L (CTA→CTG hybG ← hydrogenase 2 accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,138,4890TC 6.1% 264.7 / 7.3 98L44L (CTA→CTGhybGhydrogenase 2 accessory protein
Reads supporting (aligned to +/- strand):  ref base T (41/51);  new base C (3/3);  total (44/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACC  >  W3110S.gb/3138421‑3138558
                                                                    |                                                                     
gctttggcttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAg                                                                     <  1:860066/71‑1 (MQ=255)
     ggcttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCgg                                                                >  1:2353541/1‑71 (MQ=255)
     ggcttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCg                                                                 <  1:2857780/70‑1 (MQ=255)
      gcttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCgg                                                                >  1:3641279/1‑70 (MQ=255)
      gcttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCgg                                                                <  1:3983981/70‑1 (MQ=255)
          cATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAgg                                                            <  1:2212757/70‑1 (MQ=255)
          cATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAg                                                             <  1:1652059/69‑1 (MQ=255)
            tCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGt                                                          >  1:4105579/1‑70 (MQ=255)
            tCTTCGTCGATGATGCTCATCGCAAATCCGACGTGTACAAGCACCCACTGGCCCAGt                                                                       >  1:424443/1‑57 (MQ=37)
            tCTTCATCGATAATGCTCATGGCAAATCCGACGTGTACCAGTACCCACTGGCCCAGCAGATCGGCAGGGtt                                                         <  1:467036/71‑1 (MQ=21)
            tCTTCATCGATAATGCTCATGGCAAATCCGACGTGTACCAGTACCCACTGGCCCAGCAGATCGGCAGGGtt                                                         <  1:1355092/71‑1 (MQ=21)
            tCTTCATCGATAATGCTCATGGCAAATCCGACGTGTACCAGTACCCACTGGCCCAGCAGATCGGCAGGGtt                                                         <  1:4308966/71‑1 (MQ=21)
             cTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTa                                                        >  1:497094/1‑71 (MQ=255)
             cTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTa                                                        >  1:1599728/1‑71 (MQ=255)
             cTTCGTCGATGATGCTCATGGCAAATCAGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTa                                                        >  1:242822/1‑71 (MQ=255)
              ttCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTAc                                                       >  1:478254/1‑71 (MQ=255)
               tcgtcgATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTAcc                                                      >  1:2239612/1‑71 (MQ=255)
                cgtcgATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTa                                                                      >  1:1744304/1‑54 (MQ=255)
                cgtcgATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAg                                                                     <  1:3739337/55‑1 (MQ=255)
                cgtcgATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGtt                                                         <  1:1178168/67‑1 (MQ=255)
                cgtcgATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTAcc                                                      <  1:1223302/70‑1 (MQ=255)
                cgtcgATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTAcc                                                      >  1:3122335/1‑70 (MQ=255)
                cgtcgATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCt                                                     >  1:1947101/1‑71 (MQ=255)
                cgtcgATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCt                                                     >  1:708247/1‑71 (MQ=255)
                  tcgATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTAcc                                                      >  1:1389057/1‑68 (MQ=255)
                  tcgATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCtt                                                    <  1:3119613/70‑1 (MQ=255)
                    gatgatGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTca                                                  >  1:3249373/1‑70 (MQ=255)
                    gatgatGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTca                                                  >  1:803374/1‑70 (MQ=255)
                     atgatgCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTcaca                                                <  1:1583976/71‑1 (MQ=255)
                      tgatgCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAATAGATCGGCAGGGTTACCTTCACaa                                               <  1:1907852/71‑1 (MQ=255)
                        atgCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTAc                                                       >  1:3910619/1‑61 (MQ=255)
                        atgCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTAc                                                       >  1:3752090/1‑61 (MQ=255)
                        atgCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAAt                                             <  1:198163/71‑1 (MQ=255)
                         tgCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATc                                            <  1:3787935/71‑1 (MQ=255)
                           cTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCa                                           <  1:2113654/70‑1 (MQ=255)
                            tCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAAt                                             <  1:2025544/67‑1 (MQ=255)
                            tCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAgg                                         >  1:748226/1‑71 (MQ=255)
                             cATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAggg                                        >  1:677695/1‑71 (MQ=255)
                               tGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGc                                       >  1:1497174/1‑70 (MQ=255)
                                 gCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGtt                                                         >  1:3091837/1‑50 (MQ=255)
                                 gCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGa                                     >  1:2852832/1‑70 (MQ=255)
                                 gCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGa                                     >  1:1487199/1‑70 (MQ=255)
                                  cAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTcaca                                                >  1:3238767/1‑58 (MQ=255)
                                  cAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGa                                     <  1:1631482/69‑1 (MQ=255)
                                    aaTCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATAtt                                 <  1:536129/71‑1 (MQ=255)
                                      tCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGtt                                                         >  1:2379263/1‑45 (MQ=255)
                                      tCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTca                               >  1:100895/1‑71 (MQ=255)
                                       ccGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTca                               >  1:4279016/1‑70 (MQ=255)
                                          aCGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATc                           >  1:1027541/1‑71 (MQ=255)
                                           cGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATcg                          <  1:1847214/71‑1 (MQ=255)
                                            gtgtACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATcgc                         <  1:810628/71‑1 (MQ=255)
                                              gtACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAg                                          >  1:636083/1‑52 (MQ=255)
                                              gtACAAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCt                                                     <  1:1926305/41‑1 (MQ=255)
                                                aCCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCtt                     >  1:2979917/1‑71 (MQ=255)
                                                 ccAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTcaca                                                <  1:1341938/43‑1 (MQ=255)
                                                 ccAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTcaca                                                <  1:1956119/43‑1 (MQ=255)
                                                 ccAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCa                                           >  1:4141507/1‑48 (MQ=255)
                                                 ccAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCtt                     <  1:3760229/70‑1 (MQ=255)
                                                 ccAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCtt                     <  1:947229/70‑1 (MQ=255)
                                                 ccAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCtt                     <  1:1849442/70‑1 (MQ=255)
                                                 ccAGCACCAACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCtt                     <  1:1620254/70‑1 (MQ=255)
                                                  cAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCg                                      <  1:164295/52‑1 (MQ=255)
                                                  cAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGa                   <  1:3556885/71‑1 (MQ=255)
                                                   aGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTca                               >  1:142697/1‑58 (MQ=255)
                                                    gCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATa                 >  1:3042103/1‑71 (MQ=255)
                                                    gCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATa                 >  1:4098924/1‑71 (MQ=255)
                                                    gCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATa                 >  1:2921762/1‑71 (MQ=255)
                                                     cACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATAtt                                 <  1:2906213/54‑1 (MQ=255)
                                                     cACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATAtt                                 <  1:3430790/54‑1 (MQ=255)
                                                     cACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATAc                <  1:1759037/71‑1 (MQ=255)
                                                      aCCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATAcc               <  1:2030045/71‑1 (MQ=255)
                                                      aCCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATAcc               <  1:2462637/71‑1 (MQ=255)
                                                       cccACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACca              <  1:2687333/71‑1 (MQ=255)
                                                          aCTGGCCCAGTAGATCGGCAGGGTTACCTTCACaa                                               <  1:777451/35‑1 (MQ=255)
                                                          aCTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGc                                       <  1:876950/43‑1 (MQ=255)
                                                          aCTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCt                      <  1:3111060/60‑1 (MQ=255)
                                                          aCTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGa                   >  1:1550663/1‑63 (MQ=255)
                                                          aCTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGa                   <  1:1305178/63‑1 (MQ=255)
                                                          aCTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGAt                  <  1:3588268/64‑1 (MQ=255)
                                                          aCTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACcaca            >  1:4166119/1‑70 (MQ=255)
                                                           cTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATa          <  1:3313256/71‑1 (MQ=255)
                                                           cTGGCCCAGTAGATCGGCAGGGGTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATa          <  1:3530362/71‑1 (MQ=255)
                                                            tGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGa                                     <  1:2694340/43‑1 (MQ=255)
                                                            tGGCCCAGCAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACGTCGCGCTTGATa                 >  1:2103495/1‑63 (MQ=255)
                                                            tGGCCCAGCAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACGTCGCGCTTGATa                 >  1:1301571/1‑63 (MQ=255)
                                                            tGGCCCAGCAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACGTCGCGCTTGATa                 >  1:3645771/1‑63 (MQ=255)
                                                             ggCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATa                 >  1:768372/1‑62 (MQ=255)
                                                             ggCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACt        >  1:2606488/1‑71 (MQ=255)
                                                              gCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATAc         <  1:3773239/69‑1 (MQ=255)
                                                               cccAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTc      <  1:2225348/71‑1 (MQ=255)
                                                               cccAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTc      <  1:2792404/71‑1 (MQ=255)
                                                                ccAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCa     <  1:1594181/71‑1 (MQ=255)
                                                                ccAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCa     <  1:149769/71‑1 (MQ=255)
                                                                ccAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCa     <  1:1063022/71‑1 (MQ=255)
                                                                  aGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAAc   >  1:287050/1‑71 (MQ=255)
                                                                   gTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGa                                     <  1:3805634/36‑1 (MQ=255)
                                                                   gTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTc                                <  1:2194774/41‑1 (MQ=255)
                                                                   gTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACcaca            >  1:71840/1‑61 (MQ=255)
                                                                   gTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAAcc  <  1:1693792/71‑1 (MQ=255)
                                                                    tAGATCGGCAGGGTTACCTTCACAAATCAGGGCGa                                     >  1:3630606/1‑35 (MQ=255)
                                                                    |                                                                     
GCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACC  >  W3110S.gb/3138421‑3138558

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: