Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,138,506 C→T 5.8% G39S (GGT→AGT)  hybG ← hydrogenase 2 accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,138,5060CT 5.8% 219.0 / 3.0 69G39S (GGT→AGT) hybGhydrogenase 2 accessory protein
Reads supporting (aligned to +/- strand):  ref base C (42/23);  new base T (2/2);  total (44/25)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

ATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTG  >  W3110S.gb/3138445‑3138572
                                                             |                                                                  
atgCTCATGGCAAATCCGACGTTTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACaa                                                             <  1:882099/69‑1 (MQ=255)
  gCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACaaa                                                            <  1:1700157/68‑1 (MQ=255)
         gCAAATCCGACGTGTACCAGCACCCACTGGCCCATTAGATCGGCAGGGTTACCTTCACAAATCAGGGc                                                     >  1:1107192/1‑68 (MQ=255)
         gCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACTTTCACAGATCAGGGc                                                     >  1:1076026/1‑68 (MQ=255)
         gCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACTTTCACAAATCAGGGc                                                     >  1:913097/1‑68 (MQ=255)
         gCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGc                                                     >  1:2404356/1‑68 (MQ=255)
         gCAAATCCGACGTGTACCAGCACCAACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGc                                                     >  1:732053/1‑68 (MQ=255)
              tCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATAtt                                               >  1:1995101/1‑69 (MQ=255)
              tCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATAtt                                               >  1:1109193/1‑69 (MQ=255)
              tCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATAtt                                               >  1:1884281/1‑69 (MQ=255)
              tCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATAtt                                               >  1:1341354/1‑69 (MQ=255)
              tCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATAtt                                               >  1:452152/1‑69 (MQ=255)
              tCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATAtt                                               >  1:308085/1‑69 (MQ=255)
                       tACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCg                                                    >  1:1697330/1‑55 (MQ=255)
                       tACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGa                                                   <  1:3250484/56‑1 (MQ=255)
                       tACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGa                                                   <  1:2124845/56‑1 (MQ=255)
                       tACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGa                                                   <  1:510795/56‑1 (MQ=255)
                       tACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGa                                                   <  1:434893/56‑1 (MQ=255)
                       tACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATcgcg                                      <  1:377890/69‑1 (MQ=255)
                       tACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATcgcg                                      >  1:2136071/1‑69 (MQ=255)
                       tACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATcgcg                                      >  1:142983/1‑69 (MQ=255)
                       tACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCCATATTCACATcgcg                                      <  1:3575915/69‑1 (MQ=255)
                         ccAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCgcgc                                     >  1:3557816/1‑68 (MQ=255)
                         ccAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCgcgc                                     >  1:3384901/1‑68 (MQ=255)
                         ccAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCgcgc                                     >  1:1988623/1‑68 (MQ=255)
                         ccAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCgcgc                                     >  1:1178193/1‑68 (MQ=255)
                           aGCACCCACTGGCCCAGTAGATCGGCAGGGTTACTTTCACAAATCAGGGCGATATTCACATCGCGCtt                                   <  1:1505167/68‑1 (MQ=255)
                           aGCACCCACTGGCCCAGTAGATCGGCAGGGTTACTTTCACAAATCAGGGCGATATTCACATCGCGCtt                                   <  1:1100472/68‑1 (MQ=255)
                           aGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCtt                                   <  1:1099225/68‑1 (MQ=255)
                           aGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCtt                                   <  1:364767/68‑1 (MQ=255)
                           aGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTg                                  <  1:3606188/69‑1 (MQ=255)
                                ccACTGGCCCAGTAGATCGGCAGGGTTACCTTCACaaa                                                            >  1:221470/1‑38 (MQ=255)
                                ccACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATAc                              >  1:1175177/1‑68 (MQ=255)
                                ccACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATAc                              >  1:2946365/1‑68 (MQ=255)
                                   cTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACcaca                          <  1:187075/69‑1 (MQ=255)
                                     ggCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATa                        >  1:3208561/1‑69 (MQ=255)
                                     ggCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATa                        >  1:2367700/1‑69 (MQ=255)
                                           gTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGa                                                   >  1:64220/1‑36 (MQ=255)
                                           gTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCtt                                   <  1:942589/52‑1 (MQ=255)
                                           gTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCtt                                   <  1:1570243/52‑1 (MQ=255)
                                           gTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCtt                                   <  1:2371652/52‑1 (MQ=255)
                                              gATCGTCAGGGTTACCTTCACAAATCAGGGCGATAtt                                               >  1:1626140/1‑37 (MQ=255)
                                                 cGGCAGGGTTACCTTCACAAATCAGGGCGATATTca                                             <  1:797283/36‑1 (MQ=255)
                                                 cGGCAGGGTTACCTTCACAAATCAGGGCGATATTca                                             <  1:1072600/36‑1 (MQ=255)
                                                 cGGCAGGGTTACCTTCACAAATCAGGGCGATATTca                                             <  1:3265819/36‑1 (MQ=255)
                                                 cGGCAGGGTTACCTTCACAAATCAGGGCGATATTca                                             <  1:2043403/36‑1 (MQ=255)
                                                     aGGGTTACCTTCACAAATCAGGGCGATATTCACATcg                                        >  1:809703/1‑37 (MQ=255)
                                                     aGGGTTACCTTCACAAATCAGGGCGATATTCACATc                                         >  1:2037480/1‑36 (MQ=255)
                                                     aGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCtt                                   >  1:1659142/1‑42 (MQ=255)
                                                     aGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACAt                         >  1:1882673/1‑52 (MQ=255)
                                                     aGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACAt                         <  1:2487943/52‑1 (MQ=255)
                                                     aGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACAt                         <  1:1243595/52‑1 (MQ=255)
                                                     aGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACAt                         <  1:272636/52‑1 (MQ=255)
                                                     aGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACAt                         >  1:79714/1‑52 (MQ=255)
                                                     aGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACAt                         <  1:3479399/52‑1 (MQ=255)
                                                        gTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTg     >  1:3145579/1‑69 (MQ=255)
                                                        gTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTg     >  1:2747722/1‑69 (MQ=255)
                                                         ttACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTg     >  1:2900536/1‑68 (MQ=255)
                                                         ttACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTg     >  1:1525846/1‑68 (MQ=255)
                                                         ttACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTg     >  1:558931/1‑68 (MQ=255)
                                                         ttACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTg     >  1:1315373/1‑68 (MQ=255)
                                                          tACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTgg    >  1:628772/1‑68 (MQ=255)
                                                          tACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTgg    >  1:125850/1‑68 (MQ=255)
                                                          tACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTgg    >  1:2625471/1‑68 (MQ=255)
                                                           aCCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACca                            >  1:1816494/1‑43 (MQ=255)
                                                           aCCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACca                            >  1:2962014/1‑43 (MQ=255)
                                                           aCCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTg  >  1:22121/1‑69 (MQ=255)
                                                           aCCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTg  >  1:562974/1‑69 (MQ=255)
                                                           aCCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTg  >  1:1624126/1‑69 (MQ=255)
                                                            ccTTCACAAATCAGTGCGATATTCACATCGCGCTTGa                                 >  1:1782346/1‑37 (MQ=255)
                                                            ccTTCACAAATCAGGGCGATATTCACATCGCGCTTGa                                 >  1:1522606/1‑37 (MQ=255)
                                                             |                                                                  
ATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTG  >  W3110S.gb/3138445‑3138572

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: