Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,139,938 C→T 43.7% 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,9380CT43.7% 41.0 / 113.3 87R2H (CGT→CAT) hybDpredicted maturation element for hydrogenase 2
Reads supporting (aligned to +/- strand):  ref base C (24/25);  new base T (17/21);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 8.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

CTTCGACAATCCGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTACCGCACAGGC  >  W3110S.gb/3139870‑3140004
                                                                    |                                                                  
cTTCGACAATCCGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACg                                                                   >  1:1700250/1‑70 (MQ=255)
 ttCGACAATCCGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATATGCa                                                                 >  1:246736/1‑71 (MQ=255)
      cAATCCGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTa                                                              >  1:1129046/1‑69 (MQ=255)
       aaTCCGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAg                                                           <  1:3181229/71‑1 (MQ=255)
          ccGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCAt                                                                <  1:2958149/63‑1 (MQ=255)
          ccGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAAc                                                        >  1:4812949/1‑71 (MQ=255)
           cGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATATGCATTACAGaa                                                         >  1:1906411/1‑69 (MQ=255)
           cGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATATGCATTACAGAAcc                                                       >  1:1052741/1‑71 (MQ=255)
             cacaCCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATATGCATTACAGAACCtt                                                     <  1:4536660/71‑1 (MQ=255)
             cacaCCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACa                                                            >  1:830136/1‑64 (MQ=255)
             cacaCCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACa                                                            >  1:4420595/1‑64 (MQ=255)
               cacCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCTTAAGACTAAAATACGCATTACAGAACCTTCa                                                   >  1:2772014/1‑71 (MQ=255)
               cacCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCTTAAGACTAAAATACGCATTACAGAACCTTCa                                                   >  1:2766006/1‑71 (MQ=255)
               cacCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCa                                                   >  1:3004257/1‑71 (MQ=255)
                 ccGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATATGCATTACAGAACCtt                                                     <  1:3783067/67‑1 (MQ=255)
                 ccGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACt                                                 <  1:1758464/71‑1 (MQ=255)
                  cGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACt                                                 >  1:1321390/1‑70 (MQ=255)
                   gATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGa                                               <  1:4248945/71‑1 (MQ=255)
                      ggCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATATGCATTACAGAACCTTCACTGAAAc                                            >  1:1809202/1‑71 (MQ=255)
                       gCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATATGCATTACAGAAc                                                        >  1:3659606/1‑58 (MQ=255)
                       gCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATATGCATTACAGAAc                                                        >  1:5083360/1‑58 (MQ=255)
                          tCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCtt                                                     <  1:1111959/58‑1 (MQ=255)
                          tCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACt                                        >  1:2377691/1‑71 (MQ=255)
                           cATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCAc                                         <  1:3643531/69‑1 (MQ=255)
                           cATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCAc                                         <  1:200492/69‑1 (MQ=255)
                           cATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCAc                                         <  1:862453/69‑1 (MQ=255)
                           cATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCACt                                        <  1:4152441/70‑1 (MQ=255)
                            aTCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTc                                      <  1:2093003/71‑1 (MQ=255)
                            aTCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTc                                      >  1:3840061/1‑71 (MQ=255)
                            aTCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTc                                      >  1:240454/1‑71 (MQ=255)
                            aTCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTc                                      <  1:4747608/71‑1 (MQ=255)
                             tCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTc                                                    <  1:3908461/56‑1 (MQ=255)
                              cGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTc                                      <  1:2498966/69‑1 (MQ=255)
                               ggTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCCtt                                   <  1:1928884/71‑1 (MQ=255)
                                gTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGt                                    >  1:1716431/1‑69 (MQ=255)
                                gTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTg                                  <  1:4839980/71‑1 (MQ=255)
                                gTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTg                                  <  1:1762156/71‑1 (MQ=255)
                                 tCAGCAAAATATTGCCGACCCCTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTc                                      <  1:1238806/66‑1 (MQ=255)
                                 tCAGCAAAATATTGCCGACCCCTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTCGTTGc                                 >  1:4392688/1‑71 (MQ=255)
                                     cAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTc                                                    <  1:1055248/48‑1 (MQ=255)
                                     cAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTc                             <  1:4160151/71‑1 (MQ=255)
                                       aaaTATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAg                           >  1:28548/1‑71 (MQ=255)
                                        aaTATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGaa                                              <  1:2708941/51‑1 (MQ=255)
                                        aaTATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGc                          <  1:2737448/71‑1 (MQ=255)
                                            ttGCCGACCCCTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATc                       <  1:5136027/70‑1 (MQ=255)
                                            ttGCCGACCCCTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATc                       <  1:658631/70‑1 (MQ=255)
                                            ttGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATc                       <  1:1308103/70‑1 (MQ=255)
                                             tGCCGACCCCTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCa                     >  1:1484968/1‑71 (MQ=255)
                                               ccGACCCCTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCAc                    <  1:3971536/70‑1 (MQ=255)
                                                    cccTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACgtg              >  1:4211936/1‑71 (MQ=255)
                                                    cccTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACgtg              <  1:95582/71‑1 (MQ=255)
                                                    cccTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACgtg              <  1:3082259/71‑1 (MQ=255)
                                                    cccTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACgtg              <  1:3088000/71‑1 (MQ=255)
                                                     ccTAAGACTAAAATATGCATTAGAGAACCTTCACTGAAACCACtt                                       <  1:1831543/45‑1 (MQ=255)
                                                     ccTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCACtt                                       <  1:3461341/45‑1 (MQ=255)
                                                     ccTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCACtt                                       <  1:3457451/45‑1 (MQ=255)
                                                     ccTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACgtg              >  1:1107687/1‑70 (MQ=255)
                                                      cTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCAc                    <  1:4528869/63‑1 (MQ=255)
                                                      cTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCAc                    <  1:2045754/63‑1 (MQ=255)
                                                      cTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCAc                    <  1:2250716/63‑1 (MQ=255)
                                                      cTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACgtgt             >  1:4955275/1‑70 (MQ=255)
                                                      cTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTa            >  1:2907734/1‑71 (MQ=255)
                                                      cTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTACGTTGCCGTCAGCATCCACTACGTGTa            >  1:3716396/1‑71 (MQ=255)
                                                       tAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTc                                      <  1:1304639/44‑1 (MQ=255)
                                                       tAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTAc           >  1:350485/1‑71 (MQ=255)
                                                       tAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTAc           >  1:1853715/1‑71 (MQ=255)
                                                       tAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTAc           >  1:1285258/1‑71 (MQ=255)
                                                        aaGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTCg                                     >  1:2775337/1‑44 (MQ=255)
                                                        aaGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACgtgt             <  1:4267195/68‑1 (MQ=255)
                                                        aaGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACgtgt             <  1:1038247/68‑1 (MQ=255)
                                                        aaGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACgtgt             <  1:3240562/68‑1 (MQ=255)
                                                         aGACTAAAATATGCATTACAGAACCTTCACTGaaa                                             >  1:4171258/1‑35 (MQ=255)
                                                         aGACTAAAATATGCATTACAGAACCTTCACTGaaa                                             >  1:4769618/1‑35 (MQ=255)
                                                         aGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTa            >  1:4555215/1‑68 (MQ=255)
                                                         aGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTACCg         <  1:4271242/71‑1 (MQ=255)
                                                         aGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTACCg         <  1:3272681/71‑1 (MQ=255)
                                                          gACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTACCg         >  1:3749529/1‑70 (MQ=255)
                                                              aaaaTATGCATTACAGAACCTTCACTGAAACCACTTCGtt                                   <  1:2659287/40‑1 (MQ=255)
                                                               aaaTATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTa                  >  1:2186374/1‑56 (MQ=255)
                                                               aaaTATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTACCGCACAgg   >  1:4682143/1‑71 (MQ=255)
                                                               aaaTATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTACCGCACAgg   >  1:1150492/1‑71 (MQ=255)
                                                               aaaTATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCAACTa                  >  1:1861897/1‑56 (MQ=255)
                                                               aaaTACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCa                            >  1:3141086/1‑46 (MQ=255)
                                                               aaaTACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTACCGCACAgg   >  1:1304772/1‑71 (MQ=255)
                                                                aaTACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTACCGCACAGGc  <  1:2381343/71‑1 (MQ=255)
                                                                aaTACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTACCGCACAGGc  <  1:468296/71‑1 (MQ=255)
                                                                aaTACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTACCGCACAGGc  <  1:3338660/71‑1 (MQ=255)
                                                                    |                                                                  
CTTCGACAATCCGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTACCGCACAGGC  >  W3110S.gb/3139870‑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: