Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 987,461 G→T 17.3% L327L (CTG→CTT ycaQ → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606987,4610G→T17.3% 5.7 ‑4.2 ‑0.158981ycaQhypothetical protein
Reads supporting (aligned to +/- strand):  new base (10/1):  ref base (22/25):  total (32/26)
Fisher's exact test strand distribution p-value = 1.55e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.17e-01

CCAGCGCCGAAACGCCAGTATGGCTATTTTGTTCTGCCGTTATTACATCGTGGGCAATTAGTTGGGCGAAT  >  REL606/987426‑987496
                                   |                                   
ccAGCGCCGAAACGCCAGTATGGCTATTTTATTCTg                                     <  1:7549987/36‑1 (MQ=255)
 cAGCGCCGAAACGCCAGTATGGCTATTTTTTTCTGc                                    >  1:64011/1‑36 (MQ=255)
 cAGCGCCGAAACGCCAGTATGGCTATTTTGTTCTGc                                    <  1:2769859/36‑1 (MQ=255)
  aGCGCCGAAACGCCAGTATGGCTATTTTTTTCTGcc                                   >  1:1834732/1‑36 (MQ=255)
    cggcGAAACGCCAGTATGGCTTTTTTGTTCTGCCGt                                 <  1:5267994/33‑1 (MQ=37)
    cgcCGAAACGCCAGTATGGCTATTTTTTTCTTCCGt                                 >  1:5583968/1‑36 (MQ=37)
    cgcCGAAACGCCAGTATGGCTAATTTTTTTTGccct                                 >  1:4936454/1‑34 (MQ=21)
      ccGAAACGCCAGTATGGCTATTTTTTTCTGCCGtta                               >  1:903866/1‑36 (MQ=255)
      ccGAAACGCCAGTAT.GCTATTTTGTTCTGCCGttat                              <  1:2772710/36‑1 (MQ=38)
       cGAAACGCCAGTATGGCTATTTTTTTCTTCCGtttt                              >  1:2488127/1‑34 (MQ=25)
        gAAACGCCAGTATGGCTATTTTTTTCTTCCGttttt                             >  1:8086608/1‑36 (MQ=21)
        gAAACGCCAGTATGGCTATTTTTTTCTGCCGttttt                             >  1:3057625/1‑36 (MQ=37)
        gAAACGCCAGTATGGCTATTTTGTTCTGCCGttatt                             <  1:928657/36‑1 (MQ=255)
         aaaCGCCAGTATGGCTATTTTTTTCTTGCGtttttc                            >  1:3031538/1‑35 (MQ=16)
         aaaCGCCAGTATGGCTATTTTTTTCTTCCGttatta                            >  1:4673876/1‑36 (MQ=37)
         aaaCGCCAGTATGGCTATTTTGTTCTTCCGttatta                            >  1:5263940/1‑36 (MQ=255)
         aaaCGCCAGTATGGCTATTTTGTTCTGCCGttatta                            <  1:3474047/36‑1 (MQ=255)
          aaCGCCAGTATGGCTATTTTGTTCTGCCGTTATTAc                           <  1:3654181/36‑1 (MQ=255)
           aCGCCAGTATGGCTATTTTTTTTTTCCTTTTTTACa                          >  1:6718696/1‑36 (MQ=12)
           aCGCAAGTATGGCTATTTTAATCTTCCGtttttata                          >  1:6807062/1‑34 (MQ=9)
            cGCCAGTATGGCTATTTTTTTATTCCGTTATTACAt                         >  1:5419227/1‑36 (MQ=21)
             gCCAGTATGGCTATTTTTTTCTTCCGTTATTACATc                        >  1:775668/1‑36 (MQ=37)
             gCCAGTATGGCTATTTTTTTCTGCCGTTATTATATc                        >  1:1177199/1‑36 (MQ=37)
             gCCAGTATGGCTATTTTTTTCTGCCGTTATTACATc                        >  1:3481749/1‑36 (MQ=255)
             gCCAGTATGGCTATTTTGTTCTGCCGTTATTACAtt                        >  1:7838163/1‑35 (MQ=255)
             gCCAGTATGGCTATTTTGTTCTGCCGTTATTACATc                        <  1:5226290/36‑1 (MQ=255)
                 gTATGGCTCTTTTGTTCTGCAGTTATTACATCGTgg                    <  1:1397952/36‑1 (MQ=37)
                 gTATGGCTATTTTGTTCTGCCGTTATTACATCGTgg                    >  1:3560720/1‑36 (MQ=255)
                 gTATGGCTATTTTGTTCTGCCGTTATTACATCGTgg                    <  1:1807494/36‑1 (MQ=255)
                   aTGGCTATTTTGTTCTGCCGTTATTACATCGTGGGc                  <  1:3147749/36‑1 (MQ=255)
                    tGGCTATTTTGTTCTGCCGTTATTACATCGTGGGCa                 <  1:3685754/36‑1 (MQ=255)
                    tGGCTATTTTGTTCTGCCGTTATTACATCGTGGGCa                 >  1:1207733/1‑36 (MQ=255)
                     ggCTATTTTGTTCTGCCGTTATTACATCGTGGGCaa                >  1:1196293/1‑36 (MQ=255)
                     ggCTATTTTGTTCTGCCGTTATTACATCGTGGGCaa                <  1:1801596/36‑1 (MQ=255)
                     ggCTATTTTGTTCTGCCGTTATTACATCGTGGGCaa                <  1:4814806/36‑1 (MQ=255)
                      gCTATTTTGTTCTGCCGTTATTACATCGTGGTCAAt               >  1:2108788/1‑36 (MQ=255)
                      gCTATTTTGTTCTGCCGTTATTACATCGTGGGCAAt               <  1:1284027/36‑1 (MQ=255)
                      gCTATTTTGTTCTGCCGTTATTACATCGTGGGCAAt               >  1:8174012/1‑36 (MQ=255)
                      gCTATTTTGTTCTGCCGTTATTACATCGTGGGCAAt               >  1:4567367/1‑36 (MQ=255)
                      gCTATTTTGTTCTGCCGTTATTACATAGTGGGCAAt               >  1:6742779/1‑36 (MQ=255)
                       cTATTTTGTTCTGCCGTTATTACATCGTGGGCAAtt              <  1:742098/36‑1 (MQ=255)
                        tATTTTGTTCTGCCGTTATTACATCGTGGGCAATTa             >  1:3217149/1‑36 (MQ=255)
                          ttttGTTCTGCCGTTATTAGATCGTGGGCAATTAGt           <  1:428013/36‑1 (MQ=255)
                          ttttGTTCTGCCGTTATTACATCGTGGGCAATTAGt           <  1:6052331/36‑1 (MQ=255)
                          ttttGTTCTGCCGTTATTACATCGTGGGCAATTAGt           >  1:5546315/1‑36 (MQ=255)
                          ttttGTTCTGCCGTTATTACATCGTGGGCAATTAGt           <  1:4016970/36‑1 (MQ=255)
                          ttttGTTCTGCCGTTATTACATCGTGGGCAATTAGt           <  1:3313878/36‑1 (MQ=255)
                          ttttGTTCTGCCGTTATTACATCGTGGGCAATTAGt           >  1:1026380/1‑36 (MQ=255)
                          ttttATACTTCCTTTATTGCATCTTGGGCAATTAGt           <  1:2876440/36‑1 (MQ=12)
                            ttGTTCTGCTGTTATTACATCGTGGGCAATTAGTTg         <  1:5335981/36‑1 (MQ=255)
                            ttGTTCTGCCGTTATTACATCGTGGGCAATTAGTTg         >  1:5301541/1‑36 (MQ=255)
                            ttGTTCTGCCGTTATTACATCGTGGGCAATTAGTTg         <  1:4906625/36‑1 (MQ=255)
                             tGTTCTGCCGTTATTACATCGTGGGCAATTAGTTgg        <  1:6779980/36‑1 (MQ=255)
                              gTTCTGCCGTTATTACATCGTGGGCAATTAGTTggg       <  1:3827046/36‑1 (MQ=255)
                              gTTCTGCCGTTATTACATCGTGGGCAATTAGTTggg       >  1:3362657/1‑36 (MQ=255)
                               ttCTGCCGTTATTACATCGTGGGCAATTAGTTGGGc      <  1:1097882/36‑1 (MQ=255)
                                 cTGCCGTTATTACATCGTGGGCAATTAGTTGGGCGa    >  1:1124733/1‑36 (MQ=255)
                                  tGCCGTTATTACATCGTGGGCAATTAGTTGGGCGaa   <  1:5505563/36‑1 (MQ=255)
                                  tGCCGTTATTACATCGTGGGCAATTAGTTGGGCGaa   >  1:2192225/1‑36 (MQ=255)
                                   gCCGTTATTACATCGTGGGCAATTAGTTGGGCGAAt  >  1:2392532/1‑36 (MQ=255)
                                   |                                   
CCAGCGCCGAAACGCCAGTATGGCTATTTTGTTCTGCCGTTATTACATCGTGGGCAATTAGTTGGGCGAAT  >  REL606/987426‑987496

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG