Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,138,531 A→G 7.6% D30D (GAT→GAC hybG ← hydrogenase 2 accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,138,5310AG 7.6% 204.8 / 7.6 79D30D (GAT→GAChybGhydrogenase 2 accessory protein
Reads supporting (aligned to +/- strand):  ref base A (37/36);  new base G (4/2);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

ACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGC  >  W3110S.gb/3138463‑3138601
                                                                    |                                                                      
aCGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATc                                                                      >  1:1027541/1‑71 (MQ=255)
 cGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATcg                                                                     <  1:1847214/71‑1 (MQ=255)
  gtgtACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATcgc                                                                    <  1:810628/71‑1 (MQ=255)
      aCCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCtt                                                                >  1:2979917/1‑71 (MQ=255)
       ccAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCtt                                                                <  1:947229/70‑1 (MQ=255)
       ccAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCtt                                                                <  1:1849442/70‑1 (MQ=255)
       ccAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCtt                                                                <  1:3760229/70‑1 (MQ=255)
       ccAGCACCAACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCtt                                                                <  1:1620254/70‑1 (MQ=255)
        cAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGa                                                              <  1:3556885/71‑1 (MQ=255)
          gCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATa                                                            >  1:2921762/1‑71 (MQ=255)
          gCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATa                                                            >  1:3042103/1‑71 (MQ=255)
          gCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATa                                                            >  1:4098924/1‑71 (MQ=255)
           cACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATAc                                                           <  1:1759037/71‑1 (MQ=255)
            aCCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATAcc                                                          <  1:2462637/71‑1 (MQ=255)
            aCCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATAcc                                                          <  1:2030045/71‑1 (MQ=255)
             cccACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACca                                                         <  1:2687333/71‑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)
                  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:2792404/71‑1 (MQ=255)
                     cccAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTc                                                 <  1:2225348/71‑1 (MQ=255)
                      ccAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCa                                                <  1:1063022/71‑1 (MQ=255)
                      ccAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCa                                                <  1:1594181/71‑1 (MQ=255)
                      ccAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCa                                                <  1:149769/71‑1 (MQ=255)
                        aGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAAc                                              >  1:287050/1‑71 (MQ=255)
                         gTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACcaca                                                       >  1:71840/1‑61 (MQ=255)
                         gTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAAcc                                             <  1:1693792/71‑1 (MQ=255)
                           aGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGa                                                              >  1:3889436/1‑52 (MQ=255)
                           aGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTg                                           >  1:3436210/1‑71 (MQ=255)
                            gATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTgc                                          >  1:1589200/1‑71 (MQ=255)
                            gATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTgc                                          >  1:3006558/1‑71 (MQ=255)
                              tCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACGTCGCGCTTGATAc                                                           <  1:4358470/52‑1 (MQ=255)
                              tCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACGTCGCGCTTGATAc                                                           <  1:1984610/52‑1 (MQ=255)
                               cGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCtt                                                                >  1:1303616/1‑46 (MQ=255)
                               cGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTgcgc                                        >  1:3795777/1‑70 (MQ=255)
                               cGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTgcgc                                        >  1:4298517/1‑70 (MQ=255)
                                 gCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGAGCTTGATACAACATACTTCAACCTGCGCAAg                                     >  1:1123189/1‑71 (MQ=255)
                                   aGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCt                                   >  1:2951182/1‑71 (MQ=255)
                                     ggTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTgg                                 >  1:74415/1‑71 (MQ=255)
                                       ttACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATa                                                            >  1:911763/1‑42 (MQ=255)
                                       ttACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGt                                >  1:1419533/1‑70 (MQ=255)
                                       ttACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTg                               <  1:1435218/71‑1 (MQ=255)
                                       ttACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTg                               <  1:3196394/71‑1 (MQ=255)
                                       ttACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTg                               <  1:155523/71‑1 (MQ=255)
                                        tACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACAt                                                      >  1:3090696/1‑47 (MQ=255)
                                        tACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCaa                                      >  1:2988212/1‑63 (MQ=255)
                                         aCCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGaa                             >  1:3905552/1‑71 (MQ=255)
                                         aCCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGaa                             >  1:3200541/1‑71 (MQ=255)
                                         aCCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGaa                             >  1:2682641/1‑71 (MQ=255)
                                          ccTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAat                            <  1:2336721/71‑1 (MQ=255)
                                           cTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCa                                                <  1:457608/50‑1 (MQ=255)
                                           cTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCa                                                <  1:215396/50‑1 (MQ=255)
                                             tCACAAATCAGGGCGATATTCACATCGCGCTTGATAcc                                                          <  1:550682/38‑1 (MQ=255)
                                             tCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATc                         >  1:3622529/1‑71 (MQ=255)
                                             tCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATc                         >  1:753498/1‑71 (MQ=255)
                                              cacaAATCAGGGCGATATTCACATCGCGCTTGATACCACATa                                                     >  1:1644644/1‑42 (MQ=255)
                                              cacaAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCt                        >  1:320352/1‑71 (MQ=255)
                                              cacaAATCAGGGCGATATTCACATAGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCt                        >  1:363319/1‑71 (MQ=255)
                                              cacaAATCAGGGCGATATTAACATCGCGCTTGATACCACATa                                                     >  1:2137781/1‑42 (MQ=255)
                                                caAATCAGGGCGATATTCACGTCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATc                         >  1:4072999/1‑68 (MQ=255)
                                                caAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTc                      <  1:716120/71‑1 (MQ=255)
                                                 aaaTCAGGGCGATATTCACATCGCGCTTGATACCACATACTTc                                                 <  1:3069476/43‑1 (MQ=255)
                                                 aaaTCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCa                     <  1:2132263/71‑1 (MQ=255)
                                                   aTCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCAcc                   <  1:2105862/71‑1 (MQ=255)
                                                     cAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAAcc                                             >  1:3509651/1‑43 (MQ=255)
                                                      aGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGa                 >  1:152652/1‑70 (MQ=255)
                                                         gCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTg                               <  1:3317039/53‑1 (MQ=255)
                                                               ttCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGAc                <  1:57465/62‑1 (MQ=255)
                                                                tCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAatat                          >  1:2168989/1‑51 (MQ=255)
                                                                 cacaTCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCAcct     >  1:2377763/1‑71 (MQ=255)
                                                                 cacaTCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCAcct     >  1:3820169/1‑71 (MQ=255)
                                                                 cacaTCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCAcct     >  1:3717595/1‑71 (MQ=255)
                                                                 cacaTCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCAcc      <  1:1637892/70‑1 (MQ=255)
                                                                  acaTCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCAcctg    <  1:690219/71‑1 (MQ=255)
                                                                  acaTCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCAcctg    <  1:2478150/71‑1 (MQ=255)
                                                                    aTCGCGCTTGATACCACATACTTCAACCTTCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCAcctggc  <  1:3641506/71‑1 (MQ=255)
                                                                    aTCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTc                      <  1:1198070/51‑1 (MQ=255)
                                                                    |                                                                      
ACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGC  >  W3110S.gb/3138463‑3138601

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: