Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,871,909 G→T 35.3% 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→T35.3% 6.6 ‑4.4 0.033intergenic (‑96/‑36)yoaA/yoaBconserved protein with nucleoside triphosphate hydrolase domain/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (12/0):  ref base (11/9):  total (24/9)
Fisher's exact test strand distribution p-value = 1.17e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGGATAATGAGGACAAAATG  >  REL606/1871874‑1871944
                                   |                                   
aGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTg                                     <  1:3460803/36‑1 (MQ=255)
aGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTg                                     <  1:2747688/36‑1 (MQ=255)
 gtggtgACGGCGCCCTCACATTGTGGCAGTCTTTgc                                    <  1:5219570/36‑1 (MQ=255)
 gtggtgACGGCGACCTCACATTGTGGCAGTCTTTgc                                    <  1:4701228/36‑1 (MQ=255)
 gtggtgACGGCGACCTCACATTGTGGCAGTCTTTgc                                    <  1:1817059/36‑1 (MQ=255)
    gtgACGGCGACCTCACATTGTGGCAGTCTTTgcagc                                 <  1:2123327/36‑1 (MQ=255)
     tgACGGCGACCTCACATTGTGGCAGTCTTTgcagca                                <  1:3906393/36‑1 (MQ=255)
      gACGGCGAACTCACATTGTGGCAGTCTTTTCAGCAc                               >  1:657202/1‑36 (MQ=38)
        cGGCGACCTCACATTGTGGCAGTTTTTTTAGCTCGa                             >  1:3209786/1‑36 (MQ=19)
        cGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGa                             >  1:5891727/1‑36 (MQ=255)
          gCGACCTCACATTGTGGCAGTCTTTGCAGCACGaaa                           <  1:4016076/36‑1 (MQ=255)
          gCGACCTCACATTGTGGCAGTCTTTGCAGCACGaaa                           <  1:4908332/36‑1 (MQ=255)
           cGACCTCACATTGTGGCAGTCCTTGCAGAACGaaaa                          >  1:3009301/1‑36 (MQ=38)
           cGACCTCACATTGTGGCAGCATTTGCAGCACGaaaa                          >  1:265234/1‑36 (MQ=37)
            gACCTCACATTGTGGCAGTCTTTTCAGTACGAAAAt                         >  1:5157786/1‑36 (MQ=37)
            gACCTCACATTGTGGCAGTCTTTTCAGCACGAAAAt                         >  1:114055/1‑36 (MQ=255)
            gACCTCACATTGGGGAAGTCTGTGCTGCTCGAAAAt                         >  1:1617956/1‑36 (MQ=16)
             aCCTCACATTGTGGCAGTCTGTCCAGCACGAAAATg                        >  1:282535/1‑36 (MQ=37)
              ccTCACATTGTGGCAGTCTTTTTAGCACGAAAATgg                       >  1:5360379/1‑36 (MQ=37)
              ccTCACATTGTGGCAGTCTTTTCCGCACGCAAAttg                       >  1:4267300/1‑34 (MQ=21)
               cTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGa                      <  1:2672343/36‑1 (MQ=255)
               cTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGa                      <  1:2357692/36‑1 (MQ=255)
               cTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGa                      <  1:177182/36‑1 (MQ=255)
               cTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGa                      >  1:4194743/1‑36 (MQ=255)
                tCACATTGTGGCAGTCTTTTCAGCACGAAAATGGaa                     >  1:5633970/1‑36 (MQ=255)
                tCACATTGTGGCAGTCTTTGCAGCACGAAAATGGaa                     <  1:5611141/36‑1 (MQ=255)
                 cacaTTGTGGCAGTCTTTTTAGCACGAAAATGGAAg                    >  1:4066889/1‑36 (MQ=37)
                  acaTTGTGGCAGTCTTTTCAGCACGAAAATGGAAgg                   >  1:3312426/1‑36 (MQ=255)
                  acaTTGTGGCAGCNNTTTCAGAACGAAAATCGAAgg                   >  1:563740/1‑36 (MQ=16)
                    aTTGTGGCAGTCTTTTCAGCACGAAAATGGAAGGAt                 >  1:1630774/1‑36 (MQ=255)
                    aTTGTGGCAGGCTTTTACGCGCGAAAATTGAAgggt                 >  1:1393255/1‑34 (MQ=2)
                       gtgGCAGTCTTTGCAGCACGAAAATGGAAGGATAAt              <  1:3367975/36‑1 (MQ=255)
                         ggCAGTCTTTGCAGCACGAAAATGGAAGGATAATGa            <  1:2409059/36‑1 (MQ=255)
                         ggCAGTCTTTGCAGCACGAAAATGGAAGGATAATGa            >  1:1378902/1‑36 (MQ=255)
                            aGTCTTTGCAGCACGAAAATGGAAGGATAATGAGGa         >  1:844281/1‑36 (MQ=255)
                             gTCTTTGCAGCACGAAAATGGAAGGATAATGAGGAc        >  1:5587918/1‑36 (MQ=255)
                              tCTTTGCAGCACGAAAATGTAAGGATAATGAGGACa       <  1:4570046/36‑1 (MQ=255)
                              tCTTTGCAGCACGAAAATGGAAGGATAATGAGGACa       >  1:2412709/1‑36 (MQ=255)
                               ctgtGCAGCACGACCATGGAAGGGCAATGAGGACaa      >  1:155665/4‑36 (MQ=16)
                               cTTTGCAGCACGAAAATGGAAGGATAATGAGGACaa      >  1:2214504/1‑36 (MQ=255)
                                   gcagcaCGAAAATGGAAGGATAATGAGGACAAAatg  >  1:126929/1‑36 (MQ=255)
                                   |                                   
AGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGGATAATGAGGACAAAATG  >  REL606/1871874‑1871944

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 14 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG