Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,004,251 T→C 90.2% intergenic (‑177/+425) 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,2510T→C90.2% 2.4 ‑1.5 0.031intergenic (‑177/+425)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  new base (13/14):  ref base (3/0):  total (17/14)
Fisher's exact test strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAA  >  REL606/1004218‑1004283
                                 |                                
tcccGTGAAAGTGATACAAAGATAAAGATTTGGctt                                <  1:2173314/34‑1 (MQ=37)
tACCGTGAAAGTGCTACAAAGATAAAGATTTGGctt                                >  1:4512586/1‑36 (MQ=255)
 aCCGTGAAAGTGCTACAAAGATAAAGATTTGGCTTc                               >  1:1401959/1‑36 (MQ=255)
   cGTGAAAGTGCTACAAAGATAAAGGTTTGGCTTCaa                             <  1:5054806/36‑1 (MQ=37)
    gTGAAAGTGCTACAAAGATAAAGATTTGGCTTCaaa                            >  1:856375/1‑36 (MQ=255)
      gAAAGTGCTACAAAGATAAAGATTTGGCTTCaaaaa                          <  1:1536091/36‑1 (MQ=255)
      gAAAGTGCTACAAAGATAAAGATTTGGCTTCaaaaa                          <  1:2984700/36‑1 (MQ=255)
       aaaGTGCTACAAAGATAAAGATTTGGCTTCAAAAAg                         <  1:3007571/36‑1 (MQ=255)
        aaGTGCTACAAAGATAAAGATTTGGCTTCAAAAAGa                        <  1:8129467/36‑1 (MQ=255)
         aGTGCTACAAAGATAAAGATTTGGCTTCAAAAAGaa                       >  1:3889545/1‑36 (MQ=255)
           tGCTACAAAGATAAAGATTTGGCTTCAAAAAGaaaa                     >  1:542408/1‑36 (MQ=255)
           tGCTACAAAGATAAAGATTTGGCTTCAAAAAGaaaa                     <  1:1117082/36‑1 (MQ=255)
              tACAAAGATAAAGATTTGGCTTCAAAAAGAAAAata                  >  1:8446033/1‑36 (MQ=255)
                cAAAGATAAAGATTTGGCTTCAAAAAGAAAAATATg                <  1:1632446/36‑1 (MQ=255)
                  aagataaagatTTGGTTTCAAAAAGAAAAATATGTa              >  1:5867990/1‑36 (MQ=255)
                  aagataaagatTTGGCTTCAAAAAGAAAAATATGTa              >  1:1424064/1‑36 (MQ=255)
                  aagataaagatTTGGCTTCAAAAAGAAAAATATGTa              <  1:8365085/36‑1 (MQ=255)
                  aagataaagatTTGGCTTCAAAAAGAAAAATATGTa              >  1:7816136/1‑36 (MQ=255)
                  aagataaagatTTGGCTTCAAAAAGAAAAATATGTa              >  1:6837389/1‑36 (MQ=255)
                     ataaagatTTGGCTTCAAAAAGAAAAATATGTAAcc           >  1:6720271/1‑36 (MQ=255)
                      taaagatTTGGCTTCAAAAAGAAAAATATGTAACCa          >  1:3711722/1‑36 (MQ=255)
                      taaagatTTGGCTTCAAAAAGAAAAATATGTAACCa          <  1:1682546/36‑1 (MQ=255)
                        aggatTTGGTTTCAAAAAGAAAAATATGTAACCaaa        >  1:4409851/3‑36 (MQ=255)
                        aagatTTGGCTTCAAAAAGAAAAATATGTAACCaaa        <  1:958272/36‑1 (MQ=255)
                        aagatTTGGCTTCAAAAAGAAAAATATGTAACCaaa        <  1:1120907/36‑1 (MQ=255)
                         agatTTGGTTTCAAAAAGAAAAATATGTAACCaaaa       >  1:667111/1‑36 (MQ=255)
                         agatTTGGCTTCAAAAAGAAAAATATGTAACCaaaa       <  1:7946103/36‑1 (MQ=255)
                         agatTTGGCTTCAAAAAGAAAAATATGTAACCaaaa       >  1:6124938/1‑36 (MQ=255)
                          gatTTGGCTTCAAAAAGAAAAATATGTAACCAAAAg      >  1:3051301/1‑36 (MQ=255)
                           atTTGGCTTCAAAAAGAAAAATATGTAACCAAAAGt     <  1:6122165/36‑1 (MQ=255)
                           atTTGGATTCAAAAAGAAAAATATGTAACCAAAAGt     >  1:4832558/1‑36 (MQ=255)
                             ttGGCTTAAAAAAGAAAAATATGTAACCAAAAGTaa   <  1:5996535/36‑1 (MQ=37)
                              tGGCTTCAAAAAGAAAAATATGTAACCAAAAGTaaa  >  1:7723256/1‑36 (MQ=255)
                                 |                                
TACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAA  >  REL606/1004218‑1004283

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 37 ≤  ATCG