Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,108,495 T→G E3A (GAA→GCA)  metJ ← transcriptional repressor protein MetJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,108,4950TG96.8% 256.9 / ‑1.2 63E3A (GAA→GCA) metJtranscriptional repressor protein MetJ
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base G (35/26);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

AACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGA  >  REL606/4108422‑4108565
                                                                         |                                                                      
aaCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCa                                                                      >  2:997057/1‑76 (MQ=255)
aaCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCa                                                                      >  2:469033/1‑76 (MQ=255)
aaCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCa                                                                      >  2:469022/1‑76 (MQ=255)
aaCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCa                                                                      <  1:489574/76‑1 (MQ=255)
 aCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAg                                                                     >  1:1047480/1‑76 (MQ=255)
 aCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAg                                                                     >  1:1093357/1‑76 (MQ=255)
   cGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGcc                                                                   <  1:107650/76‑1 (MQ=255)
      aaTCTTTTTGTCTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATg                                                                <  1:603677/76‑1 (MQ=255)
      aaTCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATg                                                                >  1:26419/1‑76 (MQ=255)
      aaTCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATg                                                                >  1:1717577/1‑76 (MQ=255)
      aaTCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATg                                                                <  1:1648424/76‑1 (MQ=255)
      aaTCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATg                                                                <  1:1624418/76‑1 (MQ=255)
       aTCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATga                                                               <  1:1245658/76‑1 (MQ=255)
       aTCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATga                                                               <  1:63630/76‑1 (MQ=255)
          tttttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGAt                                                            <  1:514876/76‑1 (MQ=255)
          tttttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGAt                                                            <  1:869877/76‑1 (MQ=255)
          tttttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGAt                                                            <  1:1606566/76‑1 (MQ=255)
           ttttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATa                                                           <  1:90512/76‑1 (MQ=255)
           ttttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATa                                                           <  1:550073/76‑1 (MQ=255)
            tttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATAc                                                          >  1:293712/1‑76 (MQ=255)
              tGACTTTTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACtt                                                        <  2:3623/76‑1 (MQ=255)
              tGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACtt                                                        <  1:1621345/76‑1 (MQ=255)
                 cTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAAt                                                     <  2:1434825/76‑1 (MQ=255)
                   tGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATcc                                                   >  1:1572549/1‑76 (MQ=255)
                    gTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCct                                                  >  2:1254314/1‑76 (MQ=255)
                     ttCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCctc                                                 >  1:251769/1‑76 (MQ=255)
                     ttCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCctc                                                 >  2:1177463/1‑76 (MQ=255)
                       cACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCtt                                               >  1:1363289/1‑76 (MQ=255)
                       cACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCtt                                               >  2:228305/1‑76 (MQ=255)
                       cACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCtt                                               >  1:1226759/1‑76 (MQ=255)
                       cACTCTTCTTGCAGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCctctc                                               >  2:460068/1‑75 (MQ=255)
                           cttcttGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTc                                           >  2:857376/1‑76 (MQ=255)
                              cttGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAAt                                        >  1:1163274/1‑76 (MQ=255)
                               ttGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATa                                       >  2:826/1‑76 (MQ=255)
                               ttGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATa                                       <  1:1411301/76‑1 (MQ=255)
                               ttGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAAGCCTCTTCGTCAATa                                       >  2:1657600/1‑76 (MQ=255)
                                tGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATaa                                      >  2:1639444/1‑76 (MQ=255)
                                 gCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATaaa                                     >  2:1163529/1‑76 (MQ=255)
                                 gCCGTGCTCAGCGAATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTAAATCCTCCTCACCaaaaaa                                     >  1:534312/1‑76 (MQ=255)
                                  ccGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTTAATAAAt                                    <  1:1144723/76‑1 (MQ=255)
                                    gTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTg                                  <  2:1411538/76‑1 (MQ=255)
                                    gTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTg                                  >  1:1738211/1‑76 (MQ=255)
                                    gTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTg                                  >  2:1652530/1‑76 (MQ=255)
                                    gTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTNCGTCAATAAATTg                                  <  1:1190900/76‑1 (MQ=255)
                                      gCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTgag                                >  2:704178/1‑76 (MQ=255)
                                         cAGCGTATGGGCTGATATATTCGCCGCTACATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacc                             >  2:839066/1‑76 (MQ=255)
                                           gCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccag                           <  1:835373/76‑1 (MQ=255)
                                            cGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccaga                          >  2:1496725/1‑76 (MQ=255)
                                            cGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccaga                          >  2:456702/1‑76 (MQ=255)
                                             gTATGGGCTGATATATTCGCCGCTCCGTGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagac                         >  1:478177/1‑76 (MQ=255)
                                             gTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagac                         >  2:784258/1‑76 (MQ=255)
                                             gTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagcc                         >  2:314863/1‑74 (MQ=255)
                                             gTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagac                         <  2:594227/76‑1 (MQ=255)
                                                tGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATCAATTGAGACCAGACcac                      >  1:787401/1‑76 (MQ=255)
                                                    cTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAg                  <  1:1279691/76‑1 (MQ=255)
                                                     tGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGt                 >  1:530254/1‑76 (MQ=255)
                                                     tGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGt                 >  1:1591030/1‑76 (MQ=255)
                                                      gATATATACGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGAACACACAgat                >  2:775468/1‑74 (MQ=255)
                                                       atatatTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGCCCACACAGTTg               <  2:28739/76‑1 (MQ=255)
                                                       atatatTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTg               <  1:1242379/76‑1 (MQ=255)
                                                         atatTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAt             >  1:1674321/1‑76 (MQ=255)
                                                          tatTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAtg            >  1:37674/1‑76 (MQ=255)
                                                           atTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAtgt           <  1:1163529/76‑1 (MQ=255)
                                                              cgccgcTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTggg        <  2:55890/76‑1 (MQ=255)
                                                                  gcTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACt    >  2:24903/1‑76 (MQ=255)
                                                                    tCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGa  >  1:31578/1‑76 (MQ=255)
                                                                         |                                                                      
AACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGA  >  REL606/4108422‑4108565

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