Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,841,701 C→T 48.8% L134L (CTG→CTA hycH ← protein required for maturation of hydrogenase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,841,7010CT48.8% 6.6 / 120.9 82L134L (CTG→CTAhycHprotein required for maturation of hydrogenase 3
Reads supporting (aligned to +/- strand):  ref base C (17/25);  new base T (21/19);  total (38/44)
Fisher's exact test for biased strand distribution p-value = 3.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

GCCAGCAGCGGACCTGCGCCATCATCGCCCATCATGCTATTGCCAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATGCAGCATACTTAACAGCATTGTA  >  W3110S.gb/2841632‑2841770
                                                                     |                                                                     
gCCAGCAGCGGACCTGCGCCATCATCGCCCATCATGCTATTGCCAACACAGAGTAAAACGTCAGTCACGCa                                                                      <  1:1842003/71‑1 (MQ=255)
   agcagcGGACCTGCGCCATCATCGCCCATCATGCTATTGCCAACACAGAGTAAAACGTCAGTCACGTAGtc                                                                   >  1:2179136/1‑71 (MQ=255)
               gcgcCATCATCGCCCATCATGCTATTGCCAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATc                                                        >  1:534589/1‑70 (MQ=255)
                cgcCATCATCGCCCATCATGCTATTGCCAACACAGAGTAAAACGTCAGTCACGTAGTCTCCTCACCATCa                                                       <  1:1055145/70‑1 (MQ=255)
                cgcCATCATCGCCCATCATGCTATTGCCAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCa                                                       <  1:1323873/70‑1 (MQ=255)
                 gcCATCATCGCCCATCATGCTATTGCCAACACAGAGTAAAACGTCAGTCACGTAGTCTCCTCACCATCaaa                                                     <  1:1162027/71‑1 (MQ=255)
                 gcCATCATCGCCCATCATGCTATTGCCAACACAGAGTAAAACGTCAGTCACGTAGTCTCCTCACCATCaaa                                                     <  1:597857/71‑1 (MQ=255)
                 gcCATCATCGCCCATCATGCTATTGCCAACACAGAGTAAAACGTCAGTCACGTAGTCTCCTCACCATCaaa                                                     <  1:932542/71‑1 (MQ=255)
                   catcatCGCCCATCATGCTATTGCCAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCAAATa                                                   <  1:1178582/71‑1 (MQ=255)
                   catcatCGCCCATCATGCTATTGCCAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCAAATa                                                   <  1:960766/71‑1 (MQ=255)
                      catcGCCCATCATGCTATTGCCTACACAGAGTAAAACGTCAGTCACGTAGTCTCCTCACCATCaaa                                                     <  1:1611963/66‑1 (MQ=255)
                      catcGCCCATCATGCTATTGCCAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCAAATAGa                                                 >  1:420109/1‑70 (MQ=255)
                       atcGCCCATCATGCTATTGCCAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCAc                                                            >  1:559600/1‑58 (MQ=255)
                       atcGCCCATCATGCTATTGCCAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCAAATAGATg                                               >  1:786917/1‑71 (MQ=255)
                       atcGCCCATCATGCTATTGCCAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCAAATAGATg                                               >  1:2347433/1‑71 (MQ=255)
                       atcGCCCATCATGCTATTGCCAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCAAATAGATg                                               >  1:2725606/1‑71 (MQ=255)
                             catcatGCTATTGCTAACACAGAGTAAAACGTCAGTCACGTAGTctcc                                                                <  1:1004390/48‑1 (MQ=255)
                             catcatGCTATTGCCAACACAGAGTAAAACGTCAGTCACGTAGTctcc                                                                <  1:33338/48‑1 (MQ=255)
                             catcatGCTATTGCCAACACAGAGTAAAACGTCAGTCACGTAGTctcc                                                                <  1:1954444/48‑1 (MQ=255)
                             catcatGCTATTGCCAACACAGAGTAAAACGTCAGTCACGTAGTCTCCTCACCATCAAATAGATGGCGttt                                         >  1:77745/1‑71 (MQ=255)
                             catcatGCTATTGCCAACACAGAGTAAAACGTCAGTCACGTAGTCTCCTCACCATCAAATAGATGGCGttt                                         <  1:1879755/71‑1 (MQ=255)
                              atcatGCTATTGCCAACACAGAGTAAAACGTCAGTCACGTAGTCTCCTCACCATCAAATAGATGGCGtttt                                        <  1:1136048/71‑1 (MQ=255)
                              atcatGCTATTGCCAACACAGAGTAAAACGTCAGTCACGCAGtct                                                                  <  1:1975283/45‑1 (MQ=255)
                              atcatGCTATTGCCAACACAGAGTAAAACGTCAGTCACGCAGtct                                                                  <  1:2231660/45‑1 (MQ=255)
                              atcatGCTATTGCCAACACAGAGTAAAACGTCAGTCACGCAGtct                                                                  <  1:1182351/45‑1 (MQ=255)
                              atcatGCTATTGCCAACACAGAGTAAAACGTCAGTCACGCAGtct                                                                  <  1:1971938/45‑1 (MQ=255)
                               tcatGCTATTGCCAACACAGAGTAAAACGTCAGTCACGTAGTCTCCTCAcc                                                           <  1:2515437/51‑1 (MQ=255)
                               tcatGCTATTGCCAACACAGAGTAAAACGTCAGTCACGTAGTCTCCTCACCATCa                                                       <  1:1677653/55‑1 (MQ=255)
                               tcatGCTATTGCCAACACAGAGTAAAACGTCAGTCACGTAGTCTCCTCACCATCa                                                       <  1:773884/55‑1 (MQ=255)
                               tcatGCTATTGCCAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCAAATAg                                                  <  1:1456910/60‑1 (MQ=255)
                               tcatGCTATTGCCAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCAAATAGAt                                                <  1:1348966/62‑1 (MQ=255)
                               tcatGCTATTGCCAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCAAATAGAt                                                <  1:399207/62‑1 (MQ=255)
                                catGCTATTGCCAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCAAATAGATGGCGTTTTcc                                      >  1:1569809/1‑71 (MQ=255)
                                 atGCTATTGCCAACACAGAGTAAAACGTCAGTCACGTAGTCTCCTCACCATCaaa                                                     >  1:799888/1‑55 (MQ=255)
                                 atGCTATTGCCAACACAGAGTAAAACGTCAGTCACGTAGTCTCCTCACCATCaaa                                                     >  1:1306666/1‑55 (MQ=255)
                                 atGCTATTGCCAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCa                                                       <  1:1323067/53‑1 (MQ=255)
                                 atGCTATTGCCAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCAAATAGATGGCGTTTTCCt                                     <  1:1825333/71‑1 (MQ=255)
                                      aTTGCCAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCAAATAGATGGCGTTTTCCtgatg                                 >  1:131927/1‑70 (MQ=255)
                                       ttGCCAACACAGAGTAAAACGTCAGTCACGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCtgatg                                 <  1:1728267/69‑1 (MQ=255)
                                       ttGCCAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCaa                                                      <  1:150087/48‑1 (MQ=255)
                                       ttGCCAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCaa                                                      <  1:2470685/48‑1 (MQ=255)
                                       ttGCCAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCaa                                                      <  1:789448/48‑1 (MQ=255)
                                       ttGCCAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGaa                               >  1:746427/1‑71 (MQ=255)
                                         gCCAACACAGAGTAAAACGTCAGTCACGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCtga                                   >  1:2150415/1‑65 (MQ=255)
                                         gCCAACACAGAGTAAAACGTCAGTCACGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCtga                                   >  1:493387/1‑65 (MQ=255)
                                         gCCAACACAGAGTAAAACGTCAGTCACGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAat                              >  1:1224539/1‑70 (MQ=255)
                                           cAACACAGAGTAAAACGTCAGTCACGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATc                           <  1:51135/71‑1 (MQ=255)
                                           cAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATc                           <  1:1402991/71‑1 (MQ=255)
                                            aaCACAGAGTAAAACGTCAGTCACGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCa                          >  1:94796/1‑71 (MQ=255)
                                            aaCACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCAAATAGa                                                 >  1:761087/1‑48 (MQ=255)
                                             acacAGAGTAAAACGTCAGTCACGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCAt                         >  1:639098/1‑71 (MQ=255)
                                              cacaGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCAAATAGATGGCGTTTTCCtgatg                                 <  1:1197515/62‑1 (MQ=255)
                                               acaGAGTAAAACGTCAGTCACGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATg                        >  1:1452976/1‑70 (MQ=255)
                                                   agTAAAACGTCAGTCACGCAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATgcagca                   <  1:1426339/71‑1 (MQ=255)
                                                    gTAAAACGTCAGTCACGCAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATgcagca                   <  1:2536687/70‑1 (MQ=255)
                                                     tAAAACGTCAGTCACGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAatat                            <  1:2251000/60‑1 (MQ=255)
                                                     tAAAACGTCAGTCACGCAGTCTCCTCACCATCAAATa                                                   <  1:1780749/37‑1 (MQ=255)
                                                     tAAAACGTCAGTCACGCAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATGCAGCATa                 <  1:2252924/71‑1 (MQ=255)
                                                       aaaCGTCAGTCACGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATGCAGCATACt               >  1:487079/1‑71 (MQ=255)
                                                       aaaCGTCAGTCACGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATGCAGCATACt               >  1:2585710/1‑71 (MQ=255)
                                                       aaaCGTCAGTCACGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATGCAGCATACt               >  1:697040/1‑71 (MQ=255)
                                                       aaaCGTCAGTCACGCAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATGCAGCATACt               >  1:504042/1‑71 (MQ=255)
                                                        aaCGTCAGTCACGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATGCAGCATACtt              <  1:777695/71‑1 (MQ=255)
                                                        aaCGTCAGTCACGCAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATGCAGCATACtt              <  1:1006111/71‑1 (MQ=255)
                                                        aaCGTCAGTCACGCAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATGCAGCATACtt              <  1:2056601/71‑1 (MQ=255)
                                                             cagtcaCGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATGCAGCATACTTAACAg         >  1:2047526/1‑71 (MQ=255)
                                                             cagtcaCGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATGCAGCATACTTAACAg         >  1:1771832/1‑71 (MQ=255)
                                                             cagtcaCGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATGCAGCATACTTAACAg         >  1:2050789/1‑71 (MQ=255)
                                                             cagtcaCGCAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATGCAGCATACTTAACa          >  1:166216/1‑70 (MQ=255)
                                                                 caCGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGaa                               >  1:1832536/1‑45 (MQ=255)
                                                                 caCGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGaa                               >  1:2576933/1‑45 (MQ=255)
                                                                 caCGCAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATGCAGCATACTTAACAGCAtt     >  1:2081361/1‑71 (MQ=255)
                                                                 caCGCAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATGCAGCATACTTAACAGCAtt     >  1:201117/1‑71 (MQ=255)
                                                                  aCGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAata                             >  1:2504804/1‑46 (MQ=255)
                                                                  aCGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAata                             >  1:1288972/1‑46 (MQ=255)
                                                                  aCGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATGCAGCATACTTAACAGCATTg    <  1:2047387/71‑1 (MQ=255)
                                                                  aCGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATGCAGCATACTTAACAGCATTg    >  1:201444/1‑71 (MQ=255)
                                                                  aCGTAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATGCAGCATACTTAACAGCATTg    <  1:1285180/71‑1 (MQ=255)
                                                                  aCGCAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATGCAGCATACTTAACAGCa       >  1:2513689/1‑68 (MQ=255)
                                                                  aCGCAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATGCAGCATACTTAACAGCATTg    <  1:767812/71‑1 (MQ=255)
                                                                    gCAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGaa                               >  1:2186068/1‑42 (MQ=255)
                                                                    gCAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATGCAGCATACTTAACAGCATTGTa  >  1:972602/1‑71 (MQ=255)
                                                                     |                                                                     
GCCAGCAGCGGACCTGCGCCATCATCGCCCATCATGCTATTGCCAACACAGAGTAAAACGTCAGTCACGCAGTCTCCTCACCATCAAATAGATGGCGTTTTCCTGATGAATATCATGCAGCATACTTAACAGCATTGTA  >  W3110S.gb/2841632‑2841770

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: