Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,139,938 C→T 6.4% R2H (CGT→CAT)  hybD ← predicted maturation element for hydrogenase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,139,9380CT 6.4% 253.4 / 5.4 78R2H (CGT→CAT) hybDpredicted maturation element for hydrogenase 2
Reads supporting (aligned to +/- strand):  ref base C (39/34);  new base T (3/2);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTTCGACAATCCGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTACCGCACAGGC  >  W3110S.gb/3139868‑3140004
                                                                      |                                                                  
cGCTTCGACAATCCGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATAc                                                                    <  1:2059020/71‑1 (MQ=255)
   ttCGACAATCCGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCa                                                                 >  1:2433586/1‑71 (MQ=255)
     cGACAATCCGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCAt                                                                >  1:1375701/1‑70 (MQ=255)
     cGACAATCCGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCAt                                                                <  1:3211049/70‑1 (MQ=255)
     cGACAATCCGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCAt                                                                >  1:1519188/1‑70 (MQ=255)
       aCAATCCGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTa                                                              >  1:2465303/1‑70 (MQ=255)
       aCAATCCGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTa                                                              >  1:2852940/1‑70 (MQ=255)
        cAATCCGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTa                                                              >  1:643452/1‑69 (MQ=255)
        cAATCCGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACa                                                            <  1:2861315/71‑1 (MQ=255)
            ccGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAAc                                                        >  1:2212506/1‑71 (MQ=255)
            ccGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAAc                                                        >  1:4237391/1‑71 (MQ=255)
             cGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAAcc                                                       <  1:1796118/71‑1 (MQ=255)
              gCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATATGCATTACAGAACCt                                                      <  1:3253079/71‑1 (MQ=255)
              gCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCt                                                      >  1:3599810/1‑71 (MQ=255)
               cacaCCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCtt                                                     >  1:3486984/1‑71 (MQ=255)
                acacCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCtt                                                     <  1:1745165/70‑1 (MQ=255)
                acacCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCtt                                                     >  1:3458468/1‑70 (MQ=255)
                acacCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTc                                                    <  1:380189/71‑1 (MQ=255)
                     gATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACt                                                 >  1:1834535/1‑69 (MQ=255)
                      aTGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATATGCATTACAGAACCTTCACTGa                                               >  1:334305/1‑70 (MQ=255)
                         gCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAAc                                            >  1:346873/1‑70 (MQ=255)
                            tCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAAcc                                           >  1:2473586/1‑68 (MQ=255)
                            tCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACt                                        >  1:893563/1‑71 (MQ=255)
                            tCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACt                                        >  1:1346510/1‑71 (MQ=255)
                            tCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACt                                        >  1:4240217/1‑71 (MQ=255)
                             cATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACtt                                       >  1:234819/1‑71 (MQ=255)
                             cATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACtt                                       >  1:678589/1‑71 (MQ=255)
                             cATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACtt                                       >  1:4065556/1‑71 (MQ=255)
                              aTCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACt                                        <  1:4143233/69‑1 (MQ=255)
                              aTCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTc                                      <  1:3860026/71‑1 (MQ=255)
                               tCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTg                                                <  1:3263704/60‑1 (MQ=255)
                               tCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGa                                               <  1:2222414/61‑1 (MQ=255)
                               tCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTc                                      <  1:3123489/70‑1 (MQ=255)
                               tCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCg                                     <  1:3971269/71‑1 (MQ=255)
                               tCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCg                                     <  1:3050577/71‑1 (MQ=255)
                               tCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCg                                     <  1:382943/71‑1 (MQ=255)
                                cGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGAATTACAGAACCTTCACTGAAACCACTTCGt                                    >  1:2229295/1‑71 (MQ=255)
                                 ggTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGtt                                   <  1:3194515/71‑1 (MQ=255)
                                    cAGCAAAATATTGCCGACGCCTAAGACTTAAATACGCTTTACAGAACCTTCACt                                                 <  1:1948799/54‑1 (MQ=255)
                                      gCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGt                              <  1:2198066/71‑1 (MQ=255)
                                      gCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGt                              >  1:2911837/1‑71 (MQ=255)
                                      gCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGt                              <  1:1656013/71‑1 (MQ=255)
                                      gCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGt                              <  1:2052435/71‑1 (MQ=255)
                                        aaaaTATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTg                                                >  1:4004382/1‑51 (MQ=255)
                                         aaaTATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAg                           <  1:2666577/71‑1 (MQ=255)
                                         aaaTATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACAGaa                                              >  1:1536155/1‑52 (MQ=255)
                                          aaTATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACt                                                 >  1:860468/1‑48 (MQ=255)
                                           atatTGCCGACCCCTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCa                         >  1:1733665/1‑71 (MQ=255)
                                           atatTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTg                                                >  1:2853659/1‑48 (MQ=255)
                                           atatTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCa                         >  1:3394445/1‑71 (MQ=255)
                                            tatTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCAt                        >  1:2571429/1‑71 (MQ=255)
                                              ttGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATc                       <  1:2384215/70‑1 (MQ=255)
                                              ttGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATc                       >  1:710780/1‑70 (MQ=255)
                                              ttGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATc                       >  1:599265/1‑70 (MQ=255)
                                                 ccGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATcc                      <  1:3947834/68‑1 (MQ=255)
                                                 ccGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACt                   <  1:2484737/71‑1 (MQ=255)
                                                    aCCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACg                >  1:1489158/1‑71 (MQ=255)
                                                    aCCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACg                >  1:252603/1‑71 (MQ=255)
                                                     ccccTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACgt               >  1:759366/1‑71 (MQ=255)
                                                     ccccTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACgt               <  1:1556375/71‑1 (MQ=255)
                                                     ccccTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACgt               <  1:3040201/71‑1 (MQ=255)
                                                      cccTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACgtg              <  1:3096759/71‑1 (MQ=255)
                                                       ccTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGt                              <  1:3093450/54‑1 (MQ=255)
                                                        cTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGc                          <  1:2605800/57‑1 (MQ=255)
                                                        cTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCa                            <  1:681808/55‑1 (MQ=255)
                                                        cTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCa                         <  1:2853003/58‑1 (MQ=255)
                                                        cTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCa                     >  1:2214797/1‑62 (MQ=255)
                                                         tAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATcc                      <  1:650315/60‑1 (MQ=255)
                                                          aaGACTAAAATACGCATTACAGAACCTTCACTTaa                                              >  1:2287083/1‑35 (MQ=255)
                                                          aaGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTAcc          >  1:1738895/1‑71 (MQ=255)
                                                          aaGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTAcc          >  1:2215757/1‑71 (MQ=255)
                                                           aGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTACCg         <  1:1942545/71‑1 (MQ=255)
                                                           aGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTACCg         <  1:2318879/71‑1 (MQ=255)
                                                           aGACTAAAATACGCATTACAGAACCATCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTAccc         <  1:3327086/71‑2 (MQ=255)
                                                             aCTAAAATATGCATTACAGAACCTTCACTGAAACCACTTCGtt                                   >  1:2523409/1‑43 (MQ=255)
                                                               tAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCa                            <  1:3624154/48‑1 (MQ=255)
                                                                 aaaTACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTACCGCACAgg   >  1:3859863/1‑71 (MQ=255)
                                                                  aaTACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTa                  >  1:3157648/1‑55 (MQ=255)
                                                                  aaTACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTACCGCACAGGc  <  1:607039/71‑1 (MQ=255)
                                                                      |                                                                  
CGCTTCGACAATCCGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTACCGCACAGGC  >  W3110S.gb/3139868‑3140004

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: