Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 124,866 A→G intergenic (‑511/‑30) aroP ← / → pdhR aromatic amino acid transporter/transcriptional regulator of pyruvate dehydrogenase complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606124,8660AG98.5% 271.3 / ‑5.5 67intergenic (‑511/‑30)aroP/pdhRaromatic amino acid transporter/transcriptional regulator of pyruvate dehydrogenase complex
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (37/29);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AAAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGT  >  REL606/124793‑124936
                                                                         |                                                                      
aaaaCCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCt                                                                      >  1:1446096/1‑76 (MQ=255)
aaaaCCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCt                                                                      >  1:354076/1‑76 (MQ=255)
  aaCCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTc                                                                    >  1:670715/1‑76 (MQ=255)
    ccTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCgg                                                                  <  1:1644377/76‑1 (MQ=255)
        tATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAg                                                              >  1:824119/1‑76 (MQ=255)
         aTGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGt                                                             >  2:1232085/1‑76 (MQ=255)
          tGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTg                                                            >  1:1543615/1‑76 (MQ=255)
            gACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGc                                                          <  1:88867/76‑1 (MQ=255)
              cATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCtt                                                        >  1:1003551/1‑76 (MQ=255)
               aTAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTa                                                       <  2:1084234/76‑1 (MQ=255)
               aTAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTa                                                       >  2:911834/1‑76 (MQ=255)
                 aaGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAg                                                     >  1:215528/1‑76 (MQ=255)
                 aaGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAg                                                     >  1:270039/1‑76 (MQ=255)
                 aaGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAg                                                     >  1:800674/1‑76 (MQ=255)
                 aaGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTAAAg                                                     >  2:1018468/1‑76 (MQ=255)
                   ggTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGAc                                                   <  1:872493/76‑1 (MQ=255)
                     tGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAg                                                 <  1:1800344/76‑1 (MQ=255)
                     tGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAg                                                 <  1:460429/76‑1 (MQ=255)
                       aaTACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGa                                               >  1:884278/1‑76 (MQ=255)
                        aTACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGaa                                              >  2:353744/1‑76 (MQ=255)
                         tACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAAc                                             >  1:291865/1‑76 (MQ=255)
                          aCTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACt                                            <  1:441102/76‑1 (MQ=255)
                          aCTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACt                                            >  2:1377264/1‑76 (MQ=255)
                          aCTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACt                                            <  2:535195/76‑1 (MQ=255)
                          aCTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACt                                            <  2:726497/76‑1 (MQ=255)
                           cTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTc                                           >  1:1520529/1‑76 (MQ=255)
                            tttGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCa                                          >  1:1754836/1‑76 (MQ=255)
                             ttGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCAt                                         >  2:288047/1‑76 (MQ=255)
                             ttGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCAt                                         >  2:313995/1‑76 (MQ=255)
                             ttGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCAt                                         >  2:328544/1‑76 (MQ=255)
                             ttGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCAt                                         >  1:798931/1‑76 (MQ=255)
                             ttGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCAt                                         <  1:901842/76‑1 (MQ=255)
                              tGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATg                                        >  1:889882/1‑76 (MQ=255)
                               gTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATgg                                       <  1:1497042/76‑1 (MQ=255)
                                ttACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGc                                      >  2:248279/1‑76 (MQ=255)
                                ttACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGc                                      >  2:733889/1‑76 (MQ=255)
                                  aCTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCt                                    >  2:941621/1‑76 (MQ=255)
                                   cTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTa                                   >  2:75457/1‑76 (MQ=255)
                                   cTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTa                                   <  1:677760/76‑1 (MQ=255)
                                   cTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTa                                   <  1:549142/76‑1 (MQ=255)
                                    tttAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTAc                                  <  1:933482/76‑1 (MQ=255)
                                     ttAGCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACa                                 <  2:1828588/76‑1 (MQ=255)
                                        gCGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCa                              <  1:1613255/76‑1 (MQ=255)
                                         cGTCACAGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCaa                             <  1:1071653/76‑1 (MQ=255)
                                            cacaGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAAt                          >  2:890392/1‑76 (MQ=255)
                                            cacaGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAAt                          >  1:1801218/1‑76 (MQ=255)
                                             acaGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATc                         >  2:1865803/1‑76 (MQ=255)
                                             acaGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATc                         >  2:744782/1‑76 (MQ=255)
                                              caGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATcc                        >  2:1246235/1‑76 (MQ=255)
                                               aGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCg                       >  2:1655983/1‑76 (MQ=255)
                                               aGACATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCg                       <  2:1816802/76‑1 (MQ=255)
                                                 aCATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGcc                     <  2:659571/76‑1 (MQ=255)
                                                  cATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGcca                    <  1:865661/76‑1 (MQ=255)
                                                  cATGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGcca                    <  1:668279/76‑1 (MQ=255)
                                                  cATGAAAGTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGcca                    <  1:668258/76‑1 (MQ=255)
                                                   aTGAAATTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGccaa                   >  1:36725/1‑76 (MQ=255)
                                                      aaaTTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGccaacca                <  2:1790747/76‑1 (MQ=255)
                                                       aaTTGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGccaaccaa               >  1:1577466/1‑76 (MQ=255)
                                                         ttGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCaaaa             <  1:818895/76‑1 (MQ=255)
                                                          tGGTAAGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAAc            >  2:1413336/1‑76 (MQ=255)
                                                          tGGTAAGACCAATTGACTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAAc            <  1:963246/76‑1 (MQ=255)
                                                              aaGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACtctc        <  1:512012/76‑1 (MQ=255)
                                                              aaGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACtctc        >  1:723619/1‑76 (MQ=255)
                                                               aGACCAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTcc       <  2:1110576/76‑1 (MQ=255)
                                                                  ccAATTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGAt    <  1:287317/76‑1 (MQ=255)
                                                                    aaTTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGAtgt  <  2:1487856/76‑1 (MQ=255)
                                                                    aaTTGGCTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGAtgt  <  2:1632077/76‑1 (MQ=255)
                                                                         |                                                                      
AAAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGT  >  REL606/124793‑124936

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