Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,672,722 G→T 14.1% L199L (CTG→CTT ppnK → inorganic polyphosphate/ATP‑NAD kinase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,672,7220G→T14.1% 3.7 ‑3.1 ‑0.137597ppnKinorganic polyphosphate/ATP‑NAD kinase
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (13/18):  total (19/18)
Fisher's exact test strand distribution p-value = 4.55e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.66e-01

CGCCTATTCCCTCTCTGCAGGCGGTCCTATTCTGACCCCCTCTCTGGATGCGATTACCCTGGTGCCCAT  >  REL606/2672689‑2672757
                                 |                                   
cGCCTATTCCCTCTCTGCAGGCGGTCCTATTCTGAc                                   <  1:5577665/36‑1 (MQ=255)
 gCCTATTCCCTCTCTGCAGGCGGTCCTATTCTGAcc                                  <  1:7197637/36‑1 (MQ=255)
    tATTCCCTCTCTGCAGGCGGTCCTATTCTGAccccc                               <  1:5928887/36‑1 (MQ=255)
     aTTCCCTCTCTGCAGGCGGTCCTTTTTTTTCCCCct                              >  1:600867/1‑36 (MQ=19)
      ttCCCTCTCTGCAGGCGGTCCTATTCTGACCCCctc                             <  1:5332202/36‑1 (MQ=255)
      ttCCCTCTCTGCAGGCGGTCCTATTCTGACCCCctc                             <  1:7518066/36‑1 (MQ=255)
       tCCCTCTCTGCAGGCGGTCCTATTCTGACCCCctct                            <  1:5818064/36‑1 (MQ=255)
        cccTCTCTGCAGGCGGTCCTATTCTGACCCCctctc                           <  1:6488136/36‑1 (MQ=255)
        cccTCTCTGCAGGCGGTCCTATTCTGACCCCctctc                           <  1:7744810/36‑1 (MQ=255)
        cccTCTCTGCAGGCGGTCCCTTTCTGACCCCctctc                           >  1:4806088/1‑36 (MQ=38)
         ccgctctGCAGGCGGTCCTATTCTGACCCCctctct                          <  1:2698272/33‑1 (MQ=255)
         ccTCTCTGCAGGCGGTCCTATTTTTACCCCctctct                          >  1:1967204/1‑36 (MQ=37)
           tctctGCAGGCGGTCCTCTTTTTACCCCccctcctg                        >  1:3708807/1‑33 (MQ=255)
           tctctGCAGGCGGTCCTATTCTGACCCCCTCTCTgg                        <  1:5546429/36‑1 (MQ=255)
            ctctGCAGGCGGTCCTTTTTTTTCCCCccctctttg                       >  1:5536205/1‑33 (MQ=9)
            ctctGCAGGCGGTCCTATTCTTACCCCctctcttgg                       >  1:1961590/1‑33 (MQ=38)
             tctGCAGGCGGTCCTATTCTGACCCCCTCTCTGGAt                      <  1:1863794/36‑1 (MQ=255)
               tGCAGGCGGTCCTATTCTGACCCCCTCTCTGGATGc                    <  1:3441851/36‑1 (MQ=255)
               tGCAGGCGGTACTTTTTTCACCCCCTCCCCGGGTGc                    >  1:788750/1‑36 (MQ=0)
                  aGGCGGTCCTATTCTGACCCCCTCTCTGGATGCGAt                 <  1:4158800/36‑1 (MQ=255)
                   ggCGGTCCTATTCTGACCCCCTCTCTGGATGCGAtt                >  1:1474508/1‑36 (MQ=255)
                    gCGGTCCTATTCTGACCCCCTCTCTGGATGCGATTa               >  1:614025/1‑36 (MQ=255)
                    gCGGTCCTATTCTGACCCCCTCTCTGGATGCGATTa               <  1:409335/36‑1 (MQ=255)
                     cGGTCCTATTCTGACCCCCTCTCTGGATGCGATTAc              >  1:5550055/1‑36 (MQ=255)
                       gTCCTATTCTGACCCCCTCTCTGGATGCGATTAccc            <  1:1268770/36‑1 (MQ=255)
                         ccTATTCTGACCCCCTCTCTGGATGCGATTACCCTg          >  1:346974/1‑36 (MQ=255)
                           tATTCTGACCCCCTCTCTGGATGCGATTACCCTGGt        >  1:5112969/1‑36 (MQ=255)
                           tATTCTGACCCCCTCTCTGGATGCGATTACCCTGGt        <  1:6863947/36‑1 (MQ=255)
                            aTTCTGACCCCCTCTCTGGATGCGATTACCCTGGTg       >  1:5242731/1‑36 (MQ=255)
                            aTTCTGACCCCCTCTCTGGATGCGATTACCCTGGTg       >  1:7167853/1‑36 (MQ=255)
                            aTTCTGACCCCCTCTCTGGATGCGATTACCCTGGTg       >  1:7171800/1‑36 (MQ=255)
                              tCTGACCCCCTCTCTGGATGCGATTACCCTGGTGcc     >  1:4979204/1‑36 (MQ=255)
                              tCTGACCCCCTCTCTGGATGCGATTACCCTGGTGcc     <  1:4275395/36‑1 (MQ=255)
                              tCTGACCCCCTCTCTGGATGCGATTACCCTGGTGcc     >  1:2041844/1‑36 (MQ=255)
                              tCTGACCCCCTCTCTGGATGCGATTACCCTGGTGcc     >  1:1647465/1‑36 (MQ=255)
                               cTGACCCCCTCTCTGGATGCGATTACCCTGGTGccc    >  1:4383002/1‑36 (MQ=255)
                                tGACCCCCTCTCTGGATGCGATTACCCTGGTGCCCa   <  1:1659560/36‑1 (MQ=255)
                                 gACCCCCTCTCTGGATGCGATTACCCTGGTGCCCAt  <  1:7637624/36‑1 (MQ=255)
                                 |                                   
CGCCTATTCCCTCTCTGCAGGCGGTCCTATTCTGACCCCCTCTCTGGATGCGATTACCCTGGTGCCCAT  >  REL606/2672689‑2672757

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 19 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG