Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,135,984 A→G 20.7% S26P (TCC→CCC)  dinI ← DNA damage‑inducible protein I

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,135,9840A→G20.7% 5.7 ‑3.6 ‑0.03576dinIDNA damage‑inducible protein I
Reads supporting (aligned to +/- strand):  new base (7/0):  ref base (13/14):  total (21/14)
Fisher's exact test strand distribution p-value = 2.62e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.85e-01

TTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCT  >  REL606/1135949‑1136018
                                   |                                  
ttCATTATCAGGAAACGCATACTGAATACGGCGGGa                                    >  1:1175388/1‑36 (MQ=255)
 tCATTATCAGGAAACGCATACTGAATACGGCGGGaa                                   >  1:1401210/1‑36 (MQ=255)
  cATTATCAGGAAACGCATACTGAATACGGCGGGaaa                                  <  1:1717564/36‑1 (MQ=255)
   aTTATCAGGAAACGCATACTGAATACGGCGGGgaaa                                 >  1:6227710/1‑35 (MQ=255)
    ttATCAGGAAACGCATACTGAATACGGCGGGAAAGt                                <  1:739229/36‑1 (MQ=255)
     tATCAGGAAACGCATACTGAATACGGCGGGGAAAtt                               >  1:5504267/1‑36 (MQ=37)
      aTCAGGAAACGCATACTGAATACGGCGGGGAAGttt                              >  1:746054/1‑35 (MQ=255)
      aTCAGGAAACGCATACTGAATACGGCGGGGAAGttt                              >  1:818439/1‑35 (MQ=255)
      aTCAGGAAACGCATAATGAATACGGCGGGTAAGtta                              >  1:5839991/1‑35 (MQ=37)
        cAGGAAACGCATACTGAATACGGCGGGAAAGTTCgc                            <  1:6558285/36‑1 (MQ=255)
        cAGGAAACGCATACTGAATACGGCGGGAAAGTTCgc                            <  1:6390264/36‑1 (MQ=255)
         aGGAAACGCATACTGAATACGGCGGGGAAGGTCgcc                           >  1:7890142/1‑36 (MQ=38)
         aGGAAACGCATACTGAATACGGCGGGAAAGTTCgcc                           <  1:5322617/36‑1 (MQ=255)
          ggAAACGCATACTGAATACGGCGGGAAAGTTCgcca                          <  1:7857406/36‑1 (MQ=255)
             aaCGCATACTGAATACGGCGGTAAAGTTCgccagcc                       <  1:2403744/36‑1 (MQ=255)
              aCGCATACTGAATACGGCGGGAAAGTTCgccagcca                      <  1:4016135/36‑1 (MQ=255)
                gCATACTGAATACGGCGGGAAAGTTCGCCAGCCAgg                    <  1:339441/36‑1 (MQ=255)
                 cATACTGAATACGGCGGGGAAGATTGCCAACCAagg                   >  1:1109800/1‑36 (MQ=16)
                 cATACTGAATACGGCGGGGAAATTTTTCCGCCAggg                   >  1:7638710/1‑36 (MQ=9)
                 cATACGGAATACGGCGGGAAAGTTCGCCAGCCAggg                   <  1:345294/36‑1 (MQ=255)
                   tACTGAATACGGCGGGGAAGATTGCCAGCCAggggg                 >  1:1561642/1‑34 (MQ=21)
                   tACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCg                 <  1:4683122/36‑1 (MQ=255)
                   tACTGAATACGGCGGGAAAAATTGGCAGGCAggggg                 >  1:3344772/1‑34 (MQ=12)
                    aCTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGt                <  1:7361797/36‑1 (MQ=255)
                      tGAACAAGGCGGTAAAGTTCGCCAGCCAGGGCGTCa              <  1:7612195/36‑1 (MQ=25)
                         aTACGGCGGGAAAGTTCGCCCGCCCGGGCGGCCATa           >  1:2471207/1‑36 (MQ=16)
                           aCGACGGCAAAGTCCGCCAGCCAGGGCGTCAATAGc         <  1:4399044/36‑1 (MQ=25)
                             ggcggGAAAGTTCGCCAGCCAGGGCGTCAATAGccc       >  1:843131/1‑36 (MQ=255)
                             ggcggGAAAGTTCGCCAGCCAGGGCGTCAATAGccc       >  1:1299501/1‑36 (MQ=255)
                             ggcggGAAAGTTCGCCAGCCAGGGCGTCAATAGccc       >  1:7209655/1‑36 (MQ=255)
                             ggcggGAAAGTTCGCCAGCCAGGGCGTCAATAGccc       >  1:3022678/1‑36 (MQ=255)
                             ggcggGAAAGTTCGCCAGCCAGGGCGTCAATAGccc       >  1:5564216/1‑36 (MQ=255)
                              gcggGAAAGTTCGCCAGCCAGGGCGTCACTAGcccc      >  1:2865002/1‑36 (MQ=255)
                               cggGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCa     <  1:800455/36‑1 (MQ=255)
                                gggAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAg    >  1:5221748/1‑36 (MQ=255)
                                 ggAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGc   >  1:1600116/1‑36 (MQ=255)
                                 ggAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGc   >  1:5963191/1‑36 (MQ=255)
                                 ggAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGc   >  1:5304248/1‑36 (MQ=255)
                                  gAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCt  >  1:5631506/1‑36 (MQ=255)
                                   |                                  
TTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCT  >  REL606/1135949‑1136018

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 17 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG