Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,995,272 T→C 48.5% S524S (TCA→TCG uvrC ← excinuclease UvrABC, endonuclease subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,995,2720TC48.5% ‑6.6 / 45.8 33S524S (TCA→TCGuvrCexcinuclease UvrABC, endonuclease subunit
Reads supporting (aligned to +/- strand):  ref base T (8/9);  new base C (12/4);  total (20/13)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-01

CGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTAC  >  W3110S.gb/1995214‑1995339
                                                          |                                                                   
cGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTggcca                                                          >  1:3210038/1‑68 (MQ=255)
 gCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCa                                                        <  1:1386767/71‑1 (MQ=255)
 gCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCa                                                        <  1:1446452/71‑1 (MQ=255)
        gcgTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTaaa                                                 >  1:3020480/1‑71 (MQ=255)
        gcgTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTaaa                                                 >  1:581381/1‑71 (MQ=255)
        gcgTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTaaa                                                 >  1:3146995/1‑71 (MQ=255)
        gcgTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTaaa                                                 >  1:3153002/1‑71 (MQ=255)
        gcgTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTaaa                                                 >  1:2438752/1‑71 (MQ=255)
                aTGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTaaa                                                 <  1:3329047/63‑1 (MQ=255)
                aTGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTaaa                                                 <  1:3073252/63‑1 (MQ=255)
                aTGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTaaa                                                 <  1:2652469/63‑1 (MQ=255)
                aTGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCAc                                         >  1:262751/1‑71 (MQ=255)
                aTGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCAc                                         >  1:2513114/1‑71 (MQ=255)
                aTGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCAc                                         >  1:2408129/1‑71 (MQ=255)
                aTGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCAc                                         >  1:1823942/1‑71 (MQ=255)
                aTGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCAc                                         >  1:1345871/1‑71 (MQ=255)
                           cgcgAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTaaa                                                 <  1:3215307/52‑1 (MQ=255)
                            gcgAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCa                              >  1:2955823/1‑70 (MQ=255)
                            gcgAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCa                              >  1:577281/1‑70 (MQ=255)
                                  atGCTGGATAACATGCAGCGCGGGCGAATCTTGCGGCAAACTAAATc                                               >  1:2032114/1‑47 (MQ=39)
                                  atGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCt                                            >  1:2740619/1‑50 (MQ=255)
                                  atGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCt                                            >  1:2635264/1‑50 (MQ=255)
                                  atGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCt                                            >  1:2181252/1‑50 (MQ=255)
                                  atGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCt                                            >  1:625152/1‑50 (MQ=255)
                                  atGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCt                                            >  1:677590/1‑50 (MQ=255)
                                  atGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCt                                            >  1:902490/1‑50 (MQ=255)
                                        gATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGt                   <  1:3318286/69‑1 (MQ=255)
                                                       gggTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTAc  <  1:2898680/71‑1 (MQ=255)
                                                       gggTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTAc  <  1:245391/71‑1 (MQ=255)
                                                       gggTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTAc  <  1:3330424/71‑1 (MQ=255)
                                                       gggTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTAc  <  1:403999/71‑1 (MQ=255)
                                                       gggTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCCCAAAGAACAGCGTTTCCAGTCCAGCCTTAc  <  1:298444/71‑1 (MQ=255)
                                                       gggTGAATCTGGCGGCAAACTAAATCCATCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTAc  <  1:1321108/71‑1 (MQ=255)
                                                          |                                                                   
CGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTAC  >  W3110S.gb/1995214‑1995339

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: