Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,821,480 T→C 59.4% T316A (ACC→GCC)  recA ← DNA strand exchange and recombination protein with protease and nuclease activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,821,4800TC59.4% 21.8 / 31.2 32T316A (ACC→GCC) recADNA strand exchange and recombination protein with protease and nuclease activity
Reads supporting (aligned to +/- strand):  ref base T (5/8);  new base C (8/11);  total (13/19)
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

CGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCACCTT  >  W3110S.gb/2821421‑2821542
                                                           |                                                              
cGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGtt                                                       <  1:1739271/69‑1 (MQ=255)
cGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGtt                                                       <  1:1823905/69‑1 (MQ=255)
cGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGtt                                                       <  1:3007318/69‑1 (MQ=255)
 gTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTa                                                      <  1:429344/69‑1 (MQ=255)
    gAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATc                                                    >  1:2145275/1‑68 (MQ=255)
         cGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGtt                                                       <  1:3408118/60‑1 (MQ=255)
            gTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCa                                               >  1:3615127/1‑65 (MQ=255)
                    gcagcaACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAg                                          >  1:16059/1‑62 (MQ=255)
                    gcagcaACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCg                                   >  1:1966398/1‑69 (MQ=255)
                    gcagcaACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCg                                   >  1:2841512/1‑69 (MQ=255)
                     cagcaACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCTGGTTATCTTTCAGCCAGGCAGTCg                                   <  1:1192784/68‑1 (MQ=255)
                     cagcaACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCg                                   <  1:787499/68‑1 (MQ=255)
                     cagcaACTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCg                                   <  1:200494/68‑1 (MQ=255)
                     cagcaACTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCg                                   <  1:2352030/68‑1 (MQ=255)
                     cagcaACTCACGTACTTTCTTCTCGATCTATTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCg                                   <  1:3215392/68‑1 (MQ=255)
                      agcaACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCa                                 <  1:3212892/69‑1 (MQ=255)
                      agcaACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCa                                 <  1:3607343/69‑1 (MQ=255)
                      agcaACTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCa                                 <  1:2849476/69‑1 (MQ=255)
                        caACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCAtt                               <  1:815731/69‑1 (MQ=255)
                         aaCTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCg                                   >  1:3298157/1‑64 (MQ=255)
                         aaCTCACGTACTTTCTTCTCGATCTCTATCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTc                              >  1:3370437/1‑69 (MQ=255)
                           cTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGc                            <  1:1732669/69‑1 (MQ=255)
                            tCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGc                            >  1:707772/1‑68 (MQ=255)
                                            cGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTAGCCTGACCGa              <  1:921391/66‑1 (MQ=255)
                                               tctcTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGAtcttct        >  1:93808/1‑69 (MQ=255)
                                               tctcTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGAtcttct        >  1:665537/1‑69 (MQ=255)
                                               tctcTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGAtctt          <  1:714814/67‑1 (MQ=255)
                                                     tCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGAtcttct        <  1:2706105/63‑1 (MQ=255)
                                                     tCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTtctc       <  1:2218037/64‑1 (MQ=255)
                                                     tCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCACCtt  >  1:2885503/1‑69 (MQ=255)
                                                     tCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCACCtt  >  1:2881089/1‑69 (MQ=255)
                                                     tCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCACCtt  >  1:2058029/1‑69 (MQ=255)
                                                           |                                                              
CGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCACCTT  >  W3110S.gb/2821421‑2821542

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: