Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,821,480 T→C 76.0% 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,4800TC76.0% 94.5 / 28.4 50T316A (ACC→GCC) recADNA strand exchange and recombination protein with protease and nuclease activity
Reads supporting (aligned to +/- strand):  ref base T (8/4);  new base C (12/26);  total (20/30)
Fisher's exact test for biased strand distribution p-value = 4.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCACCTTTGTAGC  >  W3110S.gb/2821412‑2821548
                                                                    |                                                                    
gaAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGctt                                                                    >  1:2985916/1‑71 (MQ=255)
gaAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGctt                                                                    >  1:1724953/1‑71 (MQ=255)
             gAGTTCGGGTTGCTCAGCAGCAACTCACGTACTATCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCttt                                                       >  1:3306567/1‑71 (MQ=255)
                  cGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGc                                                   >  1:3525766/1‑70 (MQ=255)
                  cGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGc                                                   >  1:158038/1‑70 (MQ=255)
                  cGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGc                                                   >  1:2713232/1‑70 (MQ=255)
                       tGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAgg                                               >  1:1448143/1‑69 (MQ=255)
                        gCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGc                                                   >  1:254351/1‑64 (MQ=255)
                        gCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGc                                                   >  1:1593623/1‑64 (MQ=255)
                         cTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAg                                            >  1:3593279/1‑70 (MQ=255)
                         cTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGt                                           <  1:3465806/71‑1 (MQ=255)
                         cTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGt                                           >  1:1292612/1‑71 (MQ=255)
                            agcagcAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCg                                         <  1:1832089/70‑1 (MQ=255)
                            agcagcAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCg                                         <  1:771541/70‑1 (MQ=255)
                            agcagcAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCg                                         <  1:3504196/70‑1 (MQ=255)
                            agcagcAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCg                                         <  1:931165/70‑1 (MQ=255)
                            agcagcAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGc                                        <  1:3212327/71‑1 (MQ=255)
                             gcagcaACTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCa                                       <  1:1193140/71‑1 (MQ=255)
                             gcagcaACTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCa                                       <  1:944131/71‑1 (MQ=255)
                                gcaACTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGcc                                                  <  1:1888979/57‑1 (MQ=255)
                                gcaACTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGcc                                                  <  1:89212/57‑1 (MQ=255)
                                 caACTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCg                                   <  1:1127398/71‑1 (MQ=255)
                                  aaCTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGc                                  >  1:2484975/1‑71 (MQ=255)
                                  aaCTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGc                                  >  1:2604016/1‑71 (MQ=255)
                                  aaCTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGc                                  >  1:3003191/1‑71 (MQ=255)
                                    cTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCtt                                <  1:522637/71‑1 (MQ=255)
                                    cTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCtt                                <  1:2100480/71‑1 (MQ=255)
                                    cTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCtt                                <  1:583791/71‑1 (MQ=255)
                                    cTCACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCtt                                <  1:345526/71‑1 (MQ=255)
                                      cACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTa                              >  1:1330370/1‑71 (MQ=255)
                                      cACGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCAAGGCAGTCGCATTCGCTTTa                              >  1:2791118/1‑71 (MQ=255)
                                       aCGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTAc                             <  1:1271655/71‑1 (MQ=255)
                                       aCGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCttt                               <  1:237819/69‑1 (MQ=255)
                                       aCGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCttt                               <  1:184994/69‑1 (MQ=255)
                                       aCGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCttt                               <  1:2806008/69‑1 (MQ=255)
                                       aCGTACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTa                              <  1:3118371/70‑1 (MQ=255)
                                       aCGTACTTTCTTCTCGATCTCTTTCGCGGCTCCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCttt                               <  1:1886630/69‑1 (MQ=255)
                                          tACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCt                          >  1:3332367/1‑71 (MQ=255)
                                          tACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCt                          >  1:1474133/1‑71 (MQ=255)
                                          tACTTTCTTCTCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACACt                          >  1:2307654/1‑71 (MQ=255)
                                                cttctCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCg                     <  1:3097790/70‑1 (MQ=255)
                                                cttctCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCg                     <  1:64015/70‑1 (MQ=255)
                                                cttctCGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCg                     <  1:3101926/70‑1 (MQ=255)
                                                   ctcGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGa                    <  1:530590/68‑1 (MQ=255)
                                                   ctcGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGa                    <  1:465625/68‑1 (MQ=255)
                                                   ctcGATCTCTTTCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGa                    <  1:710813/68‑1 (MQ=255)
                                                       aTCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCg                     <  1:99054/63‑1 (MQ=255)
                                                            tttCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCAc           <  1:3465085/68‑1 (MQ=255)
                                                            tttCGCGGCTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCAc           <  1:3362752/68‑1 (MQ=255)
                                                                  ggTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCa            >  1:525276/1‑61 (MQ=255)
                                                                  ggTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCACCTTTGTAGc  >  1:722429/1‑71 (MQ=255)
                                                                  ggTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCACCTTTGTAGc  >  1:2268433/1‑71 (MQ=255)
                                                                    |                                                                    
GAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGTTATCTTTCAGCCAGGCAGTCGCATTCGCTTTACCCTGACCGATCTTCTCACCTTTGTAGC  >  W3110S.gb/2821412‑2821548

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: