Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,821,628 G→A 57.7% G266G (GGC→GGT 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,6280GA57.7% 6.3 / 29.0 26G266G (GGC→GGTrecADNA strand exchange and recombination protein with protease and nuclease activity
Reads supporting (aligned to +/- strand):  ref base G (5/6);  new base A (10/5);  total (15/11)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

CTCGATCAGCTTCTCTTTTACGCCCAGGTCAACCAGTTCGCCGTAGAAGTTGATACCTTCGCCGTAGAGGATCTGGAATTCAGCCTGTTTAAACGGCGCAGCGATTTTGTTCTTCACCACTT  >  W3110S.gb/2821568‑2821689
                                                            |                                                             
ctcGATCAGCTTCTCTTTTACGCCCAGGTCAACCAGTTCGCCGTAGAAGTTCATACCTTCACCGTAGAgg                                                      <  1:3343030/70‑1 (MQ=255)
   gATCAGCTTCTCTTTTACGCCCAGGTCAACCAGTTCGCCGTAGAAGTTGATACCTTCACCGTAGAGGa                                                     <  1:2228329/68‑1 (MQ=255)
     tCAGCTTCTCTTTTACGCCCAGGTCAACCAGTTAGCCGTAGAAGTTGATACCTTCACCGTAGAGGATCTgg                                                <  1:3446991/71‑1 (MQ=255)
       aGCTTCTCTTTTACGCCCAGGTCAACCAGTTCGCCGTAGAAGTTGATACCTTCGCCGTAGAGGATCTGGa                                               <  1:1918843/70‑1 (MQ=255)
           tctcTTTTACGCCCAGGTCAACCAGTTCGCCGTAGAAGTTGATACCTTCACCGTAGAGGATCTGGAAtt                                            >  1:1131498/1‑69 (MQ=255)
             tctTTTACGCCCAGGTCAACCAGTTCGCCGTAGAAGTTGATACCTTCACCGTAGAGGATCTGGAATTCAGc                                        >  1:1024684/1‑71 (MQ=255)
             tctTTTACGCCCAGGTCAACCAGTTCGCCGTAGAAGTTGATACCTTCACCGTAGAGGATCTGGAATTCAGc                                        >  1:1899659/1‑71 (MQ=255)
                    cGCCCAGGTCAACCAGTTCGCCGTAGAAGTTGATACCTTCGCCGTAGAGGATCTGGAATTCAGCCTGTTTa                                 <  1:2855043/71‑1 (MQ=255)
                             cAACCAGTTCGCCGTAGAAGTTGATACCTTCACCGTaga                                                        >  1:896399/1‑39 (MQ=255)
                             cAACCAGTTCGCCGTAGAAGTTGATACCTTCACCGTaga                                                        >  1:3498708/1‑39 (MQ=255)
                             cAACCAGTTCGCCGTAGAAGTTGATACCTTCACCGTaga                                                        >  1:1567502/1‑39 (MQ=255)
                                      cGCCGTAGAAGTTGATACCTTCGCCGTAGAGGATCTGGAATTCAGCCTGTTTAAACGGCGCAGCGATTTTg               >  1:2604132/1‑71 (MQ=255)
                                      cGCCGTAGAAGTTGATACCTTCGCCGTAGAGGATCTGGAATTCAGCCTGTTTAAACGGCGCAGCGATTTTg               >  1:2918582/1‑71 (MQ=255)
                                      cGCCGTAGAAGTTGATACCTTCGCCGTAGAGGATCTGGAATTCAGCCTGTTTAAACGGCGCAGCGATTTTg               >  1:418153/1‑71 (MQ=255)
                                      cGCCGTAGAAGTTGATACCTTCACCGTAGAGGATCTGGAATTCAGCCTGTTTAAACGGCGCAGCGATTTTg               >  1:2191425/1‑71 (MQ=255)
                                      cGCCGTAGAAGTTGATACCTTCACCGTAGAGGATCTGGAATTCAGCCTGTTTAAACGGCGCAGCGATTTTg               >  1:681737/1‑71 (MQ=255)
                                      cGCCGTAGAAGTTGATACCTTCACCGTAGAGGATCTGGAATTCAGCCTGTTTAAACGGCGCAGCGATTTTg               >  1:1285534/1‑71 (MQ=255)
                                        ccGTAGAAGTTGATACCTTCGCCGTAGAGGATCTGGAATTCAGCCTGTTTAAACGGCGCAGCGATTTTGtt             <  1:2303746/71‑1 (MQ=255)
                                        ccGTAGAAGTTGATACCTTCGCCGTAGAGGATCTGGAATTCAGCCTGTTTAAACGGCGCAGCGATTTTGtt             <  1:1749550/71‑1 (MQ=255)
                                        ccGTAGAAGTTGATACCTTCGCCGTAGAGGATCTGGAATTCAGCCTGTTTAAACGGCGCAGCGATTTTGtt             <  1:1470179/71‑1 (MQ=255)
                                         cGTAGAAGTTGATACCTTCACCGTAGAGGATCTGTAATTCAGCCTGTTTAAACGGCGCAGCGATTTTGttc            <  1:2967948/71‑1 (MQ=255)
                                         cGTAGAAGTTGATACCTTCACCGTAGAGGATCTGGAATTCAGCCTGTTTAAACGGCGCAGCGATTTTGttc            <  1:3042327/71‑1 (MQ=255)
                                                 ttGATACCTTCGCCGTAGAGGATCTGGAATTCAGCCTGTTTAAACGGCGCAGCGATTTTGTTCTTcaccac    <  1:182857/71‑1 (MQ=255)
                                                   gATACCTTCGCCGTAGAGGATCTGGAATTCAGCCTGTTTAAACGGCGCAGCGATTTTGTTCTTCACCACtt  >  1:2711647/1‑71 (MQ=255)
                                                   gATACCTTCGCCGTAGAGGATCTGGAATTCAGCCTGTTTAAACGGCGCAGCGATTTTGTTCTTCACCACtt  >  1:1671731/1‑71 (MQ=255)
                                                   gATACCTTCACCGTAGAGGATCTGGAATTCAGCCTGTTTAAACGGCGCAGCGAttt                 >  1:440237/1‑56 (MQ=255)
                                                            |                                                             
CTCGATCAGCTTCTCTTTTACGCCCAGGTCAACCAGTTCGCCGTAGAAGTTGATACCTTCGCCGTAGAGGATCTGGAATTCAGCCTGTTTAAACGGCGCAGCGATTTTGTTCTTCACCACTT  >  W3110S.gb/2821568‑2821689

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: