Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,871,909 G→T 19.4% intergenic (‑96/‑36) yoaA ← / → yoaB conserved protein with nucleoside triphosphate hydrolase domain/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,871,9090G→T19.4% 2.4 ‑4.3 ‑0.039intergenic (‑96/‑36)yoaA/yoaBconserved protein with nucleoside triphosphate hydrolase domain/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (8/0):  ref base (15/15):  total (24/15)
Fisher's exact test strand distribution p-value = 1.32e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.99e-01

AGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGGATAATGAGGACAAAATG  >  REL606/1871874‑1871944
                                   |                                   
aGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTg                                     <  1:1180498/36‑1 (MQ=255)
aGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTg                                     <  1:7924563/36‑1 (MQ=255)
aGTGGTGACGGCGACATCACATTGTGGCAGTCTTTg                                     <  1:6174837/36‑1 (MQ=255)
 gtggtgACGGCGACCTCACATTGTGGCAGTCTTTgc                                    <  1:979924/36‑1 (MQ=255)
 gtggtgACGGCGACCTCACATTGTGGCAGTCTTTgc                                    <  1:4981037/36‑1 (MQ=255)
   ggtgACGGCGACCTCACATTGTGGCAGTCTTTgcag                                  <  1:4602939/36‑1 (MQ=255)
    gtgACGGCGACCTCACATTGTGGCAGTCTTTgcagc                                 <  1:7445728/36‑1 (MQ=255)
      gACGGCGACCTCACATTGTGGCAGTCTTTCCAGCAc                               >  1:4691478/1‑36 (MQ=255)
       aCGGCGACCTCACATTGTGGCAGTCTTTTCAGCACg                              >  1:147675/1‑36 (MQ=255)
         ggCGACCTCACATTGTGGCAGTCTTTGCAGCACGaa                            <  1:3831747/36‑1 (MQ=255)
          gCGACCTCACATTGTGGCAGTCTTTGCAGCACGaaa                           <  1:7855294/36‑1 (MQ=255)
          gCGACCTCACATTGTGGCAGTCTTTGCAGCACGaaa                           <  1:2425017/36‑1 (MQ=255)
          gCGACCTCACATTGTGGCAGTCTTTGCAGCACGaaa                           <  1:6984621/36‑1 (MQ=255)
           cGACCTCACATTGTGGCAGTCTTTTCAGCACGaaaa                          >  1:1176700/1‑36 (MQ=255)
            gACCTCACATTGTGGCAGTCTTTTCAGCACGGAAAt                         >  1:3943070/1‑36 (MQ=37)
             aCCTCACATTGTGGCAGTTTTTTTAGCACGTAAAtt                        >  1:2301155/1‑35 (MQ=16)
             aCCTCACATTGTGGCAGTCTTTTCAGCACGgaaaag                        >  1:4033496/1‑34 (MQ=25)
             aCCTCACATTGTGGCAGTCTTTGCAGCACGAAAATg                        <  1:1802423/36‑1 (MQ=255)
             aCCTCACATTGTGGCAGTCTTTGCAGCACGAAAATg                        <  1:281521/36‑1 (MQ=255)
                tCACATTGTGGCAGTCTTTGCAGCACGAAAATggga                     >  1:4836538/1‑34 (MQ=255)
                  acaTTGTGGCAGTCTTTTCAGCACGAAAATGGAAgg                   >  1:3236453/1‑36 (MQ=255)
                  acaTTGTGGCAGTCTTTGCAGCACGAAAATGGGAgg                   >  1:4701876/1‑36 (MQ=255)
                  acaTTGTGGCAGTCTTTGCAGCACGAAAATGGAAgg                   >  1:2996348/1‑36 (MQ=255)
                   caTTGTGGCAGTCTTTTCAGCACGAAAATGGAAGGa                  >  1:564638/1‑36 (MQ=255)
                   caTTGTGGCAGTCTTTGCAGCACGAAAATGGAAGGa                  >  1:5924360/1‑36 (MQ=255)
                   caTTGTGGCAGTCTTTGCAGCACGAAAATGGAAGGa                  <  1:5248644/36‑1 (MQ=255)
                    aTTGTGGCAGTCTTTGCAGCACGAAAATGGAAGGAt                 <  1:2303658/36‑1 (MQ=255)
                     ttGTGGCAGTCTTTTCAGCACGAAAATGGAAGGGTa                >  1:2182936/1‑36 (MQ=38)
                     ttGTGGCAGTCTTTGCAGCACGAAAATTGAAGGATa                >  1:7767875/1‑36 (MQ=255)
                     ttGTGGCAGTCTTTGCAGCACGAAAATGGAAGGATa                <  1:4844734/36‑1 (MQ=255)
                      tgtgGCAGTCTTTGCAGCACGAAAATGGAAGGATaa               >  1:6986856/1‑36 (MQ=255)
                       gtgGCAGTCTTTGCAGCACGAAAATGGAAGGATAAt              >  1:1141976/1‑36 (MQ=255)
                        tgGCAGTCTTTGCAGCACGAAAATGGAAGGATAATg             >  1:6495321/1‑36 (MQ=255)
                         ggCAGTCTTTGCAGCACGAAAATGGAAGGATAATGa            <  1:4660119/36‑1 (MQ=255)
                         ggCAGTCTTTGCAGCACGAAAATGGAAGGATAATGa            <  1:1639396/36‑1 (MQ=255)
                          gCAGTCTTTGCAGCACGAAAATGGAAGGATAATGAg           >  1:5718465/1‑36 (MQ=255)
                           cAGTCTTTGCAGCACGAAAATGGAAGGATAATGAgg          >  1:4929708/1‑36 (MQ=255)
                            aGTCTTTGCAGCACGAAAATGGAAGGATAATGAGGa         <  1:3689395/36‑1 (MQ=255)
                             gTCTTTGCAGCACGGAAATGGAAGGATAATGAGGAc        >  1:5056330/1‑36 (MQ=255)
                              tCTTTGCAGCACGAAAATGGAAGGATAATGAGGACa       <  1:5250913/36‑1 (MQ=255)
                               cTTTGCAGCACGAAAATTGAAGGATAATGAGGACaa      <  1:1180942/36‑1 (MQ=255)
                               cTTTGCAGCACGAAAATGGAAGGATAATGAGGACaa      >  1:5682871/1‑36 (MQ=255)
                                tttGCAGCACGAAAATGGAAGGATAATGAGGACaaa     <  1:4996375/36‑1 (MQ=255)
                                tttGCAGCACGAAAATGGAAGGATAATGAGGACaaa     >  1:5009241/1‑36 (MQ=255)
                                 ttGCAGCACGAAAATGGAAGGATAATGAGGACaaaa    >  1:7615853/1‑36 (MQ=255)
                                 ttGCAGCACGAAAATGGAAGGATAATGAGGACaaaa    >  1:1862871/1‑36 (MQ=255)
                                   gcagcaCGAAAATGGAAGGATAATGAGGACAAAatg  >  1:6055434/1‑36 (MQ=255)
                                   gcagcaCGAAAATGGAAGGATAATGAGGACAAAatg  >  1:7590158/1‑36 (MQ=255)
                                   |                                   
AGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGGATAATGAGGACAAAATG  >  REL606/1871874‑1871944

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG