Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,004,243 A→G 20.0% intergenic (‑169/+433) ompF ← / ← asnS outer membrane porin 1a (Ia;b;F)/asparaginyl‑tRNA synthetase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,004,2430A→G20.0% 8.5 ‑2.8 ‑0.125intergenic (‑169/+433)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (10/10):  total (15/10)
Fisher's exact test strand distribution p-value = 6.13e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.35e-01

CGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAATATGTAACCAAA  >  REL606/1004210‑1004277
                                 |                                  
cGTTTCGCTACCGTGAAAGTGCTACaaagataaaga                                  <  1:486710/36‑1 (MQ=255)
cGTTTCGCTACCGTGAAAGTGCTACaaagataaaga                                  >  1:2681959/1‑36 (MQ=255)
 gTTTCGCTACCGTGAAAGTGCTACaaagataaagat                                 >  1:5012172/1‑36 (MQ=255)
   ttCGCTACCGTGAAAGTGCTACAAAGATAAAGAttt                               >  1:6255556/1‑36 (MQ=255)
    tCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTg                              >  1:3788637/1‑36 (MQ=255)
     cGCTACCGTGAAAGTGCTACAAAGATAAATATTTgg                             <  1:6044798/36‑1 (MQ=255)
     cGCTACCGTGAAAGTGCTACAAAGATAAAGATTTgg                             <  1:4558630/36‑1 (MQ=255)
       cTACCGTGAAAGTGCTACAAAGATAAAGATTTGGtt                           >  1:1730055/1‑36 (MQ=255)
       cTACCGTGAAAGTGCTACAAAGATAAAGATTTGGtt                           >  1:2114169/1‑36 (MQ=255)
        tACCGTGAAAGTGCTACAAAGATAAAGATTTGGttt                          >  1:5108582/1‑36 (MQ=255)
         aCCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTc                         <  1:2157290/36‑1 (MQ=255)
          ccGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCa                        <  1:3375348/36‑1 (MQ=255)
           cGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCaa                       >  1:1409035/1‑36 (MQ=255)
            gTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaa                      <  1:331672/36‑1 (MQ=255)
             tGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaaa                     >  1:4589056/1‑36 (MQ=255)
                aaGTGCTACAAAGATAAGGATTTGGTTTAAAAAAGa                  >  1:1740525/1‑36 (MQ=38)
                aaGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGa                  <  1:3342824/36‑1 (MQ=255)
                aaGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGa                  <  1:408890/36‑1 (MQ=255)
                aaGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGa                  <  1:1551562/36‑1 (MQ=255)
                  gTGCTACAAAGATAAGGATTTGGTTTCAAAAAGaaa                >  1:5259340/1‑36 (MQ=255)
                   tGCTACAAAGATAAAGATTTGGTTTCAAAAAGaaaa               <  1:551300/36‑1 (MQ=255)
                   tGCTACAAAGATAAAGATTTGGTTTCAAAAAGaaaa               >  1:7833442/1‑36 (MQ=255)
                     cTACAAAGATAAAGATTTGGTTTCAAAAAGAACAat             >  1:1578025/1‑36 (MQ=255)
                      tACAAAGATAAAGATTTGGTTTCAAAAAGAAAAata            <  1:5772449/36‑1 (MQ=255)
                       aCAAAGATAAGGATTTGGTTTCAAAAAGAAAAatat           >  1:3984415/1‑36 (MQ=255)
                       aCAAAGATAAAGATTTGGTTTCAAAAAGAAAAatat           >  1:7312425/1‑36 (MQ=255)
                        cAAAGATAAATATTTGGTTTCAAAAAGAAAAATATg          >  1:6986158/1‑36 (MQ=255)
                        cAAAGATAAAGATTTGGTTTCAAAAAGAAAAATATg          >  1:1786214/1‑36 (MQ=255)
                         aaagataaagatTTGGTTTCAAAAAGAAAAATATGt         >  1:4930873/1‑36 (MQ=255)
                         aaagataaagatTTGGTTTCAAAAAGAAAAATATGt         >  1:4322973/1‑36 (MQ=255)
                         aaagataaagatCTGGTTTCAAAAAGAAAAATATGt         >  1:1906397/1‑36 (MQ=255)
                          aagataaatatTTGGTTTCAAAAAGAAAAATATGTa        >  1:4569971/1‑36 (MQ=255)
                          aagataaagatTTGGTTTCAAAAAGAAAAATATGTa        >  1:1352987/1‑36 (MQ=255)
                            gataaatatTTGGTTTCAAAAAGAAAAATATGTAAc      >  1:7113437/1‑36 (MQ=255)
                            gataaagatTTGGTTTCAAAAAGAAAAATATGTAAc      >  1:600001/1‑36 (MQ=255)
                            gataaagatTTGGTTTCAAAAAGAAAAATATGTAAc      <  1:5249686/36‑1 (MQ=255)
                             ataaagatTTGGTTTCAAAAAGAAAAATATGTAAcc     <  1:4345153/36‑1 (MQ=255)
                               aaagatTTGGTTTCAAAAAGAAAAATATGTAACCaa   >  1:4991841/1‑36 (MQ=255)
                                aagatTTGGTTTCAACAAGAAAAATATGTAACCaaa  >  1:5794212/1‑36 (MQ=255)
                                aagatTTGGTTTCAAAAAGAAAAATATGTAACCaaa  >  1:5719878/1‑36 (MQ=255)
                                aagatTTGGTTTCAAAAAGAAAAATATGTAACCaaa  >  1:3328007/1‑36 (MQ=255)
                                 |                                  
CGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAATATGTAACCAAA  >  REL606/1004210‑1004277

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 37 ≤  ATCG