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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,004,2430AG94.4% 130.4 / ‑0.9 36intergenic (‑169/+433)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (14/20);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATA  >  REL606/1004170‑1004312
                                                                         |                                                                     
tgtCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataagga                                                                     <  2:355400/76‑1 (MQ=255)
tgtCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataagga                                                                     <  1:398813/76‑1 (MQ=255)
   cTTTATGTGTCTGNAGGCATCTTTCCATTAAAACTAACGTTTCGCGACCGTGAAAGTGCTACAAAGATAAGGAttt                                                                  <  2:1064251/76‑1 (MQ=255)
   cTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGAttt                                                                  >  1:523642/1‑76 (MQ=255)
    tttATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTg                                                                 <  1:182843/76‑1 (MQ=255)
       aTGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGAGCTACAAAGATAAGGTTTTGGtt                                                              >  2:836484/1‑76 (MQ=255)
        tgtgTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGttt                                                             <  2:587906/76‑1 (MQ=255)
         gtgtCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTc                                                            >  1:785713/1‑76 (MQ=255)
         gtgtCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTc                                                            <  1:274457/76‑1 (MQ=255)
         gtgtCTGCAGGCATCTTTCCATTAAAACTAACATTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTc                                                            <  1:274465/76‑1 (MQ=255)
          tgtCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCa                                                           >  1:627720/1‑76 (MQ=255)
           gtCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCaa                                                          <  1:444803/76‑1 (MQ=255)
           gtCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCaa                                                          >  1:934644/1‑76 (MQ=255)
               gCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAg                                                      <  1:846035/76‑1 (MQ=255)
               gCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAg                                                      <  1:626109/76‑1 (MQ=255)
                 aGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGaa                                                    >  2:119405/1‑76 (MQ=255)
                     gtCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGTTTTGGTTTCAAAAAGAAAAat                                                >  1:314286/2‑76 (MQ=255)
                     aTCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAat                                                >  1:314275/1‑76 (MQ=255)
                      tCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAata                                               <  2:441762/76‑1 (MQ=255)
                             aTTAAATCTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGTATTTGGTTTCAAAAAGAAAAATATGTAAcc                                        <  1:1064146/76‑1 (MQ=255)
                             aTTAAATCTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGTATTTGGTTTCAAAAAGAAAAATATGTAAcc                                        <  1:1064147/76‑1 (MQ=255)
                             aTTAAAGCTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAAcc                                        <  2:1134691/76‑1 (MQ=255)
                             aTTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAATATGTAAcc                                        >  2:261407/1‑76 (MQ=255)
                              ttAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCa                                       <  2:803064/76‑1 (MQ=255)
                               tAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCaa                                      <  1:413452/76‑1 (MQ=255)
                                aaaaCTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCaaa                                     >  2:328970/1‑76 (MQ=255)
                                     tAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGAATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTaa                                >  2:431605/1‑76 (MQ=255)
                                      aaCGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTaaa                               <  1:318476/76‑1 (MQ=255)
                                        cGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAAt                             >  1:504801/1‑76 (MQ=255)
                                               cTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGa                      <  1:167903/76‑1 (MQ=255)
                                                 aCCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAAc                    <  1:201233/76‑1 (MQ=255)
                                                  ccGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTGACCAAAAGTANAATTTAAGGagcg                   >  1:149845/1‑73 (MQ=255)
                                                  ccGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACt                   <  2:475342/76‑1 (MQ=255)
                                                   cGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACtt                  <  1:1120433/76‑1 (MQ=255)
                                                           tGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAaca          <  1:426092/76‑1 (MQ=255)
                                                           tGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAaca          >  2:751467/1‑76 (MQ=255)
                                                           tGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAaca          >  2:751480/1‑76 (MQ=255)
                                                              tACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCg       >  1:1135245/1‑76 (MQ=255)
                                                                cAAAGATAAGGAATTGGTTTCAAAAAGAAAAANATCTAACCNAAAGTCAAATTTAAGGAACNTTGTGAACACCGTc     >  2:968305/1‑76 (MQ=255)
                                                                 aaagataaggatTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCa    >  1:486731/1‑76 (MQ=255)
                                                                 aaagataaggatTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCa    >  1:486720/1‑76 (MQ=255)
                                                                   agataaggatTTGGTTTCAACAGGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCata  >  1:5521/1‑76 (MQ=255)
                                                                         |                                                                     
TGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATA  >  REL606/1004170‑1004312

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: