Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,750,741 C→A 19.1% L58F (TTG→TTT radC ← DNA repair protein RadC

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,750,7410C→A19.1% 9.1 ‑6.5 0.049174radCDNA repair protein RadC
Reads supporting (aligned to +/- strand):  new base (2/9):  ref base (28/10):  total (30/19)
Fisher's exact test strand distribution p-value = 1.48e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTC  >  REL606/3750706‑3750776
                                   |                                   
tGAACGCCACTAAATTGTTCATATTCAGAGGTTAAc                                     <  1:1705576/36‑1 (MQ=255)
tGAACGACACTAAATTGTTCATATTCAGAGGTTAAc                                     >  1:7490980/1‑36 (MQ=255)
 gAACGCCACTAAATTGTTCATATTCAGAGGTTAACa                                    <  1:2321311/36‑1 (MQ=255)
  aaCGCCACTAAATTGTTCATATTCAGAGGTTAACaa                                   <  1:5693732/36‑1 (MQ=255)
  aaCGCCACTAAATTGTTCATATTCAGAGGTTAACaa                                   <  1:3662173/36‑1 (MQ=255)
  aaCGCCACTAAATTGTTAATATTCAGAGGTTAACaa                                   <  1:3240189/36‑1 (MQ=255)
    cGCCACTAAATTGTTCATATTCAGAGGTTAACAAGc                                 <  1:7388676/36‑1 (MQ=255)
     gCCACTAAATTGTTCATATTCAGAGGTTAACAAGcc                                >  1:6548722/1‑36 (MQ=255)
     gCCACTAAATTGTTCATATTCAGAGGTTAACAAGcc                                >  1:6309431/1‑36 (MQ=255)
      ccACTAAATTGTTCATATTCAGAGGTTAACAAGCCa                               <  1:4805461/36‑1 (MQ=255)
      ccACTAAATTGTTCATATTCAGAGGTTAACAAGCCa                               <  1:4639496/36‑1 (MQ=255)
       cACTAAATTGTTCATATTCAGAGGTTAACAAGCCAt                              <  1:1772232/36‑1 (MQ=255)
       cACTAAATTGTTCATATTCAGAGGTTAACAAGCCAt                              >  1:599086/1‑36 (MQ=255)
        aCTAAATTGTTAATATTAAGAGGTTAAAAAGCCATa                             <  1:46752/36‑1 (MQ=37)
         cTAAATTGTTCATATTCAGAGGTTAACAAGCCATaa                            >  1:7050846/1‑36 (MQ=255)
         cTAAATTGGTCATATTCAGAGGTTAACAAGCCATaa                            >  1:5955674/1‑36 (MQ=255)
          tAAATTGTTCATATTCAGAGGTTAACAAGCCATaaa                           >  1:6020228/1‑36 (MQ=255)
          tAAATTGTTCATATTCAGAGGTTAACAAGCCATaaa                           >  1:3587595/1‑36 (MQ=255)
          tAAATTGGTCATATTCAGAGGTTAAAAAGCCATaaa                           >  1:3113382/1‑36 (MQ=37)
            aaTTGTTCATATTCAGAGGTTAAAAAGCCATAAaga                         >  1:11791/1‑36 (MQ=255)
            aaTTGTTCACATTCAGAGGTTAACAAGCCATAAaga                         >  1:5444162/1‑36 (MQ=255)
             atggTTCATATTAAGAGGTAAAAAAGCCATAAagag                        <  1:2044701/33‑1 (MQ=16)
             aTTGTTCATATTCAGAGGTTAACAAGCCATAAagag                        >  1:7219205/1‑36 (MQ=255)
              ttGTTCATATTCAGAGGTTAACAAGCCATAAagaga                       >  1:6409834/1‑36 (MQ=255)
              ttGTTCATATTCAGAGGTTAACAAGCCATAAagaga                       >  1:2534808/1‑36 (MQ=255)
              ttGTTCATATTCAGAGGTTAACAAGCCATAAagaga                       >  1:5446606/1‑36 (MQ=255)
               tGTTCATATTCAGAGGTTAACAAGCCATAAagagag                      >  1:2961429/1‑36 (MQ=255)
                gTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGc                     >  1:1927305/1‑36 (MQ=255)
                gTTCATATTCAGAGGTTAAAAAGCCATAAAGAGAGc                     <  1:5420879/36‑1 (MQ=255)
                gTTCATATTAAGAGGTAAAAAAGCCATAAAGAGAGc                     <  1:4718833/36‑1 (MQ=21)
                gTTAATATTAAGAGGTTAACAAGCCATAAAGAGAGc                     <  1:2296683/36‑1 (MQ=38)
                gTTAATATTAAGAGGTTAAAAAGCCATAAAGAGAGc                     <  1:221170/36‑1 (MQ=37)
                 ttCATATTCAGAGGTTAACAAGCCATAAAGAGAGcc                    >  1:3302345/1‑36 (MQ=255)
                 ttCATATTCAGAGGTTAACAAGCCATAAAGAGAGcc                    >  1:3999726/1‑36 (MQ=255)
                 ttCATATTAAGAGGTTAACAAGCCATAAAGAGAGcc                    <  1:7325129/36‑1 (MQ=255)
                  taatatTAAGAGGTAAAAAAGCCATAAAGAGAGCCg                   <  1:1090706/34‑1 (MQ=18)
                   cATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGa                  >  1:1954536/1‑36 (MQ=255)
                    atatTCAGAGGTTAACAAGCCATAAAGAGAGCCGaa                 >  1:7158134/1‑36 (MQ=255)
                    atatTCAGAGGTTAACAAGCCATAAAGAGAGCCGaa                 >  1:208264/1‑36 (MQ=255)
                     tatTCAGAGGTTAACAAGCCATAAAGAGAGCCGaaa                >  1:3754094/1‑36 (MQ=255)
                      atTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAAt               >  1:1191839/1‑36 (MQ=255)
                      atTAAGAGGTTAAAAAGCCATAAAGAGAGCCGAAAt               <  1:5507034/36‑1 (MQ=38)
                        taagagGTAAAAAAGCCATAAAGAGAGCCGAAATtc             <  1:2288952/34‑1 (MQ=25)
                          agagGTTAACAAGCCATAAAGAGAGCCGAAATtctc           >  1:6063831/1‑36 (MQ=255)
                          agagGTTAACAAGCCATAAAGAGAGCCGAAATtctc           >  1:900898/1‑36 (MQ=255)
                           gagGTTAACAAGCCATAAAGAGAGCCGAAATTCTcc          >  1:6304657/1‑36 (MQ=255)
                            agGTTAAAAAGCCATAAAGAGAGCCGAAATTCTCCa         <  1:6575180/36‑1 (MQ=255)
                            agGGTAACAAGCCATAAAGAGAGCCGAAATTCTCCa         >  1:1680918/1‑36 (MQ=255)
                              gTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGc       >  1:5191400/1‑36 (MQ=255)
                                tAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCAt     >  1:6406197/1‑36 (MQ=255)
                                 aaCAAGCCATAAAGAGAGCCGAAATTCTCCAGCAtc    >  1:3285391/1‑36 (MQ=255)
                                   cAAGCCATAAAGAGAGCCGCAATTCTCCAGCAtctc  >  1:2687894/1‑36 (MQ=255)
                                   |                                   
TGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTC  >  REL606/3750706‑3750776

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 22 ≤  ATCG  < 31 ≤  ATCG  < 37 ≤  ATCG