Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,139,938 C→T 41.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,9380CT41.7% 11.0 / 25.6 24R2H (CGT→CAT) hybDpredicted maturation element for hydrogenase 2
Reads supporting (aligned to +/- strand):  ref base C (11/3);  new base T (3/7);  total (14/10)
Fisher's exact test for biased strand distribution p-value = 3.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

CGACAATCCGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTACCGCACAGGCCAT  >  W3110S.gb/3139873‑3140007
                                                                 |                                                                     
cGACAATCCGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATATGCAt                                                                   >  1:409842/1‑70 (MQ=255)
 gACAATCCGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATATGCAtt                                                                  <  1:445872/70‑1 (MQ=255)
        cGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAAcc                                                          >  1:223020/1‑71 (MQ=255)
        cGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAAcc                                                          <  1:458139/71‑1 (MQ=255)
                 aTGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATATGCATTACAGAACCTTCACTGaa                                                 <  1:1080774/71‑1 (MQ=255)
                 aTGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGATCCTTCACTGa                                                  >  1:1699688/1‑70 (MQ=255)
                       tCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACt                                           >  1:1383915/1‑71 (MQ=255)
                          tCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGaa                                                            >  1:944266/1‑51 (MQ=255)
                                   aaaaTATTGCCGACCCCTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGt                                 >  1:1318251/1‑69 (MQ=255)
                                   aaaaTATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGt                                 >  1:277233/1‑69 (MQ=255)
                                    aaaTATTGCCGACCCCTAAGACTAAAATATGCATTGCAGAACCTTCACTGAAACCACTTCGTTGCCGTCAg                              >  1:1365637/1‑71 (MQ=255)
                                    aaaTATTGCCGACCCCTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAg                              <  1:1133730/71‑1 (MQ=255)
                                    aaaTATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGtt                                      >  1:4097/1‑63 (MQ=255)
                                    aaaTATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGt                                       >  1:709441/1‑62 (MQ=255)
                                       tatTGCCGACCCCTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCAt                           <  1:988482/71‑1 (MQ=255)
                                         ttGCCGACCCCTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATcc                         <  1:399313/71‑1 (MQ=255)
                                             cGACCCCTAAGACTAAAATATGCATTACAGAACCtt                                                        <  1:314851/36‑1 (MQ=255)
                                             cGACCCCTAAGACTAAAATATGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACt                      <  1:927723/70‑1 (MQ=255)
                                                   cTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTa               >  1:1658795/1‑71 (MQ=255)
                                                            aaaTACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACgtgt                >  1:521001/1‑61 (MQ=255)
                                                            aaaTACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACgtgt                >  1:952258/1‑61 (MQ=255)
                                                             actACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTACCGc           <  1:344586/63‑1 (MQ=255)
                                                               tACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTACCGCACAGGCCa   <  1:24073/71‑1 (MQ=255)
                                                                aCGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTACCGCACAGGCCAt  >  1:1368434/1‑71 (MQ=255)
                                                                 |                                                                     
CGACAATCCGCACACCGATGGCTTCATCGGTCAGCAAAATATTGCCGACCCCTAAGACTAAAATACGCATTACAGAACCTTCACTGAAACCACTTCGTTGCCGTCAGCATCCACTACGTGTACCGCACAGGCCAT  >  W3110S.gb/3139873‑3140007

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: